计量经济学实验四 序列相关的检验与修正.doc

实验四 序列相关的检验与修正实验目的1、理解序列相关的含义后果、2、学会序列相关的检验与消除方法实验内容利用下表资料，试建立我国城乡居民储蓄存款模型，并检验模型的自相关性。表3 我国城乡居民储蓄存款与GDP统计资料（1978年100）年份存款余额YGDP指数X年份存款余额YGDP指数X1978210.60100.019895146.90271.31979281.00107.619907034.20281.71980399.50116.019919107.00307.61981523.70122.1199211545.40351.41982675.40133.1199314762.39398.81983892.50147.6199421518.80449.319841214.70170.0199529662.25496.519851622.60192.9199638520.84544.119862237.60210.0199746279.80592.019873073.30234.0199853407.47638.219883801.50260.7一、模型的估计0、准备工作。建立工作文件，并输入数据。1、相关图分析SCAT X Y相关图表明，GDP指数与居民储蓄存款二者的曲线相关关系较为明显。现将函数初步设定为线性、双对数等不同形式，进而加以比较分析。2、估计模型，利用LS命令分别建立以下模型线性模型： LS Y C X (-6.706) (13.862)0.9100 F192.145 S.E5030.809双对数模型：GENR LNY=LOG(Y) GENR LNX=LOG(X) LS LNY C LNX (-31.604) (64.189)0.9954 F4120.223 S.E0.12213、选择模型比较以上模型，可见各模型回归系数的符号及数值较为合理。各解释变量及常数项都通过了检验，模型都较为显著。比较各模型的残差分布表。线性模型的残差在较长时期内呈连续递减趋势而后又转为连续递增趋势，残差先呈连续递增趋势而后又转为连续递减趋势，因此，可以初步判断这种函数形式设置是不当的。而且，这个模型的拟合优度也较双对数模型低，所以又可舍弃线性模型。双对数模型具有很高的拟合优度，因而初步选定回归模型为双对数回归模型。二、模型自相关的检验1.图示法其一，残差序列et的变动趋势图。菜单：QuickGraphline，在对话框中输入resid；或者用命令操作，直接在命令行输入：line X。其二，作et-1和et之间的散点图。菜单：QuickGraphScatter，在对话框中输入resid(-1) resid；或者用命令操作，直接在命令行输入：scat resid(-1) resid。2.DW检验因为n21，k1，取显著性水平0.05时，查表得1.22，1.42，而00.7062DW，所以存在（正）自相关。3.LM(BG)检验在方程窗口中点击View/Residual Test/Series Correlation LM Test，并选择滞后期为2，则会得到如图4-1所示的信息。图4-1 双对数模型的BG检验图中，=11.31531，临界概率P=0.0034，因此辅助回归模型是显著的，即存在自相关性。又因为，的回归系数均显著地不为0，说明双对数模型存在一阶和二阶自相关性。三、自相关的修正（1）自相关系数的估计主要的方法有：A. 根据和DW统计量之间的近似关系，取的估计为：1-DW/2B. 直接取=1C. 采用杜宾两步法估计。LS Y C Y(-1) X X(-1)，Y(-1)的系数估计即为的估计D. 科克伦-奥科特迭代法。首先产生残差序列，命名为e，然后e对其滞后1阶回归（无常数项），LS e e(-1)，e(-1) 的系数估计作为的估计（2）加入AR项在LS命令中加上AR(1)和AR(2)，使用迭代估计法估计模型。键入命令：LS LNY C LNX AR(1) AR(2)则估计结果如图4-2所示。图4-2 加入AR项的双对数模型估计结果 图4-2表明，调整后模型的DW1.6445，n19，k1，取显著性水平0.05时，查表得1.18，1.40，而1.6445DW4，说明模型不存在一阶自相关性；再BG检验（图4-3），也表明不存在高阶自相关性，因此，中国城乡居民储蓄存款的双对数模型为： (-25.263) (52.683)0.9982 F2709.985 S.E0.0744 DW1.6445图4-3习题1. 下表给出了美国1960-1995年36年间个人实际可支配收入X和个人实际消费支出Y的数据。美国个人实际可支配收入和个人实际消费支出单位：100亿美元年份个人实际可支配收入X个人实际消费支出Y年份个人实际可支配收入X个人实际消费支出Y196019611962196319641965196619671968196919701971197219731974197519761977157162169176188200211220230237247256268287285290301311143146153160169180190196207215220228242253251257271283197819791980198119821983198419851986198719881989199019911992199319941995326335337345348358384396409415432440448449461467478493295302301305308324341357371382397406413411422434447458注：资料来源于Economic Report of the President，数据为1992年价格。要求：（1）用普通最小二乘法估计收入消费模型；（2）检验收入消费模型的自相关状况（5%显著水平）；（3）用适当的方法消除模型中存在的问题。2.下表是北京市连续19年城镇居民家庭人均收入与人均支出的数据。 北京市19年来城镇居民家庭收入与支出数据表（单位：元）年份顺序人均收入(元)人均生活消费支出(元)商品零售物价指数(%)人均实际收入(元)人均实际支出(元)12345678910111213141516171819450.18 491.54 599.40 619.57 668.06 716.60 837.65 1158.84 1317.33 1413.24 1767.67 1899.57 2067.33 2359.88 2813.10 3935.39 5585.88 6748.68 7945.78359.86 408.66 490.44 511.43 534.82 574.06 666.75 923.32 1067.38 1147.60 1455.55 1520.41 1646.05 1860.17 2134.65 2939.60 4134.12 5019.76 5729.45100.00 101.50 108.60 110.20 112.30 113.00 115.40 136.80 145.90 158.60 193.30 229.10 238.50 258.80 280.30 327.70 386.40 435.10 466.90450.18 484.28 551.93 562.22 594.89 634.16 725.87 847.11 902.90 891.07 914.47 829.14 866.81 911.85 1003.60 1200.91 1445.62 1551.06 1701.82359.86 402.62 451.60 464.09 476.24 508.02 577.77 674.94 731.58 723.58 753.00 663.64 690.17 718.77 761.56 897.04 1069.91 1153.70 1227.13要求：（1）建立居民收入消费函数； （2）检验模型中存在的问题，并采取适当的补救措施预以处理； （3）对模型结果进行经济解释。3.中国19802007年全社会固定资产投资总额X与工业总产值Y的统计资料如下表所示。单位：亿元年份全社会固定资产投资（X）工业增加值（Y）年份全社会固定资产投资（X）工业增加值（Y）1980910.91996.5199417042.119480.719819612048.4199520019.324950.619821230.42162.3199622913.529447.619831430.12375.6199724941.132921.419841832.92789.0199828406.234018.419852543.23448.7199929854.735861.519863120.63967.0200032917.740033.619873791.74585.8200137213.543580.619884753.85777.2200243499.947431.319894410.46484.0200355566.654945.5199045176858.0200470477.465210.019915594.58087.1200588773.677230.819928080.110284.52006109998.291310.9199313072.314188.02007137323.9107367.2试问：(1) 当设定模型为时，是否存在序列相关性？(2) 若按照一阶自相关假设试用广义最小二乘法估计原模型。(3) 采用差分形式与作为新数据，估计模型该模型是否存在序列相关？实验五 多重共线性的检验和修正实验目的1、理解多重共线性的含义与后果、2、学会序多重共线性的修正实验内容1、例表4是19781997年我国钢材产量（万吨）、生铁产量（万吨）、发电量（亿千瓦时）、固定资产投资（亿元）、国内生产总值（亿元）、铁路运输量（万吨）的统计资料。表4 我国钢材产量及其它相关经济变量统计资料年份钢材产量Y生铁产量X1发电量X2固定资产投资X3国内生产总值X4铁路运输量X51978220834792566668.7232641101191979249736732820699.3640381118931980271638023006746.945181112791981267034173093638.2148621076731982292035513277805.952951134951983307237383514885.26593511878419843372400137701052.43717112407419853693438441071523.51896413070919864058506444951795.321020213563519874386550349732101.691196314065319884689570454522554.861492814494819894859582058482340.5216909151489199051536238621225341854815068119915638676567753139.032161815289319926697758975394473.762663815762719937716895683956811.353463416266319948428974192819355.354675916309319958980105291007010702.975847816585519969338107231081312185.796788516880319979979115111135613838.96744631697342、多重共线性的检验（1）综合统计检验法若 在OLS法下：R2与F值较大，但t检验值较小，则可能存在多重共线性。（2）简单相关系数检验利用相关系数可以分析解释变量之间的两两相关情况。在Eviews软件中可以直接计算相关系数矩阵。本例中，在Eviews软件命令窗口中键入：COR X1 X2 X3 X4 X5或在包含所有解释变量的数组窗口中点击ViewCorrelations，其结果如图1所示。由相关系数矩阵可以看出，解释变量之间的相关系数均为0.93以上，即解释变量之间是高度相关的。图5-1（3）判定系数检验法当解释变量多余两个且变量之间呈现出较复杂的相关关系时，可以通过建立辅助回归模型来检验多重共线性。本例中，在Eviews软件命令窗口中键入：LS X1 C X2 X3 X4 X5LS X2 C X1 X3 X4 X5LS X3 C X1 X2 X4 X5LS X4 C X1 X2 X3 X5LS X5 C X1 X2 X3 X4得到相应的回归结果，分析每个方程对应的F值和T值，来检验这些变量间是否相关以及相关联程度。对应的回归结果如下图所示。图5-2图5-3图5-4图5-5图5-6上述每个回归方程的F检验值都非常显著，方程回归系数的T检验值表明：X1与X5、X2与X3、X3与X5、X4与X、X5与X1、X3、X4的T检验值较小，这些变量之间可能不相关或相关程度较小。3、多重共线性的克服逐步回归（一）建立基本的一元线性回归方程（1）被解释变量对每一个解释变量进行初始回归，选取拟合优度最高的首先进入方程；根据经济理论分析和回归结果，可知钢材产量和生铁产量关联度最大，所以建立基本的一元回归方程：（2）然后把其余解释变量逐步引入模型，根据拟合优度选出最优方程。表5 钢材产量预测模型逐步回归结果模型X1X2X3X4X5Y=f(X1)0.9214 (56.807)0.9949 0.9941 Y=f(X1,X2)0.4159 (3.5394)0.4872 (4.3234)0.9974 0.9970 Y=f(X1,X3)0.959 (14.185)0.0249 (-0.5738)0.9950 0.9940 Y=f(X1,X4)0.9414 (13.025)-0.0025 (-0.2846)0.9945 0.9938 Y=f(X1,X5)0.8578 (20.229)0.0084 (-0.2846)0.9919 0.9910 Y=f(X1,X2,X3)0.405 (2.835)0.491 (4.1225)0.0046 (0.1424)0.99690.9974Y=f(X1,X2,X4)0.4433 (3.4857)0.4911 (4.2748)-0.0039 (-0.6347)0.99690.9974Y=f(X1,X2,X5)0.4073 (3.1797)0.5025 (3.6357)-0.001 (-0.2041)0.99690.9974所以，建立的多元回归模型为：Y = -287.68669 + 0.4159*X1 + 0.4872*X2习题1.下表给出了中国商品进口额Y、国内生产总值GDP、消费者价格指数CPI。年份商品进口额(亿元)国内生产总值(亿元)居民消费价格指数(1985=100)19851257.88964.410019861498.310202.2106.519871614.211962.5114.319882055.114928.3135.819892199.916909.2160.219902574.318547.9165.219913398.721617.8170.819924443.326638.1181.719935986.234634.4208.419949960.146759.4258.6199511048.158478.1302.8199611557.467884.6327.9199711806.574462.6337.1199811626.178345.2334.4199913736.482067.5329.7200018638.889468.1331.0200120159.297314.8333.3200224430.3105172.3330.6200334195.6117251.9334.6资料来源：中国统计年鉴，中国统计出版社2000年、2004年。请考虑下列模型：(1)利用表中数据估计此模型的参数。(2)你认为数据中有多重共线性吗？ (3)进行以下回归：根据这些回归你能对数据中多重共线性的性质说些什么？(4)假设数据有多重共线性，但在5%水平上个别地显著，并且总的F检验也是显著的。对这样的情形，我们是否应考虑共线性的问题？2. 理论上认为影响能源消费需求总量的因素主要有经济发展水平、收入水平、产业发展、人民生活水平提高、能源转换技术等因素。为此，收集了中国能源消费总量Y (万吨标准煤)、国内生产总值(亿元)X1(代表经济发展水平)、国民总收入(亿元)X2(代表收入水平)、工业增加值(亿元)X3、建筑业增加值(亿元)X4、交通运输邮电业增加值(亿元)X5(代表产业发展水平及产业结构)、人均生活电力消费 (千瓦小时)X6(代表人民生活水平提高)、能源加工转换效率(%)X7(代表能源转换技术)等在1985-2002年期间的统计数据，具体如下：年份能源消费国民总收入工业建筑业交通运输邮电人均生活电力消费能源加工转换效率yX1X2X3X4X5X6X71985766828989.18964.43448.7417.9406.921.368.2919868085010201.410202.23967.0525.7475.623.268.3219878663211954.511962.54585.8665.8544.926.467.4819889299714922.314928.35777.2810.0661.031.266.5419899693416917.816909.26484.0794.0786.035.366.5119909870318598.418547.96858.0859.41147.542.467.2199110378321662.521617.88087.11015.11409.746.965.9199210917026651.926638.110284.51415.01681.854.666199311599334560.534634.414143.82284.721