□張家澍

（河南大學歐亞國際學院，河南 開封 475004）

1 研究背景

糧食安全問題關(guān)乎國計民生。我國人口基數(shù)較大且正處于發(fā)展中國家階段，糧食問題顯得尤為突出。隨著社會經(jīng)濟轉(zhuǎn)型和糧食種植面積減少，我國糧食產(chǎn)量曾經(jīng)連續(xù)5 年減產(chǎn)。2003 年以后，我國糧食產(chǎn)量穩(wěn)步增長，截至2021 年，我國糧食產(chǎn)量連續(xù)7 年保持在6.5 億t 以上。隨著社會主義現(xiàn)代化建設(shè)進程加快對社會發(fā)展的影響，房地產(chǎn)企業(yè)大量開墾土地進行住宅建設(shè)，耕地面積逐年減少對我國糧食安全構(gòu)成了明顯威脅。

糧食安全生產(chǎn)在任何時候都不能放松。必須把滿足人們對基礎(chǔ)物質(zhì)資料的需求放在首位，穩(wěn)定糧食產(chǎn)量，抓好糧食生產(chǎn)。在此背景下，考察糧食產(chǎn)量影響因素對糧食產(chǎn)量的影響效果及程度十分必要。通過對糧食產(chǎn)量影響因素的考察，準確科學預測糧食產(chǎn)量，是保障糧食安全和農(nóng)業(yè)現(xiàn)代化發(fā)展的重要任務(wù)之一。查明糧食產(chǎn)量的影響因素，準確預測糧食產(chǎn)量，有利于開展包括國民生產(chǎn)總值、國民經(jīng)濟結(jié)構(gòu)等在內(nèi)的相關(guān)領(lǐng)域研究，對現(xiàn)代經(jīng)濟持續(xù)健康發(fā)展和人們生活消費質(zhì)量不斷提高具有重要意義。

糧食產(chǎn)量直接關(guān)系到國家和社會長期穩(wěn)定。在一定程度上，國家穩(wěn)定與糧食產(chǎn)量有著密切的關(guān)系?，F(xiàn)實生活中，除增加農(nóng)作物灌溉面積、施肥量外，收購價較低、收益不理想、流通渠道不順等都會影響我國農(nóng)戶的生產(chǎn)積極性，進而可能影響全國糧食作物的耕種面積。應(yīng)以保障農(nóng)戶利益為基本出發(fā)點，頒布實施一系列惠民政策，有效保障全國糧食總產(chǎn)量，加強對糧食的宏觀調(diào)控，切實保障農(nóng)民收入。

2 糧食產(chǎn)量影響因素的文獻綜述

目前國內(nèi)學者對糧食產(chǎn)量影響因素的分析主要采用回歸分析、分位數(shù)回歸、灰色關(guān)聯(lián)分析等方法。巫琦玲和張葵（2017）在近年有關(guān)糧食產(chǎn)量文獻的基礎(chǔ)上，歸納出使用頻率最高的影響因素與準確率最高的模型，建立以糧食產(chǎn)量為因變量，以糧食播種面積、有效灌溉面積、農(nóng)業(yè)機械總動力、化肥施用折純量、受災面積、成災面積為自變量的多元線性回歸模型，對2000—2015 年樣本數(shù)據(jù)進行了詳細分析。

李玉虎和陳俊美（2019）選取2000—2017 年河南省統(tǒng)計數(shù)據(jù)，選用灰色關(guān)聯(lián)分析模型，對影響河南省糧食產(chǎn)量的8 個因子進行建模分析，得出農(nóng)藥施用實物量、化肥施用折純量和耕地面積是影響糧食產(chǎn)量的重要因素，農(nóng)村用電量對糧食產(chǎn)量的影響程度最低，灌溉面積、鄉(xiāng)村從業(yè)人員、農(nóng)用機械總動力和農(nóng)用塑料薄膜使用量對糧食產(chǎn)量的影響程度依次降低等結(jié)論。

毛偉（2012）利用HP 濾波處理1979—2009 年湖北省糧食產(chǎn)量數(shù)據(jù)，得出其波動特征，并將各時間序列的長期趨勢剔除，使序列僅含波動項以進行分位數(shù)分析。其研究結(jié)果表明，在糧食產(chǎn)量的低分位點上，各因素對糧食增產(chǎn)的作用明顯；在其他分位點處，作用不明顯。由此提出，要確保糧食播種面積，科學投放化肥，保障農(nóng)業(yè)生產(chǎn)用電，鼓勵糧食生產(chǎn)規(guī)?；晟萍Z食生產(chǎn)的風險管理工具和制度，以保證糧食產(chǎn)量和保障糧食安全。

3 我國糧食產(chǎn)量影響因素的實證模型構(gòu)建、檢驗及估計

文章研究有效灌溉面積、肥料施用量、大型拖拉機數(shù)量、農(nóng)用柴油機數(shù)量4 個因素對糧食產(chǎn)量的影響，從回歸分析的角度對糧食產(chǎn)量影響因素進行分析。對通過回歸分析建立的模型加以統(tǒng)計檢驗、計量經(jīng)濟學等相關(guān)檢驗，并將模型加以實際運用，結(jié)合其他學者的研究方法進行驗證，為糧食生產(chǎn)提供一些切實可行的建議。

3.1 數(shù)據(jù)來源、變量設(shè)定及模型設(shè)計

數(shù)據(jù)來源于《中國統(tǒng)計年鑒》收集的31 個省（自治區(qū)、直轄市）的糧食產(chǎn)量及相關(guān)影響因素數(shù)據(jù)。相關(guān)數(shù)據(jù)準確、詳細，可以用于研究糧食產(chǎn)量的影響因素。根據(jù)相關(guān)數(shù)據(jù)與信息，建立如下多元回歸模型。

式中：代表糧食產(chǎn)量；代表有效灌溉面積；代表肥料施用量；代表大型拖拉機數(shù)量；代表農(nóng)用柴油機數(shù)量。

將收集的樣本數(shù)據(jù)，采取Rstudio 軟件進行操作運算，由回歸結(jié)果得出初步線性回歸模型如下。

在實證研究結(jié)果中，=0.873 6，=44.94。

3.2 樣本數(shù)據(jù)的描述性統(tǒng)計

從表1 可以看出，每個變量的標準差和極差都比較大，而二者均可以用來衡量數(shù)據(jù)的離散程度，故可以得出該數(shù)據(jù)的離散程度相對較大。由此說明，各地區(qū)糧食產(chǎn)量存在較大差異，各地區(qū)之間農(nóng)業(yè)發(fā)展水平存在較大差異，各地區(qū)之間發(fā)展極不平衡。

表1 變量間的描述性統(tǒng)計

3.3 多重共線性檢驗與逐步回歸

為了確保模型設(shè)計的科學性，文章采取多重共線性檢驗來觀察變量之間的相關(guān)關(guān)系。由表2 可知，變量與之間的相關(guān)系數(shù)為0.816 9，同時，與間的相關(guān)系數(shù)也比較大，達0.704 4。當解釋變量之間相關(guān)性很強時，即有理由說明模型存在明顯的多重共線性。因此，必須對模型進行相關(guān)統(tǒng)計檢驗和計量經(jīng)濟學檢驗，并予以修正，消除模型的多重共線性。下面采用逐步回歸法，解決模型出現(xiàn)的多重共線性問題。

表2 相關(guān)系數(shù)矩陣

采用逐步回歸的方法，依次建立因變量與自變量的線性回歸模型?；谥鸩交貧w命令所提取的3 個模型分別進行線性回歸，由回歸結(jié)果可知，建立的關(guān)于f（，，，）模型調(diào)整后的判別系數(shù)值為0.854 2，但僅有變量和的檢驗顯著；關(guān)于f（，，）模型調(diào)整后的判別系數(shù)值為0.853 5，并且仍只有變量和通過了關(guān)于回歸系數(shù)顯著性判斷的檢驗；關(guān)于f（，）模型調(diào)整后的判別系數(shù)為0.849 8，盡管相比于前兩次回歸，模型擬合優(yōu)度稍有下降，但此時模型中包含的全部變量和均通過了顯著性檢驗。因此，確定最終模型為剔除了變量和后，關(guān)于和的回歸分析，最終模型以=f（，）為最優(yōu)。

為了提升研究的科學性和嚴謹性，采用vif 方差膨脹因子的大小判斷是否存在共線性，默認膨脹因子小于10 為可容忍的范圍，即方差膨脹因子vif 小于10 就認為模型不存在多重共線性；若vif 大于10，則有理由認為模型存在多重共線性。變量、的方差膨脹因子均為1.299 718＜10，故通過逐步回歸對模型的修正后，此時模型不存在多重共線性，見表3。

表3 變量vif 值

與此同時，各回歸系數(shù)在顯著性水平下均顯著，故最終認定此模型合理，模型如下。

根據(jù)模型（3）可知，肥料施用量每增加1 萬t，糧食產(chǎn)量平均提高6.633 1 萬t；大型拖拉機數(shù)量每增加1 000 臺，糧食產(chǎn)量平均提高3.525 6 萬t。可見，肥料施用量是影響我國糧食產(chǎn)量的顯著性因素。從現(xiàn)實意義看，過度施用肥料會在一定程度上對糧食產(chǎn)量產(chǎn)生負面影響。從回歸結(jié)果得知，隨機干擾項的數(shù)值為77.855 4，說明仍有其他因素影響糧食產(chǎn)量。

3.4 統(tǒng)計檢驗與異方差檢驗

3.4.1 統(tǒng)計檢驗

從前文研究可以發(fā)現(xiàn)，兩個影響因素變量所對應(yīng)的參數(shù)估計值均為正值，即肥料施用量的適當增加和大型拖拉機數(shù)量農(nóng)業(yè)機械數(shù)量的增加對糧食產(chǎn)量均有正向促進作用。=0.859 8，調(diào)整=0.849 8，仍比較靠近1。因此，文章認為該模型的擬合是可行的且擬合效果比較好，變化的80%以上均可由解釋變量進行解釋。

考察方差整體顯著性，=85.84，計算相關(guān)臨界值，值仍小于顯著性水平，表明選擇線性模型是合理的，回歸方程是顯著的，即和聯(lián)合起來對模型因變量的作用效果是顯著的。

考察單個變量的顯著性，此時變量和的檢驗值均大于0.05 的置信水平，故此時參數(shù)均是顯著的。

3.4.2 異方差檢驗

通過對模型進行BP 異方差檢驗，其中原假設(shè)H為模型不存在異方差，備擇假設(shè)H為模型存在異方差。由程序運行結(jié)果可知，值為0.226 6，高于顯著性水平0.05，認為回歸模型中不會出現(xiàn)異方差現(xiàn)象。

4 我國糧食產(chǎn)量影響因素的實證結(jié)果及對策

4.1 我國糧食產(chǎn)量影響因素實證模型的估計結(jié)果

由模型（2）可知，影響產(chǎn)量的4 個因素均對糧食產(chǎn)量具有正向促進作用。當其他變量保持不變時，有效灌溉面積每增多1 單位，糧食產(chǎn)量相應(yīng)平均增多0.173 5 萬t；肥料施用量每施加1 萬t，糧食產(chǎn)量相應(yīng)平均增多4.575 2 萬t。當其他變量保持不變時，大型拖拉機數(shù)量每投入1 000 臺，糧食產(chǎn)量相應(yīng)平均增多3.170 8 萬t；農(nóng)用柴油機數(shù)量每投入1 000 臺，糧食產(chǎn)量相應(yīng)平均增多0.444 4 萬t。=0.873 6，調(diào)整=0.854 2，比較靠近1，認為該模型擬合是可行的且擬合效果比較好，變化的80%以上可有解釋變量進行解釋。=44.94，計算相關(guān)臨界值，同時值小于顯著性水平，即回歸方程是顯著的。變量和的t 檢驗值大于0.05 的置信水平，故此時參數(shù)均是顯著的。

由模型（3）可知，當其他變量保持不變時，肥料施用量每提高1 萬t，糧食產(chǎn)量相應(yīng)提高6.633 1 萬t；保持肥料施用量不變，大型拖拉機每增多1 000 臺，糧食產(chǎn)量相應(yīng)增多3.525 6 萬t。利用模型（3）進行預測分析，設(shè)置解釋變量值，預測相應(yīng)的糧食產(chǎn)量。由表4可得出在給定具體值的基礎(chǔ)上對應(yīng)的95%的取值范圍（上界和下界）。其中，fit 為糧食產(chǎn)量的預測均值產(chǎn)量，lwr 和upr 分別為取300、取400 時對應(yīng)的糧食產(chǎn)量區(qū)間估計的上下界，即設(shè)定肥料施用量為300 萬t，設(shè)定大型拖拉機數(shù)量為400 000 臺，通過predict 函數(shù)預測出糧食產(chǎn)量的概率為0.95 的相應(yīng)區(qū)間。

表4 我國糧食產(chǎn)量預測模型實證結(jié)果

4.2 提升我國糧食產(chǎn)量的建議

雖然目前我國糧食產(chǎn)量處于自給自足并有剩余的狀態(tài)，但隨著人口增加和耕地面積減少，保障糧食安全的任務(wù)依然艱巨。各地區(qū)之間農(nóng)業(yè)發(fā)展水平存在很大的不平衡性，未來應(yīng)當注重協(xié)調(diào)地區(qū)間農(nóng)業(yè)發(fā)展，充分發(fā)揮各地區(qū)優(yōu)勢，實現(xiàn)協(xié)調(diào)發(fā)展?；趯ξ覈Z食產(chǎn)量影響因素的分析結(jié)果，在最終模型設(shè)定的基礎(chǔ)上，提出提升我國糧食產(chǎn)量的政策建議。

第一，合理施用化肥，增加有機肥施用量。合理施用化肥對提升糧食產(chǎn)量有一定的促進作用，但過度施用化肥會帶來一系列問題，因此應(yīng)該注重對化肥的合理施用。施用有機肥可以培肥地力，有助于提高產(chǎn)量和土地資源利用率。應(yīng)以科學的態(tài)度合理施肥，增加有機肥施用量。

第二，提升農(nóng)業(yè)機械化水平，推廣對大型拖拉機的使用。大型拖拉機數(shù)量增加能夠使耕地效率得到顯著提升，故應(yīng)進一步加強基礎(chǔ)設(shè)施和配套設(shè)施建設(shè)，加大對大型拖拉機的推廣力度。同時，政府應(yīng)當鼓勵科研人員研發(fā)更高效、先進的農(nóng)用機械，提高農(nóng)業(yè)生產(chǎn)效率。

第三，政府應(yīng)當審時度勢，加大財政對糧食生產(chǎn)的支持力度。在糧食產(chǎn)量不穩(wěn)定時，應(yīng)當對農(nóng)產(chǎn)品以正常價格回購，避免農(nóng)民生活水平大幅下降。在農(nóng)業(yè)健康發(fā)展時，注意穩(wěn)定市價，并適度提升糧價，促進農(nóng)民收入增加。除模型中所呈現(xiàn)的糧食產(chǎn)量影響因素外，還有很多其他因素影響糧食產(chǎn)量。政府部門有必要增加政策投入，自覺承擔起對糧食產(chǎn)品市場價格合理干預和管理的責任，避免谷賤農(nóng)傷，并采取糧食補助等優(yōu)惠政策幫助農(nóng)戶增收，提升農(nóng)民生產(chǎn)糧食的積極性，從而提高糧食產(chǎn)量和質(zhì)量。