回归分析MATLAB实现.pptVIP

图3.21 残差示意图 % 删除异常点程序并建模 [b1,bint1,r1,rint1,stats1]=regress([y(1:12);y(14:m)],[ones(m-1,1),[x(1:12,:);x(14:m,:)]]) rcoplot(r1,rint1); 删除异常点后，残差示意图如图2-21所示，此时没有异常点，改进回归模型的系数、系数置信区间与统计量参见表3.10 表3.10改进回归模型的系数、系数置信区间与统计量 回归系数 回归系数估计值 回归系数置信区间 ?0 328.4616 [290.6145 ,366.3087] ?1 -4.3594 [-4.8880 ,-3.8308] ?2 0.0010 [0.00073,0.0012] R2= 0.9931，F= 858.5846，p0.0001，s2 =0.0150 我们将表3.9与表3.10加以比较，可以发现：可决系数从0.9863提高到0.9931，F统计量从468.4118提高到858.5846，由此可知改进后的模型显著性提高。 图 3-22删除异常点后残差示意图 图3.21 残差示意图 3.2.3 多元线性回归实例 例3.2.3现代服务业是社会分工不断深化的产物，随着经济的发展，科学技术的进步，现代服务业的发展受到多种因素和条件的影响。不仅受到经济总体发展水平的影响，还受到第二产业、就业、投入等因素的影响，从这几个主要方面出发，利用江苏省统计年鉴的有关数据，通过建立多元线性回归模型对1990-2008年各种因素对现代服务业的影响进行回归分析。假如构建如下江苏省服务业增长模型： Y代表江苏省服务业的增加值(单位：亿元)，反映了江苏省服务业发展的总体水平。 x1~x4表示影响江苏省服务业发展的四种主要因素和影响，其中 x1代表江苏省人均GDP(单位：元)，说明江苏省总体经济发展水平对服务业的影响； x2代表江苏省第二产业的增加值(单位：亿元)，说明了工业发展对服务业的影响，体现了生产性服务业的需求规模； x3表示江苏省服务业的就业人数(单位：万人)； x4表示江苏省服务业资本形成总额(单位：亿元)，主要体现服务业投资的经济效应。 表3.11 江苏省关于服务业发展及各影响因素相关数据 年份 Y服务业增加值 省人均GDP 第二产业增加值 服务业就业人数 服务业资本形成总额 1989 37.76 2038 70.24 589.74 252.01 1990 28.13 2109 35.53 623.19 275.82 1991 93.58 2353 101.33 640.95 330.71 1992 160.62 3106 325.34 706.39 439.32 1993 286.58 4321 478.79 786.37 620.97 1994 277.12 5801 588.72 855.97 858.91 1995 387.11 7319 528.49 920.45 1102.71 1996 367.16 8471 358.86 975.66 1293.43 1997 291.77 9371 337.74 1025.22 1370.21 1998 280.01 10049 228.24 1102.31 1624.74 1999 227.61 10695 280.05 1151.68 1773.37 2000 329.16 11765 515.74 1192.02 1903.37 2001 385.44 12882 471.57 1263.77 2131.87 2002 437.02 14396 697.03 1341.86 2189.78 2003 601.39 16830 1182.62 1407.63 2686.57 2004 704.72 20223 1650.88 1443.37 3362.19 2005 1291.11 24560 1917.05 1542.46 3930.56 2006 1360.09 28814 1895.8 1625.06 4628.59 2007 1769.28 33928 2055.56 1713.33 5287.91 解：%输入各影响因素的数据 x0=[2038 70.24 589.74 252.01 2109 35.53 623.19 275.82 ……………………………………………………………………… 33928 2055.56 1713.33 5287.91]; y=[37.76,28.13,93.58,160.62,286.58,277.12,387.11,367.16,291.77,280.01,227.61,329.16,385.4