張曉東 張宇 宋繼華 王欣蘭

摘要 利用2009年1月—2017年12月的月度數(shù)據(jù)，建立了VAR（2）模型，分析了活豬存欄量、仔豬價格、能繁母豬量、玉米價格、豆粕價格對活豬價格傳導(dǎo)效應(yīng)。結(jié)果表明：存欄量、能繁母豬量、仔豬價格與活豬價格存在雙向的Granger因果關(guān)系，活豬價格是玉米價格單向Granger原因，豆粕價格與活豬價格不存在Granger 因果關(guān)系；脈沖響應(yīng)分析表明，對活豬價格波動沖擊較大的因素有活豬價格自身、豆粕價格、活豬存欄量；方差分解結(jié)果表明，生豬價格波動貢獻(xiàn)率由大到小的順序是活豬自身價格、豆粕價格、活豬存欄量、仔豬價格、玉米價格和能繁母豬量波動?；趯?shí)證分析結(jié)果，提出幾點(diǎn)政策建議。

關(guān)鍵詞 VAR模型；活豬；價格波動；傳導(dǎo)效應(yīng)

中圖分類號 S-9文獻(xiàn)標(biāo)識碼 A文章編號 0517-6611（2018）15-0208-04

Abstract An VAR（2）model is established， based on the monthly statistical data from Jan. 2009 to Dec. 2017， to analyze the conductive effect of the stock quantity of hogs， the price of piglets， the quantity of breeding sows， the price of corn and the price of soybean on hogs price fluctuation. The results showed that there is a two-way Granger causality between hogs price and the factors of stock quantity， breeding sows and price of piglets；the price of hogs is the one-way Granger cause of corn price；and there is no Granger causality between the price of soybean meal and the price of hogs. The analysis of impulse response shows that the factors that contribute greater impulse are the price of hogs itself， the price of soybean and the stock quantity of hogs. The result of variance decomposition gives the descending order of contribution rate towards hogs price： the price of hogs itself， the price of soybean， the stock quantity of hogs， the price of piglets， the price of corn and the quantity fluctuation of breeding sows. Finally some political suggestions are put forward based on the results of empirical analysis.

Key words VAR model；Hogs；Price fluctuation；Conductive effect

我國養(yǎng)豬業(yè)自20世紀(jì)80年代中期市場化改革后，活豬價格開始大幅度波動。2010年4月—2011年11月出現(xiàn)一次價格大漲，漲幅高達(dá)107%，其后出現(xiàn)多次波動。2015年3月—2016年5月，又出現(xiàn)一次價格的大漲，漲幅高達(dá)66.67%?；钬i價格的大起大落不僅會影響生豬養(yǎng)殖者收入，還影響到其他農(nóng)產(chǎn)品價格的平穩(wěn)。

關(guān)于豬價的研究在波動周期[1-2]、波動原因[3-5]等方面取得了很多有價值的研究成果，對該研究有重要的啟示，但也存在一定的缺陷。他們在分析的養(yǎng)豬業(yè)價格波動時多采用豬肉價格為研究對象，在影響因素的選取上，多以玉米價格、仔豬價格為主，活豬存欄量、能繁母量多采用統(tǒng)計(jì)年鑒的年底存欄量數(shù)據(jù)，不能直接反映養(yǎng)豬業(yè)的波動規(guī)律，把豆粕價格作為影響因素較少，將活豬價格、玉米價格、豆粕價格、仔豬價格、活豬存欄量、能繁母量作為一個整體系統(tǒng)來考慮則更少。筆者采用系統(tǒng)的思維，采用權(quán)威網(wǎng)站發(fā)布的月度數(shù)據(jù)，運(yùn)用VAR模型實(shí)證分析活豬價格的波動原因，分析影響因素的傳導(dǎo)作用，為我國政府部門推進(jìn)農(nóng)業(yè)供給側(cè)結(jié)構(gòu)性改革，穩(wěn)定生豬生產(chǎn)提供一定的參考。

1 變量選取、數(shù)據(jù)來源與研究方法

1.1 變量選取與數(shù)據(jù)來源

養(yǎng)豬業(yè)一個生產(chǎn)周期至少需要1年，在生產(chǎn)過程中活豬的存欄量直接決定了活豬的供給量，而活豬供給量取決于能繁母豬的數(shù)量。因此，從活豬自身生產(chǎn)系統(tǒng)考慮，選取活豬存欄量、能繁母豬量作為影響變量，分別用NPIG、SOW表示。養(yǎng)豬業(yè)生產(chǎn)成本中，玉米、豆粕在飼料成本中所占比重較高，仔豬的成本也是養(yǎng)豬業(yè)生產(chǎn)成本的構(gòu)成部分。因此，從養(yǎng)豬業(yè)生產(chǎn)成本考慮，選取玉米價格、豆粕價格和仔豬價格作為影響變量，分別用CORN、BEAN、PIGLET表示。活豬價格作為被解釋變量，用HOG表示。

該研究選取2009年1月—2017年12月月度數(shù)據(jù)，共108個樣本，活豬存欄量、能繁母豬量數(shù)據(jù)來源農(nóng)業(yè)部畜牧司公布數(shù)據(jù)，玉米價格、豆粕價格和仔豬價格數(shù)據(jù)來源于中國畜牧信息網(wǎng)。為消除時間序列里的異方差現(xiàn)象，使變量的趨勢線性化，將變量取自然對數(shù)。活豬價格、活豬存欄量、能繁母豬量、玉米價格、豆粕價格和仔豬價格分別記為lnHOG、lnNPIG、lnSOW、lnCORN、lnBEAN、lnPIGLET。

1.2 研究方法

該研究建立VAR模型對活豬價格波動傳導(dǎo)效應(yīng)進(jìn)行研究。VAR（向量自回歸模型）模型是包含多個方程的非結(jié)構(gòu)化模型，即變量之間的關(guān)系不是以經(jīng)濟(jì)理論來為基礎(chǔ)的。其建模思想是把每個外所有生變量作為內(nèi)生變量滯后值的函數(shù)來構(gòu)造模型，其實(shí)質(zhì)上是考察多個變量之間的動態(tài)互動關(guān)系[6-8]。VAR模型一般表達(dá)式為：

2 結(jié)果與分析

2.1 數(shù)據(jù)平穩(wěn)性檢驗(yàn)

在回歸分析之前，為避免非平穩(wěn)時間序列產(chǎn)生的偽回歸，需對數(shù)據(jù)進(jìn)行平穩(wěn)性檢驗(yàn)。該研究采用常見的ADF檢驗(yàn)，結(jié)果見表1。在0.05的顯著性水平條件下，lnHOG、lnNPIG、lnSOW、lnCORN、lnBEAN、lnPIGLET都是非平穩(wěn)的。

2.2 協(xié)整關(guān)系檢驗(yàn)

為進(jìn)一步檢驗(yàn)變量間是否存在長期穩(wěn)定的關(guān)系，該研究采用 Johansen 協(xié)整檢驗(yàn)，檢驗(yàn)變量間的協(xié)整關(guān)系，檢驗(yàn)結(jié)果見表2。跡檢驗(yàn)結(jié)果顯示，拒絕原假設(shè)“存在至多一個協(xié)整關(guān)系”，即變量存在2個協(xié)整關(guān)系，說明這幾個變量間存在協(xié)整關(guān)系，可以通過取對數(shù)序列建立 VAR 模型。

2.3 VAR模型分析

VAR模型的不同的滯后期會導(dǎo)致模型估計(jì)結(jié)果的顯著不同，滯后期的選擇依據(jù)根據(jù)AIC或SC值最小原則來確定。AIC和SC值均是在滯后2期為最小，確定出滯后階數(shù)為VAR（2）模型。VAR（2）模型參數(shù)結(jié)果見表3。

2.4 VAR（2）模型穩(wěn)定性檢驗(yàn)

對于滯后長度為2的VAR（2）模型而言，AR特征多項(xiàng)式有10個復(fù)數(shù)根，2個實(shí)數(shù)根。檢驗(yàn)結(jié)果見圖1。AR特征多項(xiàng)式的根的倒數(shù)小于1，既AR特征多項(xiàng)式復(fù)數(shù)根和實(shí)數(shù)根都位于單位圓內(nèi)，由此可以得出VAR（ 2）模型是穩(wěn)定的。

2.5 Granger 因果關(guān)系檢驗(yàn)

為判斷活豬價格波動與其影響因素之間是否存在因果關(guān)系，以及它們的影響方向，故運(yùn)用 Granger 因果關(guān)系檢驗(yàn)進(jìn)行分析。由于取對數(shù)后的數(shù)據(jù)均為非平穩(wěn)，因此，采用H-P濾波，分解出其波動項(xiàng)，采用ADF檢驗(yàn)進(jìn)行檢驗(yàn)，波動項(xiàng)數(shù)據(jù)均平穩(wěn)。Granger 因果關(guān)系檢驗(yàn)結(jié)果見表4。

由表4可知，在0.05顯著水平下，存欄量與活豬價格存在雙向的Granger 因果關(guān)系，即存欄量的波動會引起豆粕價格的波動，活豬價格的波動也會引起存欄量的波動；活豬價格與能繁母豬量存在雙向的Granger 因果關(guān)系，即能繁母豬量的波動會引起活豬價格的波動，活豬價格的波動也會引起能繁母豬量的波動；活豬價格與仔豬價格存在雙向的Granger 因果關(guān)系，即仔豬價格的波動會引起活豬價格的波動，活豬價格的波動也會引起仔豬價格的波動；活豬價格僅是玉米價格的單向的Granger 因果關(guān)系，即活豬價格的波動會引起玉米價格的波動，玉米價格的波動不會引起活豬價格的波動；豆粕價格與活豬價格不存在Granger 因果關(guān)系，即豆粕價格的波動不會引起活豬價格的波動，活豬價格的波動也不會引起豆粕價格的波動。

由上分析得出，活豬存欄量、能繁母豬量、仔豬價格是活豬價格的重要的因素，活豬存欄量、能繁母豬量、仔豬價格與活豬價格構(gòu)成一個長期相互影響動態(tài)的系統(tǒng)。豆粕價格、玉米價格并不是活豬價格的重要影響因素，原因是我國玉米的生產(chǎn)量過剩，其用途也非常廣泛，養(yǎng)豬業(yè)消耗的玉米僅是玉米產(chǎn)量的一小部分，豆粕是大豆產(chǎn)品的副產(chǎn)品，養(yǎng)豬業(yè)消耗掉的豆粕也只是產(chǎn)量的一部分。

2.6 脈沖響應(yīng)分析

該研究使用脈沖響應(yīng)函數(shù)分析活豬存欄量（lnNPIG）、仔豬價格（lnPIGLET）、能繁母豬量（lnSOW）、玉米價格（lnCORN）、豆粕價格（lnBEAN）對活豬價格（lnHOG）的動態(tài)影響過程，結(jié)果見圖2。

由圖2可以看出，活豬價格受到一個標(biāo)準(zhǔn)差大小沖擊后，活豬價格對自身的沖擊有正向影響，在第3期到達(dá)最大值，其后迅速減小，第9期后沖擊影響趨于平緩，這種效果將持續(xù)10個月以上?；钬i存欄量對活豬價格在第1～4期有正向影響，第4期以后有負(fù)向影響，第7期后沖擊趨于平緩?；钬i存欄量增加，在4～5個月后，出欄數(shù)量增加，活豬價格降低，說明活豬存欄量對豬肉價格的促進(jìn)作用有一定的滯后性，符合活豬生產(chǎn)規(guī)律。仔豬價格對活豬價格的沖擊在第1～4期有負(fù)向影響，沖擊影響較弱，第4期后有正向影響，沖擊影響較弱。能繁母豬量對活豬價格的沖擊在1～5期有負(fù)向影響，第5期后有正向影響，但第1～10期沖擊影響均較弱。玉米價格對活豬價格的沖擊與能繁母豬量對活豬價格沖擊效果類似。豆粕價格對活豬價格的沖擊為負(fù)向影響，第1期為0，其后沖擊影響變強(qiáng)，在第4期后沖擊影響趨于平緩，豆粕價格的沖擊對活豬價格上漲有一定的影響，且影響是持續(xù)的。

2.7 方差分解

在VAR模型的基礎(chǔ)上對lnHOG進(jìn)行方差分解，可以得出：lnHOG變動主要由自身沖擊產(chǎn)生的，當(dāng)其穩(wěn)定時貢獻(xiàn)率約占60%左右，活豬價格波動也是由于自身的供求關(guān)系造成的；其次是豆粕價格，當(dāng)其穩(wěn)定時貢獻(xiàn)率約占25%左右；再次是活豬存欄量，當(dāng)其穩(wěn)定時貢獻(xiàn)率約占10%左右，能繁母豬量對活豬價格貢獻(xiàn)率最低，僅為0.6%左右。對生豬價格波動貢獻(xiàn)率由大到小的順序是活豬自身價格、豆粕價格、活豬存欄量、仔豬價格、玉米價格和能繁母豬量波動。

3 結(jié)論與建議

3.1 結(jié)論

利用2009年1月—2017年12月的月度數(shù)據(jù)，分析了活豬存欄量、仔豬價格、能繁母豬量、玉米價格、豆粕價格對活豬價格的傳導(dǎo)效應(yīng)，得出如下結(jié)論：

（1）確定出滯后階數(shù)為VAR（2）模型，且VAR（2）模型是穩(wěn)定的。

（2）Granger檢驗(yàn)表明：存欄量與活豬價格存在雙向Granger原因，即存欄量的波動會導(dǎo)致活豬價格的波動，活豬價格的波動也會導(dǎo)致存欄量的波動；活豬價格與能繁母豬量存在雙向Granger原因，即活豬價格的波動會導(dǎo)致能繁母豬量的波動，能繁母豬量的波動也會導(dǎo)致活豬價格的波動；活豬價格與仔豬價格存在雙向Granger原因，即活豬價格的波動會導(dǎo)致仔豬價格的波動，仔豬價格也會導(dǎo)致活豬價格的波動；活豬價格是玉米價格雙單向Granger原因，即活豬價格的波動會導(dǎo)致玉米價格的波動，而玉米價格也不會導(dǎo)致活豬價格的波動；豆粕價格與活豬價格不存在Granger 因果關(guān)系，即活豬價格的波動不會導(dǎo)致豆粕價格的波動，豆粕價格也不會導(dǎo)致活豬價格的波動。

（3）脈沖響應(yīng)分析表明：對活豬價格波動沖擊較大的因素有活豬價格自身、豆粕價格，活豬存欄量。

（4）方差分解結(jié)果表明：生豬價格波動60%左右是由其自身波動引起的，25%左右是由豆粕價格波動引起的，10%左右是由活豬存欄量波動引起的，5%左右是由仔豬價格、玉米價格和能繁母豬量波動引起的。

3.2 政策建議

結(jié)合前文分析，提出以下幾點(diǎn)政策建議：

（1）存欄量波動與活豬價格波動互為因果關(guān)系，活豬價格波動與能繁母豬量波動互為因果關(guān)系，活豬價格波動與仔豬價格波動存互為因果關(guān)系，活豬價格波動僅是玉米價格波動的原因。由此可以看出，活豬價格波動與活豬存欄量、仔豬價格、能繁母豬量波動是一個動態(tài)的系統(tǒng)，所以應(yīng)該建立一個生豬價格波動的中長期動態(tài)監(jiān)測系統(tǒng)，在動態(tài)監(jiān)測中，應(yīng)當(dāng)把活豬存欄量、仔豬價格、能繁母豬量作為監(jiān)測重點(diǎn)，及時發(fā)布準(zhǔn)確信息，合理引導(dǎo)生產(chǎn)和流通，促使活豬價格平穩(wěn)運(yùn)行。

（2）活豬價格的波動會引起玉米價格的波動，這是由于我國的養(yǎng)豬業(yè)規(guī)?；潭冗€偏低，到2016年我國規(guī)?；B(yǎng)殖為50%左右，還存在大量的散養(yǎng)戶，散養(yǎng)戶在活豬價格上漲時，大量囤積玉米，引起了玉米價格的波動。因此要積極推進(jìn)農(nóng)業(yè)供給側(cè)結(jié)構(gòu)性改革，發(fā)展規(guī)模高效養(yǎng)豬業(yè)，實(shí)行規(guī)模化、現(xiàn)代化和集約化養(yǎng)豬，穩(wěn)定生豬生產(chǎn)。

（3）雖然玉米是養(yǎng)豬業(yè)的主要飼料，但玉米價格波動不會引起活豬價格的波動。這主要是由于我國的玉米產(chǎn)量大，養(yǎng)豬業(yè)消耗掉的玉米僅是我國玉米產(chǎn)量的一部分，另外我國養(yǎng)豬業(yè)主要分布在四川、河南、湖南、江西、福建、山東等幾個主要省份，而玉米主產(chǎn)區(qū)分布在黑龍江、吉林、遼寧、內(nèi)蒙古、山西、河北等幾個省份，二者的分布區(qū)域不一致，養(yǎng)豬業(yè)不能充分消化玉米的產(chǎn)量。因此，要優(yōu)化全國生豬養(yǎng)殖區(qū)域布局，引導(dǎo)產(chǎn)能向玉米主產(chǎn)區(qū)以及環(huán)境容量大的地區(qū)轉(zhuǎn)移。

參考文獻(xiàn)

[1] 陳蓉.我國生豬生產(chǎn)波動周期分析[J].農(nóng)業(yè)技術(shù)經(jīng)濟(jì)，2009（3）：77-86.

[2] 岳冬冬，王征兵.我國生豬生產(chǎn)波動周期研究[J].農(nóng)業(yè)技術(shù)經(jīng)濟(jì)，2010（10）：18-25.

[3] 徐萬剛.當(dāng)前我國新一輪豬價波動的成因與對策分析[J].天府新論，2011（6）：74-79.

[4] 吳英.基于蛛網(wǎng)理論的我國近10年生豬價格波動分析[J].廣東農(nóng)業(yè)科學(xué)，2011（3）：106-108.

[5] 寧攸涼，喬娟.中國生豬價格波動的影響與成因探究[J].中國畜牧業(yè)，2010（2）：52-56.

[6] 張曉東，王超群，王寧，等.基于VAR模型的生豬價格波動的實(shí)證研究[J].安徽農(nóng)業(yè)科學(xué)，2014，42（8）：2475-2477，2480.

[7] 黃海峰，嵇炳繽，任亮，等.農(nóng)產(chǎn)品價格對CPI傳導(dǎo)效應(yīng)研究：以秈稻、生豬、棉花價格為例[J].價格理論與實(shí)踐，2016（12）：94-97.

[8] 吳振信，薛冰，王書平.基于VAR模型的油價波動對我國經(jīng)濟(jì)影響分析[J].中國管理科學(xué)，2011，19（1）：21-28.