甘林針 郭 敬 張立中

(北京林業(yè)大學(xué)經(jīng)濟管理學(xué)院，北京，100083)

毛皮產(chǎn)業(yè)是中國北方重要的特色農(nóng)業(yè)產(chǎn)業(yè)，近幾年發(fā)展迅速。根據(jù)中國皮革協(xié)會2014～2016年發(fā)布的《中國貂、狐、貉取皮數(shù)量統(tǒng)計報告》和《中國貂、狐、貉存欄數(shù)量統(tǒng)計報告》，中國是全球最大的毛皮動物養(yǎng)殖國家，其三大主要毛皮動物的存欄總量達到1億只左右，年取皮量達到7500萬張左右，占全球皮張產(chǎn)量的一半以上。改革開放以來，毛皮市場從計劃經(jīng)濟轉(zhuǎn)為市場經(jīng)濟，轉(zhuǎn)型過程以及轉(zhuǎn)型后期，由于錯綜復(fù)雜的影響因素，毛皮價格“急漲暴跌”，超出合理的價格波動區(qū)間，以及養(yǎng)殖戶預(yù)期價格與市場價格之間出現(xiàn)偏差，都導(dǎo)致了中國農(nóng)業(yè)“惠農(nóng)惠民”的目標(biāo)在毛皮產(chǎn)業(yè)上發(fā)生錯位[1]。毛皮價格的異常波動影響了毛皮產(chǎn)業(yè)相關(guān)從業(yè)人員的收入，毛皮加工企業(yè)惡意壓縮加工成本致使毛皮質(zhì)量降低，高市場風(fēng)險甚至導(dǎo)致部分毛皮動物養(yǎng)殖戶退出生產(chǎn)等一系列問題[2]，嚴重制約著毛皮產(chǎn)業(yè)的健康發(fā)展。

毛皮作為特色農(nóng)產(chǎn)品，其價格波動特征及影響因素不僅受到相關(guān)從業(yè)人員的關(guān)注，近幾年也引起學(xué)術(shù)界的重視，部分研究組織和專家學(xué)者通過不同方法對國內(nèi)外毛皮價格問題進行研究，取得了豐碩成果。中國養(yǎng)殖業(yè)可持續(xù)發(fā)展戰(zhàn)略研究項目組[3]發(fā)現(xiàn)毛皮價格波動十分頻繁，周期不定，波幅較大，中國畜產(chǎn)品流通協(xié)會[4]指出整個裘皮市場行情、養(yǎng)殖數(shù)量、流行趨勢、其他因素等都對毛皮價格波動有影響。張立垚等[5]對比分析國內(nèi)外毛皮價格，在ARIMA模型預(yù)測的基礎(chǔ)上發(fā)現(xiàn)國際毛皮價格波動頻繁，劉鵬杰等[6]也發(fā)現(xiàn)國內(nèi)毛皮價格波動起伏較大，楊錫濤等[7]通過Comperz曲線對國內(nèi)毛皮均價進行生命周期分析，發(fā)現(xiàn)水貂皮價格波動幅度較大，狐皮價格波動周期不穩(wěn)定。綜合來看造成毛皮價格波動的影響因素主要包括經(jīng)濟因素和非經(jīng)濟因素，孫晨希[8]、李光玉等[9]、甘林針等[10]分析發(fā)現(xiàn)經(jīng)濟因素主要包括毛皮動物養(yǎng)殖和皮張加工成本、毛皮動物養(yǎng)殖數(shù)量、毛皮市場需求等，張旭等[11]、姜東暉等[12]、于超等[13]、李永欣[14]指出非經(jīng)濟因素主要包括毛皮質(zhì)量、國內(nèi)外產(chǎn)業(yè)政策影響、交易方式以及流行趨勢等。

總結(jié)已有的文獻可知，在內(nèi)容上學(xué)者們大多發(fā)現(xiàn)毛皮價格異常波動較為明顯，波動特征無規(guī)律性，并且影響毛皮價格波動的因素錯綜復(fù)雜，為本文的研究提供了一定的基礎(chǔ)。在研究方法上，學(xué)者們大多還是使用簡單的對比分析和定性解釋。而在其他農(nóng)產(chǎn)品價格波動的研究上，羅冬暉等[15]發(fā)現(xiàn)大量學(xué)者使用H-P濾波法對農(nóng)產(chǎn)品價格的波動進行趨勢分析和波動分析，并采用VAR模型、誤差修正模型等計量經(jīng)濟模型進行農(nóng)產(chǎn)品價格波動原因和影響因素的分析，得出豐富的研究成果。然而以上用于其他農(nóng)產(chǎn)品價格波動研究作出重要貢獻的方法卻鮮少用在毛皮價格波動的研究上。

鑒于此，本文選取因色澤艷麗、堅韌耐磨、輕柔暖和[16]而在毛皮市場上深受消費者喜愛的藍狐(北極狐Alopexlagopus)皮為例。首先以農(nóng)產(chǎn)品價格形成理論為基礎(chǔ)，結(jié)合文獻成果，梳理影響藍狐皮價格波動的影響因素并選取相應(yīng)的量化指標(biāo)；接著運用H-P濾波法考察毛皮價格波動周期長度、周期幅度等波動特征；最后采用向量誤差修正模型(VEC)探索毛皮價格波動各影響因素對毛皮價格的影響程度及動態(tài)影響機制，得出結(jié)論可為增加養(yǎng)殖戶及相關(guān)從業(yè)人員的收入、穩(wěn)定毛皮價格和推進毛皮產(chǎn)業(yè)健康發(fā)展提供理論參考依據(jù)。

1 研究方法

1.1 研究方法

1.1.1 H-P濾波法

Hodrick-Prescott濾波是一種被廣泛用于分析序列組成成分中長期趨勢的方法，Hodrick和Prescott(1980)在分析戰(zhàn)后美國經(jīng)濟周期的論文中首次使用該方法，該方法的原理如下：

(1)

(2)

1.1.2 向量誤差修正模型(VEC)

Engle和Granger將協(xié)整與誤差修正模型結(jié)合起來，建立了向量誤差修正模型(VEC)，可以稱之為含有協(xié)整約束的VAR模型，前提條件是變量間存在協(xié)整關(guān)系。變量間的長期和短期信息皆包含在內(nèi)，該模型先是通過協(xié)整檢驗進行長期動態(tài)分析，再通過脈沖響應(yīng)函數(shù)和方差分解進行短期波動分析，綜合得出數(shù)據(jù)本身存在的規(guī)律[18]，研究各影響因素對藍狐皮價格波動的影響程度及動態(tài)影響機制。

1.1.3 數(shù)據(jù)處理軟件

利用Excel、Stata13.0進行數(shù)據(jù)分析。

1.2 變量選擇

根據(jù)農(nóng)產(chǎn)品價格形成理論，農(nóng)產(chǎn)品價格是由供求關(guān)系及價值規(guī)律共同作用的結(jié)果，其影響因素包含供給因素、需求因素及外部沖擊[18]。結(jié)合相關(guān)文獻的研究成果，考慮到數(shù)據(jù)可獲得性，本文選取藍狐皮價格(Y)為被解釋變量。解釋變量包含供給因素、需求因素和外部沖擊，其中供給因素包括藍狐養(yǎng)殖成本(X1)，用漁業(yè)產(chǎn)品生產(chǎn)價格指數(shù)(藍狐的主要日糧是海雜魚)表示；藍狐皮加工成本(X2)，用皮革和毛皮產(chǎn)業(yè)固定總投資表示。需求因素為消費者需求能力(X3)，用高收入群體的收入水平表示。外部沖擊包括國際毛皮市場沖擊(X4)，用毛皮和人造毛皮及其制品進出口額之比表示；藍狐皮替代品價格(X5)，用同為長毛裘皮的貉子(貉Nyctereutesprocyonoides)皮價格表示。

6組數(shù)據(jù)的時間段是1995～2015年，其中藍狐皮價格和貉子皮價格的數(shù)據(jù)來自中國皮毛信息網(wǎng)和相關(guān)文獻整理；毛皮和人造毛皮及其制品進出口額之比通過《中國統(tǒng)計年鑒》中毛皮和人造毛皮及其制品進口額和毛皮和人造毛皮及其制品出口額計算得到；漁業(yè)產(chǎn)品生產(chǎn)價格指數(shù)、高收入群體的收入水平、皮革和毛皮產(chǎn)業(yè)固定總投資3個變量的數(shù)據(jù)來自《中國統(tǒng)計年鑒》和國家統(tǒng)計局。

2 實證分析

2.1 藍狐皮價格波動特征分析

通過H-P濾波法，得到1995～2015年藍狐皮價格波動結(jié)果圖1，橫軸表示年份，左縱軸表示藍狐皮價格絕對波動，右縱軸表示藍狐皮價格及趨勢值。按照“波谷-波谷”的周期劃分方法，1995～2015年間中國藍狐皮價格大體經(jīng)歷了5次波動，1995～1998年為第1個周期，1998～2003年為第2個周期，2003～2007年為第3個周期，2007～2012年為第4個周期，2012年至今為第5個周期。

從總體上看藍狐皮價格波動呈現(xiàn)以下特征：一是藍狐皮價格波動頻繁，周期較短。藍狐皮價格波動的5個周期中，波長分別是4年、6年、5年、6年和4年，均屬于中短周期，顯然藍狐皮價格波動的影響因素較多，供需因素和外部沖擊都對藍狐皮價格產(chǎn)生牽制。二是藍狐皮價格擴張期長，收縮期短。在前4個完整波動周期中，擴張期分別為1年、3年、3年和3年，平均擴張期為2.5年；收縮期分別為2年、2年、1年和2年，平均收縮期為1.75年。顯然藍狐皮價格的下跌非常迅速，而恢復(fù)卻需要更長的時間，表明藍狐皮具有畜牧業(yè)產(chǎn)品養(yǎng)殖周期長的特點，市場反應(yīng)具有一定滯后性。三是藍狐皮價格波動幅度較大。藍狐皮價格波動的5個周期中，波幅分別是500、380、220、335和277，價格起伏十分明顯。1995～2015年間，藍狐皮價格最低為200元/張，最高為800元/張，上漲幅度高達300%，顯然藍狐皮價格對各影響因素的變動非常敏感。四是樣本期內(nèi)藍狐皮價格呈現(xiàn)先下跌后上漲的趨勢。根據(jù)H-P濾波法分解的趨勢成分可以直觀看出，雖然部分年份藍狐皮價格存在一定的波動，但整體上呈現(xiàn)先下跌后上漲的趨勢。

圖1 H-P濾波分析結(jié)果Fig.1 Results of H-P filter processing method

2.2 藍狐皮價格波動影響因素分析

為了降低異方差和共線性問題，本文對數(shù)據(jù)進行取對數(shù)處理，用LY、LX1、LX2、LX3、LX4、LX5分別代表藍狐皮價格、藍狐養(yǎng)殖成本、藍狐皮加工成本、消費者需求能力、國際毛皮市場沖擊、藍狐皮替代品價格的對數(shù)。

2.2.1 變量平穩(wěn)性檢驗及協(xié)整檢驗

首先采用ADF檢驗方法對序列進行單位根平穩(wěn)性檢驗，結(jié)果如表1所示，在10%顯著性水平下，樣本期6個變量均存在單位根，為非平穩(wěn)序列，而各變量的一階差分是平穩(wěn)序列。因此，可以判定樣本期6個變量均為非平穩(wěn)的一階單整序列，這6個一階單整序列滿足協(xié)整檢驗平穩(wěn)性的要求。

在單位根檢驗的基礎(chǔ)上，進一步采用Johansen檢驗對6個變量進行協(xié)整檢驗，結(jié)果如表2所示。跡檢驗和最大特征根檢驗都表明在5%的顯著性水平下，無法拒絕“協(xié)整向量個數(shù)為2”的原假設(shè)，證明6個變量之間存在長期均衡關(guān)系。

表1 單位根檢驗結(jié)果

Tab.1 Results of the unit root test

注：(1)DLY、DLX1、DLX2、DLX3、DLX4、DLX5為變量對數(shù)的一階差分；(2)最優(yōu)滯后階數(shù)根據(jù)SIC信息準(zhǔn)則選取

Note：(1)DLY、DLX1、DLX2、DLX3、DLX4andDLX5are the first order difference of the logarithm of variables.(2)The optimal delay order is selected according to SIC information criterion

表2 協(xié)整檢驗結(jié)果

Tab.2 Results of cointegration test

2.2.2 格蘭杰因果分析

協(xié)整檢驗結(jié)果只能解釋變量之間是否存在長期的均衡關(guān)系，但是這種關(guān)系是否構(gòu)成因果關(guān)系還需要進一步驗證，格蘭杰因果檢驗是考察變量間因果關(guān)系較為便捷的工具，對變量做格蘭杰因果檢驗，結(jié)果見表3。觀察各項假設(shè)的P值，若P<0.05，則在5%的顯著性水平下拒絕原假設(shè)，由此可知，藍狐皮加工成本、消費者需求能力、國際毛皮市場沖擊3個變量是藍狐皮價格的格蘭杰原因，而藍狐皮養(yǎng)殖成本和藍狐皮替代品價格2個變量不是藍狐皮價格的格蘭杰原因。

藍狐養(yǎng)殖成本不是藍狐皮價格的格蘭杰原因，從毛皮動物養(yǎng)殖環(huán)節(jié)來看，這反映了海雜魚等生鮮日糧運輸、貯藏成本居高不下，而毛皮動物飼料業(yè)迅速發(fā)展，養(yǎng)殖者減少了生鮮日糧，增加了飼料的使用，漁業(yè)產(chǎn)品價格指數(shù)所代表的藍狐養(yǎng)殖成本影響逐漸減小。從整個毛皮產(chǎn)業(yè)來看，說明目前制約毛皮產(chǎn)業(yè)發(fā)展的主要原因是加工、經(jīng)營等環(huán)節(jié)，藍狐皮加工成本是藍狐皮價格的格蘭杰原因也證明了這一點，藍狐皮加工環(huán)節(jié)對藍狐皮價格的影響較大。貉子皮價格不是藍狐皮價格的格蘭杰原因，反映了雖然貉子皮和藍狐皮同屬長毛裘皮，但貉子皮大多數(shù)是中低端毛皮，使用較為大眾，而藍狐皮大多數(shù)是中高端毛皮，屬于奢侈品消費，二者消費者有所差異，以高收入群體的收入水平所表示的消費者需求能力是藍狐皮價格的格蘭杰原因也證明了這一點，藍狐皮的消費者主要是收入較高的群體，說明毛皮產(chǎn)業(yè)內(nèi)部需要協(xié)調(diào)發(fā)展，根據(jù)不同種類毛皮的特點，拓寬銷路，積極開發(fā)多元市場。國際毛皮市場沖擊是藍狐皮價格的格蘭杰原因，反映了毛皮市場極強的外向型特點，說明中國毛皮銷售在重視國內(nèi)貿(mào)易市場的同時，國際市場的沖擊也不容忽視。

2.2.3 脈沖響應(yīng)分析及預(yù)測方差分解

脈沖響應(yīng)分析描述的是一個內(nèi)生變量的沖擊給其他內(nèi)生變量所帶來的影響，可以在一定程度上反映變量之間的影響程度，圖2為藍狐皮價格的脈沖響應(yīng)結(jié)果，其中，橫軸代表沖擊作用滯后期數(shù)(單位：年)，縱軸代表脈沖響應(yīng)水平。圖2中第1行3個圖依次為藍狐養(yǎng)殖成本、藍狐皮加工成本和消費者需求能力對藍狐皮價格波動的影響，第2行3個圖依次為國際毛皮市場沖擊、藍狐皮替代品價格和藍狐皮價格自身對藍狐皮價格波動的影響。可以明顯看出6個圖都出現(xiàn)正負交替的波動現(xiàn)象，其中藍狐皮價格自身的影響最大，其次是代表供給因素的藍狐養(yǎng)殖成本和藍狐皮加工成本對藍狐皮價格波動的沖擊幅度，并且幅度隨時間的延長逐漸減?。唤又谴硇枨笠蛩氐南M者需求能力對藍狐皮價格的沖擊幅度，也隨時間的延長逐漸減??；最后是代表外部沖擊的國際毛皮市場沖擊和藍狐皮替代品價格對藍狐皮價格的沖擊幅度一直較小，隨時間的延長逐漸趨于平緩。

表3 格蘭杰因果關(guān)系檢驗結(jié)果

Tab.3 Results of granger causality test

圖2 脈沖響應(yīng)結(jié)果分析Fig.2 Results of the impulse response

為進一步分析各變量和藍狐皮價格之間的影響程度和動態(tài)影響機制，對藍狐皮價格進行預(yù)測方差分解，由表4結(jié)果可以看出，在預(yù)測的第1期，藍狐皮價格波動不能由除了自身以外的其余5個變量解釋，第2期開始，藍狐皮價格自身引起的波動程度逐漸降低，從第1期的100%逐漸下降至第30期的41.19%，同時，隨著預(yù)測期的延長，其他5個變量對藍狐皮價格波動的影響程度逐漸增加，進一步印證了協(xié)整檢驗的結(jié)果，藍狐皮價格與其他5個變量存在長期關(guān)系。到第30期時，藍狐皮價格波動自身的影響程度雖有所下降，但仍然在6個變量中影響程度最大，占到41.19%，依次是藍狐養(yǎng)殖成本的25.99%，藍狐皮加工成本的24.00%，消費者需求能力的6.28%，藍狐皮替代品價格的1.40%，國際毛皮市場沖擊的1.13%，比例大小與脈沖響應(yīng)圖的結(jié)果相一致。

表4 預(yù)測方差分解結(jié)果

Tab.4 Results of forecast variance decomposition

3 討論

本文采用H-P濾波法對藍狐皮價格波動進行趨勢分析和波動分析，并采用向量誤差修正模型(VEC)就各影響因素對藍狐皮價格波動的影響程度及動態(tài)影響機制進行研究，以期為目前大多數(shù)通過對毛皮價格定性和靜態(tài)分析所得出的結(jié)論做補充。其結(jié)果表明藍狐皮價格波動頻繁、幅度較大、周期較短，在波動周期內(nèi)擴張期長、收縮期短，總體呈現(xiàn)先下跌后上漲的趨勢。市場供需因素和外部沖擊都對藍狐皮價格波動產(chǎn)生影響，但各類因素的影響程度不同，長期來看，代表供給因素的藍狐養(yǎng)殖成本和藍狐皮加工成本對藍狐皮價格產(chǎn)生了較大的沖擊，代表需求因素的消費者需求能力對藍狐皮價格波動的影響程度較?。徊㈦S著預(yù)測期的延長呈現(xiàn)小幅度的增長趨勢，代表外部沖擊的國際毛皮市場沖擊和藍狐皮替代品價格對藍狐皮價格波動的影響程度較小。并且除了藍狐皮價格自身的影響，其他因素對藍狐皮價格波動的影響都具有滯后期。

根據(jù)以上發(fā)現(xiàn)可知，我國毛皮價格波動劇烈，穩(wěn)定性較弱，產(chǎn)業(yè)發(fā)展還不成熟，產(chǎn)業(yè)鏈各環(huán)節(jié)都亟待改善。為了穩(wěn)定毛皮價格、促進毛皮產(chǎn)業(yè)發(fā)展，本文結(jié)合研究結(jié)果提出以下建議：一是發(fā)展毛皮動物規(guī)?；B(yǎng)殖，實行統(tǒng)一化、標(biāo)準(zhǔn)化的生產(chǎn)以降低養(yǎng)殖成本；二是提升毛皮加工水平，提高毛皮質(zhì)量，制定優(yōu)質(zhì)毛皮商標(biāo)，增強我國毛皮市場競爭力；三是組建毛皮拍賣行，設(shè)立毛皮質(zhì)量標(biāo)準(zhǔn)，適度學(xué)習(xí)國際拍賣行的先進經(jīng)驗，優(yōu)化我國毛皮貿(mào)易環(huán)境；四是建立行業(yè)組織，及時掌握信息不對稱的市場行情，宏觀調(diào)控毛皮產(chǎn)業(yè)各環(huán)節(jié)的供求關(guān)系，科學(xué)制定毛皮產(chǎn)業(yè)發(fā)展目標(biāo)，提高從業(yè)人員收入，增強我國毛皮貿(mào)易的總實力。