張曉東 張宇 宋繼華 王欣蘭

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

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

中圖分類(lèi)號(hào) S-9文獻(xiàn)標(biāo)識(shí)碼 A文章編號(hào) 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

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

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

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

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

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

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

1.2 研究方法

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

2 結(jié)果與分析

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

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

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

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

2.3 VAR模型分析

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

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

對(duì)于滯后長(zhǎng)度為2的VAR（2）模型而言，AR特征多項(xiàng)式有10個(gè)復(fù)數(shù)根，2個(gè)實(shí)數(shù)根。檢驗(yàn)結(jié)果見(jiàn)圖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)

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

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

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

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

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

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

2.7 方差分解

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

3 結(jié)論與建議

3.1 結(jié)論

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

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

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

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

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

3.2 政策建議

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

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

（2）活豬價(jià)格的波動(dòng)會(huì)引起玉米價(jià)格的波動(dòng)，這是由于我國(guó)的養(yǎng)豬業(yè)規(guī)?；潭冗€偏低，到2016年我國(guó)規(guī)?；B(yǎng)殖為50%左右，還存在大量的散養(yǎng)戶，散養(yǎng)戶在活豬價(jià)格上漲時(shí)，大量囤積玉米，引起了玉米價(jià)格的波動(dò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è)的主要飼料，但玉米價(jià)格波動(dòng)不會(huì)引起活豬價(jià)格的波動(dòng)。這主要是由于我國(guó)的玉米產(chǎn)量大，養(yǎng)豬業(yè)消耗掉的玉米僅是我國(guó)玉米產(chǎn)量的一部分，另外我國(guó)養(yǎng)豬業(yè)主要分布在四川、河南、湖南、江西、福建、山東等幾個(gè)主要省份，而玉米主產(chǎn)區(qū)分布在黑龍江、吉林、遼寧、內(nèi)蒙古、山西、河北等幾個(gè)省份，二者的分布區(qū)域不一致，養(yǎng)豬業(yè)不能充分消化玉米的產(chǎn)量。因此，要優(yōu)化全國(guó)生豬養(yǎng)殖區(qū)域布局，引導(dǎo)產(chǎn)能向玉米主產(chǎn)區(qū)以及環(huán)境容量大的地區(qū)轉(zhuǎn)移。

參考文獻(xiàn)

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

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

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

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

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

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

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

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