多重共线性的检验与处理.doc

实验名称：多重共线性的检验与处理实验时间：2011.12.10实验要求：主要是学习多重共线性的检验与处理，主要是研究解释变量与其余解释变量之间有严重多重共线性的模型，分析变量之间的相关系数。通过具体案例建立 模型，然后估计参数，求出相关的数据。再对模型进行检验，看数据之间是否存在多重共线性。最后利用所求出的模型来进行修正。实验内容：实例：我国钢材供应量分析通过分析我国改革开放以来（19781997）钢材供应量的历史资料，可以建立一个单一方程模型。根据理论及对现实情况的认识，影响我国钢材供应量 Y（万吨）的主要因素有：原油产量X1（万吨），生铁产量X2（万吨），原煤产量X3（万吨），电力产量X4（亿千瓦小时），固定资产投资X5(亿元)， 国内生产总值X6(亿元)，铁路运输量X7（万吨）。（一）建立我国钢材供应量的计量经济模型：（二）估计模型参数，结果为：Dependent Variable: YMethod: Least SquaresDate: 11/02/09 Time: 16:09Sample: 1978 1997Included observations: 20VariableCoefficientStd. Errort-StatisticProb.C139.2362718.24930.1938550.8495X1-0.0519540.090753-0.5724830.5776X20.1275320.1324660.9627510.3547X3-24.2942797.48792-0.2492030.8074X40.8632830.1867984.6214750.0006X50.3309140.1055923.1338890.0086X6-0.0700150.025490-2.7467550.0177X70.0023050.0190870.1207800.9059R-squared0.999222 Mean dependent var5153.350Adjusted R-squared0.998768 S.D. dependent var2511.950S.E. of regression88.17626 Akaike info criterion12.08573Sum squared resid93300.63 Schwarz criterion12.48402Log likelihood-112.8573 F-statistic2201.081Durbin-Watson stat1.703427 Prob(F-statistic)0.000000由此可见，该模型 可绝系数很高，F检验值2201.081，明显显著。但当 ， 系数的t检验不显著，而且系数的符号与预期的相反，这表明很可能存在严重的多重共线性。（三）计算各解释变量的相关系数，选择 数据，得相关系数矩阵（表3.1）。表3.1 相关系数矩阵X2X3X4X5X6X7X2 1.000000 0.964400 0.994921 0.969686 0.972530 0.931689X3 0.964400 1.000000 0.974809 0.894963 0.913344 0.982943X4 0.994921 0.974809 1.000000 0.959613 0.969105 0.945444X5 0.969686 0.894963 0.959613 1.000000 0.996169 0.827643X6 0.972530 0.913344 0.969105 0.996169 1.000000 0.846079X7 0.931689 0.982943 0.945444 0.827643 0.846079 1.000000由相关系数矩阵可以看出，各解释变量相互之间的相关系数较高，证实确实存在严重多重共线性。（四）修正多重共线性采用逐步回归的办法，去检验和解决多重共线性问题。分别做Y对 的一元回归，结果如表3.2所示表3.2 一元回归结果变量参数估计量1.1817840.926212926.71780.8840470.5724510.1086650.106826t统计量10.1062957.8201715.8724362.4938115.4789216.5453511.455240.8501710.9946450.9333170.9954120.9301230.9383030.8793750.8418470.9943470.9296120.9951570.9262410.9348750.872673其中，加入 的方程 最大，以 为基础，顺次加入其他变量逐步回归，结果如表3.3所示。表3.3 加入新变量的回归结果（一）变量X1X2X3X4X5X6X7X4、X1-0.074944（-1.359593）0.932329（24.46548）0.995375X4、X20.429201（3.895858）0.476623（4.534407）0.997291X4、X3-125.0949（-1.960538）0.996689（16.91041）0.995818X4X50.808630（16.81022）0.052646（1.633786）0.995568X4、X60.858914（14.63712）0.003283（0.441961）0.994931X4、X70.92763921.44406-0.005928（-1.065860）0.995194经比较，新加入 的方程 =0.997291，改进最大，而且各参数的t检验显著，选择保留 ,再加入其他新变量逐步回归，结果如表3.4所示：表3.4 加入新变量的回归结果（二）X1X2X3X4X5X6X7X4、X2、X1-0.052631-1.2556220.4096523.7425050.5290874.7452150.997380X4、X2、X30.3868193.559933-81.25478-1.6042270.5900214.8000560.997521在 、 基础上加入 后的方程 增大，但 参数的t检验不显著，甚至 的符号也变得不合理。加入后的方程 增大，但 参数的t检验不显著，甚至 的符号也变得不合理。加入 后的方程 下降，而且 参数的t检验不显著。加入 后的方程下降，而且 参数的t检验不显著，甚至 的符号也变得不合理。加入 后的方程 下降，而且 参数的t检验不显著，甚至的符号也变得不合理。这说明 、 、 、 、 引起严重多重共线性，应予剔除。