我国粮食产量的影响因素分析计量经济学模型.doc

-我国粮食产量的影响因素分析一.研究背景：改革开放以来，中国经济迅速开展，人口增长迅猛，对粮食的需求日益增加粮食产量无疑成了影响中国经济开展的重大因素同时，粮食的产量直接关系到农业劳动力的生活水平，因此，"三农〞问题成为中国经济研究的热点问题，提高粮食产量，关注农村居民收入迫在眉睫为此，本文将就粮食产量影响因素进展分析，希望从中发现一些对粮食产量关键作用的因素二.研究方案与数据的搜集统计：影响粮食总产量的因素有很多，包括粮食作物耕种面积、粮食面积单产、有效灌溉面积、化肥用量、农药用量、农业机械总动力、农用塑料薄膜用量、受灾面积、成灾面积等，根据实际情况及模型建立需要选取其中五个作为研究对象，分别农业化肥施用量〔*1〕，粮食播种面积(*2)，成灾面积(*3)，农业机械总动力(*4)，农业劳动力(*5)。

表中列出了中国粮食生产的相关数据，拟建立中国粮食生产函数：表1 中国粮食生产与相关投入资料年份粮食产量(万吨)农业化肥施用量〔万公斤〕粮食播种面积〔千公顷〕成灾面积〔公顷〕农业机械总动力〔万千瓦〕农业劳动力〔万人〕1983 38728 1660 114047 16209 18022 31151 1984 40731 1740 112884 15264 19497 30868 1985 37911 1776 108845 22705 20913 31130 1986 39151 1931 110933 23656 22950 31254 1987 40208 1999 111268 20393 24836 31663 1988 39408 2142 110123 23945 26575 32249 1989 40755 2357 112205 24449 28067 33225 1990 44624 2590 113466 17819 28708 38914 1991 43529 2806 112314 27814 29389 39098 1992 44264 2930 110560 25895 30308 38699 1993 45649 3152 110509 23133 31817 37680 1994 44510 3318 109544 31383 33802 36628 1995 46662 3594 110060 22267 36118 35530 1996 50454 3828 112548 21233 38547 34820 1997 49417 3981 112912 30309 42021 34840 1998 51230 4084 113787 25181 45208 35177 1999 50839 4124 113161 26731 48996 35768 2000 46218 4146 108463 34374 52574 36043 2001 45264 4254 106080 31793 55172 36513 2002 45706 4339 103891 27319 57930 36870 2003 43070 4412 99410 32516 60387 36546 2004 46947 4637 101606 16297 64028 35269 2005 48402 4766 104278 19966 68398 33970 2006 49804 4928 104958 24632 72522 32561 2007 50160 5108 105638 25064 76590 31444 注：这里由于没有从事粮食生产的农业劳动数据，用第一产业劳动力替代。

资料来源："中国统计年鉴"〔1995，2021〕研究假设：农业化肥施用量〔*1〕与粮食产量正相关粮食播种面积(*2) 与粮食产量正相关成灾面积(*3) 与粮食产量负相关农业机械总动力(*4) 与粮食产量正相关农业劳动力(*5) 与粮食产量正相关三、模型的估计、检验、确认1.画散点图由于点较分散，将他们取对数，使其更集中设A1=log(1) A2=log(2) A3=log(3) A4=log(4) A5=log(5) Z=log(y), 做散点图如左侧由图可以看出，log〔y〕和log(*1),log(*4)有较为明显的线性关系，建立多元回归模型2.用OLS估计模型：LogY=b0+b1log(*1)+b2log(*2)+b3log(*3)+b4log(*4)+b5log(*5)Logy=-4.173+0.381log(*1)+1.222log(*2)-0.081log(*3)-log(*4)-0.101log(*5)从模型可以看出，*1,*2,*3均通过了显著性检验，且估计量的系数符合经济含义，*4和*5未通过显著性检验，且系数为负，不符合经济含义模型整体R^2为0.981587，F-statistic为202.6826>2.74 (0.05水平下的F统计量值)，DW=1.79 ，可以看出模型整体较优，但个别解释变量没有通过显著性检验，具有多重共线性。

由于我们更关心多重共线性的程度，所以运用KLEIN判别法：图中可以看出log(*4)和log(*1)存在高度相关性，但并没有超过R^2，不是有害的3.运用逐步回归法抑制多重共线性：用每个*对y进展简单回归，按R^2排序：①　 Log(y)=8.902+0.224log(*1)T=0.000 0.000R^2=0.1101 DW=0.939②　 Log(y)=15.1574-0.3834log(*2)T=0.0174 0.4595R^2=0.02 DW=0.33③　 Log(y)=9.619+0.108log(*3)T=0.0000 0.2177R^2=0.0652 DW=0.597④　 Log(y)=8.9490+0.16697log(*4)T=0.000 0.0000R^2=0.602 DW=0.62⑤　 Log(y)=5.6007+0488731log(*5)T=0.0319 0.0485R^2=0.158 DW=0.32排序后：R1^2=77% R4^2=60.2% R5^2=15.8% R3^2=6.5% R2^2=2.4% 由此可见，粮食生产受农业化肥施用量的影响最大，与经历相符合，选Log(y)=8.092+0.224log(*1)为初始回归模型，依次引入Log(*4) log(*5) log(*3) log(*2)进展回归，寻找最正确回归方程〔见下表〕〔Y=log(y)〕表2 逐步回归结果如表中所示：Log(*4)因为经济含义不符合，剔除Log(*5)未通过显著性检验，剔除Log(*3)未通过显著性检验，剔除Log(*2)通过显著性检验，经济意义符合，R^2=0.94,均优于前面的，AIC，SC都有所降低，故保存。

Log(y)=-6.2856+0.2978log(*1)+1.2586log(*2)R^2=0.9452 F+statistic=189.9002 DW=1.59由上表可以看出，该方程为最优模型：各解释变量均通过显著性检验，R^2 较优，F检验通过, 变量的系数均符合经济含义，已剔除多重共线性4.受线性约束回归的F检验 Wald检验：通过此方法进一步验证是否应剔除解释变量log(*3) Log(*4)和log〔*5〕:可以看出，F检验和Wald检验的伴生概率都<0.05,说明原假设不成立，约束条件b3=b4=b5=0不成立，进一步检验b3=0是否成立可以看出，F检验和Wald检验的伴生概率都<0.05,说明原假设不成立，约束条件b3=0不成立，不应剔除解释变量log(*3),在此根底上检验是否应剔除log(*4)和log(*5).可以看出，F检验和Wald检验的伴生概率都>0.05,说明原假设成立，约束条件成立，b4=b5=0,可以剔除log(*4)和log(*5).结合以上分析，在之前得出的模型中参加解释变量Log(*3)，建立多元回归模型：由模型可知，所有解释变量均通过显著性检验，系数符号符合经济含义，R^2=0.978616，高于前面得出的模型，所以该模型为最优模型。

由以下列图可知，真实值和估计值拟合的很好最优模型为：logy=-5.999638+0.323385log(*1)+1.290729log(*2)-0.086754log(*3)R^2=0.978616 DW=1.41 F=320.34385.异方差检验：怀特检验：Obs*R-squared 的伴生概率为0.2128，说明原假设成立，模型不存在异方差性6. 自相关性检验：DW=1.413, dl=1.12 du=1.66,落在不确定区域，无法判断LM法检验自相关性：据检验结果可知，obs*R-squared 的伴生概率>0.05，承受原假设，该模型不存在自相关性7、格兰杰〔Granger〕因果检验：运用该检验可以判断两个变量在时间上的先导-滞后关系，由于该因果关系与哲学意义上的因果关系还是有区别的，如果可以证明"*是y的格兰杰原因〞，只是说明"*中包含了预测y的有效信息〞检验结果见下表:表3 格兰杰因果检验结果分析F检验〔prob〕格兰杰因果检验滞后1期滞后2期滞后4期滞后6期log(*1)不是log(y)的格兰杰原因0.05150.00480.02980.1133log(*2)不是log(y)的格兰杰原因0.18220.05160.20570.1628log(*3)不是log(y)的格兰杰原因0.28340.31930.05340.1342log(*4)不是log(y)的格兰杰原因0.21750.43320.28070.0342log(*5)不是log(y)的格兰杰原因0.70230.34010.26090.0004以上结果显示：在滞后2期和滞后期为4的时候，log(*1)是log(y)的格兰杰原因。

在滞后期为6的时候，log〔4〕和log〔5〕是log〔y〕的格兰杰原因其余的不管在滞后几期都不是log(y)的格兰杰原因以上结果说明：农业化肥的施用量包含了较多的对粮食产量的有。