关于影响粮食产量因素的回归分析

1、.关于影响粮食产量因素的回归分析摘要：我国土地资源稀缺，人口多而粮食需求量大，因此粮食产量的稳定增长，直接影响着人民生活和社会的稳定与发展。粮食生产的不稳定性对国民经济的影响是不可忽略的，主要体现在:粮食生产不稳定会引发粮食供求关系的变动,尤其当国家粮食储备不足的时候,很容易导致粮价上涨,从而影响整个宏观经济。通过回归分析粮食产量波动的原因，并据此提出相应的对策，对保障粮食生产持续稳定发展，具有重要意义。关键词:线性回归 回归分析 粮食产量 宏观经济 稳定发展一、引言本文按照计量经济分析方法,以19932012年中国粮食产量及其重要因素的时间序列数据为样本,对影响中国粮食生产的多种因素进行了分

2、析。选用了粮食产量、受灾面积、化肥施用量、粮食作物播种面积 农机动力、农村用电量，以粮食产量作为因变量,其它5个指标作为解释变量进行回归分析。（一）建立模型通过对中国粮食生产及影响因素的初步定性分析后假设,粮食产量与其它3个指标之间存在多元线性关系, 即粮食受灾面积，化肥施用量，粮食作物播种面积，存在着线性关系,也即可以把粮食产量的线性回归模型初步设定为：,其中,y:粮食产量（chanliang）, 1 受灾面积（szmj），2化肥施用量(hfsyl)，3粮食作物播种面积(bzmj)，x4农机动力，x5农村用电量，然后利用已有的数据进行模型拟合,以便发现这些因素之间存在的数量关系。（2） 数据

3、搜集和来源根据相应年度的中国统计年鉴、中国农村统计年鉴、中国农业发展报告,选用了粮食产量、受灾面积，化肥施用量，粮食作物播种面积，农机动力，农村用电量这6个指标,把这6个指标的1993精品.2012年20年间的时间序列数据进行回归分析,来分析这些因素与粮食产量的关系。以粮食产量作为因变量,其它3个指标作为解释变量进行回归分析。按照计量经济分析方法对影响中国粮食生产的多种因素进行了分析。数据如下：粮食产量受灾面积化肥施用量粮食作物播种面积农机动力农村用电量年份(万吨)y(万公顷)x1(万吨)x2(千公顷)x3（万千瓦）x4（亿千瓦时）x5199345648.84882.93151.9110509

4、318171244.8199444510.15504.33317.9109544338031473.9199546661.84582.13593.7110060361181655.7199650453.54698.93827.9112548385471812.7199749417.15342.93980.7112912420161980.1199851229.55014.54083.7113787452082042.1199950838.64998.14124.3113161489962173.4200046217.55468.84146.4108463525742421.3200145263

5、.75221.54253.8106080551722610.8200245705.84711.94339.4103891579302993.4200343069.55450.64411.699410603873432.9200446946.93710.64636.6101606640283933.0200548402.23881.84766.2104278683984375.3200649746.14109.14927.7104957725224895.8200750150.24899.25107.8105638765895509.9200852850.53999.05239.01067928

6、21905713.2200953082.04721.45404.4108986874966104.4201054641.03742.65561.7109876927806632.4201157121.13247.05704.2110573977357139.6201258957.02496.05838.91112051025597508.5本次采用的估计模型为精品.其中是回归系数 。二、模型的参数估计与分析线性利用sas对模型进行拟合，参数估计和检验，用最小二乘法得到线性回归方程的形式如下：y=-30416-1.24359*x1+5.28523*x2+0.61924*x3-0.14305*x4

7、+0.95538*x5 其中x4为农机动力（万千瓦），但是得到的参数为负值，与其经济学意义不相符，故此推测存在严重的多重共线性。此外，该模型的f值221.48，可决系数r方为0.9875，修正可决系数为0.9831说明该模型的显著性成立，方程解释变量x1,x2,x3,x系数的t检验伴随概率小于5%，即t检验显著性成立，但是解释变量x5系数的t检验伴随概率大于5%，接受系数为0的假设.因此该方程需要调整，该模型可能存在多重共线性和序列相关、异方差等问题。精品.三、模型的多重共线性检验所谓多重共线性是指解释变量之间违背了相互独立的假设，及某两个或多个解释变量之间出现了相关性。多重共线性的检验主要应

8、用逐步回归法。首先计算出各个变量之间的简单相关系数矩阵，初步判断相关性。通过sas程序运行，得到各解释变量的相关系数表如下：1）第一步：由于统计量x2的f值最大，为24.47，说明它最显著，此外r 方 =0.5761和 c(p) = 459.3327，所以x2为最先选入的统计量。精品.2）第二步：其他变量继续进行统计分析，出最显著的统计量x3，此时r方=0.9548，比上一次更加显著，而c（p）=36.7162,和上一次比下降的非常明显，说明x3这个统计量是有效果的。精品.3）第三步：继续进行上述步骤，得出最显著的统计量x1，此时r2为0.9831，而c(p)=6.9133，下降的非常明显，说

9、明x1是有效果的。4）第四步：继续上述步骤，发现在水平为0.05下，其他各变量没有通过显著性检验，所以可以剔除x4，x5。回归方程为：y=-29857-1.29658x1+3.72461x2+0.32648x3精品.4、 异方差性检验精品.以上各统计量似乎没有那个参数的t检验是显著的，且可决系数比较小。但怀特统计量nr2=20*0.1563=3.126,该值小于5%显著性水平下自由度为9的分布的相应临界值16.92.因此，接受同方差性的检验。5、 序列相关性由上表得知，d.w.=1.708，查询d.w.分布表得知，n=20,k=4时，=1.00,=1.68。d.w.4-,则无自相关。6、 统计

10、学检验选取2012年数据来进行统计检验：y=-29857-1.29658*2496.0+3.72461*5838.9+0.32648*111205=58989.300049与实际值误差为32.300049，误差为0.0543%，与实际情况吻合的很好。精品.7、 总结和建议中国的粮食生产问题,不仅是中国经济界的重要研究课题,而且也越来越受到世界经济学家的重视。要提高粮食产量，必须积极稳妥地推进农业机械化的发展：要把主要农产品生产过程机械化和产业化经营有机结合起来；对农业机械化进行结构性调整；因地制宜，有重点的推荐地区农业机械化；大力促进农业技术进步，重视农村的基础教育；建立与农业机械化相适应的农

11、村经济体制。纵观中国农村现状，与其他产业相比，农业的发展一直比较缓慢。扩大耕作面积，提高单产，实现机械化、规模化生产是我国农业健康发展的必由之路。8、 参考文献吴玉鸣.中国粮食生产主要影响因素的多因素动态关联分析j,农业经济问题,1998(1)戚世均等.中国粮食生产潜力及未来粮食生产研究j.郑州粮食学院学报,2000(3)庞皓，计量经济学m，西南财经大学出版社，2001年8月第一版周四军，对我国粮食生产影响因素的计量分析，统计与决策m，2003.赵慧江，基于回归分析的粮食产量影响因素分析，怀化学院学报m，2009. 相关程序：参数估计：1）data grain1;input year y x1

12、 x2 x3 x4 x5 ;cards;精品.199345648.84882.93151.9110509318171244.8199444510.15504.33317.9109544338031473.9199546661.84582.13593.7110060361181655.7199650453.54698.93827.9112548385471812.7199749417.15342.93980.7112912420161980.1199851229.55014.54083.7113787452082042.1199950838.64998.14124.311316148996217

13、3.4200046217.55468.84146.4108463525742421.3200145263.75221.54253.8106080551722610.8200245705.84711.94339.4103891579302993.4200343069.55450.64411.699410603873432.9200446946.93710.64636.6101606640283933.0200548402.23881.84766.2104278683984375.3200649746.14109.14927.7104957725224895.8200750150.24899.25

14、107.8105638765895509.9200852850.53999.05239.0106792821905713.2200953082.04721.45404.4108986874966104.4201054641.03742.65561.7109876927806632.4201157121.13247.05704.2110573977357139.6201258957.02496.05838.91112051025597508.5run; proc print data=grain1; title 原始样本值观测; run; proc reg data=grain1; model

15、y=x1 x2 x3 x4 x5/dw;run;2）共线性检验：data grain1;input year y x1 x2 x3 x4 x5;cards2;199345648.84882.93151.9110509318171244.8199444510.15504.33317.9109544338031473.9199546661.84582.13593.7110060361181655.7199650453.54698.93827.9112548385471812.7199749417.15342.93980.7112912420161980.1199851229.55014.54083

16、.7113787452082042.1199950838.64998.14124.3113161489962173.4精品.200046217.55468.84146.4108463525742421.3200145263.75221.54253.8106080551722610.8200245705.84711.94339.4103891579302993.4200343069.55450.64411.699410603873432.9200446946.93710.64636.6101606640283933.0200548402.23881.84766.2104278683984375.

17、3200649746.14109.14927.7104957725224895.8200750150.24899.25107.8105638765895509.9200852850.53999.05239.0106792821905713.2200953082.04721.45404.4108986874966104.4201054641.03742.65561.7109876927806632.4201157121.13247.05704.2110573977357139.6201258957.02496.05838.91112051025597508.5run; proc print da

18、ta=grain1; title 粮食生产原始样本值观测; run; proc reg data=grain1; model y=x1 x2 x3 x4 x5/dw;run;proc corr;/*求相关系数矩阵*/var x1 x2 x3 x4 x5;run;proc reg data=grain1; model y=x1 /dw;run;proc reg data=grain1;model y=x2/dw;run;proc reg data=grain1;model y=x3/dw;run;proc reg data=grain1;model y=x4/dw;run;proc reg da

19、ta=grain1;model y=x5/dw;run;精品./*以下是逐步回归过程*/proc reg data=grain1;model y=x1 x2 x3 x4 x5/ dw selection=stepwise details=all slentry=0.05 slstay=0.05;title 粮食模型逐步回归;run;data grain;set grain1 ;yguji=y=-30416-1.24359*x1+5.28523*x2+0.61924*x3-0.14305*x4+0.95538*x5;/*加波浪线的数值是前面估计得来的参数值，需要从前面的回归结果中得到。*/wuc

20、ha=y-yguji;xdwucha=wucha/y;keep year y et wucha_1 wucha y x1 x2 x3 x4 x5guji xdwucha;et=abs(wucha); wucha_1=lag(wucha);proc print data=grain;run;/*proc reg data=grain;model wucha=wucha_1;title 自相关检验;run;proc reg data=grain;model et=x1 x2 x3 x4 x5;title 用戈里瑟法检验异方差;run;*/data grain2;merge grain1 grain

21、;/*合并数据集*/keep year y yguji wucha xdwucha;proc print data=grain2;title 粮食模型模拟结果;run;3）异方差检验：data liangshi;input y x1 x2 x3 ;cards3;45648.84882.93151.9110509精品.44510.15504.33317.910954446661.84582.13593.711006050453.54698.93827.911254849417.15342.93980.711291251229.55014.54083.711378750838.64998.1412

22、4.311316146217.55468.84146.410846345263.75221.54253.810608045705.84711.94339.410389143069.55450.64411.69941046946.93710.64636.610160648402.23881.84766.210427849746.14109.14927.710495750150.24899.25107.810563852850.53999.05239.010679253082.04721.45404.410898654641.03742.65561.710987657121.13247.05704

23、.211057358957.02496.05838.9111205;run;proc print data=nongcun;title 中国粮食产量;run;data nongcun2;set liangshi;run;proc reg data=liangshi2;model y=x1 x2 x3;title liangshi2;run;data lianghsi3;set liangshi2;e=y-(-29857-1.29658*x1+3.72461*x2+0.63248*x3);e1=abs(e);e2=e*e;run;proc gplot data=liangshi3;精品.symb

24、ol v=plus i=jion;plot e2*x2;run;data liangshi4;set liangshi3;x6=x1*x1;x7=x2*x2;x8=x3*x3;x9=x1*x2;x10=x2*x3;x11=x1*x3;proc reg data=liangshi4;model e2=x1 x2 x3 x6 x7 x8 x9 x10 x11;title 怀特检验法检验异方差;run;proc reg data=liangshi4;model e2=x1 x2 x3 x6 x7 x8 ;title 没有交叉项的怀特检验;run;4）序列相关性检验：data my;input year y x1 x2 x3 ;cards;1993 45648.84882.93151.9110509199444510.15504.33317.9109544199546661.84582.13593.7110060199650453