糧食產(chǎn)量分析

1、計(jì)量經(jīng)濟(jì)學(xué)課程論文Y E U YANG-EN UNIVERSITY計(jì)量經(jīng)濟(jì)學(xué)論文題 目：中國(guó)糧食生產(chǎn)影響因素研究姓 名：陳好妹 院 系： 財(cái)政金融 專 業(yè)： 財(cái)政學(xué)年 級(jí)： 07級(jí) 學(xué) 號(hào)： 011420070048指導(dǎo)老師： 舒輔琪副教授 完 成日期：二O一O年五月三十日目 錄摘 要 3 引 言4影響我國(guó)糧食產(chǎn)量的實(shí)證分析4結(jié)論與建議24參考文獻(xiàn) 25摘 要糧食是人類最基本的生活消費(fèi)品，一個(gè)國(guó)家的糧食問(wèn)題是關(guān)系到本國(guó)的國(guó)計(jì)民生的頭等大事。人們都知道，農(nóng)業(yè)是國(guó)民經(jīng)濟(jì)發(fā)展的基礎(chǔ)，糧食是基礎(chǔ)的基礎(chǔ)，因此糧食生產(chǎn)是關(guān)系到一個(gè)國(guó)家生產(chǎn)與發(fā)展的一個(gè)永恒的主題。建國(guó)以來(lái)我國(guó)的糧食產(chǎn)量出現(xiàn)了多次的波動(dòng)，給消

2、費(fèi)者和生產(chǎn)者帶來(lái)了很不利的影響，所以了解影響糧食生產(chǎn)因素很重要。通過(guò)計(jì)量經(jīng)濟(jì)學(xué)方法創(chuàng)建我國(guó)糧食生產(chǎn)函數(shù)，我們會(huì)發(fā)現(xiàn)糧食播種、化肥施用量、受災(zāi)面積是影響糧食生產(chǎn)的三大因素，其中糧食播種面積的影響最大。關(guān)鍵詞: 糧食產(chǎn)量、播種面積、化肥施用量、受災(zāi)面積 引言總所周知，糧食是我們?nèi)祟惿靡匝永m(xù)的最基礎(chǔ)的物質(zhì)條件，沒(méi)有糧食這個(gè)物質(zhì)基礎(chǔ)，人類將無(wú)法生存?；仡櫸覈?guó)糧食的生產(chǎn)情況，我們會(huì)發(fā)現(xiàn)，隨著社會(huì)的發(fā)展，技術(shù)水平的提高，從整體來(lái)講我過(guò)糧食產(chǎn)量呈上升的趨勢(shì)。在改革開(kāi)放（1978年）以前我國(guó)糧食產(chǎn)量緩慢增長(zhǎng)，一直都存于30000萬(wàn)噸以下。改革開(kāi)放后，我國(guó)糧食產(chǎn)量從30000萬(wàn)噸一路震蕩走高，糧食生產(chǎn)得到快

3、速發(fā)展，但波動(dòng)也更頻繁復(fù)雜。在1996年總產(chǎn)量首次跨上50000萬(wàn)噸的大臺(tái)階，達(dá)到了50453萬(wàn)噸，增長(zhǎng)率為8.13%。但在2000年開(kāi)始出現(xiàn)了幾年的連續(xù)減產(chǎn)的現(xiàn)象，曾一路降到43069萬(wàn)噸的局面，一下子退回到十年前的水平，讓人擔(dān)憂。從2004年以來(lái)的5年里，我國(guó)糧食產(chǎn)量連續(xù)5年增產(chǎn)。在2008年糧食產(chǎn)量達(dá)到52870萬(wàn)噸，據(jù)中國(guó)農(nóng)業(yè)部稱，中國(guó)糧食產(chǎn)量可能達(dá)到歷史最高水平。從歷史的發(fā)展趨勢(shì)中，不難看出糧食產(chǎn)量的波動(dòng)性。因此，對(duì)我國(guó)糧食生產(chǎn)影響因素的實(shí)證研究就顯得十分有必要，以此尋找我國(guó)糧食穩(wěn)定增產(chǎn)的有效途徑。影響我國(guó)糧食產(chǎn)量因素的實(shí)證分析（一）模型的設(shè)定影響糧食生產(chǎn)的因素很多，有勞動(dòng)力、物質(zhì)投

4、入、土地、生產(chǎn)方式、技術(shù)進(jìn)步、生產(chǎn)結(jié)構(gòu)、制度因素、氣候變化和自然災(zāi)害等等因素都影響著糧食產(chǎn)量。為了基本涵蓋這些基本因素，本文選擇了以農(nóng)業(yè)化肥施用量、糧食播種面積、成災(zāi)面積、有效灌溉面積、農(nóng)業(yè)機(jī)械總動(dòng)力等為解釋變量，以糧食產(chǎn)量為被解釋變量。對(duì)于這些影響因素，我找到了1995到2008年的關(guān)于這些因素的數(shù)據(jù)，借此來(lái)分析中國(guó)糧食產(chǎn)量的影響因素和它們是如何來(lái)影響糧食產(chǎn)量的，以下是我對(duì)所找的數(shù)據(jù)做的一些說(shuō)明。1、糧食產(chǎn)量。作為被解釋變量，從表中我們可以知道它是波動(dòng)不定的，但總體趨勢(shì)還是增長(zhǎng)的。2、播種面積。隨著播種面積的減少，糧食產(chǎn)量也會(huì)相應(yīng)的減產(chǎn)，二者成正相關(guān)的關(guān)系。3、成災(zāi)面積。成災(zāi)面積的增加會(huì)使糧

5、食產(chǎn)量減少，它們是負(fù)相關(guān)的關(guān)系。4、有效灌溉面積。指具有一定的水源，地塊比較平整，灌溉工程或設(shè)備已經(jīng)配套，在一般年景下，當(dāng)年能夠進(jìn)行正常灌溉的耕地面積。在一般情況下，有效灌溉面積應(yīng)等于灌溉工程或設(shè)備已經(jīng)配備，能夠進(jìn)行正常灌溉的水田和水澆地面積之和。它是反映我國(guó)耕地抗旱能力的一個(gè)重要指標(biāo)。與糧食產(chǎn)量成正相關(guān)。5、農(nóng)業(yè)機(jī)械總動(dòng)力。包括耕作機(jī)械、排灌機(jī)械、收獲機(jī)械、農(nóng)用運(yùn)輸機(jī)械、植物保護(hù)機(jī)械、牧業(yè)機(jī)械、林業(yè)機(jī)械、漁業(yè)機(jī)械和其他農(nóng)業(yè)機(jī)械等。它的增加也會(huì)使糧食產(chǎn)量增加。（二）變量的定義假設(shè)糧食產(chǎn)量與農(nóng)業(yè)化肥施用量、糧食播種面積、成災(zāi)面積、有效灌溉面積、農(nóng)業(yè)機(jī)械總動(dòng)力之間存在線性關(guān)系，其中Y表示糧食產(chǎn)量X

6、1表示播種面積X2表示成災(zāi)面積X3表示農(nóng)用機(jī)械總動(dòng)力X4表示有效灌溉面積X5化肥施用量樣本時(shí)間從1995年到2008年樣本大?。簄=14(三) 樣本數(shù)據(jù)的選擇以2009年版的中國(guó)統(tǒng)計(jì)年鑒2009為資料來(lái)源，使用了1995年到2008年的中國(guó)糧食產(chǎn)量和播種面積、成災(zāi)面積、有效灌溉面積、農(nóng)業(yè)機(jī)械總動(dòng)力和農(nóng)業(yè)化肥施用量的數(shù)據(jù)，對(duì)影響糧食產(chǎn)量的因素進(jìn)行實(shí)證分析。表一：中國(guó)1995-2008年的糧食產(chǎn)量及其相關(guān)影響因素統(tǒng)計(jì)表年份糧食產(chǎn)量（萬(wàn)噸）/Y播種面積（千公頃）/X1成災(zāi)面積（千公頃）/X2農(nóng)用機(jī)械總動(dòng)力(萬(wàn)千瓦）/X3有效灌溉面積(千公頃）/X4化肥施用量(萬(wàn)噸）/X51995 46662 110

7、060 45824 36118 49281 3594 1996 50454 112548 46991 38547 50381 3828 1997 49417 112912 53427 42016 51239 3981 1998 51230 113787 50145 45208 52296 4084 1999 50839 113161 49980 48996 53158 4124 2000 46218 108463 54688 52574 53820 4146 2001 45264 106080 52215 55172 54249 4254 2002 45706 103891 46946 579

8、30 54355 4339 2003 43070 99410 54506 60387 54014 4412 2004 46947 101606 37106 64028 54478 4637 2005 48402 104278 38818 68398 55029 4766 2006 49804 104958 41091 72522 55750 4928 2007 50160 105638 48992 76590 56518 5108 2008 52871 106793 39990 82190 58472 5239 數(shù)據(jù)來(lái)源：2009年版的中國(guó)統(tǒng)計(jì)年鑒2009（四）經(jīng)濟(jì)計(jì)量模型的建立建立糧食產(chǎn)量與

9、農(nóng)業(yè)化肥施用量、糧食播種面積、成災(zāi)面積、有效灌溉面積、農(nóng)業(yè)機(jī)械總動(dòng)力的一個(gè)5元線性回歸模型: 其中：0 、1 、2 、3 、4、5是待定參數(shù).是隨機(jī)誤差項(xiàng) 樣本大小: n=14 總平方和的自由度: n-1回歸平方和的自由度：k=5 (解釋變量的個(gè)數(shù))殘差平方和的自由度: n-k-1 待定參數(shù)個(gè)數(shù)：k+1=6利用Eviews軟件，采用已搜集的數(shù)據(jù)對(duì)模型進(jìn)行OLS回歸，結(jié)果如下表所示：Dependent Variable: YMethod: Least SquaresDate: 05/29/10 Time: 12:43Sample: 1995 2008Included observations:

10、14CoefficientStd. Errort-StatisticProb.  C-36622.5116184.84-2.2627670.0535X10.6181520.0774927.9769390.0000X2-0.1191200.039245-3.0353100.0162X3-0.1494520.172520-0.8662880.4116X4-0.1321930.379116-0.3486890.7363X59.0857783.6519022.4879580.0376R-squared0.968455    Mean depe

11、ndent var48360.29Adjusted R-squared0.948740    S.D. dependent var2768.524S.E. of regression626.8128    Akaike info criterion16.01670Sum squared resid3143154.    Schwarz criterion16.29058Log likelihood-106.1169    Hannan-

12、Quinn criter.15.99135F-statistic49.12175    Durbin-Watson stat2.060123Prob(F-statistic)0.000009由此可見(jiàn)，該模型=0.968455，調(diào)整的=0.948740很高，F(xiàn)=49.12175明顯顯著，但是當(dāng)a=0.05時(shí)，ta/2（14-5）= t0.05/2（14-5）=2.262，此時(shí)，X3，X4的T檢驗(yàn)不顯著，這表明可能存在很嚴(yán)重的多重共線性。（五）多重共線性的的檢驗(yàn)計(jì)算各解釋變量的相關(guān)系數(shù)，選擇X1,X2,X3,X4,X5數(shù)據(jù)，通過(guò)Eviews軟件，得到相關(guān)系數(shù)

13、矩陣如下：X1X2X3X4X5X110.3108225895494686-0.6486044282501619-0.5654178356088598-0.5724658205807768X20.31082258954946861-0.4914505578111244-0.3698118956560495-0.5079255069214623X3-0.6486044282501619-0.491450557811124410.96764785817280910.9898645265333025X4-0.5654178356088598-0.36981189565604950.9676478581

14、72809110.9539102051937531X5-0.5724658205807768-0.50792550692146230.98986452653330250.95391020519375311由相關(guān)系數(shù)矩陣可以看出，各解釋變量相互之間的相關(guān)系數(shù)較高，證實(shí)確實(shí)存在多重共線性.（六）修正多重共線性 表6-1 Y對(duì)X1的一元回歸結(jié)果Dependent Variable: YMethod: Least SquaresDate: 05/29/10 Time: 13:20Sample: 1995 2008Included observations: 14CoefficientStd. Erro

15、rt-StatisticProb.  C9851.44715228.340.6469150.5299X10.3585590.1416742.5308720.0264R-squared0.348014    Mean dependent var48360.29Adjusted R-squared0.293682    S.D. dependent var2768.524S.E. of regression2326.743    Akaike info cri

16、terion18.47389Sum squared resid64964769    Schwarz criterion18.56518Log likelihood-127.3172    Hannan-Quinn criter.18.46544F-statistic6.405314    Durbin-Watson stat0.610605Prob(F-statistic)0.026377 表6-2 Y對(duì)X2的一元回歸結(jié)果Dependent Variable: YMetho

17、d: Least SquaresDate: 05/29/10 Time: 13:23Sample: 1995 2008Included observations: 14CoefficientStd. Errort-StatisticProb.  C55590.326057.3779.1772930.0000X2-0.1531970.127421-1.2022960.2524R-squared0.107509    Mean dependent var48360.29Adjusted R-squared0.033135 

18、60;  S.D. dependent var2768.524S.E. of regression2722.270    Akaike info criterion18.78788Sum squared resid88929059    Schwarz criterion18.87918Log likelihood-129.5152    Hannan-Quinn criter.18.77943F-statistic1.445515  

19、  Durbin-Watson stat0.547277Prob(F-statistic)0.252437 表6-3 Y對(duì)X3的一元回歸結(jié)果Dependent Variable: YMethod: Least SquaresDate: 05/29/10 Time: 13:24Sample: 1995 2008Included observations: 14CoefficientStd. Errort-StatisticProb.  C46568.113233.14814.403340.0000X30.0313370.0549490.5702870.57

20、90R-squared0.026387    Mean dependent var48360.29Adjusted R-squared-0.054747    S.D. dependent var2768.524S.E. of regression2843.299    Akaike info criterion18.87488Sum squared resid97012161    Schwarz criterion18.96617L

21、og likelihood-130.1242    Hannan-Quinn criter.18.86643F-statistic0.325228    Durbin-Watson stat0.755302Prob(F-statistic)0.579004 表6-4 Y對(duì)X4的一元回歸結(jié)果Dependent Variable: YMethod: Least SquaresDate: 05/29/10 Time: 13:26Sample: 1995 2008Included observations: 14Coeff

22、icientStd. Errort-StatisticProb.  C37349.9817404.922.1459430.0530X40.2046960.3232740.6331970.5385R-squared0.032331    Mean dependent var48360.29Adjusted R-squared-0.048308    S.D. dependent var2768.524S.E. of regression2834.606    

23、;Akaike info criterion18.86876Sum squared resid96419874    Schwarz criterion18.96005Log likelihood-130.0813    Hannan-Quinn criter.18.86031F-statistic0.400939    Durbin-Watson stat0.732162Prob(F-statistic)0.538485表6-5 Y對(duì)X5的一元回歸結(jié)果 Dependent

24、Variable: YMethod: Least SquaresDate: 05/29/10 Time: 13:27Sample: 1995 2008Included observations: 14CoefficientStd. Errort-StatisticProb.  C41638.486966.5455.9769200.0001X51.5316611.5784090.9703830.3510R-squared0.072761    Mean dependent var48360.29Adjusted R-squared-0.

25、004509    S.D. dependent var2768.524S.E. of regression2774.759    Akaike info criterion18.82608Sum squared resid92391436    Schwarz criterion18.91737Log likelihood-129.7826    Hannan-Quinn criter.18.81763F-statistic0.941

26、643    Durbin-Watson stat0.742663Prob(F-statistic)0.351009 表6-6 一元回歸結(jié)果 變量X1X2X3X4X5參數(shù)估計(jì)量0.358559-0.1531970.0313370.2046961.531661T2.530872-1.2022960.5702870.6331970.970383R20.3480140.1075090.0263870.0323310.8653140.2936820.033135-0.054747-0.048308-0.004509排序12543可見(jiàn)Y和X1最好以X1為基礎(chǔ)，順次

27、加入其他變量逐步回歸，做4個(gè)二元回歸，結(jié)果如下表 表6-7 Y對(duì)X1,X2二元回歸結(jié)果Dependent Variable: YMethod: Least SquaresDate: 05/29/10 Time: 13:49Sample: 1995 2008Included observations: 14CoefficientStd. Errort-StatisticProb.  C10847.9411867.280.9141040.3803X10.4654720.1161124.0088220.0021X2-0.2644150.089257-2.9623890.0129R

28、-squared0.637342    Mean dependent var48360.29Adjusted R-squared0.571404    S.D. dependent var2768.524S.E. of regression1812.477    Akaike info criterion18.03019Sum squared resid36135802    Schwarz criterion18.16713Log l

29、ikelihood-123.2113    Hannan-Quinn criter.18.01751F-statistic9.665784    Durbin-Watson stat0.661049Prob(F-statistic)0.003778表6-8 Y對(duì)X1,X3二元回歸結(jié)果Dependent Variable: YMethod: Least SquaresDate: 05/29/10 Time: 13:51Sample: 1995 2008Included observations: 14Coeffici

30、entStd. Errort-StatisticProb.  C-40365.7510778.69-3.7449580.0032X10.7294800.0898108.1225110.0000X30.1815080.0285056.3676440.0001R-squared0.860868    Mean dependent var48360.29Adjusted R-squared0.835571    S.D. dependent var2768.524S.E. of regression1

31、122.633    Akaike info criterion17.07215Sum squared resid13863349    Schwarz criterion17.20909Log likelihood-116.5051    Hannan-Quinn criter.17.05947F-statistic34.03069    Durbin-Watson stat1.335633Prob(F-statistic)0.000

32、019表6-9 Y對(duì)X1,X4二元回歸結(jié)果Dependent Variable: YMethod: Least SquaresDate: 05/29/10 Time: 13:52Sample: 1995 2008Included observations: 14CoefficientStd. Errort-StatisticProb.  C-64207.8621053.06-3.0498120.0111X10.6178900.1142895.4063810.0002X40.8590560.2140624.0131090.0020R-squared0.735406 

33、   Mean dependent var48360.29Adjusted R-squared0.687298    S.D. dependent var2768.524S.E. of regression1548.153    Akaike info criterion17.71492Sum squared resid26364557    Schwarz criterion17.85186Log likelihood-121.0045

34、0;   Hannan-Quinn criter.17.70225F-statistic15.28653    Durbin-Watson stat0.775275Prob(F-statistic)0.000667表6-10 Y對(duì)X1,X5二元回歸結(jié)果Dependent Variable: YMethod: Least SquaresDate: 05/29/10 Time: 13:53Sample: 1995 2008Included observations: 14CoefficientStd. Errort-Statis

35、ticProb.  C-46430.529700.650-4.7863310.0006X10.6729520.0717689.3767260.0000X55.1306990.6704807.6522810.0000R-squared0.896893    Mean dependent var48360.29Adjusted R-squared0.878147    S.D. dependent var2768.524S.E. of regression966.4230  &#

36、160; Akaike info criterion16.77249Sum squared resid10273707    Schwarz criterion16.90943Log likelihood-114.4074    Hannan-Quinn criter.16.75981F-statistic47.84275    Durbin-Watson stat1.509845Prob(F-statistic)0.000004表6-11 加入新變量的回歸結(jié)果（一

37、）變量X1X2X3X4X5X1,X20.465472(4.008822)-0.264415(-2.962389)0.571404X1,X30.729480(8.122511)0.181508(6.367644)0.835571X1,X40.617890(5.406381)0.859056(4.013109) 0.687298X1,X50.672952(9.376726)5.130699(7.652281)0.878147可見(jiàn)：Y與X1,X5組合最好選擇保留X5，再加入其他新變量逐步回歸，結(jié)果如下表：表6-12 Y對(duì)X1,X5,X2三元回歸結(jié)果Dependent Variable: YMetho

38、d: Least SquaresDate: 05/29/10 Time: 14:40Sample: 1995 2008Included observations: 14CoefficientStd. Errort-StatisticProb.  C-37137.566923.738-5.3638020.0003X10.6781630.04798514.132780.0000X54.3296790.4946428.7531590.0000X2-0.1342790.035111-3.8244150.0033R-squared0.958131   &

39、#160;Mean dependent var48360.29Adjusted R-squared0.945571    S.D. dependent var2768.524S.E. of regression645.9000    Akaike info criterion16.01412Sum squared resid4171869.    Schwarz criterion16.19671Log likelihood-108.0989   

40、; Hannan-Quinn criter.15.99722F-statistic76.28040    Durbin-Watson stat1.750457Prob(F-statistic)0.000000表6-13 Y對(duì)X1,X5,X3三元回歸結(jié)果Dependent Variable: YMethod: Least SquaresDate: 05/29/10 Time: 14:42Sample: 1995 2008Included observations: 14CoefficientStd. Errort-StatisticProb.&#

41、160; C-51786.319872.338-5.2455970.0004X10.5722270.0957865.9739980.0001X512.234244.7959882.5509330.0288X3-0.2622920.175525-1.4943320.1660R-squared0.915714    Mean dependent var48360.29Adjusted R-squared0.890429    S.D. dependent var2768.524S.E. of regressi

42、on916.4243    Akaike info criterion16.71379Sum squared resid8398335.    Schwarz criterion16.89638Log likelihood-112.9965    Hannan-Quinn criter.16.69689F-statistic36.21475    Durbin-Watson stat1.850221Prob(F-statistic)0.

43、000011表6-14 Y對(duì)X1,X5,X4三元回歸結(jié)果Dependent Variable: YMethod: Least SquaresDate: 05/29/10 Time: 14:43Sample: 1995 2008Included observations: 14CoefficientStd. Errort-StatisticProb.  C-21518.2613718.30-1.5685800.1478X10.6619960.06129710.799850.0000X58.4667491.5739335.3793570.0003X4-0.7134610.313

44、686-2.2744470.0462R-squared0.932046    Mean dependent var48360.29Adjusted R-squared0.911660    S.D. dependent var2768.524S.E. of regression822.8607    Akaike info criterion16.49841Sum squared resid6770997.    Schwarz cri

45、terion16.68099Log likelihood-111.4889    Hannan-Quinn criter.16.48151F-statistic45.71972    Durbin-Watson stat2.554920Prob(F-statistic)0.000004表6-15 加入新變量的回歸結(jié)果(二)變量X1X2X3X4X5X1,X5,X20.67816314.13278-0.134279-3.8244154.3296798.7531590.945571X1,X5,X30.5722275.97

46、3998-0.262292-1.49433212.234242.5509330.890429X1,X5,X40.66199610.79985-0.713461-2.2744478.4667495.3793570.911660由此看出，Y對(duì)X1,X2,X5的組合最好最后結(jié)果：Y = -37137.5600712 + 0.678162714038*X1 - 0.134279330758*X2 + 4.32967873942*X5=0.958131， F=76.28040， DW=1.750457=0.945571ta/2(n-k)= t0.05/2（14-5）=2.262Fa（k-1，n-k）= F0.05（5-1，14-5）= F0.05（4，9）=6經(jīng)濟(jì)解釋： 播種面積每增加1千公頃，糧食產(chǎn)量就增加0.678萬(wàn)噸。 成災(zāi)面積每增加1千公頃，糧食產(chǎn)量就減少0.134萬(wàn)噸。 化肥施用量每增加1萬(wàn)噸，糧食產(chǎn)量就