计量经济学 上机实验手册.doc

《计量经济学 上机实验手册.doc》由会员分享，可在线阅读，更多相关《计量经济学 上机实验手册.doc（18页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、【精品文档】如有侵权，请联系网站删除，仅供学习与交流计量经济学 上机实验手册.精品文档.实验三 异方差性实验目的：在理解异方差性概念和异方差对OLS回归结果影响的基础上，掌握进行异方差检验和处理的方法。熟练掌握和运用Eviews软件的图示检验、G-Q检验、怀特（White）检验等异方差检验方法和处理异方差的方法加权最小二乘法。实验内容：书P116例4.1.4：中国农村居民人均消费函数中国农村居民民人均消费支出主要由人均纯收入来决定。农村人均纯收入除从事农业经营的收入外，还包括从事其他产业的经营性收入以及工资性收入、财产收入和转移支付收入等。为了考察从事农业经营的收入和其他收入对中国农村居民消费

2、支出增长的影响，建立双对数模型：其中，Y表示农村家庭人均消费支出，X1表示从事农业经营的纯收入，X2表示其他来源的纯收入。表4.1.1列出了中国内地2006年各地区农村居民家庭人均纯收入及消费支出的相关数据。表4.1.1 中国2006年各地区农村居民家庭人均纯收入与消费支出（单位：元）地区人均消费支出Y从事农业经营的纯收入X1其他来源纯收入X2地区人均消费支出Y从事农业经营的纯收入X1其他来源纯收入X2 北 京 5724.5958.37317.2 湖 北 2732.51934.61484.8 天 津 3341.11738.94489.0 湖 南 3013.31342.62047.0 河 北 2

3、495.31607.12194.7 广 东 3886.01313.93765.9 山 西 2253.31188.21992.7 广 西 2413.91596.91173.6 内蒙古 2772.02560.8781.1 海 南 2232.22213.21042.3 辽 宁 3066.92026.12064.3 重 庆 2205.21234.11639.7 吉 林 2700.72623.21017.9 四 川 2395.014051597.4 黑龙江 2618.22622.9929.5 贵 州 1627.1961.41023.2 上 海 8006.05328606.7 云 南 2195.61570

4、.3680.2 江 苏 4135.21497.94315.3 西 藏 2002.21399.11035.9 浙 江 6057.21403.15931.7 陕 西 2181.01070.41189.8 安 徽 2420.91472.81496.3 甘 肃 1855.51167.9966.2 福 建 3591.41691.43143.4 青 海 2179.01274.31084.1 江 西 2676.61609.21850.3 宁 夏 2247.01535.71224.4 山 东 3143.81948.22420.1 新 疆 2032.42267.4469.9 河 南 2229.31844.614

5、16.4 注：从事农业经营的纯收入由从事第一产业的经营总收入与从事第一产业的经营支出之差计算，其他来源的纯收入由总纯收入减去从事农业经营的纯收入后得到。 资料来源：中国农村住户调查年鉴（2007）、中国统计年鉴（2007）。实验步骤：一、创建文件1.建立工作文件CREATE U 1 31 【其中的“U”表示非时序数据】2.录入与编辑数据Data Y X1 X2 【意思是：同时录入Y、X1和X2的数据】3.保存文件单击主菜单栏中FileSave或Save as输入文件名、路径保存。二、数据分析1.散点图Scat X1 Y 从散点图可看出，农民农业经营的纯收入与农民人均消费支出呈现一定程度的正相关

6、。Scat X2 Y 从散点图可看出，农民其他来源纯收入与农民人均消费支出呈现较高程度的正相关。2.数据取对数处理Genr LY=LOG(Y)Genr LX1=LOG(X1)Genr LX2=LOG(X2)三、模型OLS参数估计与统计检验LS LY C LX1 LX2得到模型OLS参数估计和统计检验结果：Dependent Variable: LYMethod: Least SquaresSample: 1 31Included observations: 31VariableCoefficientStd. Errort-StatisticProb. C3.2660681.0415913.13

7、56530.0040LX10.1502140.1085381.3839750.1773LX20.4774530.0515959.2538530.0000R-squared0.779878 Mean dependent var7.928613Adjusted R-squared0.764155 S.D. dependent var0.355750S.E. of regression0.172766 Akaike info criterion-0.581995Sum squared resid0.835744 Schwarz criterion-0.443222Log likelihood12.0

8、2092 F-statistic49.60117Durbin-Watson stat1.780981 Prob(F-statistic)0.000000【注意：在学术文献中一般以这种形式给出回归方程的输出结果，而不是把上面的软件输出结果直接粘贴到文章中】可决系数0.779878，调整可决系数0.764155，显示模型拟合程度较高;同时,F检验统计量49.60117，在5%的显著性水平下通过方程总体显著性检验。可认为农民农业经营的收入和其他收入整体与农村居民消费支出的线性关系显著成立。变量X2和截距项均在5%的显著性水平下通过变量显著性检验，但X1在10%的显著水平下仍不能通过检验。四、异方差检

9、验对于双对数模型，由于【二者均为弹性系数】，可认为其他来源的纯收入而不是从事农业经营的纯收入的增长，对农户人均消费的增长更有刺激作用。也就是说，不同地区农村人均消费支出的差别主要来源于非农经营收入及工资收入、财产收入等其他来源收入的差别，因此，如果模型存在异方差性，则可能是X2引起的。1.图示检验法观察残差的平方与LX2的散点图。残差（resid）残差（resid）变量数据是模型参数估计命令完成后由Eviews软件自动生成（在Workfile框里可找到），无需人工操作获得。注意，resid保留的是最近一次估计模型的残差数据。残差的平方与LX2的散点图Scat LX2 resid2从上图可大体判

10、断出模型存在递增型异方差性。2G-Q法检验异方差补充：先定义一个变量T，取值为1、2、31【分别代表各省市】，用于在做完G-Q检验之后，再按T排序，使数据顺序还原。Data T 【提示：输入1、2、31】将所有原始数据按照X2升序排列。Sort X2Show Y X1 X2 LY LX1 LX2显示各个变量数据的目的是查看一下，所有变量数据是否按X2升序排列好了。将31对样本数据，去掉中间的7对，形成两个容量均为12的子样本，即1-12和20-31。对1-12的子样本做普通最小二乘估计，并记录残差平方和。Smpl 1 12 【意思是：将样本区间由1-31，改为1-12】Ls LY C LX1

11、LX2Dependent Variable: LYMethod: Least SquaresSample: 1 12Included observations: 12VariableCoefficientStd. Errort-StatisticProb. C3.1412081.1223582.7987570.0208LX10.3983850.0787915.0562340.0007LX20.2347510.1097472.1390090.0611R-squared0.739693 Mean dependent var7.700532Adjusted R-squared0.681847 S.D

12、. dependent var0.156574S.E. of regression0.088316 Akaike info criterion-1.803481Sum squared resid0.070197 Schwarz criterion-1.682255Log likelihood13.82089 F-statistic12.78726Durbin-Watson stat1.298449 Prob(F-statistic)0.002343子样本1： =0.070197对20-31的子样本做普通最小二乘估计，并记录残差平方和。Smpl 20 31 【意思是：将样本区间由1-12，改为2

13、0-31】Ls LY C LX1 LX2Dependent Variable: LYMethod: Least SquaresSample: 20 31Included observations: 12VariableCoefficientStd. Errort-StatisticProb. C3.9936441.8840542.1197080.0631LX1-0.1137660.159962-0.7112050.4950LX20.6201680.1116545.5543800.0004R-squared0.876931 Mean dependent var8.239746Adjusted R

14、-squared0.849582 S.D. dependent var0.375812S.E. of regression0.145754 Akaike info criterion-0.801478Sum squared resid0.191197 Schwarz criterion-0.680251Log likelihood7.808868 F-statistic32.06485Durbin-Watson stat2.565362 Prob(F-statistic)0.000080子样本2：=0.191197异方差检验在5%与10%的显著性水平下，自由度为（9,9）的F分布临界值分别为与

15、。因此5%显著性水平下不能拒绝同方差假设，但在10%的显著性水平下拒绝。补充：怀特检验软件操作：在原始模型的OLS方程对象窗口中，选择view/Residual test/White Heteroskedasticity。Eviews提供了包含交叉项的怀特检验“White Heteroskedasticity（cross terms）”和没有交叉项的怀特检验“White Heteroskedasticity（no cross terms）”这样两个选择。【问题：如果是刚做完上面的G-Q检验，如何得到原始模型？答案：先恢复成全样本，再按T排序，然后做OLS回归。】SMPL 1 31 【意思是：将

16、样本区间恢复到1-31】补充：将样本数据按T升序排列，使数据顺序还原。Sort T 【意思是：将数据顺序还原】Ls LY C LX1 LX2下面是在原始模型的OLS方程对象窗口中，选择view/Residual test/White Heteroskedasticity，然后进行包含交叉项的怀特检验“White Heteroskedasticity（cross terms）”所得到的输出结果【最上方显示了两个检验统计量：F统计量和White统计量nR2；下方显示的是以OLS的残差平方为被解释变量的辅助回归方程的回归结果】：White Heteroskedasticity Test:F-stat

17、istic9.833740 Probability0.000027Obs*R-squared20.55085 Probability0.000985Test Equation:Dependent Variable: RESID2Method: Least SquaresDate: 05/03/11 Time: 17:21Sample: 1 31Included observations: 31VariableCoefficientStd. Errort-StatisticProb. C10.243285.4745221.8710820.0731LNX1-2.3290701.116442-2.0

18、861530.0473LNX120.1491140.0581072.5661950.0167LNX1*LNX20.0193330.0412650.4685070.6435LNX2-0.4573070.454020-1.0072380.3235LNX220.0211010.0133571.5796940.1267R-squared0.662931 Mean dependent var0.026959Adjusted R-squared0.595517 S.D. dependent var0.042129S.E. of regression0.026794 Akaike info criterio

19、n-4.229312Sum squared resid0.017948 Schwarz criterion-3.951766Log likelihood71.55434 F-statistic9.833740Durbin-Watson stat1.462377 Prob(F-statistic)0.000027可见，怀特统计量nR2=20.55085【=310.662931】，大于自由度n【也即辅助回归方程中解释变量的个数】为5的c2分布临界值，因此，在5%的显著性水平下拒绝同方差的原假设。五、采用加权最小二乘法处理异方差【以下内容和教材P118-120不一样，但是我们必须掌握的重点以原始模型

20、的OLS回归残差的绝对值的倒数为权数，手工完成加权最小二乘估计】LS LY C LX1 LX2 Genr E=resid 【意思是：记录双对数模型OLS估计的残差】用残差的绝对值的倒数对LY、LX1、LX2做加权：Genr LYE=LY/abs(E)Genr LX1E=LX1/abs(E)Genr LX2E=LX2/abs(E)Genr CE=1/abs(E) LS LYE CE LX1E LX2EDependent Variable: LYEMethod: Least SquaresSample: 1 31Included observations: 31VariableCoefficien

21、tStd. Errort-StatisticProb. CE3.3265760.17357219.165430.0000LX1E0.1509450.0248196.0818320.0000LX2E0.4678640.00978247.830950.0000R-squared0.999984 Mean dependent var188.5956Adjusted R-squared0.999983 S.D. dependent var245.8904S.E. of regression1.014851 Akaike info criterion2.959126Sum squared resid28

22、.83782 Schwarz criterion3.097899Log likelihood-42.86646 Durbin-Watson stat1.522872可以看出，lnX1参数的t统计量有了显著改进，这表明在1%显著性水平下，都不能拒绝从事农业生产带来的纯收入对农户人均消费支出有着显著影响的假设。六、检验加权的回归模型是否还存在异方差1.检验是否由LX1E引起异方差Sort LX1E 【意思是：将原始数据按LX1E升序排列】子样本1的回归：Smpl 1 12LS LYE CE LX1E LX2EDependent Variable: LYEMethod: Least SquaresS

23、ample: 1 12Included observations: 12VariableCoefficientStd. Errort-StatisticProb. CE1.8316892.0877290.8773600.4031LX1E0.3574790.2180921.6391210.1356LX2E0.4621270.1060424.3579420.0018R-squared0.989665 Mean dependent var35.95002Adjusted R-squared0.987369 S.D. dependent var9.562672S.E. of regression1.0

24、74733 Akaike info criterion3.194340Sum squared resid10.39546 Schwarz criterion3.315566Log likelihood-16.16604 Durbin-Watson stat2.429720子样本1：=10.39546子样本2的回归：Smpl 20 31LS LYE CE LX1E LX2EDependent Variable: LYEMethod: Least SquaresDate: 05/01/11 Time: 23:23Sample: 20 31Included observations: 12Varia

25、bleCoefficientStd. Errort-StatisticProb. CE3.3225350.18672317.793940.0000LX1E0.1529940.0270385.6584170.0003LX2E0.4665590.01048344.506230.0000R-squared0.999989 Mean dependent var397.9217Adjusted R-squared0.999987 S.D. dependent var292.1047S.E. of regression1.057514 Akaike info criterion3.162037Sum sq

26、uared resid10.06502 Schwarz criterion3.283264Log likelihood-15.97222 Durbin-Watson stat1.471565子样本2：=10.06502异方差检验 【注意做题的步骤】提出假设 计算检验统计量：在5%的显著性水平下，自由度为（9,9）的F分布临界值分别为。因此5%显著性水平下不能拒绝同方差假设。2.检验是否由LX2E引起异方差Smpl 1 31 【意思是：将样本区间复原】Sort lx2e 【意思是：将原始数据按LX2E升序排列】子样本1的回归：Smpl 1 12LS LYE CE LX1E LX2EDepende

27、nt Variable: LYEMethod: Least SquaresSample: 1 12Included observations: 12VariableCoefficientStd. Errort-StatisticProb. CE1.4746501.8155630.8122270.4376LX1E0.4032710.1814662.2222990.0534LX2E0.4693070.0983074.7739070.0010R-squared0.991889 Mean dependent var36.19042Adjusted R-squared0.990086 S.D. depe

28、ndent var10.00683S.E. of regression0.996372 Akaike info criterion3.042927Sum squared resid8.934823 Schwarz criterion3.164153Log likelihood-15.25756 Durbin-Watson stat2.031453子样本1：=8.934823子样本2的回归：Smpl 20 31LS LYE CE LX1E LX2EDependent Variable: LYEMethod: Least SquaresSample: 20 31Included observati

29、ons: 12VariableCoefficientStd. Errort-StatisticProb. CE3.3225350.18672317.793940.0000LX1E0.1529940.0270385.6584170.0003LX2E0.4665590.01048344.506230.0000R-squared0.999989 Mean dependent var397.9217Adjusted R-squared0.999987 S.D. dependent var292.1047S.E. of regression1.057514 Akaike info criterion3.

30、162037Sum squared resid10.06502 Schwarz criterion3.283264Log likelihood-15.97222 Durbin-Watson stat1.264061子样本2：=10.06502异方差检验 【注意做题的步骤】提出假设 计算检验统计量：在5%的显著性水平下，自由度为（9,9）的F分布临界值分别为。因此5%显著性水平下不能拒绝同方差假设。【结论：用OLS估计的残差绝对值的倒数作为权数，对存在异方差的模型加权，然后采用OLS估计，则一定会消除异方差。】最终通过异方差检验的估计方程为：实验四 序列相关性实验目的：在理解序列相关性的基本概念

31、、序列相关的严重后果的基础上，掌握进行序列相关检验和处理的方法。熟练掌握Eviews软件的图示检验、DW检验、拉格朗日乘数（LM）检验等序列相关性检验方法和处理序列相关性的方法广义差分法。实验内容：书P132例4.2.1（例2.6.2）：中国居民总量消费函数建立总量消费函数是进行宏观经济管理的重要手段。为了从总体上考察中国居民收入与消费的关系，P56表2.6.3给出了中国名义支出法国内生产总值GDP、名义居民总消费CONS以及表示宏观税负的税收总额TAX、表示价格变化的居民消费价格指数CPI（1990=100），并由这些数据整理出实际支出法国内生产总值GDPC=GDP/CPI、居民实际消费总支

32、出Y=CONS/CPI，以及实际可支配收入X=（GDP-TAX）/CPI。表2.6.3 中国居民总量消费支出与收入资料单位：亿元年份GDPCONSCPITAXGDPCXY19783605.61759.146.21519.287802.56678.83806.719794092.62011.547.07537.828694.27551.64273.219804592.92331.250.62571.709073.77944.24605.519815008.82627.951.90629.899651.88438.05063.919825590.02902.952.95700.0210557.39

33、235.25482.419836216.23231.154.00775.5911510.810074.65983.219847362.73742.055.47947.3513272.811565.06745.719859076.74687.460.652040.7914966.811601.77729.2198610508.55302.164.572090.3716273.713036.58210.9198712277.46126.169.302140.3617716.314627.78840.0198815388.67868.182.302390.4718698.715794.09560.5

34、198917311.38812.697.002727.4017847.415035.59085.5199019347.89450.9100.002821.8619347.816525.99450.9199122577.410730.6103.422990.1721830.918939.610375.8199227565.213000.1110.033296.9125053.022056.511815.3199336938.116412.1126.204255.3029269.125897.313004.7199450217.421844.2156.655126.8832056.228783.4

35、13944.2199563216.928369.7183.416038.0434467.531175.415467.9199674163.633955.9198.666909.8237331.933853.717092.5199781658.536921.5204.218234.0439988.535956.218080.6199886531.639229.3202.599262.8042713.138140.919364.1199991125.041920.4199.7210682.5845625.840277.020989.3200098749.045854.6200.5512581.51

36、49238.042964.622863.92001108972.449213.2201.9415301.3853962.546385.424370.12002120350.352571.3200.3217636.4560078.051274.026243.22003136398.856834.4202.7320017.3167282.257408.128035.02004160280.463833.5210.6324165.6876096.364623.130306.22005188692.171217.5214.4228778.5488002.174580.433214.4200622117

37、0.580120.5217.6534809.72101616.385623.136811.2资料来源：根据中国统计年鉴（2001，2007）整理。实验步骤：一、创建文件1.建立工作文件CREATE A 1978 2006 【其中的“A”表示年度数据】2.录入与编辑数据Data X Y3.保存文件单击主菜单栏中FileSave或Save as输入文件名、路径保存。二、数据分析:趋势图Plot X Y 【意思是：同时画出Y和X的趋势图】从X和Y的趋势图中可看出它们存在共同变动趋势。三、OLS参数估计与统计检验LS Y C XDependent Variable: YMethod: Least Sq

38、uaresSample: 1978 2006Included observations: 29VariableCoefficientStd. Errort-StatisticProb. C2091.295334.98696.2429140.0000X0.4375270.00929747.059500.0000R-squared0.987955 Mean dependent var14855.72Adjusted R-squared0.987509 S.D. dependent var9472.076S.E. of regression1058.633 Akaike info criterion

39、16.83382Sum squared resid30259014 Schwarz criterion16.92811Log likelihood-242.0903 F-statistic2214.596Durbin-Watson stat0.277155 Prob(F-statistic)0.000000从OLS估计的结果看，模型拟合较好：可决系数，截距项和斜率项的t检验值均大于5%显著性水平下自由度为n-2=27的临界值。而且，斜率项符合经济理论中边际消费倾向在0与1之间的绝对收入假说。斜率项0.438表明，在19782006年间，以1990年价计算的中国居民可支配总收入每增加1亿元，居民

40、消费支出平均增加0.438亿元。四、序列相关性检验1.图示检验法残差与时间t的关系图（趋势图）Plot resid相邻两期残差之间的关系图Scat resid(-1) resid从两个关系图看出，随机误差项呈正序列相关性。2.D.W.检验法D.W.检验值为0.277155，表明在5%显著性水平下，n=29，k=2（包括常数项），查表得，由于D.W.=0.277155，故存在正序列相关。五、处理序列相关1.修正模型设定偏误（剔除虚假序列相关）首先面临的问题是，模型的序列相关是纯序列相关，还是由于模型设定有偏误而导致的虚假序列相关。从X和Y的趋势图中看到它们表现出共同的变动趋势，因此有理由怀疑较高

41、的0.987955部分地是由这一共同的变化趋势带来的。为了排除时间序列模型中这种随时间变动而具有的共同变化趋势的影响，一种解决方案是在模型中引入时间趋势项，将这种影响分离出来。由于本例中可支配收入X与消费支出Y均呈非线性变化态势，因此引入的时间变量T（T=1,2，29）以平方的形式出现，回归模型变化为：编辑变量Tdata T 在数据表中输入1-29。做如下的回归Ls Y C X T2 Dependent Variable: YMethod: Least SquaresSample: 1978 2006Included observations: 29VariableCoefficientStd. Errort-StatisticProb. C3328.191195.032617.064790.0000X0.1761520.0259866.7787880.0000T221.655822.12418310.194890.0000R-squared0.997590 Mean dependent