多重共线性处理经典例题.doc

理论上认为影响能源消费需求总量的因素主要有经济发展水平、收入水平、产业发展、人民生活水平提高、能源转换技术等因素。为此，收集了中国能源消费总量Y (万吨标准煤)、国民总收入(亿元)X1(代表收入水平)、国内生产总值 (亿元)X2(代表经济发展水平)、工业增加值(亿元)X3、建筑业增加值(亿元)X4、交通运输邮电业增加值(亿元)X5(代表产业发展水平及产业结构)、人均生活电力消费 (千瓦小时)X6(代表人民生活水平提高)、能源加工转换效率(%)X7(代表能源转换技术)等在1985-2007年期间的统计数据，具体如表4.2所示。表4.12 19852007年统计数据年份能源消费国民总收入国内生产总值工业增加值建筑业增加值交通运输邮电增加值人均生活电力消费能源加工转换效率yX1X2X3X4X5X6X71985766829040.790163448.7417.9406.921.368.2919868085010274.410275.23967525.7475.623.268.3219878663212050.612058.64585.8665.8544.926.467.4819889299715036.815042.85777.281066131.266.5419899693417000.916992.3648479478635.366.5119909870318718.318667.86858859.41147.542.467.2199110378321826.221781.58087.11015.11409.746.965.9199210917026937.326923.510284.514151681.854.666.0019931159933526035333.9141882266.52205.661.267.32199412273748108.548197.919480.72964.72898.372.765.2199513117659810.560793.724950.63728.83424.183.571.05199613894870142.571176.629447.64387.44068.593.171.519971377987765344621.64593101.869.23199813221483024.384402.334018.44985.85178.4106.669.4419991338318818989677.135861.55172.15821.8118.269.19200013855398000.599214.64003.65522.37333.4132.469.042001143199108068.2109655.243580.65931.78406.1144.669.032002151797119095.7120332.747431.36465.59393.4156.369.042003174990135174135822.854945.57490.810098.4173.769.42004203227159586.7159878.3652108694.312147.6190.270.712005223319183956.1183084.876912.910133.810526.1216.771.082006246270213131.7211923.591310.911851.112481.1249.471.242007265583251483.2249529.9107367.214014.114604.1274.971.25资料来源：中国统计年鉴，中国统计出版社2000、2008年版。要求：1)建立对数多元线性回归模型，分析回归结果。2)如果决定用表中全部变量作为解释变量，你预料会遇到多重共线性的问题吗？为什么？3)如果有多重共线性，你准备怎样解决这个问题？明确你的假设并说明全部计算。参考解答：(1)建立对数线性多元回归模型，引入全部变量建立对数线性多元回归模型如下:生成: lny=log(y), 同样方法生成: lnx1,lnx2,lnx3,lnx4,lnx5,lnx6,lnx7.作全部变量对数线性多元回归,结果为: 从修正的可决系数和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：修正理论假设，在高度相关的变量中选择相关程度最高的变量进行回归建立模型：而对变量取对数后，能源消费总量的对数与人均生活电力消费的对数相关程度最高，可建立这两者之间的回归模型。如B：进行逐步回归，直至模型符合需要研究的问题，具有实际的经济意义和统计意义。采用逐步回归的办法，去检验和解决多重共线性问题。分别作对的一元回归，结果如下:一元回归结果：变量lnX1lnX2lnX3lnX4lnX5lnX6lnX7参数估计值0.3160.3150.2770.2970.2730.4218.73t统计量14.98514.629.71813.2211.71716.1734.648可决系数0.9140.9110.8180.8930.8670.9260.507调整可决系数0.9100.9060.8090.8880.8610.9220.484其中加入lnX6的方程调整的可决系数最大, 以lnX6为基础, 顺次加入其他变量逐步回归。结果如下表: 变量lnX1lnX2lnX3lnX4lnX5lnX6lnX7lnX6lnX1-0.186(-0.698)0.666(1.891)0.920lnX6lnX2-0.251(-1.021)0.753(2.308)0.922lnX6lnX30.061(1.548)0.341(5.901)0.927lnX6lnX4-0.119(-0.897)0.585(3.167)0.921lnX6lnX5-0.623(-7.127)1.344(10.314)0.977lnX6lnX70.391(11.071)0.924经比较，新加入lnX5的方程调整可决系数改进最大, 各参数的t检验也都显著，但是lnX5参数的符号与经济意义不符合。若再加入其他变量后的逐步回归,若剔除不显著的变量和无经济意义的变量后, 仍为第一步所建只包含lnX6的一元回归模型。如果需要建立多元线性回归模型，则需寻找新的变量或改变模型形式。例如, 不取对数作全部变量多元线性回归,结果为:可以看出还是有严重多重共线性。作逐步回归:分别作一元回归得到：变量 X1 X2X3X4X5X6X7参数估计值0.73330.73531.665513.190910.8980678.005819332.30t 统计量26.469825.362718.025725.963613.514722.42294.70240.97090.96840.93930.96970.89690.95990.51290.96950.96690.93640.96830.89200.95800.4897以X1为基础加入其他变量, 结果为:X1X2X3X4X5X6X7X1，X26.6399（0.0022）-5.9308(0.0054)0.9785X1，X30.5512(0.0000)0.4349(0.0821)0.9726X1，X40.5040(0.3356)4.1326(0.6580)0.9683X1，X51.0516(0.0000)-5.0269(0.013)0.9766X1，X61.0075(0.0088)-255.80