计量实验报告研究人均粮食产量的影响因素.docx

1、计量实验报告研究人均粮食产量的影响因素计 量 经 济 学 实 习 报 告班级： 学号： ： 【摘 要】 本报告通过统计分析1983年至2000年我国粮食的相关数据，研究人均粮食产量的影响因素，选取亩均施肥量 、人均播种面积 、人均受灾面积 、亩均机械动力 四个因素为解释变量，利用利用Eviews软件，建立回归模型进行回归分析、参数检验和模型修正从而得出最终模型。由模型可知：人均粮食产量与亩均施肥量、人均播种面积呈正相关关系，与人均受灾面积呈负相关关系。【关键字】 粮食产量 多元回归分析 检验和修正1.文献综述中国是世界上最大的粮食生产国之一，同时也是粮食的消费大国。一直以来各种农业科技迅速的发

2、展，带动了我国经济社会的发展。随着人口的增长和贸易全球化的进程不断深化，粮食安全问题已渐渐为人们所关注。关于影响粮食产量的因素，很多前人对其做过了分析，现有文献中也出现了许许多多的粮食生产模型。如通过粮食总产量、粮食播种面积、化肥费用、其他物质费用、粮食成灾面积、时间虚变量，建立柯布道格拉斯生产函数。例如张素文，李晓青等主要运用多元回归模型的方法分析了湖南省近50年来粮食播种面积，粮食产量的总体变化趋势1,王伏虎2分别从时间空间角度，粮食价格角度，粮食资源属性和资源供给结构等方面建立了粮食供需平衡函数，并确立了粮食安全模式。总结下来，影响粮食产量的主要因素有：粮食播种面积、有效灌溉面积、农业机

3、械总动力、粮食作物受灾面积、农用化肥施用量、粮食单产、种粮劳动力数量等。经研究分析，其中一些因素已被认为对粮食产量影响不显著，各因素之间也还存在着相关性。现有文献在某些变量上也达成了一致，如种植面积、施肥量等对粮食产量的影响，但某些因素的影响仍然寻在分歧。2.前期准备首先从众多的影响因素之中，选择出对因变量影响最大的四个解释变量：亩均施肥量、人均播种面积、人均受灾面积、亩均机械动力；然后通过计量经济学对模型进行多元回归分析、多重共线性检验和修正、异方差性检验和修正、自相关性检验和修正，从而得出一个拟合程度较优、估计参数显著的最终模型。为了考察这些因素对人均粮食产量的影响，构造如下模型：其中，表

4、示人均粮食产量，表示亩均施肥量，表示人均播种面积，表示人均受灾面积，表示亩均机械动力。下表列出了从1983年到2000年18年期间人均粮食产量和亩均施肥量、人均播种面积、人均受灾面积、亩均机械动力的统计数据：年份人均粮食产量亩均施肥量人均播种面积人均受灾面积亩均机械动力19831.2238229620.014553653.6039386820.5122214810.56950365119841.285497870.015412283.5626952820.4817421490.61533848819851.249065120.016314943.5861489550.7480783490.68

5、902690119861.2850297040.01740333.6410870780.7764466470.75327403419871.3024943310.017968333.6044055720.6605992870.80453514719881.2528095070.019446443.500891730.7612197470.84483893219891.2562999950.021007093.4587937920.7536474470.86518395219901.3388378180.022828873.4042795770.5346260470.86131579619911

6、.2732878380.024984423.2853511490.813600770.85967185719921.3004671390.026503263.2482298680.7607809150.89044275319931.3725637590.028521663.32275950.6955577870.95666632619941.3615659690.030288293.3509634360.9600095441.0340070319951.4431025680.032652193.4037947080.6886452551.11701124219961.5639607690.03

7、4011273.4887354160.6581753481.19487049119971.5235749150.035254893.4811884730.9344564031.29539477519981.5702008190.035888993.487574480.7717982981.38562636419991.5447043310.03644633.4383108790.8122010951.48870617820001.4091927740.03822873.3070508421.0480676881.6029880333.回归模型建立与检验根据表中数据，运用Eviews3.1软件建

8、立回归模型进行多元回归分析，OLS法的估计结果如下：（2.304216） （7.690759） (6.828676) (3.782882) (1.338807)=0.969793，=0.960498，D.W.= 1.708077，F=104.3393可决系数检验：此模型的可决系数和修正后的可决系数分别为=0.969793，=0.960498，表明人均粮食产量的变化中，可由各个解释变量的水平和变化解释的比重占到了96%以上，模型在整体上拟合得很好。参数t检验：由于n-k-1=18-4-1=13,所以t检验的自由度为13，从而在5%的显著性水平下t分布临界值为=2.160。以上数据显示，截距项、所

9、对应的t值分别为=-2.304216，=7.690759，=6.828676，=-3.782882，=-1.338807。通过比较可知，系数的t值绝对值小于临界值，所以该系数与0没有显著差异，其余4个t值都通过了显著性检验。F检验：模型的F值为104.3393，而临界值=3.18，模型F值远远大于临界值，说明在5%的显著性水平下，模型在总体上是高度显著的。4.多重共线性检验及修正4.1.多重共线性检验4.1.1.相关系数检验：变量的相关系数矩阵YX1X2X3X4Y 1.000000 0.8731000.172571 0.303278 0.843149X1 0.873100 1.0000000.

10、563081 0.605531 0.950784X20.1725710.563081 1.0000000.4302990.414674X3 0.303278 0.6055310.430299 1.000000 0.652859X4 0.843149 0.9507840.414674 0.652859 1.000000从上表可知，和相关系数高达0.950784，两者高度正相关。4.1.2.辅助回归判定系数检验：将亩均施肥量和亩均机械动力进行回归，OLS法的估计结果如下：（0.117237）（12.27389）=0.903989，=0.897989，D.W.= 0.218306，F=150.648

11、4。 =0.903989，辅助模型总体高度显著，前参数的t值12.27389=2.120，可认为显著不为0。以上数据说明此模型拟合程度很好，因此，亩均施肥量和亩均机械动力之间存在显著的线性关系。4.1.3.方差膨胀因子检验：方差膨胀因子VIF10，因此，模型存在较严重的多重共线性。4.2.多重共线性修正(1)运用OLS方法逐一求Y对各个解释变量的回归，结果如下：（21.72368）（7.163276）=0.762303，=0.747447，D.W.= 1.253261，F=51.31252。（2.321127）（0.700797）=0.029781，=-0.030858，D.W.= 0.385

12、200，F=0.491116。（8.292303）（1.273071）=0.091978，=0.035226，D.W.= 0.792076，F=1.620711。（18.94393）（6.272487）=0.710900，=0.692831，D.W.= 1.130649，F=39.34410。通过比较分析，人均粮食产量和亩均施肥量的线性关系较强，拟合程度较好。(2)在第一步选出的最优回归模型的基础上，分别代入、，结果如下：（1.940879）（12.21398） (5.022285)=0.911358，=0.899540，D.W.= 1.933729，F=77.11045。（18.94989）

13、（8.455421） (2.764412)=0.842529，=0.821533，D.W.= 0.757341，F=40.12776。（20.24167）（1.838573） (0.335044)=0.764068，=0.732611，D.W.= 1.274303，F=24.28887。通过比较分析，第一个模型可决系数有明显提高，且比其他模型高，各个解释变量的系数也都通过显著性检验，因此，人均粮食产量和亩均施肥量、人均播种面积、人均受灾面积的线性关系较强，拟合程度较好。在代入后，可决系数已无明显提高，且的系数为负，没有经济意义，所以将删除。在删除后，模型的统计检验均有较大改善，经过上述逐步回归

14、分析，表明和、的回归模型为较优，最终模型回归结果如下：（1.826592）（18.84028） (7.080862) (4.701493)=0.965628，=0.958262，D.W.= 1.558951，F=131.1013。处理后的模型的可决系数和修正可决系数分别为=0.969793，=0.960498，表明模型在整体上拟合得很好。临界值=2.145，通过对比，虽然截距项没有通过显著性检验，但、所对应的系数都是显著的。模型的F=131.1013=3.34，模型F值远远大于临界值，说明在5%的显著性水平下，模型在总体上是高度显著的。5.异方差检验和修正5.1.图示法残差平方散点图通过回归模型的残差平方散点图-可以判断，残差平方差异很小，没有出现明显差异趋势，因此初步认为模型不存在异方差性。5.2.White检验法建立辅助回归方程 通过多元回归分析，OLS法的回归方程如下：其中回归方程的=11.65357，由于服从自由度为9的卡方分布，查表可得，在5%的显著性水平下，卡方分布的临界值=16.92，所以则应当接受原假设；查表得=2.306，由回归数据得到的各参数t值得绝对值均小于临界值，即模型参数都不显著。因此，模型不存在异方差性。综上所述，模型不存在异方差性，因此不需要修正。6.自相关性检验