计量经济学实验三
计量经济学上机实验手册

实验三异方差性实验目的:在理解异方差性概念和异方差对OLS回归结果影响的基础上,掌握进行异方差检验和处理的方法;熟练掌握和运用Eviews软件的图示检验、G-Q检验、怀特White 检验等异方差检验方法和处理异方差的方法——加权最小二乘法;实验内容:书P116例4.1.4:中国农村居民人均消费函数中国农村居民民人均消费支出主要由人均纯收入来决定;农村人均纯收入除从事农业经营的收入外,还包括从事其他产业的经营性收入以及工资性收入、财产收入和转移支付收入等;为了考察从事农业经营的收入和其他收入对中国农村居民消费支出增长的影响,建立双对数模型:其中,Y表示农村家庭人均消费支出,X1表示从事农业经营的纯收入,X2表示其他来源的纯收入;表4.1.1列出了中国内地2006年各地区农村居民家庭人均纯收入及消费支出的相关数据;表4.1.1 中国2006年各地区农村居民家庭人均纯收入与消费支出单位:元注:从事农业经营的纯收入由从事第一产业的经营总收入与从事第一产业的经营支出之差计算,其他来源的纯收入由总纯收入减去从事农业经营的纯收入后得到;资料来源:中国农村住户调查年鉴2007、中国统计年鉴2007;实验步骤:一、创建文件1.建立工作文件CREATE U 1 31 其中的“U”表示非时序数据2.录入与编辑数据Data Y X1 X2 意思是:同时录入Y、X1和X2的数据3.保存文件单击主菜单栏中File→Save或Save as→输入文件名、路径→保存;二、数据分析1.散点图①Scat X1 Y从散点图可看出,农民农业经营的纯收入与农民人均消费支出呈现一定程度的正相关;②Scat X2 Y从散点图可看出,农民其他来源纯收入与农民人均消费支出呈现较高程度的正相关;2.数据取对数处理Genr LY=LOG YGenr LX1=LOG X1Genr LX2=LOG X2三、模型OLS 参数估计与统计检验 LS LY C LX1 LX2得到模型OLS 参数估计和统计检验结果:Dependent Variable: LY Method: Least Squares Sample: 1 31Variable CoefficientStd. Errort-StatisticProb.C LX1 R-squaredMean dependent var Adjusted R-squared . dependent var . of regression Akaike info criterion Sum squared resid Schwarz criterion Log likelihood F-statistic 注意:在学术文献中一般以这种形式给出回归方程的输出结果,而不是把上面的软件输出结果直接粘贴到文章中可决系数,调整可决系数,显示模型拟合程度较高;同时,F 检验统计量,在5%的显着性水平下通过方程总体显着性检验;可认为农民农业经营的收入和其他收入整体与农村居民消费支出的线性关系显着成立;变量X2和截距项均在5%的显着性水平下通过变量显着性检验,但X1在10%的显着水平下仍不能通过检验;四、异方差检验对于双对数模型,由于12(0.150214)(0.477453)ββ=<=二者均为弹性系数,可认为其他来源的纯收入而不是从事农业经营的纯收入的增长,对农户人均消费的增长更有刺激作用;也就是说,不同地区农村人均消费支出的差别主要来源于非农经营收入及工资收入、财产收入等其他来源收入的差别,因此,如果模型存在异方差性,则可能是X2引起的;1.图示检验法观察残差的平方与LX2的散点图;①残差resid残差resid变量数据是模型参数估计命令完成后由Eviews软件自动生成在Workfile 框里可找到,无需人工操作获得;注意,resid保留的是最近一次估计模型的残差数据;②残差的平方与LX2的散点图Scat LX2 resid^2从上图可大体判断出模型存在递增型异方差性;2.G-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的子样本做普通最小二乘估计,并记录残差平方和RSS;1Smpl 1 12 意思是:将样本区间由1-31,改为1-12Ls LY C LX1 LX2Dependent Variable: LYMethod: Least Squares Sample: 1 12C LX1 LX2R-squaredMean dependent var Adjusted R-squared . dependent var . of regression Akaike info criterion Sum squared resid Schwarz criterion Log likelihood F-statistic Durbin-Watson statProbF-statistic子样本1:12ln 3.1412080.398385ln 0.234751ln Y X X e =+++1RSS =④对20-31的子样本做普通最小二乘估计,并记录残差平方和2RSS ; Smpl 20 31 意思是:将样本区间由1-12,改为20-31 Ls LY C LX1 LX2Dependent Variable: LY Method: Least Squares Sample: 20 31Included observations: 12C LX1 R-squaredMean dependent var Adjusted R-squared . dependent var . of regression Akaike info criterion Sum squared resid Schwarz criterion Log likelihood F-statistic Durbin-Watson statProbF-statistic子样本2:12ln 3.9936440.113766ln 0.6201681ln Y X X e =-++2RSS =⑤异方差检验在5%与10%的显着性水平下,自由度为9,9的F分布临界值分别为0.05(9,9) 3.18F=与0.10(9,9) 2.44F=;因此5%显着性水平下不能拒绝同方差假设,但在10%的显着性水平下拒绝;补充:怀特检验软件操作:在原始模型的OLS方程对象窗口中,选择view/Residual test/White Heteroskedasticity;Eviews提供了包含交叉项的怀特检验“White Heteroskedasticitycross terms”和没有交叉项的怀特检验“White Heteroskedasticityno cross terms”这样两个选择;问题:如果是刚做完上面的G-Q检验,如何得到原始模型答案:先恢复成全样本,再按T排序,然后做OLS回归;SMPL 1 31 意思是:将样本区间恢复到1-31补充:将样本数据按T升序排列,使数据顺序还原;Sort T 意思是:将数据顺序还原Ls LY C LX1 LX2下面是在原始模型的OLS方程对象窗口中,选择view/Residual test/White Heteroskedasticity,然后进行包含交叉项的怀特检验“White Heteroskedasticitycross terms”所得到的输出结果最上方显示了两个检验统计量:F统计量和White统计量nR2;下方显示的是以OLS的残差平方为被解释变量的辅助回归方程的回归结果:F-statistic ProbabilityTest Equation:Dependent Variable: RESID^2Method: Least SquaresDate: 05/03/11 Time: 17:21Sample: 1 31C LNX1 LNX1^2 LNX1LNX2 LNX2 R-squaredMean dependent var Adjusted R-squared . dependent var . of regression Akaike info criterion Sum squared resid Schwarz criterion Log likelihood F-statistic 可见,怀特统计量nR 2==31×,大于自由度也即辅助回归方程中解释变量的个数为5的2分布临界值07.115205.0=)(χ,因此,在5%的显着性水平下拒绝同方差的原假设; 五、采用加权最小二乘法处理异方差以下内容和教材P118-120不一样,但是我们必须掌握的重点——以原始模型的OLS 回归残差的绝对值的倒数为权数,手工完成加权最小二乘估计LS LY C LX1 LX2Genr 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: LYE Method: Least Squares Sample: 1 31CELX1ER-squared Mean dependent varAdjusted R-squared . dependent var. of regression Akaike info criterionSum squared resid Schwarz criterionLog likelihood Durbin-Watson stat可以看出,lnX1参数的t统计量有了显着改进,这表明在1%显着性水平下,都不能拒绝从事农业生产带来的纯收入对农户人均消费支出有着显着影响的假设;六、检验加权的回归模型是否还存在异方差1.检验是否由LX1E引起异方差Sort LX1E 意思是:将原始数据按LX1E升序排列①子样本1的回归:Smpl 1 12LS LYE CE LX1E LX2EDependent Variable: LYEMethod: Least SquaresSample: 1 12Variable Coefficient Std. Error t-Statistic Prob.CELX1ER-squared Mean dependent varAdjusted R-squared . dependent var. of regression Akaike info criterionSum squared resid Schwarz criterionLog likelihood Durbin-Watson stat子样本1:RSS=1②子样本2的回归:Smpl 20 31LS LYE CE LX1E LX2EDependent Variable: LYE Method: Least Squares Date: 05/01/11 Time: 23:23 Sample: 20 31Variable CoefficientStd. Errort-StatisticProb.CE LX1E R-squaredMean dependent var Adjusted R-squared . dependent var . of regression Akaike info criterion Sum squared resid Schwarz criterion Log likelihoodDurbin-Watson stat子样本2:2RSS =③异方差检验 注意做题的步骤提出假设 22012:H σσ= 22112:H σσ≠ 计算检验统计量:在5%的显着性水平下,自由度为9,9的F 分布临界值分别为0.05(9,9) 3.18F =;因此5%显着性水平下不能拒绝同方差假设;2.检验是否由LX2E 引起异方差Smpl 1 31 意思是:将样本区间复原Sort lx2e 意思是:将原始数据按LX2E 升序排列 ①子样本1的回归: Smpl 1 12LS LYE CE LX1E LX2EDependent Variable: LYE Method: Least Squares Sample: 1 12CE LX1E R-squaredMean dependent var Adjusted R-squared . dependent var . of regression Akaike info criterion Sum squared resid Schwarz criterion Log likelihoodDurbin-Watson stat子样本1:1RSS = ②子样本2的回归: Smpl 20 31LS LYE CE LX1E LX2EDependent Variable: LYE Method: Least Squares Sample: 20 31Included observations: 12CE LX1E R-squaredMean dependent var Adjusted R-squared . dependent var . of regression Akaike info criterion Sum squared resid Schwarz criterion Log likelihoodDurbin-Watson stat子样本2:2RSS =③异方差检验 注意做题的步骤提出假设 22012:H σσ= 22112:H σσ≠ 计算检验统计量:在5%的显着性水平下,自由度为9,9的F 分布临界值分别为0.05(9,9) 3.18F =;因此5%显着性水平下不能拒绝同方差假设;结论:用OLS 估计的残差绝对值的倒数作为权数,对存在异方差的模型加权,然后采用OLS估计,则一定会消除异方差;最终通过异方差检验的估计方程为:实验四序列相关性实验目的:在理解序列相关性的基本概念、序列相关的严重后果的基础上,掌握进行序列相关检验和处理的方法;熟练掌握Eviews软件的图示检验、DW检验、拉格朗日乘数LM检验等序列相关性检验方法和处理序列相关性的方法——广义差分法;实验内容:书P132例4.2.1:中国居民总量消费函数建立总量消费函数是进行宏观经济管理的重要手段;为了从总体上考察中国居民收入与消费的关系,P56表2.6.3给出了中国名义支出法国内生产总值GDP、名义居民总消费CONS以及表示宏观税负的税收总额TAX、表示价格变化的居民消费价格指数CPI1990=100,并由这些数据整理出实际支出法国内生产总值GDPC=GDP/CPI、居民实际消费总支出Y=CONS/CPI,以及实际可支配收入X=GDP-TAX/CPI;表2.6.3 中国居民总量消费支出与收入资料单位:亿元年份GDP CONS CPI TAX GDPC X Y19781979198019811982198319841985198619871988198919901991199219931994199519961997199819992000200120022003200420052006资料来源:根据中国统计年鉴2001,2007整理;实验步骤:一、创建文件1.建立工作文件CREATE A 1978 2006 其中的“A”表示年度数据2.录入与编辑数据Data X Y3.保存文件单击主菜单栏中File→Save或Save as→输入文件名、路径→保存;二、数据分析:趋势图Plot X Y 意思是:同时画出Y和X的趋势图从X和Y的趋势图中可看出它们存在共同变动趋势;三、OLS参数估计与统计检验LS Y C XDependent Variable: YMethod: Least Squares Sample: 1978 2006C R-squaredMean dependent var Adjusted R-squared . dependent var . of regression Akaike info criterion Sum squared residSchwarz criterion Log likelihood F-statistic Durbin-Watson statProbF-statistic从OLS 估计的结果看,模型拟合较好:可决系数20.9880R =,截距项和斜率项的t 检验值均大于5%显着性水平下自由度为n-2=27的临界值0.025(27) 2.05t =;而且,斜率项符合经济理论中边际消费倾向在0与1之间的绝对收入假说;斜率项表明,在1978—2006年间,以1990年价计算的中国居民可支配总收入每增加1亿元,居民消费支出平均增加亿元;四、序列相关性检验 1.图示检验法①残差与时间t 的关系图趋势图 Plot resid②相邻两期残差之间的关系图 Scat resid-1 resid从两个关系图看出,随机误差项呈正序列相关性;.检验值为,表明在5%显着性水平下,n=29,k=2包括常数项,查表得1.34L d =, 1.48U d =,由于.= 1.34L d <=,故存在正序列相关;五、处理序列相关1.修正模型设定偏误剔除虚假序列相关首先面临的问题是,模型的序列相关是纯序列相关,还是由于模型设定有偏误而导致的虚假序列相关;从X 和Y 的趋势图中看到它们表现出共同的变动趋势,因此有理由怀疑较高的2R =部分地是由这一共同的变化趋势带来的;为了排除时间序列模型中这种随时间变动而具有的共同变化趋势的影响,一种解决方案是在模型中引入时间趋势项,将这种影响分离出来;由于本例中可支配收入X 与消费支出Y 均呈非线性变化态势,因此引入的时间变量TT=1,2,……,29以平方的形式出现,回归模型变化为:①编辑变量T data T在数据表中输入1-29; ②做如下的回归 Ls Y C X T^2Dependent Variable: Y Method: Least Squares Sample: 1978 2006 Included observations: 29C X T ^2R-squaredMean dependent var Adjusted R-squared . dependent var . of regression Akaike info criterionSum squared resid 6054792. Schwarz criterionLog likelihood F-statistic 得到如下的修正模型:可见,T 2的t 统计量显着;但是,修正的模型.值仍然较低,没有通过5%显着性水平下的.检验n=29,k=3时,27.1=L D ,56.1=U D ,因此该模型仍存在正序列相关性;补充:序列相关性的拉格朗日乘数检验LM检验在EViews软件中,如果在上面的OLS回归方程界面直接做残差序列的LM检验,那么得到的是如下结果,和书上P133结果不一致:原因:EViews在做LM检验时,为了不损失样本,把滞后残差序列的“前样本”缺失值设定为0Presample missing value lagged residuals set to zero.;这样,它的样本容量仍然是n,而不是n-p;回归结果和书上也有不同;解决办法:要使软件的LM检验结果和教材P133结果一致,办法是进行OLS估计之后,先把残差序列resid用genr生成另一序列e,再做辅助回归,即:genr e=resid先做含1阶滞后残差的辅助回归:ls e c x t^2 e-1Dependent Variable: EMethod: Least SquaresDate: 04/26/13 Time: 07:08Sample adjusted: 1979 2006Included observations: 28 after adjustmentsVariable Coefficient Std. Error t-Statistic Prob.CXT^2E-1R-squared Mean dependent varAdjusted R-squared . dependent var. of regression Akaike info criterionSum squared resid 2103016. Schwarz criterionLog likelihood Hannan-Quinn criter.F-statistic Durbin-Watson statProbF-statisticLM检验统计量必须自己算:LM=n-pR2=29-1=由于该值大于显着性水平为5%、自由度为1的2分布临界值84.31205.0=)(χ,由此判断原模型存在1阶序列相关;再做含2阶滞后残差的辅助回归: ls e c x t^2 e-1 e-2Dependent Variable: E Method: Least Squares Date: 04/26/13 Time: 07:32 Sample adjusted: 1980 2006Included observations: 27 after adjustmentsVariable Coefficient Std. Error t-Statistic Prob. C X T^2 E-1 E-2R-squaredMean dependent var Adjusted R-squared . dependent var . of regressionAkaike info criterion Sum squared resid 1806465. Schwarz criterion Log likelihood Hannan-Quinn criter. F-statistic Durbin-Watson statProbF-statisticLM 检验统计量必须自己算:LM=n-pR 2=29-2=由于该值大于显着性水平为5%、自由度为2的2分布临界值99.52205.0=)(χ,由此判断原模型存在序列相关;但2~-t e 的系数未通过5%的显着性检验,表明在5%的显着性水平下不存在2阶序列相关性;所以,结合前面含1阶、2阶滞后残差的辅助回归结果,可以判断在5%的显着性水平下仅存在1阶序列相关性;2.广义差分法处理序列相关①Ls Y C X T^2 AR1Dependent Variable: Y Method: Least Squares Sampleadjusted: 1979 2006Included observations: 28 after adjusting endpoints Variable CoefficientStd. Errort-StatisticProb.C X T^2 AR1R-squaredMean dependent var Adjusted R-squared . dependent var . of regression Akaike info criterionSum squared resid 2164144. Schwarz criterionLog likelihood F-statistic AR1前的参数值即为随机扰动项的1阶序列相关系数,在5%的显着性水平下显着;.= ,在5%显着性水平下,1.18.. 1.65L U d DWd =<<=样本容量为28,无法判断广义差分变换后模型是否已不存在序列相关;②继续引入AR2以下内容和教材P133-134的做法不同,但是我们必须掌握的基本做法Ls Y C X T ^2 AR1 AR2Dependent Variable: Y Method: Least Squares Sampleadjusted: 1980 2006Included observations: 27 after adjusting endpointsC X T^2 AR1 AR2R-squaredMean dependent var Adjusted R-squared. dependent var. of regression Akaike info criterionSum squared resid 1834086. Schwarz criterionLog likelihood F-statisticInverted AR Roots .53 .53+.32iAR2前的参数在10%的显着性水平下显着不为0;且.= ,接近于2,认为在10%显着性水平下,已不存在序列相关;但是,在5%的显着性水平下,则没必要引入AR2;注意:教材P133用LM检验的结果是,引入AR1 的回归方程在5%的显着性水平下已不存在序列相关性,因而不需要引入AR2;补充:下面是针对引入AR1的回归方程式的LM检验的命令操作和检验结果:首先,采用上面得到的1阶自回归系数1也即AR1的系数,做如下的1阶广义差分变量的OLS回归注:与式等价:Ls y-1 c x-1 t^t-1^2Dependent Variable: Y-1Method: Least SquaresDate: 06/02/13 Time: 11:07Sample adjusted: 1979 2006Included observations: 28 after adjustmentsVariable Coefficient Std. Error t-Statistic Prob.CX-1T^T-1^2R-squared M ean dependent varAdjusted R-squared . dependent var. of regression A kaike info criterionSum squared resid 2164144. S chwarz criterionLog likelihood H annan-Quinn criter.F-statistic D urbin-Watson statProbF-statistic然后,将上述1阶广义差分方程的残差序列resid 记为e :genr e=resid 最后,做如下的辅助回归:ls e c x-1 t^t-1^2 e-1Dependent Variable: E Method: Least Squares Date: 06/02/13 Time: 11:16 Sample adjusted: 1980 2006Included observations: 27 after adjustmentsVariable CoefficientStd. Errort-StatisticProb.C X-1 T^T-1^2 E-1R-squaredM ean dependent var Adjusted R-squared . dependent var . of regression A kaike info criterionSum squared resid 1965048. S chwarz criterionLog likelihood H annan-Quinn criter. F-statistic D urbin-Watson statProbF-statistic于是,LM 检验统计量:LM=27=;查表,当显着性水平为5%时,自由度为1的2的临界值)(1205.0χ为;上述LM <)(1205.0χ,表明模型的随机误差项已不存在序列相关;。
计量经济学上机实验

计量经济学上机实验上机实验一:一元线性回归模型实验目的:EViews软件的基本操作实验内容:对线性回归模型进行参数估计并进行检验上机步骤:中国内地2011年中国各地区城镇居民每百户计算机拥有量和人均总收入一.建立工作文件:1.在主菜单上点击File\New\Workfile;2.选择时间频率,A3.键入起始期和终止期,然后点击OK;二.输入数据:1.键入命令:DATA Y X2.输入每个变量的统计数据;3.关闭数组窗口(回答Yes);三.图形分析:1.趋势图:键入命令PLOT Y X2.相关图:键入命令 SCAT Y X 散点图:趋势图:上机结果:Yˆ11.958+0.003X=s (βˆ) 5.6228 0.0002t (βˆ) 2.1267 11.9826prob 0.0421 0.00002=0.831 R2=0.826 FR=143.584 prob(F)=0.0000上机实验二:多元线性回归模型实验目的:多元回归模型的建立、比较与筛选,掌握基本的操作要求并能根据理论对分析结果进行解释实验内容:对线性回归模型进行参数估计并进行检验上机步骤:商品的需求量与商品价格和消费者平均收入趋势图:散点图:上机结果:i Yˆ=132.5802-8.878007X1-0.038888X2s (βˆ) 57.118 4.291 0.419t (βˆ) 2.321 -2.069 -0.093prob 0.0533 0.0773 0.9286 R2=0.79 R2=0.73 F =13.14 prob(F)=0.00427三:非线性回归模型实验目的:EViews软件的基本操作实验内容:对线性回归模型进行参上机步骤:我国国有独立核算工业企业统计资料一.建立工作文件:1.在主菜单上点击File\New\Workfile;2.选择时间频率,A3.键入起始期和终止期,然后点击OK;二.输入数据:1.键入命令:DATA Y L K2.输入每个变量的统计数据;3.关闭数组窗口(回答Yes);三.图形分析:1.趋势图:键入命令PLOT Y K L2.相关图:键入命令 SCAT Y K L四.估计回归模型:键入命令LS Y C K L上机结果:Y =4047.866K1.262204L-1.227157s (βˆ) 17694.18 0232593 0.759696t (βˆ) 0.228768 5.426669 -1.615325prob 0.8242 0.0004 0.1407R2=0.989758 R2=0.987482 F=434.8689 prob(F)=0.0000上机实验四:异方差实验目的::掌握异方差的检验与调整方法的上机实现实验内容:我国制造工业利润函数行业销售销售行业销售销售实验步骤:一.检验异方差性1.图形分析检验:1) 观察Y、X相关图:SCAT Y X2) 残差分析:观察回归方程的残差图LS Y C X在方程窗口上点击Residual按钮;2. Goldfeld-Quant检验:SORT XSMPL 1 10LS Y C X(计算第一组残差平方和)SMPL 19 28LS Y C X(计算第二组残差平方和)计算F统计量,判断异方差性3.White检验:SMPL 1 28LS Y C X在方程窗口上点击:View\Residual\Test\White Heteroskedastcity 由概率值判断异方差性。
计量经济学实训报告

计量经济学实训报告一、实验设计:本次实验是基于计量经济学的理论知识和方法,通过对已有的数据进行回归分析,验证理论假设的可行性。
实验的目的是了解计量经济学在实际应用中的重要性,以及掌握回归分析等基本方法。
二、实验过程:1.数据收集:我们选择了一个包含多个变量的数据集,包括自变量和因变量,旨在通过回归模型来预测因变量的取值。
2.数据清洗:对收集到的数据进行清洗和预处理,包括处理缺失值、异常值等。
3.变量选择:根据计量经济学的原理和假设,选择适合的自变量和因变量,并对其进行初步的分析。
4.模型建立:根据选择的自变量和因变量,建立回归模型,并假设一些条件。
5.模型估计:利用统计软件对建立的回归模型进行估计和拟合,获得回归系数和拟合度等相关参数。
6.模型诊断与检验:对建立的回归模型进行诊断和检验,检查模型的拟合度和有效性。
7.结果分析:根据模型估计和检验结果,分析自变量对因变量的影响程度和显著性等,并解读模型。
三、实验结果:经过以上的实验过程和分析,我们得到了以下结论:1.自变量X对因变量Y的影响具有统计显著性;2.自变量X1对因变量Y的影响程度较大,而自变量X2的影响相对较小;3.拟合度较高,模型的解释能力较强。
四、实验感想:通过本次实验,我们深刻认识到计量经济学在实际问题中的重要性。
通过建立回归模型,我们可以对研究对象的变量关系进行实证分析,从而对问题进行解释和预测。
同时,我们也了解到了回归分析中的一些注意事项,如数据的选择和处理、模型的建立和检验等。
在今后的学习中,我们将进一步掌握和应用计量经济学的方法,提高对实际问题的分析和解决能力。
同时,我们也意识到计量经济学的方法和理论需要结合实际问题来进行应用,只有在实际问题中进行实践和应用,才能更好地理解和掌握计量经济学的知识。
计量经济实验报告多元(3篇)

第1篇一、实验目的本次实验旨在通过多元线性回归模型,分析多个自变量与因变量之间的关系,掌握多元线性回归模型的基本原理、建模方法、参数估计以及模型检验等技能,提高运用计量经济学方法解决实际问题的能力。
二、实验背景随着经济的发展和社会的进步,影响一个变量的因素越来越多。
在经济学、管理学等领域,多元线性回归模型被广泛应用于分析多个变量之间的关系。
本实验以某地区居民消费支出为例,探讨影响居民消费支出的因素。
三、实验数据本实验数据来源于某地区统计局,包括以下变量:1. 消费支出(Y):表示居民年消费支出,单位为元;2. 家庭收入(X1):表示居民家庭年收入,单位为元;3. 房产价值(X2):表示居民家庭房产价值,单位为万元;4. 教育水平(X3):表示居民受教育程度,分为小学、初中、高中、大专及以上四个等级;5. 通货膨胀率(X4):表示居民消费价格指数,单位为百分比。
四、实验步骤1. 数据预处理:对数据进行清洗、缺失值处理和异常值处理,确保数据质量。
2. 模型设定:根据理论知识和实际情况,建立多元线性回归模型:Y = β0 + β1X1 + β2X2 + β3X3 + β4X4 + ε其中,Y为因变量,X1、X2、X3、X4为自变量,β0为截距项,β1、β2、β3、β4为回归系数,ε为误差项。
3. 模型估计:利用统计软件(如SPSS、R等)对模型进行参数估计,得到回归系数的估计值。
4. 模型检验:对估计得到的模型进行检验,包括以下内容:(1)拟合优度检验:通过计算R²、F统计量等指标,判断模型的整体拟合效果;(2)t检验:对回归系数进行显著性检验,判断各变量对因变量的影响是否显著;(3)方差膨胀因子(VIF)检验:检验模型是否存在多重共线性问题。
5. 结果分析:根据模型检验结果,分析各变量对因变量的影响程度和显著性,得出结论。
五、实验结果与分析1. 拟合优度检验:根据计算结果,R²为0.812,F统计量为30.456,P值为0.000,说明模型整体拟合效果较好。
计量经济学实验报告

计量经济学实验报告:马艺菡学号:4班级:9141070302任课教师:静文实验题目简单线性回归模型分析一实验目的与要求目的:影响财政收入的因素可能有很多,比如国生产总值,经济增长,零售物价指数,居民收入,消费等。
为研究国生产总值对财政收入是否有影响,二者有何关系。
要求:为研究国生产总值变动与财政收入关系,需要做具体分析。
二实验容根据1978-1997年中国国生产总值X和财政收入Y数据,运用EV软件,做简单线性回归分析,包括模型设定,模型检验,模型检验,得出回归结果。
三实验过程:(实践过程,实践所有参数与指标,理论依据说明等)简单线性回归分析,包括模型设定,估计参数,模型检验,模型应用。
(一)模型设定为研究中国国生产总值对财政收入是否有影响,根据1978-1997年中国国生产总值X和财政收入Y,如图11978-1997年中国国生产总值和财政收入(单位:亿元)1996 66850.5 7407.991997 73452.5 8651.14根据以上数据作财政收入Y 和国生产总值X的散点图,如图2从散点图可以看出,财政收入Y和国生产总值X大体呈现为线性关系,所以建立的计量经济模型为以下线性模型:(二)估计参数1、双击“Eviews”,进入主页。
输入数据:点击主菜单中的File/Open/EV Workfile—Excel—GDP.xls;2、在EV主页界面点击“Quick”菜单,点击“Estimate Equation”,出现“Equation Specification”对话框,选择OLS估计,输入““y c x”,点击“OK”。
即出现回归结果图3;参数估计结果为:Y=857.8375+0.100036iX(67.12578)(0.002172)t=(12.77955)(46.04910)2r=0.991583F=2120.520S.E.=208.5553DW=0.864 0323、在“Equation”框中,点击“Resids”,出现回归结果的图形(图4):剩余值(Residual)、实际值(actual),拟合值(fitted)4、.(三)模型检验1.经济意义检验回归模型为:Y=857.8375+0.100036*X(其中Y为财政收入,iX为国生产总值;)所估计的参数=0.100036,说明国生产总值每增加1亿元,财政收入平均增加0.100036亿元。
EVIEWS心得

计量经济学作业(3)eviews软件学习心得姓名:林君泓班级:1008106 学号:1100800130 学院:机电工程学院(二学位)eviews软件学习心得实验中,我完成模型的参数估计,模型的统计检验,建立了一元线性回归模型和多元线性回归模型的经济计量模型,并对模型进行了异方差和自相关性检验以及对模型的修正,使得模型更加的合理。
实验过程使我对经济计量建模过程有一个直观感性的认识,并比较熟悉了现代计量经济分析软件的实际操作流程。
在整个操作过程中,我们体会和获取到用eviews软件对经济原理进行验证的乐趣与经验,通过eviews软件的应用,免去了大量的运算过程,使得我们分析问题更加的方便快捷,而且比自己计算时更加准确。
虽然在实验过程中,由于对软件不熟悉,上机操作时不可避免的遇到一些问题,但这些经验却锤炼了我发现问题的眼光,丰富了我们分析问题的思路。
而且在老师和同学的帮助下,我能够顺利的运用eviews软件对一些经济数据进行分析。
实验中,老师结合案例,现场的演示,细心的对我们进行指导,使我对eviews软件有了更深层的了解,学会了对软件进行简单的操作,对实际的经济问题进行分析与检验。
使原本枯燥、繁琐、难懂的课本知识变得简洁化,跨越理论和实践的鸿沟。
当然,在使用软件的同时虽然有时会遇到步骤和结果不同的情况,但我们可以对模型进行检验和修正,使之更能准确的分析经济问题。
通过本次实验,我也深刻体会到,eviews是一门十分实用的软件,对以后的学习有着很大的作用。
而如何正确和合理的使用便是当前最重要的任务。
实习中,我们能够直观而充分地体会到老师课堂讲授内容的精华之所在,这提高了手动操作软件、数量化分析与解决问题的能力,还可以培养我在处理实验经济问题的严谨的科学的态度,并且避免了课堂知识与实际应用的脱节。
本次实验的收获、体会、经验、问题和教训,使我初步投身于计量经济学,通过利用eviews软件将所学到的计量知识进行实践,让我加深了对理论的理解和掌握,直观而充分地体会到老师课堂讲授内容的精华之所在。
实验报告计量经济学

计量经济学实验报告书实验二、实验开设对象本实验的开设对象为《计量经济学》课程的学习者,实验为必修内容、实验目的实验二、掌握计量经济学多元模型的建立,模型形式的设定,模型拟合度、t检验和F 检验判断过程;三、实验环境微型计算机(要求必须能够连接In ternet,且安装有Eviews6.0软件。
)四、实验成果根据所给定的范例数据和要求,利用Eviews6.0软件对其进行分析和处理,并撰写实验报告。
Workflle U*mTLEDViaw | Prc-c d Oku"^1 | P-ranl N HIWH rra«x«DW BL *▼ | I Sia-r^ Tranap-aiiB E-drlI3T ■3MB ■工:xi 沁4b3-¥ XtX2IP 阴rn 丁也电niSb0.6534101985175 479724.11729 0.057131inn. IH ^I I :史Nfl 昭却* n 1*寻 1SB7 壬 g B2S£I7-2-4.13-112 ” D 日皿N 10BS2J.17 3J9 0.74200<5 I 總HP 71. 1 HURT口 TTiflHR?23:7j2S3:21.7S-103 D.7487B6-IB9-I £55 5541 2a.344*ie 0.7300821R>Ri77nn )npeii in 口 口 丁7■口sji-4 鬧 13 S1437 D.76B2&71^94 3&3 E7&& 17-93 17^ 0 61320BTRiR ■刍Hon R,Df»ri in :1:7口□ 口 riAHH433 03:2H 1占:&也斗出-IBBT 4眄 44&Z is.33333 0.9171051DEII1SiD 1 HUA ia CHI 孑pp □ 071斗口 Tis.ess«e 1.006117ZDDD &丁口 48TS 1庁方"5昌 1.069^627DD1 & 1 U 74+4 13 U7Q3Q 1 了曰□斗12002 67& 4-3^2 Ifi 12>D€2 1.^845072QQ3 T33 0&54is1.5301963DD4"iiI 葩 I Grc-up: LflN RJ I LE J D WcdJil*: (JNTTTLE&rLinfcrtiaKT'. J |optic-rii Jupdata Ad-dTri^L ・・l <oraph: UMTTT L ED Wnrkfii ■:: <jNTin"LED::Urrtrt:l«d i,i PtCTc|obj«ct j|^!Print|HMnBCarjpK Opliion-Si—Grap*! typ«-OetalwiSrapH dat-a:Fit Ihnesi!Axi^i tKJV iJdrr :1^1^ ■|~s l li^«■C^K£i[U¥|X1O[k&*朗X21333137?D146 |23 -IBD-IS ft fii 122-41^3-4 1S6.773324 OB&^D0 *£^41Qi^as175.470724 317230EW134fosei laa.teaa24 2D&&1 C €441251537 206.SJ9724 13-112 G1>QS8226.273224 1734&G.742<XM1339 231 aes?22 3G7B40 73511-321>E190237.2836S1.751D30.74^76619912S5.!ifiJ12D 3G4-SB0.73OTB21992286.390613 9DB3D0 7707171393 32i 90531E 519BT0 TBAZUj?363.27C517.BB174 O.S132tlS1995390.SO9S-IE 32DDE.0W7M11995433.932515.BZ244Q WB Mfi19S7ilGgjdiSS15.233BE0/9171Miggg50 1.385 J15.DG7B90 97H4A1199953J.9-392 1 CMW1172000 575.-337915.3E55412001 fiig n7ddldi B7-D59 1 2W4152002 570.J12215.12953 1 ML4W72003 733.CJC5d!15.424BD% si^iggzog 4in* _ b回Groupi UrrrnLED Worwila UNrTTTLEDiiUrt4iecr>. . 5 X[vfcaw] [ Ptlnt] M«n・]rriMM_| [ifWi. F J [ WDrt[Tkiam口■■[lE曰5M(I IL'L;. Grnun: UNTril l O Warlcf ik< UNTITI. ri? IJnfcrilwiA,「召斫i凶。
计量经济学课程

《计量经济学》课程实验教学手册上海立信会计学院一、实验教学大纲注:类型指:操作性、设计性、综合性。
实验形式:指实验室集中、课外分散等实验形式。
三、实验项目实验一 多元线性回归[实验目的] 使学生掌握数据的输入,作散点图,计算有关统计量(均值、方差、标准差、相关系数矩阵及方差协方差矩阵等)、掌握线性回归分析方法。
[实验内容] ①数据输入、输出及处理;②作散点图;③建立线性回归模型,对结果加以分析。
[实验步骤] 数据输入处理—作散点图—建立模型—假设检验—分析预测或点击File New Program ,打开程序输入窗口,同时执行多条命令 create a 1980 1998read D:\tao\PPT\Econometrics\zdata\P42.xls 2scat x yequation eq1.ls y c xshow eq1pagestruct(end=@last+2) *scalar n=eq1.@regobsscalar k=eq1.@ncoefx(n+1)=1763x(n+2)=1863eq1.forecast yf y_seplot y yfgenr ypl=yf(n)-@qtdist(0.975,n-k)*y_se(n)genr ypu=yf(n)+@qtdist(0.975,n-k)*y_se(n)genr ycu=yf(n)+@qtdist(0.975,n-k)*@sqrt(y_se(n)^2-@se^2)genr ycl=yf(n)-@qtdist(0.975,n-k)*@sqrt(y_se(n)^2-@se^2)show @coefs(2)show @stderr(2)genr beta1l=@coefs(2)-@qtdist(0.975,n-k)*@stderr(2)genr beta1u=@coefs(2)+@qtdist(0.975,n-k)*@stderr(2)genr yf1=yf(n,n+1)group P42 yf1 ypl ypu ycl ycushow P42[实验方法]上机[实验条件] 利用统计计量软件EViews[实验指导]1.打开EViews软件,双击桌面EViews快捷方式图标。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
实 验 三:多元回归模型与非线性回归模型【实验目的】掌握多元回归模型参数估计,特别是非线性回归模型的转化、参数估计及检验方法。
【实验内容】一、多元回归模型参数估计;二、生成序列以及可线性化模型的参数估计;三、不可线性化模型的迭代估计法的Eviews 软件的实现方式。
【实验数据】建立我国国有独立核算工业企业生产函数。
根据生产函数理论,生产函数的基本形式为:()ε,,,K L t f Y =。
其中,L 、K 分别为生产过程中投入的劳动与资金,时间变量t 反映技术进步的影响。
表3-1列出了我国1978-1994年期间国有独立核算工业企业的有关统计资料;其中产出Y 为工业总产值(可比价),L 、K 分别为年末职工人数和固定资产净值(可比价)。
资料来源:根据《中国统计年鉴-1995》和《中国工业经济年鉴-1995》计算整理【实验步骤】Y=AK一、建立多元线性回归模型㈠建立包括时间变量的三元线性回归模型; μββββ++++=L K T Y 3210在命令窗口依次键入以下命令即可: ⒈建立工作文件: CREATE A 78 94 ⒉输入统计资料: DATA Y L K ⒊生成时间变量t : GENR T=@TREND(77) ⒋建立回归模型: LS Y C T L K 则生产函数的估计结果及有关信息如图3-1所示。
图3-1 我国国有独立核算工业企业生产函数的估计结果 因此,我国国有独立工业企业的生产函数为:K L t y 7764.06667.06789.7732.675ˆ+++-= (模型1)t =(-0.252) (0.672) (0.781) (7.433)9958.02=R 9948.02=R 551.1018=F 模型的计算结果表明,我国国有独立核算工业企业的劳动力边际产出为0.6667,资金的边际产出为0.7764,技术进步的影响使工业总产值平均每年递增77.68亿元。
回归系数的符号和数值是较为合理的。
9958.02=R ,说明模型有很高的拟合优度,F 检验也是高度显著的,说明职工人数L 、资金K 和时间变量t 对工业总产值的总影响是显著的。
从图3-1看出,解释变量资金K 的t 统计量值为7.433,表明资金对企业产出的影响是显著的。
但是,模型中其他变量(包括常数项)的t 统计量值都较小,未通过检验。
因此,需要对以上三元线性回归模型做适当的调整,按照统计检验程序,一般应先剔除t 统计量最小的变量(即时间变量)而重新建立模型。
㈡建立剔除时间变量的二元线性回归模型; 命令:LS Y C L K则生产函数的估计结果及有关信息如图3-2所示。
图3-2 剔除时间变量后的估计结果因此,我国国有独立工业企业的生产函数为:K L y8345.02085.127.2387ˆ++-= (模型2) t =(-2.922) (4.427) (14.533)9956.02=R 9950.02=R 953.1589=F 从图3-2的结果看出,回归系数的符号和数值也是合理的。
劳动力边际产出为1.2085,资金的边际产出为0.8345,表明这段时期劳动力投入的增加对我国国有独立核算工业企业的产出的影响最为明显。
模型2的拟合优度较模型1并无多大变化,F 检验也是高度显著的。
这里,解释变量、常数项的t 检验值都比较大,显著性概率都小于0.05,因此模型2较模型1更为合理。
㈢建立非线性回归模型——C-D 生产函数。
C-D 生产函数为:εβαe K AL Y =,对于此类非线性函数,可以采用以下两种方式建立模型。
方式1:转化成线性模型进行估计; 在模型两端同时取对数,得:εβα+++=K L A Y ln ln ln ln在EViews 软件的命令窗口中依次键入以下命令:GENR LNY=log (Y ) GENR LNL=log (L ) GENR LNK=log (K ) LS LNY C LNL LNK 则估计结果如图3-3所示。
图3-3 线性变换后的C-D 生产函数估计结果即可得到C-D 生产函数的估计式为:K L yln 6737.0ln 6045.09513.1ˆln ++-= (模型3) t = (-1.172) (2.217) (9.310)9958.02=R 9951.02=R 407.1641=F 即:6737.06045.01424.0ˆK L y= 从模型3中看出,资本与劳动的产出弹性都是在0到1之间,模型的经济意义合理,而且拟合优度较模型2还略有提高,解释变量都通过了显著性检验。
方式2:迭代估计非线性模型,迭代过程中可以作如下控制:对不可线性化的非线性回归模型,适用于使用迭代估计法直接进行估计,例如,对于非线性CD 生产函数回归模型:βαK AL Y =,估计命令为:【命令方式】格式: NLS 变量=表达式NLS Y =C(1)*L^C(2)*K^C(3)其中,C (1),C (2),C (3)分别表示待估计的回归系数; 【菜单方式】在主菜单上点击Quick\ Estimate Equation ,并在弹出的方程描述对话框中直接输入方程的数学形式:Y =C(1)*L^C(2)*K^C(3)点击OK,系统将采用迭代估计法求解参数估计值。
说明:1.利用迭代法估计非线性回归模型时,需要事先确定待估参数的初始值。
可以用PARAM 命令直接设定,也可以在工作文件窗口中双击序列C ,然后在序列窗口中依次输入C(1)、C(2)、C(3)等参数的初始值。
2.迭代估计是一种近似估计,并且估计结果还会受到参数初始值和误差精度设定的影响。
因此,对于可线性化的非线性模型,最好还是将其转化成线性模型进行估计。
3.在方程描述窗口中是点击按纽Options ,可以设置迭代估计的最大迭代次数(Max Iterations )和误差精度(Convergence ),以便控制迭代估计的收敛过程(如图所示)。
操作:⑴在工作文件窗口中双击序列C ,输入参数的初始值; ⑵在方程描述框中点击Options ,输入精度控制值。
控制过程:①参数初值:0,0,0;迭代精度:10-3; 则生产函数的估计结果如图3-4所示。
图3-4 生产函数估计结果此时,函数表达式为:0317.101161.197.4721ˆK L y -= (模型4)t =(0.313)(-2.023)(8.647)9840.02=R 9817.02=R可以看出,模型4中劳动力弹性α=-1.01161,资金的产出弹性β=1.0317,很显然模型的经济意义不合理,因此,该模型不能用来描述经济变量间的关系。
而且模型的拟合优度也有所下降,解释变量L 的显著性检验也未通过,所以应舍弃该模型。
②参数初值:0,0,0;迭代精度:10-5;图3-5 生产函数估计结果从图3-5看出,将收敛的误差精度改为10-5后,迭代100次后仍报告不收敛,说明在使用迭代估计法时参数的初始值与误差精度或迭代次数设置不当,会直接影响模型的估计结果。
③参数初值:0,0,0;迭代精度:10-5,迭代次数1000;图3-6 生产函数估计结果此时,迭代953次后收敛,函数表达式为:6649.06110.01450.0ˆK L y = (模型5)t =(0.581)(2.267)(10.486)9957.02=R 9950.02=R从模型5中看出,资本与劳动的产出弹性都是在0到1之间,模型的经济意义合理,9957.02=R ,具有很高的拟合优度,解释变量都通过了显著性检验。
将模型5与通过方式1所估计的模型3比较,可见两者是相当接近的。
④参数初值:1,1,1;迭代精度:10-5,迭代次数100;图3-7 生产函数估计结果此时,迭代14次后收敛,估计结果与模型5相同。
比较方式2的不同控制过程可见,迭代估计过程的收敛性及收敛速度与参数初始值的选取密切相关。
若选取的初始值与参数真值比较接近,则收敛速度快;反之,则收敛速度慢甚至发散。
因此,估计模型时最好依据参数的经济意义和有关先验信息,设定好参数的初始值。
二、比较、选择最佳模型估计过程中,对每个模型检验以下内容,以便选择出一个最佳模型:㈠回归系数的符号及数值是否合理;㈡模型的更改是否提高了拟合优度;㈢模型中各个解释变量是否显著;㈣残差分布情况以上比较模型的㈠、㈡、㈢步在步骤一中已有阐述,现分析步骤一中5个不同模型的残差分布情况。
分别在模型1~模型5的各方程窗口中点击View/Actual, Fitted, Residual/ Actual, Fitted, Residual Table(图3-8),可以得到各个模型相应的残差分布表(图3-9至图3-13)。
可以看出,模型4的残差在前段时期内连续取负值且不断增大,在接下来的一段时期又连续取正值,说明模型设定形式不当,估计过程出现了较大的偏差。
而且,模型4的表达式也说明了模型的经济意义不合理,不能用于描述我国国有工业企业的生产情况,应舍弃此模型。
模型1的各期残差中大多数都落在σˆ±的虚线框内,且残差分别不存在明显的规律性。
但是,由步骤一中的分析可知,模型1中除了解释变量K之外,其余变量均为通过变量显著性检验,因此,该模型也应舍弃。
模型2、模型3、模型5都具有合理的经济意义,都通过了t检验和F检验,拟合优度非常接近,理论上讲都可以描述资本、劳动的投入与产出的关系。
但从图3-13看出,模型5的近期误差较大,因此也可以舍弃该模型。
最后将模型2与模型3比较发现,模型3的近期预测误差略小,拟合优度比模型2略有提高,因此可以选择模型2为我国国有工业企业生产函数。
图3-8 回归方程的残差分析图3-9 模型1的残差分布图3-10 模型2的残差分布图3-11 模型3的残差分布图3-12 模型4的残差分布图3-13 模型5的残差分布总结:以我国生产函数为例一、建立多元线性回归模型。
μβββ+++=L K Y 210 1.建立工作文件: CREATE A 1980 2007 2.输入统计资料: DATA GDP K L 3.建立回归模型: LS GDP C K L 二、经过直接代换可线性化的非线性回归模型的参数估计 1.命令方式: 分别建立以下模型:双对数模型:LS log(GDP) C log(K)对数模型: LS GDP C log(K) 指数模型: LS log(GDP) C K 二次多项式模型:LS GDP C K K ∧2 2.菜单方式:在方程窗口点击Estimate 按钮,并在弹出的方程描述框内依次输入上述模型。
3.迭代估计:在方程设定框内依次输入: GDP=C (1)*K^C(2)估计幂函数模型 GDP=C (1)*EXP (C (2)*K )估计指数函数模型GDP=C (1)+C (2)*K+C (3)*K^2 估计二次多项式模型 估计过程中,可以点击Option 按钮设定迭代次数和估计精度。