计量经济学实验3多元回归模型

1、精品文档精品文档目录1建立多元线性回归模型3（一）建立包括时间变量的三元线性回归模型；31. 建立工作文件：2. 输入统计资料：3. 生成时间变量t ：4. 建立回归模型：3333CREATEA7894DATAYLKGENRT=TREND(77)LSYCTLK（二）建立剔除时间变量的二元线性回归模型；4（三）建立非线性回归模型C-D生产函数。5比较、选择最佳模型8（一）回归系数的符号及数值是否合理；8（二）模型的更改是否提高了拟合优度；8（三）模型中各个解释变量是否显著；8（四）残差分布情况8实验三多元回归模型【实验目的】掌握建立多元回归模型和比较、筛选模型的方法。【实验内容】建立我国国有独立

2、核算工业企业生产函数。根据生产函数理论，生产函数的基本形式为：Y=f（t,L,K，n。其中，L、K分别为生产过程中投入的劳动与资金，时间变量t反映技术进步的影响。表3-1列出了我国1978-1994年期间国有独立核算工业企业的有关统计资料；其中产出Y为工业总产值（可比价），L、K分别为年末职工人数和固定资产净值（可比价）。表3-1我国国有独立核算工业企业统计资料年份时间t工业总产值Y（亿元）职工人数L（万人）固定资产K（亿元）197813289.1831392225.70197923581.2632082376.34198033782.1733342522.81198143877.863488

3、2700.90198254151.2535822902.19198364541.0536323141.76198474946.1136693350.95198585586.1438153835.79198695931.3639554302.251987106601.6040864786.051988117434.0642295251.901989127721.0142735808.711990137949.5543646365.791991148634.8044727071.351992159705.5245217757.2519931610261.6544988628.77199417109

4、28.6645459374.34资料来源：根据中国统计年鉴1995和中国工业经济年鉴-1995»计算整理【实验步骤】建立多元线性回归模型（一）建立包括时间变量的三元线性回归模型;在命令窗口依次键入以下命令即可:1 .建立工作文件:CREATE A 78942 .输入统计资料:DATA Y L K3 .生成时间变量t： GENRT=TREND(77)4 .建立回归模型:LS Y C T L K则生产函数的估计结果及有关信息如图3-1所示图3-1我国国有独立核算工业企业生产函数的估计结果因此，我国国有独立工业企业的生产函数为：y? - -675.32 77.6789t 0.6667 L

5、0.7764 K（模型1）t=(-0.252)(0.672)(0.781)(7.433)R2=0.9958R2=0.9948F=1018.551模型的计算结果表明，我国国有独立核算工业企业的劳动力边际产出为0.6667,资金的边际产出为0.7764,技术进步的影响使工业总产值平均每年递增77.68亿元。回归系数的符号和数值是较为合理的。R2=0.9958,说明模型有很高的拟合优度，F检验也是高度显著的，说明职工人数L、资金K和时间变量t对工业总产值的总影响是显著的。从图3-1看出，解释变量资金K的t统计量值为7.433,表明资金对企业产出的影响是显著的。但是，模型中其他变量(包括常数项)的t统

6、计量值都较小，未通过检验。因此，需要对以上三元线性回归模型做适当的调整，按照统计检验程序，一般应先剔除t统计量最小的变量(即时间变量)而重新建立模型。(二)建立剔除时间变量的二元线性回归模型;命令：LSYCLK则生产函数的估计结果及有关信息如信I3-2所小。VariableCoefficientStd.Errort-StatisticProbC-2387.269816.8895-2.9223900.0111L1.2095320.2730204.42G5280.0006K08344960.05742114.532070.0000R-squiared0.995617Meandependentvar

7、6407.249AdjustedR*squared0.994990S,D.dependerrtvar2466742S.E.ofregression176.0069Akaikeinfocriterion13,33771Sumsquaredresid433697.8Schwar；criterion1348475Loglikelihood-110.3705F-statistic1509.953Durbin-Watsonstat1.481994Prob(F-statistic)0.000000图3-2剔除时间变量后的估计结果因此，我国国有独立工业企业的生产函数为：夕=2387.27+1.2085L+0

8、.8345K(模型2)t=(-2.922)(4.427)(14.533)R2=0.9956R2=0.9950F=1589.953从图3-2的结果看出，回归系数的符号和数值也是合理的。劳动力边际产出为1.2085,资金的边际产出为0.8345,表明这段时期劳动力投入的增加对我国国有独立核算工业企业的产出的影响最为明显。模型2的拟合优度(99.56%)较模型1(99.58%)并无多大变化，F检验也是高度显著的(F0.05(2,14)=3.74<1589.953)。这里，解释变量、常数项的t(t0.025(14)=2.145)检验值都比较大，显著性概率都小于0.05,因此模型2较模型1更为合理

9、。(三)建立非线性回归模型一一C-D生产函数C-D生产函数为：Y=ALaKpes,对于此类非线性函数，可以采用以下两种方式建立模型。方式1:转化成线性模型进行估计；在模型两端同时取对数，得：lny=lnA二二InL：InK；在EViews软件的命令窗口中依次键入以下命令：GENRLNY=log(Y)GENRLNL=log(L)GENRLNK=log(K)LSLNYCLNLLNK则估计结果如图3-3所示。VariableCoefficientStSErrorbStatiisticProbC-1.9512531.665320-1.171698D.26D9LNL0.6044670.2726972.2

10、166250.0437LNk06736580.0723579.31013100000R-squared0.995753Meandependent恻8E92837AdjustedR-squared0.995147S.D.dependentvar0.394921S.E.ofregression0.027512Akaikeinfccriterion-4189602Sumsquaredresid0.010597Schwarzcriterion-4.042564Loglikelihood38,61162F-startistic1641.407Durbin-Watsorstat1.338201Pirob(

11、F-statistic)O.ODDDOO图3-3线性变换后的C-D生产函数估计结果即可得到C-D生产函数的估计式为：ln?=1.9513+0.60451nL+0.67371nK(模型3)t=(-1.172)(2.217)(9.310)R2=0.9958R2=0.9951F=1641.407即：?=0.1424L0.6045K0.6737从模型3中看出，资本与劳动的产出弹性都是在0到1之间，模型的经济意义合理，而且(99.58%)拟合优度较模型2(99.56%)还略有提高，解释变量都通过了显著性检验。方式2:迭代估计非线性模型，迭代过程中可以作如下控制:在工作文件窗口中双击序列C,输入参数的初始

12、值；在方程描述框中点击Options,输入精度控制值控制过程：参数初值：0,0,0;迭代精度：IO：则生产函数的估计结果如图3-4所示。Convergenceachievedafter8iterations¥=C(1)*L?C为电C0)CoefficientStdErrort-StatisticPrab.0(1)4721.97115109.750.31251107693C-1.0161470.5023022.0229800.06260(3)10316880.1193078.6473610.0000R-squared0984002Meandependentvair6407.249Adj

13、ustedR-squared0.981717S.D.dependentwar2486742S,E.ofregression336.2467Akaikeinfocriteinon14,63235Sumsquaredresid1582866.Schwarzcriterion14,77939Loglikelihood-121.3760Durbin-Watsonstat0.652402图3-4生产函数估计结果此时，函数表达式为：力4721.97L,.°1161K1.0317(模型4)t=(0.313)(2.023)(8.647)_22一一R=0.9840R=0.9817可以看出，模型4中劳动

14、力弹性«=-1.01161,资金的产出弹性P=1.0317,很显然模型的经济意义不合理，因此，该模型不能用来描述经济变量间的关系。而且模型的拟合优度也有所下降，解释变量L的显著性检验也未通过，所以应舍弃该模型。参数初值：0,0,0;迭代精度：105;Convergencenotachievedafter100iterationsY=C(1*LACp)*KACp)CoefficientStdErrort-StatisticPrab.C(1)2195,28066742250.32B9190.7471c-0.8959000476646-1,8795920.0811咐1.0043770.11

15、29228.S944200.0000R-squared0985678Meandependent向6407.249AdjustedR-squared0983632S.D.dependentvar2486.742S.E.ofregression318.1453Akaikeinfocriterion14.52168Sumsquaredresid1417030Schwarzcriterion1466672Loglikelihood-120.4343Durtin-Watsonstat0.687845图3-5生产函数估计结果从图3-5看出，将收敛的误差精度改为105后，迭代100次后仍报告不收敛，说明在使

16、用迭代估计法时参数的初始值与误差精度或迭代次数设置不当，会直接影响模型的估计结果。参数初值：0,0,0;迭代精度：105,迭代次数1000;Convergenceachievedafter953iterationsY=C仅(2)*KyCoeffiicientStd.Errort-StatisticProb.c(1)0.U5D440.2496630.5809610.5706C00.611034026959122665210.0398C(3)0.664899006340810.486100.0000R-squared0995551Meandeperdentvar6407249AdjustedR-s

17、quared0.995030S.D,dependentvar2486742S,E.ofregression175.3106Akaikeinfocriterion13,32978Sumsquaredresid430273JSchwarzcriterion13,47682Loglikelihood1103031Durbin-Watsonstat1.442927图3-6生产函数估计结果此时，迭代953次后收敛，函数表达式为：?=0.1450L°.6110K0.6649(模型5)t=(0.581)(2.267)(10.486)_22R=0.9957R=0.9950从模型5中看出，资本与劳动的

18、产出弹性都是在0到1之间，模型的经济意义合理，R2=0.9957,具有很高的拟合优度，解释变量都通过了显著性检验。将模型5与通过方式1所估计的模型3比较，可见两者是相当接近的。参数初值：1,1,1；迭代精度：105,迭代次数100;Convergenceachievedafter14iterationsY=C*L?CQWCCoefficientStd.Errort-StatisticProb0(1)01450440.2496620.5B09610.5705c0.6110340.2696912.2665210.0398C0)0.6648990.06340610486100.0000R-squar

19、ed0995651Meandependentvar6407.249AdjustedR-squared0.995030S.Ddependentvar24B6.742S.E.ofregression175.3106Akaikeinfocriterion13.32978Sumsquaredresid430273.1Schwarzcriterion13.47602Loglikelihood-110.3031Durbin-Watsonstat1.442927图3-7生产函数估计结果此时，迭代14次后收敛，估计结果与模型5相同。比较方式2的不同控制过程可见，迭代估计过程的收敛性及收敛速度与参数初始值的选取

20、密切相关。若选取的初始值与参数真值比较接近，则收敛速度快；反之，则收敛速度慢甚至发散。因此，估计模型时最好依据参数的经济意义和有关先验信息，设定好参数的初始值。二、比较、选择最佳模型估计过程中，对每个模型检验以下内容，以便选择出一个最佳模型：（一）回归系数的符号及数值是否合理；（二）模型的更改是否提高了拟合优度；（三）模型中各个解释变量是否显著；（四）残差分布情况以上比较模型的、步在步骤一中已有阐述，现分析步骤一中5个不同模型的残差分布情况。分别在模型1模型5的各方程窗口中点击View/Actual,Fitted,Residual/Actual,Fitted,ResidualTable（图3-

21、8）,可以得到各个模型相应的残差分布表（图3-9至图3-13）。可以看出，模型4的残差在前段时期内连续取负值且不断增大，在接下来的一段时期又连续取正值，说明模型设定形式不当，估计过程出现了较大的偏差。而且，模型4的表达式也说明了模型的经济意义不合理，不能用于描述我国国有工业企业的生产情况，应舍弃此模型。模型1的各期残差中大多数都落在土夕的虚线框内，且残差分别不存在明显的规律性。但是，由步骤一中的分析可知，模型1中除了解释变量K之外，其余变量均为通过变量显著性检验，因此，该模型也应舍弃。模型2、模型3、模型5都具有合理的经济意义，都通过了t检验和F检验，拟合优度非常接近，理论上讲都可以描述资本、

22、劳动的投入与产出的关系。但从图3-13看出，模型5的近期误差较大，因此也可以舍弃该模型。最后将模型2与模型3比较发现，模型3的近期预测误差略小，拟合优度比模型2略有提高，因此可以选择模型2为我国国有工业企业生产函数。Equation:EQO1Torkfile:DUOYUANHUIGUI画回工Rejree岂ntationsma|Fhex*|尸口4匚.工七|Staatx|R.fid*|EstimationOutputAIActualFittedj.Residual”Actual,Fitted#Refidu£TableCvariincftMitrisActual,FitteLResidua

23、l全即hCoe££icifentTests)Residu«lTests*ResidualGraphStandardizedResidualGraphStabilityTestsficientStdErrort-StatisticProbknR7RR9nm-nX103051T77.678931156731067153SL0.666£650.8536260.780980K077641701044597432745051360.44800.0000v图3-8回归方程的残差分析ohsActualFitted|ResidualResidualPlot1978328

24、9.183223.096609123581263463.73117.5333782.173739134304263877.864057.74-1798854161.264354.38-2031254541.054651.39-1103434946.114916.1629.95255586145467.611185335931.36600079-6942666601.606541.4360.17107434.067076133579257721.017615461056477949.558196.3436.7858634.80888382-249.0229705.52952671178.3031

25、0261.61026574.0703610928710953.624.9455'<仁A1979198019811982198319841985198619871988198919901991199219931994图3-9模型1的残差分布图3-10模型2的残差分布图3-11模型3的残差分布obsActualFittedResidualResidualPlot19783269.183753.28-464.1041119793581.263927.89-346628198037S2174017.51-235339->119B13877.66411711-239.2461119824151.254315.82-164566