4.1 假设在模型Yi12X2i3X3iui中,X2与X3之间的相关系数为零,于是有人建议你进行如下回归:
Yi12X2iu1iY
i13X3iu2i(1)是否存在ˆ2ˆ2且ˆ3ˆ3?为什么? (2)ˆ1会等于ˆ1或ˆ1或两者的某个线性组合吗? (3)是否有varˆ2varˆ2且varˆ3varˆ3?
练习题4.1参考解答:
(1) 存在ˆ2ˆ2且ˆ3ˆ3。 2因为ˆix2i3ii3i2ix3i2yxyxxx2x2xx2
2i3i2i3i当X2与X3之间的相关系数为零时,离差形式的
x2ix3i0
2有ˆi2i3ii2i2yxxyxx2x22ˆ2 2i3ix2i同理有:ˆ3ˆ3 (2) ˆ1会等于ˆ1或ˆ1的某个线性组合 因为 ˆ1Yˆ2X2ˆ3X,且3ˆ1Yˆ2X2,ˆ1Yˆ3X3 由于ˆ2ˆ2且ˆ3ˆ3,则 ˆ1Yˆ2X2Yˆ2X2ˆYˆ12X 2 ˆYˆ11Yˆ3X3Yˆ3X3ˆ3X 3则 ˆYˆXˆXYYˆ1Yˆ112233XX2X3ˆ1ˆ1Y2X3(3) 存在varˆ2varˆ2且varˆ3varˆ3。
1
ˆ因为var2x1r22i2232
ˆ当r230时,var2x1rx22i2232222iˆ2 varˆvarˆ3 同理,有var3
4.2在决定一个回归模型的“最优”解释变量集时人们常用逐步回归的方法。在逐步回归中既可采取每次引进一个解释变量的程序(逐步向前回归),也可以先把所有可能的解释变量都放在一个多元回归中,然后逐一地将它们剔除(逐步向后回归)。加进或剔除一个变量,通常是根据F检验看其对ESS的贡献而作出决定的。根据你现在对多重共线性的认识,你赞成任何一种逐步回归的程序吗?为什么?
练习题4.2参考解答:
根据对多重共线性的理解,逐步向前和逐步向后回归的程序都存在不足。逐步向前法不能反映引进新的解释变量后的变化情况,即一旦引入就保留在方程中;逐步向后法则一旦某个解释变量被剔出就再也没有机会重新进入方程。而解释变量之间及其与被解释变量的相关关系与引入的变量个数及同时引入哪些变量而呈现出不同,所以要寻找到“最优”变量子集则采用逐步回归较好,它吸收了逐步向前和逐步向后的优点。
4.3 下表给出了中国商品进口额Y、国内生产总值GDP、居民消费价格指数CPI。
表4.11 中国商品进口额、国内生产总值、居民消费价格指数 年份 1985 1986 1987 1988 1989 1990 1991 1992 1993
商品进口额 (亿元) 1257.8 1498.3 1614.2 2055.1 2199.9 2574.3 3398.7 4443.3 5986.2 国内生产总值 (亿元) 9016.0 10275.2 12058.6 15042.8 16992.3 18667.8 21781.5 26923.5 35333.9 居民消费价格指数(1985=100) 100.0 106.5 114.3 135.8 160.2 165.2 170.8 181.7 208.4 2
1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 9960.1 11048.1 11557.4 11806.5 11626.1 13736.4 18638.8 20159.2 24430.3 34195.6 46435.8 54273.7 63376.9 73284.6 48197.9 60793.7 71176.6 78973.0 84402.3 89677.1 99214.6 109655.2 120332.7 135822.8 159878.3 183084.8 211923.5 249529.9 258.6 302.8 327.9 337.1 334.4 329.7 331.0 333.3 330.6 334.6 347.7 353.9 359.2 376.5 资料来源:《中国统计年鉴》,中国统计出版社2000年、2008年。
请考虑下列模型:lnYt1+2lnGDPt3lnCPItui 1)利用表中数据估计此模型的参数。 2)你认为数据中有多重共线性吗? 3)进行以下回归:
lnYtA1+A2lnGDPtv1ilnYtB1+B2lnCPItv2ilnGDPtC1C2lnCPItv3i根据这些回归你能对数据中多重共线性的性质说些什么?
ˆ和ˆ在5%水平上个别地显著,并且总的F检验也是显4)假设数据有多重共线性,但23著的。对这样的情形,我们是否应考虑共线性的问题?
练习题4.3参考解答: (1) 参数估计结果如下
3
ln(进口)3.0601.657ln(GDP)1.057ln(CPI) (0.337) (0.092) (0.215)R20.992 R20.991 F1275.093(括号内为标准误)
(2)居民消费价格指数的回归系数的符号不能进行合理的经济意义解释,且且CPI与进口之间的简单相关系数呈现正向变动。可能数据中有多重共线性。 计算相关系数:
(3)最大的CI=108.812,表明GDP与CPI之间存在较高的线性相关。
(4)分别拟合的回归模型如下:
lnY4.09071.2186ln(GDP) t= (-10.6458) (34.6222)R20.9828 R20.9820 F1198.698
4
lnY5.44242.6637ln(CPI) t= (-4.3412) (11.6809)R20.8666 R20.8603 F136.4437
ln(GDP)1.43802.2460ln(CPI) t=(-1.9582) (16.8140)R20.9309 R20.9276 F282.7107单方程拟合效果都很好,回归系数显著,可决系数较高,GDP和CPI对进口分别有显著的单一影响,在这两个变量同时引入模型时影响方向发生了改变,这只有通过相关系数的分析才能发现。
(5)如果仅仅是作预测,可以不在意这种多重共线性,但如果是进行结构分析,还是应该引起注意。
4.4 自己找一个经济问题来建立多元线性回归模型,怎样选择变量和构造解释变量数据矩阵X才可能避免多重共线性的出现?
练习题4.4参考解答:
本题很灵活,主要应注意以下问题:
5
(1)选择变量时要有理论支持,即理论预期或假设;变量的数据要足够长,被解释变量与解释变量之间要有因果关系,并高度相关。
(2)建模时尽量使解释变量之间不高度相关,或解释变量的线性组合不高度相关。
4.5 克莱因与戈德伯格曾用1921-1950年(1942-1944年战争期间略去)美国国内消费Y和工资收入X1、非工资—非农业收入X2、农业收入X3的时间序列资料,利用OLSE估计得出了下列回归方程:
ˆ8.1331.059X10.452X20.121X3Y (8.92) (0.17) (0.66) (1.09) R20.95 F107.37括号中的数据为相应参数估计量的标准误差。试对上述模型进行评析,指出其中存在的问
题。
练习题4.5参考解答:
从模型拟合结果可知,样本观测个数为27,消费模型的判定系数R0.95,F统计量为107.37,在0.05置信水平下查分子自由度为3,分母自由度为23的F临界值为3.028,计算的F值远大于临界值,表明回归方程是显著的。模型整体拟合程度较高。 依据参数估计量及其标准误,可计算出各回归系数估计量的t统计量值:
2t08.1330.91,8.92t11.0596.10,0.17t20.4520.69,0.66t30.1210.111.09除t1外,其余的tj值都很小。工资收入X1的系数的t检验值虽然显著,但该系数的估计值过大,该值为工资收入对消费边际效应,因为它为1.059,意味着工资收入每增加一美元,消费支出的增长平均将超过一美元,这与经济理论和常识不符。
另外,理论上非工资—非农业收入与农业收入也是消费行为的重要解释变量,但两者的t检验都没有通过。这些迹象表明,模型中存在严重的多重共线性,不同收入部分之间的相互关系,掩盖了各个部分对解释消费行为的单独影响。
4.6 理论上认为影响能源消费需求总量的因素主要有经济发展水平、收入水平、产业发展、人民生活水平提高、能源转换技术等因素。为此,收集了中国能源消费总量Y (万吨标准煤)、国民总收入(亿元)X1(代表收入水平)、国内生产总值 (亿元)X2(代表经济发展水平)、工业增加值(亿元)X3、建筑业增加值(亿元)X4、交通运输邮电业增加值(亿元)X5(代表产业发
6
展水平及产业结构)、人均生活电力消费 (千瓦小时)X6(代表人民生活水平提高)、能源加工转换效率(%)X7(代表能源转换技术)等在1985-2007年期间的统计数据,具体如表4.2所示。
表4.12 1985~2007年统计数据
年份 能源消费 y 1985 1986 1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 76682 80850 86632 92997 96934 98703 103783 109170 115993 122737 131176 138948 137798 132214 133831 138553 143199 151797 174990 203227 223319 246270 265583 国民 国内生 工业 建筑业 交通运输邮电 人均生活 能源加工 总收入 产总值 增加值 增加值 增加值 电力消费 转换效率 X1 X2 X3 X4 X5 X6 X7 9040.7 10274.4 12050.6 15036.8 17000.9 18718.3 21826.2 26937.3 35260 48108.5 59810.5 70142.5 77653.1 83024.3 88189 98000.5 108068.2 119095.7 135174 159586.7 183956.1 9016 10275.2 12058.6 15042.8 16992.3 18667.8 21781.5 26923.5 35333.9 48197.9 60793.7 71176.6 78973 84402.3 89677.1 99214.6 109655.2 120332.7 135822.8 159878.3 183084.8 3448.7 417.9 3967 525.7 4585.8 665.8 5777.2 810 6484 794 6858 859.4 8087.1 1015.1 10284.5 1415 14188 2266.5 19480.7 2964.7 24950.6 3728.8 29447.6 4387.4 32921.4 4621.6 34018.4 4985.8 35861.5 5172.1 4003.6 5522.3 43580.6 5931.7 47431.3 6465.5 54945.5 7490.8 65210 8694.3 76912.9 10133.8 406.9 475.6 544.9 661 786 1147.5 1409.7 1681.8 2205.6 2898.3 3424.1 4068.5 4593 5178.4 5821.8 7333.4 8406.1 9393.4 10098.4 12147.6 10526.1 12481.1 14604.1 21.3 23.2 26.4 31.2 35.3 42.4 46.9 54.6 61.2 72.7 83.5 93.1 101.8 106.6 118.2 132.4 144.6 156.3 173.7 190.2 216.7 249.4 274.9 68.29 68.32 67.48 66.54 66.51 67.2 65.9 66.00 67.32 65.2 71.05 71.5 69.23 69.44 69.19 69.04 69.03 69.04 69.4 70.71 71.08 71.24 71.25 213131.7 211923.5 91310.9 11851.1 251483.2 249529.9 107367.2 14014.1 资料来源:《中国统计年鉴》,中国统计出版社2000、2008年版。
要求:
1)建立对数多元线性回归模型,分析回归结果。
2)如果决定用表中全部变量作为解释变量,你预料会遇到多重共线性的问题吗?为什么? 3)如果有多重共线性,你准备怎样解决这个问题?明确你的假设并说明全部计算。
练习题4.6参考解答:
(1)建立对数线性多元回归模型,引入全部变量建立对数线性多元回归模型如下: 生成: lny=log(y), 同样方法生成: lnx1,lnx2,lnx3,lnx4,lnx5,lnx6,lnx7.
7
作全部变量对数线性多元回归,结果为:
从修正的可决系数和F统计量可以看出,全部变量对数线性多元回归整体对样本拟合很好,,各变量联合起来对能源消费影响显著。可是其中的lnX3、lnX4、lnX6对lnY影响不显著,而且lnX2、lnX5的参数为负值,在经济意义上不合理。所以这样的回归结果并不理想。
(2) 预料此回归模型会遇到多重共线性问题, 因为国民总收入与GDP本来就是一对关联指标;而工业增加值、建筑业增加值、交通运输邮电业增加值则是GDP的组成部分。这两组指标必定存在高度相关。
解释变量国民总收入(亿元)X1(代表收入水平)、国内生产总值(亿元)X2(代表经济发展水平)、工业增加值(亿元)X3、建筑业增加值(亿元)X4、交通运输邮电业增加值(亿元)X5(代表产业发展水平及产业结构)、人均生活电力消费 (千瓦小时)X6(代表人民生活水平提高)、能源加工转换效率(%)X7(代表能源转换技术)等很可能线性相关,计算相关系数如下:
可以看出lnx1与lnx2、lnx3、lnx4、lnx5、lnx6之间高度相关,许多相关系数高于0.900以上。如果决定用表中全部变量作为解释变量,很可能会出现严重多重共线性问题。 (3)因为存在多重共线性,解决方法如下:
A:修正理论假设,在高度相关的变量中选择相关程度最高的变量进行回归建立模型:而对变量取对数后,能源消费总量的对数与人均生活电力消费的对数相关程度最高,可建立
8
这两者之间的回归模型。如
lny9.9320.421lnx6 (0.116) (0.026)R20.926 R20.922 F261.551B:进行逐步回归,直至模型符合需要研究的问题,具有实际的经济意义和统计意义。采用逐步回归的办法,去检验和解决多重共线性问题。分别作lnY对
lnX1,lnX2,lnX3,lnX4,lnX5,lnX6,lnX7的一元回归,结果如下:
一元回归结果:
变量 lnX1 lnX2 0.315 14.62 0.911 0.906 lnX3 0.277 9.718 0.818 0.809 lnX4 0.297 13.22 0.893 0.888 lnX5 0.273 11.717 0.867 0.861 lnX6 0.421 16.173 0.926 0.922 lnX7 8.73 4.648 0.507 0.484 参数估计值 0.316 t统计量 可决系数 调整可决系数 14.985 0.914 0.910 其中加入lnX6的方程调整的可决系数最大, 以lnX6为基础, 顺次加入其他变量逐步回归。结果如下表:
变量 lnX6 lnX1 lnX6 lnX2 lnX6 lnX3 lnX6 lnX4 lnX6 lnX5 lnX6 lnX7 lnX1 -0.186 (-0.698) -0.251 (-1.021) 0.061 (1.548) -0.119 (-0.897) -0.623 (-7.127) lnX2 lnX3 lnX4 lnX5 lnX6 0.666 (1.891) 0.753 (2.308) 0.341 (5.901) 0.585 (3.167) 1.344 (10.314) 0.391 (11.071) 0.924 0.977 0.921 0.927 0.922 lnX7 R2 0.920 经比较,新加入lnX5的方程调整可决系数改进最大, 各参数的t检验也都显著,但是lnX5参数
的符号与经济意义不符合。若再加入其他变量后的逐步回归,若剔除不显著的变量和无经济意义的变量后, 仍为第一步所建只包含lnX6的一元回归模型。
如果需要建立多元线性回归模型,则需寻找新的变量或改变模型形式。 例如, 不取对数作全部变量多元线性回归,结果为:
9
可以看出还是有严重多重共线性。作逐步回归: 分别作一元回归得到: 变量 参数估计值 t 统计量 X1 0.7333 26.4698 0.9709 0.9695 X1 6.6399(0.0022) 0.5512 (0.0000) 0.5040 (0.3356) 1.0516 (0.0000) 1.0075 (0.0088) 0.7499 (0.0000) X2 0.7353 0.9684 0.9669 X2 -5.9308 (0.0054) X3 1.6655 0.9393 0.9364 X3 0.4349 (0.0821) X4 X4 13.1909 25.9636 0.9697 0.9683 X5 X5 10.8980 13.5147 0.8969 0.8920 X6 -255.80 (0.438) X6 678.0058 22.4229 0.9599 0.9580 X7 -813.44 (0.5988) X7 19332.30 4.7024 0.5129 0.4897 25.3627 18.0257 R2 R2 X1,X2 X1,X3 X1,X4 X1,X5 X1,X6 X1,X7 以X1为基础加入其他变量, 结果为: R2 0.9785 0.9726 0.9683 0.9766 0.9690 0.9684 4.1326 (0.6580) -5.0269 (0.013) 注: 括号中为p值. 可以发现加入X2、X5、X6、X7后参数的符号不合理,加入X4后并不显著。只有加入X3后修正的可决系数有所提高,而且参数符号的经济意义合理, X3参数估计值的p值为0.0821,在10%的显著性水平下是显著的。所以相对较为合理的模型估计结果可以为:
10
4.7 在本章开始的“引子”提出的“农业的发展反而会减少财政收入吗?”的例子中,如果所采用的数据如下表所示
表4.13 1978-2007年财政收入及其影响因素数据
建筑业增加受灾面积财政收入农业增加值工业增加值总人口(万最终消费年份 值(亿(千公(亿元)CS (亿元)NZ (亿元)GZ 人)TPOP (亿元)CUM 元)JZZ 顷)SZM 1978 1132.3 1027.5 1607 138.2 96259 2239.1 50790 1979 1146.4 1270.2 1769.7 143.8 97542 2633.7 39370 1980 1159.9 1371.6 1996.5 195.5 98705 3007.9 44526 1981 1175.8 1559.5 2048.4 207.1 100072 3361.5 39790 1982 1212.3 1777.4 2162.3 220.7 101654 37148 33130 1983 1367 1978.4 2375.6 270.6 103008 4126.4 34710 1984 1642.9 2316.1 2789 316.7 104357 4846.3 31890 1985 2004.8 2564.4 3448.7 417.9 105851 5986.3 44365 1986 2122 2788.7 3967 525.7 107507 6821.8 47140 1987 2199.4 3233.0 4585.8 665.8 109300 7804.6 42090 1988 2357.2 3865.4 5777.2 810 111026 9839.5 50870 1989 2664.9 4265.9 6484 794 112704 11164.2 46991 1990 2937.1 5062.0 6858 859.4 114333 12090.5 38474 1991 3149.48 5342.2 8087.1 1015.1 115823 14091.9 55472 1992 3483.37 5866.6 10284.5 1415 117171 17203.3 51333 1993 4348.95 6963.8 14188 2266.5 118517 21899.9 48829 1994 5218.1 9572.7 19480.7 2964.7 119850 29242.2 55043 1995 6242.2 12135.8 24950.6 3728.8 121121 36748.2 45821 1996 7407.99 14015.4 29447.6 4387.4 122389 43919.5 46989 1997 8651.14 14441.9 32921.4 4621.6 123626 48140.6 53429 1998 9875.95 14817.6 34018.4 4985.8 124761 51588.2 50145 1999 11444.08 14770.0 35861.5 5172.1 125786 55636.9 49981 2000 13395.23 14944.7 4003.6 5522.3 126743 61516 54688 2001 16386.04 15781.3 43580.6 5931.7 127627 66878.3 52215 2002 18903.64 16537.0 47431.3 6465.5 128453 71691.2 47119 2003 21715.25 17381.7 54945.5 7490.8 129227 77449.5 54506 2004 26396.47 21412.7 65210 8694.3 129988 87032.9 37106 2005 31649.29 22420.0 76912.9 10133.8 130756 96918.1 38818 2006 38760.20 24040.0 91310.9 11851.1 131448 110595.3 41091 2007 51321.78 28095.0 107367.2 14014.1 132129 128444.6 48992 (资料来源:《中国统计年鉴2008》,中国统计出版社2008年版) 试分析:为什么会出现本章开始时所得到的异常结果?怎样解决所出现的问题?
练习题4.7参考解答:
(1)根据样本数据得到各解释变量的样本相关系数矩阵如下: 样本相关系数矩阵
11
解释变量之间相关系数较高,特别是农业增加值、工业增加值、建筑业增加值、最终消费之间,相关系数都在0.9以上。这显然与第三章对模型的无多重共线性假定不符合。(2)解决方案:
采用逐步回归的方式,可以得到没有共线性的回归模型,但可能存在设定偏误。 合并工业增加值与建筑业增加值,得到财政收入与第二产业的回归。 取对数再回归,可以减低共线性。
12
因篇幅问题不能全部显示,请点此查看更多更全内容