王春鳳 沈曼君

（1．江蘇開放大學(xué)，江蘇 南京 210036；2．河海大學(xué)，江蘇 南京 210036）

青貯飼料由高含水量植物性飼料發(fā)酵而成，比新鮮飼料更耐儲(chǔ)存，營養(yǎng)更豐富。青貯玉米的廣泛種植，為發(fā)展肉牛、肉羊等草食畜牧業(yè)提供了支撐[1-2]。2015 年，青貯玉米飼料價(jià)格受到國外低價(jià)玉米沖擊[3]，國內(nèi)青貯玉米價(jià)格連年下降?！凹Z改飼”政策出臺(tái)，緩解了國內(nèi)玉米市場需求。青貯玉米飼料作為草食畜牧業(yè)養(yǎng)殖的重要支撐，價(jià)格的頻繁波動(dòng)不利于草食畜牧業(yè)的發(fā)展，對國家糧食安全產(chǎn)生影響。因此，研究我國青貯玉米飼料價(jià)格波動(dòng)的因素，制定應(yīng)對策略極為重要。從成本視角分析，張玲等[4]以河北省為樣本，分析不同模式下青貯玉米飼料的成本效益，得出種養(yǎng)結(jié)合模式下的青貯玉米飼料成本效益要優(yōu)于外購模式。從需求視角分析，謝長城[5]以養(yǎng)殖產(chǎn)能恢復(fù)性增長為切入點(diǎn)，提出“去庫存、降成本、補(bǔ)短板”為核心的供給側(cè)結(jié)構(gòu)性改革，使玉米及相關(guān)替代品進(jìn)口減少，飼料原料成本下降。從供給視角分析，呂曉英等[6]構(gòu)建了系統(tǒng)動(dòng)力學(xué)模型和多方案模擬，得出未來我國玉米供需關(guān)系將長期處于偏緊格局。從經(jīng)濟(jì)效益視角分析，石子墨等[7]對青貯玉米飼料展開探究，得出在全株玉米青貯的過程中，二次發(fā)酵、有氧腐敗嚴(yán)重影響玉米青貯發(fā)酵的品質(zhì)與價(jià)格，給草食畜牧業(yè)養(yǎng)殖者帶來經(jīng)濟(jì)損失。張紅燕等[8]指出，與單一青貯飼料相比，混合青貯飼料營養(yǎng)物質(zhì)更豐富，可以明顯提高養(yǎng)殖效益。目前，關(guān)于青貯玉米飼料價(jià)格的研究多停留在成本、供需以及經(jīng)濟(jì)方面，對價(jià)格波動(dòng)的研究較少。青貯玉米飼料價(jià)格波動(dòng)的影響因素較廣，單一視角并不具有普適性。本文借助VAR模型實(shí)證，分析青貯玉米飼料價(jià)格波動(dòng)的自身因素、供給因素以及需求因素，篩選出我國青貯玉米飼料價(jià)格波動(dòng)的主要因素，助力我國草食畜牧業(yè)穩(wěn)步發(fā)展。

1 實(shí)證設(shè)計(jì)

1.1 實(shí)證方法

VAR 模型可借助當(dāng)期變量對滯后變量進(jìn)行回歸[9-10]。VAR 模型主要用于分析經(jīng)濟(jì)沖擊對經(jīng)濟(jì)變量的影響。對于青貯玉米飼料價(jià)格，影響因素可能源于內(nèi)部、外部，數(shù)量較多，適用于VAR模型。VAR模型一般形式為：

式中：n=1,2,…,N；Yn為內(nèi)生變量向量；Xn為外生變量向量；p為滯后期數(shù)；N為樣本數(shù)；A、B為系數(shù)矩陣；εn為擾動(dòng)向量。

1.2 數(shù)據(jù)來源與處理

2022 年3 月—2023 年2 月青貯玉米飼料價(jià)格波動(dòng)見圖1。由圖1 可知，青貯玉米飼料價(jià)格存在波動(dòng)性。以2005 年1 月—2022 年12 月的204 個(gè)月數(shù)據(jù)為樣本，分析影響青貯玉米飼料價(jià)格變動(dòng)的因素。本文數(shù)據(jù)來源為歷年《中國統(tǒng)計(jì)年鑒》《中國農(nóng)產(chǎn)品價(jià)格調(diào)查年鑒》、中國飼料行業(yè)信息網(wǎng)以及艾媒數(shù)據(jù)中心網(wǎng)站。將豆粕價(jià)格、肉牛存欄量、肉羊存欄量、生豬存欄量、牛肉價(jià)格、羊肉價(jià)格、豬肉價(jià)格、城鎮(zhèn)人均可支配收入、鄉(xiāng)村人均可支配收入為自變量，青貯玉米價(jià)格為因變量。為避免變量之間異方差導(dǎo)致結(jié)果不準(zhǔn)確，對變量數(shù)值進(jìn)行對數(shù)處理，即青貯玉米價(jià)格（lnSCP）、豆粕價(jià)格（lnSMP）、肉牛存欄量（lnBCS）、肉羊存欄量（lnMSS）、生豬存欄量（lnLHS）、牛肉價(jià)格（lnBP）、羊肉價(jià)格（lnMP）、豬肉價(jià)格（lnPP）、城鎮(zhèn)人均可支配收入（lnCCDI）、鄉(xiāng)村人均可支配收入（lnRPCDI）。

圖1 2022年3月—2023年2月青貯玉米飼料價(jià)格波動(dòng)

2 實(shí)證分析

2.1 單位根檢驗(yàn)（見表1、表2）

表1 原始變量平穩(wěn)性檢驗(yàn)結(jié)果

表2 一階差分變量平穩(wěn)性檢驗(yàn)結(jié)果

為避免出現(xiàn)偽回歸問題，在模型回歸之前選用ADF檢驗(yàn)法檢驗(yàn)各變量序列平穩(wěn)性[11]。統(tǒng)計(jì)量低于顯著性水平5%臨界值，說明變量平穩(wěn)，拒絕原假設(shè)。統(tǒng)計(jì)量高于顯著性水平5%臨界值，需要對各變量進(jìn)行差分處理，直至所得結(jié)果拒絕原假設(shè)，得到最終平穩(wěn)階數(shù)。Stata17.0軟件進(jìn)行模型檢驗(yàn)。

由表1 可知，ADF 各變量統(tǒng)計(jì)量均大于顯著性水平5%臨界值，說明原始變量為非平穩(wěn)序列，需對其一階差分序列進(jìn)行單位根檢驗(yàn)。

由表2 可知，一階差分處理后，變量具備平穩(wěn)性，可進(jìn)行下一步分析。

2.2 格蘭杰因果檢驗(yàn)（見表3）

表3 格蘭杰因果檢驗(yàn)結(jié)果

VAR 模型檢驗(yàn)前需進(jìn)行變量間格蘭杰因果關(guān)系檢驗(yàn)[12-13]。P值低于5%說明拒絕原假設(shè)，認(rèn)為格蘭杰因果關(guān)系成立。由表3 可知，我國青貯玉米飼料價(jià)格受豆粕價(jià)格、肉牛存欄量、肉牛價(jià)格、肉羊存欄量、肉羊價(jià)格鄉(xiāng)村人均可支配收入的影響。豆粕價(jià)格、牛肉價(jià)格、羊肉價(jià)格受青貯玉米飼料價(jià)格的影響。

2.3 VAR模型構(gòu)建及檢驗(yàn)（見表4）

表4 最優(yōu)滯后階數(shù)選擇

為保證VAR模型計(jì)算結(jié)果精確，需選擇最優(yōu)的滯后階數(shù)[14-15]。由表4 可知，經(jīng)過滯后期準(zhǔn)則檢驗(yàn)，F(xiàn)PE、AIC、SC和HQ準(zhǔn)則均認(rèn)為應(yīng)建立VAR（2）模型。為確保模型精確性，最優(yōu)滯后階數(shù)應(yīng)選擇二階。

VAR 模型是非結(jié)構(gòu)化的多方程模型[16-17]，將其兩邊各項(xiàng)指標(biāo)取對數(shù)，得到式（2）：

式中：n=1,2,…,216；αi為滯后項(xiàng)系數(shù)；εn為方程誤差項(xiàng)。

借助AR 根圖檢驗(yàn)VAR 模型穩(wěn)定性。單位根均落于圓內(nèi)，VAR 模型穩(wěn)定。反之，單位根均落于圓外情況，該模型可能為無效估計(jì)，無法進(jìn)行下一步研究[18]。

AR特征多項(xiàng)式逆根圖見圖2。

圖2 AR特征多項(xiàng)式逆根圖

由圖2 可知，所有單位根均落于圓內(nèi)，表明各指標(biāo)系統(tǒng)穩(wěn)定，VAR模型有效。

2.4 方差分解

方差分解刻畫出不同變量沖擊貢獻(xiàn)在總貢獻(xiàn)中的比重[19-20]，對比出不同變量效應(yīng)的相對大小。為探究我國青貯玉米飼料價(jià)格波動(dòng)的影響因素，對VAR模型各變量進(jìn)行方差分解，見表5。

表5 方差分解

由表5 可知，玉米飼料價(jià)格波動(dòng)的因素分為3 個(gè)方面。第一方面，從飼料價(jià)格視覺分析，青貯玉米飼料自身價(jià)格對價(jià)格波動(dòng)的影響最大，第1 期貢獻(xiàn)率最高為100%，第10 期下降至52.71%。第二方面，從供給視角分析，豆粕價(jià)格對青貯玉米飼料價(jià)格波動(dòng)的影響最大，第10 期最高貢獻(xiàn)率為8.92%；肉牛存欄量最高貢獻(xiàn)率次之，為7.49%；豬肉價(jià)格貢獻(xiàn)率最小，為1.74%。第三方面，從需求視角分析，城市人均可支配收入、鄉(xiāng)村人均可支配收入對青貯玉米飼料價(jià)格波動(dòng)的最高貢獻(xiàn)率分別為4.93%、8.01%。整體分析，除青貯玉米飼料自身價(jià)格對青貯玉米飼料價(jià)格變化的貢獻(xiàn)率最大，豆粕價(jià)格在第10 期對青貯玉米飼料價(jià)格變化的影響最顯著，鄉(xiāng)村人均可支配收入貢獻(xiàn)率次之，豬肉價(jià)格貢獻(xiàn)率最低。

3 結(jié)論與建議

3.1 結(jié)論

本文以2005—2021 年青貯玉米飼料價(jià)格月度數(shù)據(jù)為樣本，借助VAR模型探究青貯玉米飼料價(jià)格本身、豆粕價(jià)格、肉牛存欄量等10個(gè)因素對青貯玉米飼料價(jià)格波動(dòng)的影響。

從直接影響視角分析，青貯玉米飼料自身價(jià)格對青貯玉米飼料價(jià)格波動(dòng)影響最大。短期內(nèi)對價(jià)格波動(dòng)的影響最大。隨著時(shí)間推移，影響程度有所下降。

從供給視角分析，豆粕價(jià)格對青貯玉米飼料價(jià)格波動(dòng)的影響最大，第10 期貢獻(xiàn)率為8.92%。肉牛存欄量、肉羊存欄量、生豬存欄量在第10期的最高貢獻(xiàn)率依次為7.49%、4.33%、2.03%；羊肉價(jià)格對青貯玉米飼料價(jià)格波動(dòng)的影響最大，貢獻(xiàn)率為5.88%，牛肉價(jià)格貢獻(xiàn)率為3.96%，豬肉價(jià)格貢獻(xiàn)率為1.74%。

從需求視角分析，城市人均可支配收入、鄉(xiāng)村人均可支配收入對青貯玉米飼料價(jià)格波動(dòng)的最高貢獻(xiàn)率分別為4.93%、8.01%。

因此，青貯玉米飼料價(jià)格波動(dòng)的主要因素為豆粕價(jià)格、鄉(xiāng)村人均可支配收入、肉牛存欄量以及羊肉價(jià)格。

3.2 建議

第一，完善動(dòng)態(tài)監(jiān)測預(yù)警系統(tǒng)。政府相關(guān)部門需綜合青貯玉米飼料價(jià)格的主要因素，從多個(gè)維度構(gòu)建綜合性“價(jià)格波動(dòng)預(yù)警體系”，以便青貯玉米飼料價(jià)格產(chǎn)生異常波動(dòng)時(shí)，作出響應(yīng)。首先，政府對價(jià)格異常波動(dòng)進(jìn)行判別，確定預(yù)警指標(biāo)臨界值，設(shè)置預(yù)警區(qū)間。其次，政府需定期發(fā)布玉米、豆粕、小麥等青貯玉米飼料相關(guān)的市場信息，防止價(jià)格大幅波動(dòng)。最后，地方部門結(jié)合當(dāng)?shù)仫暳袭a(chǎn)業(yè)與養(yǎng)殖產(chǎn)業(yè)的實(shí)際情況，以青貯玉米飼料成本收益指數(shù)等指標(biāo)，預(yù)測青貯玉米飼料價(jià)格的趨勢。

第二，推進(jìn)市場結(jié)構(gòu)性改革。在供給端層面，政府需持續(xù)優(yōu)化糧食結(jié)構(gòu)，擴(kuò)大青貯玉米等飼料作物種植面積，強(qiáng)化青貯玉米飼料供給。按照“政府推動(dòng)、企業(yè)帶動(dòng)、農(nóng)戶聯(lián)動(dòng)、市場互動(dòng)”的思路，實(shí)施訂單式種植、合同式種植等模式，提高青貯玉米飼料原料供給的穩(wěn)定性。政府還應(yīng)加速推進(jìn)農(nóng)村閑置土地流轉(zhuǎn)進(jìn)程，建立集中規(guī)劃與大片作業(yè)的青貯玉米飼料原料示范區(qū)，借助政府優(yōu)惠補(bǔ)貼政策，培育一批高度規(guī)?；?、專業(yè)化的現(xiàn)代飼料產(chǎn)業(yè)群，提高青貯玉米飼料供給能力。政府根據(jù)農(nóng)牧結(jié)合、循環(huán)發(fā)展的要求，統(tǒng)籌優(yōu)化區(qū)域養(yǎng)殖布局。如，對于牛、羊養(yǎng)殖量超過環(huán)境承載能力的地區(qū)，政府相關(guān)部門嚴(yán)格控制養(yǎng)殖總量，引導(dǎo)畜牧養(yǎng)殖向非超載區(qū)轉(zhuǎn)移。

第三，強(qiáng)化產(chǎn)業(yè)鏈抗風(fēng)險(xiǎn)能力。對上游生產(chǎn)環(huán)節(jié)，政府部門應(yīng)支持青貯玉米飼料原料的種植、運(yùn)輸、加工、存儲(chǔ)等環(huán)節(jié)的基礎(chǔ)設(shè)施建設(shè)，提高青貯玉米飼料生產(chǎn)率。對中游流通環(huán)節(jié)，需推動(dòng)主產(chǎn)區(qū)與主銷區(qū)之間的定向調(diào)運(yùn)，提升運(yùn)輸效率，避免流通不暢引發(fā)的產(chǎn)業(yè)鏈斷裂風(fēng)險(xiǎn)。對下游消費(fèi)環(huán)節(jié)，從多元化青貯玉米飼料銷售端完善。