許壘晶 張競文 王妍

[提要] 本文從國家統(tǒng)計局網(wǎng)站選取2004～2018年數(shù)據(jù)，以糧食產(chǎn)量的時間序列及其三個因素為樣本，建立多元線性回歸模型。針對河北省糧食產(chǎn)量對糧食播種面積、農(nóng)用化肥施用折純量和糧食受災(zāi)面積的依賴性進行研究，并提出提高河北省糧食產(chǎn)量的相應(yīng)建議。

關(guān)鍵詞：河北;糧食產(chǎn)量;回歸分析;對策建議

中圖分類號：F32 文獻標識碼：A

收錄日期：2020年2月16日

一、引言

雖然前幾年我國可以自行解決溫飽問題，但耕地面積減少、大量耕地浪費和農(nóng)業(yè)機械化布局不均對糧食生產(chǎn)造成的制約仍然不容忽視。

河北省可用耕地9，780.7萬畝，僅占河北省總面積即18.88萬平方千米的34.54%，可用耕地較少。本文綜合分析各種因素對糧食生產(chǎn)的影響，為提高糧食生產(chǎn)量，提出針對性的建設(shè)方案。

二、國內(nèi)研究現(xiàn)狀

李福道等以1992～2012年貴州省糧食總產(chǎn)量、糧食播種面積和有效灌溉面積為數(shù)據(jù)源，利用LOS分析貴州省糧食產(chǎn)量的影響因素，得出貴州省糧食生產(chǎn)的關(guān)鍵因素是糧食播種面積和有效灌溉面積的結(jié)論。

梁雨薇等選用1983～2012年數(shù)據(jù)，采用嶺回歸的方法對影響中國糧食產(chǎn)量的因素進行分析，得出化肥施用量、勞動力數(shù)量和糧食播種面積對糧食總產(chǎn)量有正相關(guān)影響，受影響面積對糧食總產(chǎn)量有負相關(guān)影響。

毛偉選用分位數(shù)回歸方法對湖北省1979～2009年糧食總產(chǎn)量影響因素進行分析，得出了糧食產(chǎn)量受糧食生產(chǎn)技術(shù)創(chuàng)新、農(nóng)業(yè)生產(chǎn)資料價格不斷攀升、投入糧食生產(chǎn)的人力等方面的影響。

趙俊偉等根據(jù)1994～2011年統(tǒng)計局網(wǎng)站上的相關(guān)資料、選取了糧食播種面積和農(nóng)機總動力等五個指標，對影響江西省糧食產(chǎn)量的因素進行了回歸分析，得到了影響江西省糧食產(chǎn)量的主要因素是糧食播種面積和農(nóng)機總動力的結(jié)論。

三、河北省糧食生產(chǎn)現(xiàn)狀

（一）糧食總產(chǎn)量現(xiàn)狀。從統(tǒng)計數(shù)據(jù)可以看出，河北省糧食生產(chǎn)總體呈上升趨勢，但增速逐年放緩的趨勢。雖然2018年河北省糧食總產(chǎn)量3，708.6萬噸，居全國第五位。但與2017年相比，糧食產(chǎn)量甚至出現(xiàn)了負增長。夏季糧食產(chǎn)量1，466.47萬噸，秋季糧食產(chǎn)量2，234.39萬噸。其中，玉米產(chǎn)量1，941.15萬噸，小麥產(chǎn)量1，450.73萬噸。糧食結(jié)構(gòu)單一，以玉米、小麥為主。

（二）糧食播種面積現(xiàn)狀。糧食播種面積大體上呈上升趨勢和較小的波動性，耕地面積增長速度緩慢。2018年河北省糧食播種面積為6，538.68千公頃，居全國第四位。糧食收獲季節(jié)分為夏秋兩季。

（三）機械化水平現(xiàn)狀。河北省農(nóng)作物種植綜合機械化程度達到80%，種植業(yè)機械化水平的山麓達到62%，小麥玉米生產(chǎn)過程中基本實現(xiàn)機械化，區(qū)域優(yōu)勢作物弱勢部位機械化水平明顯提高。但主要經(jīng)濟作物機械化薄弱環(huán)節(jié)未取得良好突破。農(nóng)業(yè)機械化和農(nóng)機裝備產(chǎn)業(yè)升級還面臨著一定的挑戰(zhàn)。

四、模型設(shè)定、數(shù)據(jù)搜集整理與參數(shù)估計

（一）模型構(gòu)建。以河北省2004～2018年近15年糧食產(chǎn)量、糧食播種面積、農(nóng)用化肥施用折純量和受災(zāi)面積為基礎(chǔ)，采用相關(guān)性分析Eviews10軟件，輸入數(shù)據(jù)后，通過相關(guān)圖、相關(guān)表可以看出每一個因素與糧食產(chǎn)量具有某種程度上的線性相關(guān)關(guān)系，進而將三大因素作為自變量，將糧食產(chǎn)量作為一個因變量，建立了一個多變量線性回歸方程：

Y=β0+β1X1+β2X2+β3X3+μ

上式中，Y代表河北省糧食產(chǎn)量，X1代表河北省糧食播種面積，X2代表化肥施用折純量，X3代表糧食受災(zāi)面積，μ代表隨機擾動項。

選用指標如灌溉面積、機械總動力等會導(dǎo)致高度的多重共線性，所以放棄了這些指標。由于分析工具的局限性與自身能力的不足遺漏了影響糧食產(chǎn)量的其他重要的解釋變量。

（二）樣本數(shù)據(jù)搜集與整理。本文數(shù)據(jù)來源于國家統(tǒng)計局網(wǎng)站，選用了河北省糧食產(chǎn)量（Y）、糧食播種面積（X1）、農(nóng)用化肥施用折純量（X2）、糧食受災(zāi)面積（X3）四個指標，對這些指標2004～2018年的15年時間序列數(shù)據(jù)的回歸分析。

（三）參數(shù)估計。將收集的數(shù)據(jù)用Eviews10軟件進行操作運算，得到上述模型的參數(shù)值。

首先分別對每個因素進行線性回歸分析結(jié)果如下：

Y=-9538.123+2.008584X1

（-5.103139）（6.781580）

Y=-4856.424+24.82491X2

（-5.661002）（9.321785）

Y=3800.728-0.413750X3

（14.39442）（-2.668803）

將以上三個因素結(jié)合分析對河北省糧食產(chǎn)量的影響結(jié)果如下：

Y=-7705.478+0.886623X1+16.43553X2-0.027382X3

（-5.435568）（3.158290）（5.611760）（-0.424090）

（四）時間序列的平穩(wěn)性檢驗。Y的ADF檢驗中，從Y的參數(shù)值來看，其t統(tǒng)計量為-3.933364（在0.05置信度下），小于測試臨界值-3.882778，拒絕存在單位根的零假設(shè)，所以可斷定河北省的糧食產(chǎn)量Y時間序列是穩(wěn)定的。

X1、X2、X3的ADF檢驗中，從參數(shù)值來看，其t統(tǒng)計量在置信度為0.05下分別為-3.791172、-3.933364、-3.098896分別小于其臨界值-2.442196、1.227849、-1.758476，均拒絕存在單位根的零假設(shè)，可斷定河北省糧食播種面積、化肥施用折純量、糧食受災(zāi)面積X1、X2、X3時間序列都是穩(wěn)定的。

五、模型的檢驗與修正

（一）經(jīng)濟意義檢驗。通過回歸處理，可以得出X1、X2的系數(shù)為正值，表明在其他情況不變的前提下Y（河北省糧食產(chǎn)量）與X1（河北省糧食播種面積）和X2（農(nóng)用化肥施用折純量）呈正比，與現(xiàn)實相符。X3的系數(shù)為負值，表明在其他條件不變的前提下Y（河北省糧食產(chǎn)量）與X3（糧食受災(zāi)面積）呈反比。與現(xiàn)實相符，因而所制定的模型是符合經(jīng)濟意義的。

（二）統(tǒng)計意義檢驗

1、擬合優(yōu)度檢驗。樣本可決定系數(shù)R2=0.94，修正的樣本可決定系數(shù)R2=0.92系數(shù)較高，說明模型擬合程度較好，河北省糧食產(chǎn)量可以被糧食播種面積等三個要素所解釋。

2、t檢驗。參數(shù)的t檢驗原假設(shè)H0：βj=0（j=0，1，2，3），在顯著水平α=0.05情況下查表得tα/2（n-k-1）=2.201，由分析結(jié)果得β0、β1、β2、β3對應(yīng)的t統(tǒng)計量分別為-5.4355683、158290、5.611760、-0.424090。β0、β1、β2的絕對值大于2.201，均拒絕H0，當(dāng)假設(shè)其他變量不變的情況下，X1、X2對解釋變量有顯著影響。X3的t檢驗未通過，分析第一個原因可能是由于此因素對河北省糧食產(chǎn)量的貢獻率不高，第二個原因是由于存在多重共線性，導(dǎo)致模型未通過檢驗。

3、F檢驗。方程的整體性檢驗原假設(shè)H0：β1=β2=β3，在顯著性α=0.05情況下查表得F（2，13）=6.7，F(xiàn)^F（2，13），拒絕原假設(shè)H0，即X1、X2、X3結(jié)合起來對Y有顯著性影響。

（三）計量經(jīng)濟學(xué)檢驗與修正

1、序列相關(guān)性檢驗。此模型為時間序列，所以對模型進行序列相關(guān)檢驗DW=2.054713，接受H0，隨機誤差項μt之間不存在序列相關(guān)性。

2、異方差檢驗。因為是時間序列所以不存在異方差。

3、多重共線性檢驗。rX1X2=0.69，rX1X3=-0.63，rX2X3=-0.49基本相關(guān)系數(shù)表明解釋變量不存在高度相關(guān)，也就是不存在多重線性。

通過處理，模型得到了修正，并且通過了相應(yīng)的t、F檢驗等，最終確定如下模型：

Y=-8104.577+0.940657X1+16.47913X2

（-7.919951）（3.895841）（2.825167）

最終通過相應(yīng)檢驗，得出最終結(jié)論：河北省糧食播種面積、河北省農(nóng)用化肥施用折純量是影響河北省糧食產(chǎn)量的主要因素。

六、建議

（一）擴大有效耕地面積。加大土地托管和土地流轉(zhuǎn)力度。總體來看，河北糧食生產(chǎn)穩(wěn)中有增，糧食播種面積在一定程度上也不斷地擴大，為保持糧食產(chǎn)量持續(xù)增長，河北省應(yīng)該加大土地托管和土地流轉(zhuǎn)力度，使閑置和廢棄耕地得到再利用，加大土地資源整合力度，整平土地和山脊，擴大有效耕地面積。但也不能一味的擴大耕地面積，要保證耕地質(zhì)量、改善土質(zhì)結(jié)構(gòu)、提高土壤肥力。

（二）合理控制化肥施用量。河北省化肥施用折純量大致呈逐年上升的趨勢。然而，過度使用化肥會導(dǎo)致土壤性質(zhì)惡化、產(chǎn)品質(zhì)量惡化和環(huán)境污染。為減少以上問題并且對河北省糧食產(chǎn)量的增加起促進作用，要合理控制化肥施用量，鼓勵使用經(jīng)過處理的農(nóng)家肥等有機肥，提高化肥質(zhì)量。

主要參考文獻：

[1]李福奪，楊興洪.貴州省糧食產(chǎn)量影響因素的偏最小二乘回歸分析[J].四川師范大學(xué)學(xué)報（自然科學(xué)版），2016.39.

[2]丁華，楊耀旭.河南省糧食產(chǎn)量影響因素的實證分析——基于多元線性回歸計量經(jīng)濟模型[J].糧食科技與經(jīng)濟，2015.40.

[3]梁雨薇，倪萍，王哲.中國糧食產(chǎn)量影響因素的嶺回歸分析[J].決策與信息，2015.

[4]毛偉.基于分位數(shù)回歸的糧食產(chǎn)量影響因素分析——以湖北省為例[J].湖南財政經(jīng)濟學(xué)院學(xué)報，2012.28.

[5]趙俊偉，王帥，王杰，李順德.基于回歸分析的糧食產(chǎn)量影響因素分析——以江西省為例[J].經(jīng)濟研究導(dǎo)刊，2013.

[6]劉宇鵬，賈健.河北省糧食綜合生產(chǎn)能力的現(xiàn)狀、存在問題及對策選擇[J].北方經(jīng)濟，2010.226.

[7]巫琦玲，張葵.基于回歸分析的糧食產(chǎn)量影響因素分析[J].糧食科技與經(jīng)濟，2017.