付蓮蓮,馮家璇,趙一恒
(江西農(nóng)業(yè)大學(xué)計算機(jī)與信息工程學(xué)院,南昌 330045)
近期,非洲豬瘟喋喋不休,供給、消費(fèi)和疫情等因素交織,使得生豬價格波動更為復(fù)雜,價格空間聯(lián)動性及差異加大、網(wǎng)絡(luò)特征更為顯著。來自經(jīng)濟(jì)環(huán)境、宏觀政策、動物疫病等不確定性的響應(yīng)會改變生豬生產(chǎn)者和消費(fèi)者的反應(yīng),引起不必要的市場失靈[1],價格波動呈現(xiàn)復(fù)雜性。在價格波動特征方面,國內(nèi)外學(xué)者主要基于變異系數(shù)、HP濾波分解、正弦周期序列、ARCH模型、Hilbert-Huang變換等分析工具研究了各國家、各區(qū)域的生豬價格波動周期[2-4]。近期,有學(xué)者分析了生豬價格波動的非線性、非對稱和機(jī)構(gòu)突變性,生豬價格月度時間序列發(fā)生了結(jié)構(gòu)突變,價格波動的機(jī)理具有異質(zhì)性[5-6]。國外學(xué)者Tweeten(1972)最早提出農(nóng)產(chǎn)品價格非對稱傳導(dǎo)理論,之后學(xué)者們研究了各種食品價格的傳導(dǎo)[7-8]、畜產(chǎn)品價格的傳導(dǎo)機(jī)制[9-10]。目前中國國內(nèi)學(xué)術(shù)界對生豬產(chǎn)業(yè)的研究中探究其縱向傳導(dǎo)是主要,其方法有以下:一是非對稱誤差修正模型[11-14];二是VAR或SVAR模型[15];三是GARCH模型[16,17];四是基于協(xié)整、格蘭杰檢驗、同步系數(shù)研究省域價格的空間傳導(dǎo)[11,18]。
已有研究主要運(yùn)用HP濾波分解、季節(jié)調(diào)整變異系數(shù)和ARCH模型族分析了價格周期、非對稱等時間序列特征,這些是從時間維度描述價格的波動特征,從統(tǒng)計學(xué)角度對價格的時間周期和聚集特征進(jìn)行了刻畫。但許多價格波動的狀態(tài)具有網(wǎng)絡(luò)、模態(tài)、混沌等內(nèi)部動力學(xué)特征[19],價格的升降更多地呈現(xiàn)內(nèi)部網(wǎng)絡(luò)特征,已有研究方法對價格的網(wǎng)絡(luò)特征無法處理,而復(fù)雜網(wǎng)絡(luò)具有小世界、無標(biāo)度特性,對常規(guī)的沖擊具有魯棒性。通過對價格進(jìn)行粗粒化處理,轉(zhuǎn)變成字符串,根據(jù)字符串的走向繪制節(jié)點(diǎn)和邊,基于圖論建立生豬價格的復(fù)雜網(wǎng)絡(luò),通過分析網(wǎng)絡(luò)節(jié)點(diǎn)的強(qiáng)度、加權(quán)集聚系數(shù)、模態(tài)特征、核心模態(tài)特征,明晰價格波動的網(wǎng)絡(luò)特征和不同模態(tài)之間的傳導(dǎo)。有學(xué)者已將它成功應(yīng)用到原油價格與股票價格、房地產(chǎn)、外匯等金融市場價格的波動及其傳導(dǎo)方面[20-21]。為此,本文對生豬價格波動網(wǎng)絡(luò)中生豬價格波動模態(tài)統(tǒng)計特征、生豬價格波動核心模態(tài)分析以及生豬價格波動模態(tài)傳導(dǎo)等特征進(jìn)行分析,探討生豬價格在不同時間的波動規(guī)律,找出價格上升或下降的時間段,有利于生豬市場在非洲豬瘟等不良沖擊下迅速、有效地調(diào)整回正常狀態(tài),將不確定因素對生豬市場的負(fù)面沖擊減少到最小。
運(yùn)用復(fù)雜網(wǎng)絡(luò)的方法對生豬價格波動轉(zhuǎn)化后的模態(tài)進(jìn)行多方面分析,尋找其波動特性。首先把原始月度數(shù)據(jù)按照“3σ原則”進(jìn)行粗?;幚?,把生豬價格的漲跌數(shù)據(jù)轉(zhuǎn)化為符號序列并構(gòu)建有向加權(quán)復(fù)雜網(wǎng)絡(luò)。通過尋找其與時間序列相關(guān)的動力學(xué)統(tǒng)計特征量來分析生豬價格波動復(fù)雜網(wǎng)絡(luò)的內(nèi)在規(guī)律性。
選取待宰活豬價格原始數(shù)據(jù)(元/千克),根據(jù)統(tǒng)計學(xué)的方法把原始數(shù)據(jù)轉(zhuǎn)化為環(huán)比價格指數(shù),繪制出價格波動圖(見圖1),數(shù)據(jù)來源于中國畜牧信息網(wǎng),數(shù)據(jù)區(qū)間為:2000年1月至2018年9月。
圖1 生豬價格的波動趨勢
(1)
表1 生豬價格數(shù)據(jù)粗粒化處理過程
生豬價格波動序列本該有125(53)種模態(tài),但經(jīng)過統(tǒng)計發(fā)現(xiàn)生豬價格波動的模態(tài)共出現(xiàn)34種,其余91種模態(tài)并未出現(xiàn)。出現(xiàn)的34種模態(tài)可初步得出連續(xù)3個月中生豬價格波動較為平緩,連續(xù)3個月中沒有出現(xiàn)異常上漲“B”和異常下跌“b”同時發(fā)生的現(xiàn)象,部分模態(tài)示意圖如下。
圖2 價格網(wǎng)絡(luò)的部分模態(tài)圖
統(tǒng)計發(fā)現(xiàn)29組模態(tài)間轉(zhuǎn)換次數(shù)均為1次,模態(tài)之間轉(zhuǎn)換次數(shù)最多的是從aaA模態(tài)轉(zhuǎn)化到aAA模態(tài),轉(zhuǎn)換了15次,其次有3組模態(tài)轉(zhuǎn)換均轉(zhuǎn)換了13次。以每個模態(tài)作為圖的頂點(diǎn),模態(tài)之間的轉(zhuǎn)換作為邊,構(gòu)建一個復(fù)雜網(wǎng)絡(luò)來描述生豬價格波動的狀態(tài)(見圖3)。
Edge={開始模態(tài),結(jié)束模態(tài),權(quán)重(轉(zhuǎn)換次數(shù))}
圖3中連線粗細(xì)代表了兩節(jié)點(diǎn)(即兩模態(tài))之間轉(zhuǎn)化次數(shù)的多少。出度最大的6個節(jié)點(diǎn)已標(biāo)黑,關(guān)聯(lián)程度較大的6種模態(tài)形成一個閉環(huán),即aaA>aAA>AAA>AAa>Aaa>aaa>aaA,可初步推斷生豬價格波動是個漸進(jìn)的傳導(dǎo)過程。
在有向復(fù)雜網(wǎng)絡(luò)中,加權(quán)度值越大即權(quán)值越大說明在復(fù)雜網(wǎng)絡(luò)中的重要性越大,其中加權(quán)出度的分布計算公式如下:
(2)
WIi代表節(jié)點(diǎn)i的加權(quán)出度,N代表所有節(jié)點(diǎn)加權(quán)出度值的和,Wij表示節(jié)點(diǎn)到節(jié)點(diǎn)的距離。根據(jù)式(2)計算出價格波動復(fù)雜網(wǎng)絡(luò)模態(tài)的加權(quán)出度,如表2所示。
表2 34種模態(tài)的加權(quán)出度
從表2發(fā)現(xiàn)AAA,aaa,Aaa,aaA,AAa,aAA六個節(jié)點(diǎn)的加權(quán)出度分布的和高達(dá)79.26%,整個生豬價格波動復(fù)雜網(wǎng)絡(luò)具有較強(qiáng)的傳導(dǎo)作用,即生豬價格波動模態(tài)中這6種模態(tài)向其它模態(tài)轉(zhuǎn)化較多。生豬價格為連續(xù)3個月正常上漲的模態(tài)(AAA)的次數(shù)最多,由此模態(tài)共轉(zhuǎn)出50次,其次是連續(xù)3個月正常下跌的模態(tài)(aaa)的次數(shù),由此模態(tài)共轉(zhuǎn)出46次,6個節(jié)點(diǎn)中剩下也是“正常上漲”和“正常下跌”相組合,所有出現(xiàn)模態(tài)中,代表“異常上漲”和“異常下跌”的字符“B”和“b”沒有出現(xiàn)在主要的6個節(jié)點(diǎn)中,同時在所有出現(xiàn)的34種模態(tài)中出現(xiàn)的很少,不存在“B”和“b”即“異常上漲”和“異常下跌”同時出現(xiàn)在同一個模態(tài)的極端狀況,表明3個月的時間段內(nèi)生豬價格波動較為平穩(wěn),但34種模態(tài)中并未出現(xiàn)兩個及兩個以上的“持平”狀態(tài),說明在3個月的時間段內(nèi)生豬價格在持續(xù)波動。
在統(tǒng)計粗?;蟮姆栃蛄小癆,a,B,b,o”(分別代表“正常漲跌”,“異常漲跌”和“持平”)中,代表“正常上漲”的符號“A”出現(xiàn)了103次,頻率最高。代表“正常下跌”的符號“a”出現(xiàn)了102次,剩下的“B,b,o”分別出現(xiàn)了9次,4次和2次,即“異常漲跌”和“持平”出現(xiàn)的頻率很低。表2中連續(xù)正常上升模態(tài)(AAA)轉(zhuǎn)出50次,連續(xù)3個月正常下跌的模態(tài)(aaa)轉(zhuǎn)出46次,而其余4個主要的轉(zhuǎn)出模態(tài)中均是正常范圍的漲跌,且上漲的正常上漲的趨勢較強(qiáng)。綜上可知價格以不同的頻率呈現(xiàn)34種模態(tài),但以“AAA,aaa,Aaa,aaA,AAa,aAA”這6種模態(tài)為主,其中價格以連續(xù)3次呈同樣的狀態(tài)為主,除去意外疫情等,短期內(nèi)價格可以預(yù)測。價格大多在正常范圍內(nèi)波動,不太會出現(xiàn)暴漲暴跌的狀態(tài),長期來看價格呈上漲趨勢。
網(wǎng)絡(luò)集聚系數(shù)反映該節(jié)點(diǎn)與相鄰節(jié)點(diǎn)之間關(guān)系的緊密程度,其計算公式如下:
(3)
其中,ki、wi、wij表示節(jié)點(diǎn)的度、節(jié)點(diǎn)加權(quán)度、權(quán)重,aijajkaki表示i、j、k能否構(gòu)成三角形。
CCi為加權(quán)集聚系數(shù),CCi越大,說明該節(jié)點(diǎn)的影響程度越大。根據(jù)式(3)求得34個節(jié)點(diǎn)的CCi如表3。
表3 網(wǎng)絡(luò)節(jié)點(diǎn)加權(quán)集聚系數(shù)
由表3可知網(wǎng)絡(luò)中34個節(jié)點(diǎn)只有AAA,aAa,aAA,AaA,aaa,Aaa,aaA,AAa的加權(quán)集聚系數(shù)非零,其順序為:aAa>aAA>AAA>AaA>aaa>Aaa>aaA>AAa,即此8個節(jié)點(diǎn)的近鄰節(jié)點(diǎn)之間可以形成三角形以上的關(guān)聯(lián),為尋找主要節(jié)點(diǎn)加權(quán)集聚系數(shù)與主要節(jié)點(diǎn)加權(quán)出度的關(guān)聯(lián)性,將這8個節(jié)點(diǎn)的加權(quán)集聚系數(shù)與6個主要節(jié)點(diǎn)加權(quán)出度結(jié)合繪制條形圖。
圖4 加權(quán)出度和加權(quán)集聚系數(shù)
從圖4可判斷加權(quán)集聚系數(shù)和加權(quán)出度無強(qiáng)相關(guān)性,AAA模態(tài)和aaa模態(tài)的加權(quán)集聚系數(shù)和加權(quán)出度值均較大,說明當(dāng)生豬價格連續(xù)3個月“正常上漲-正常上漲-正常上漲”或“正常下跌-正常下跌-正常下跌”時,市場的生豬價格會進(jìn)行整頓,模態(tài)轉(zhuǎn)換的概率較大即連漲和連跌轉(zhuǎn)化為其它模態(tài);另外,aAa和AaA模態(tài)都是加權(quán)出度較低但加權(quán)集聚系數(shù)較大,說明在生豬市場中生豬價格上下連續(xù)波動的頻率不高,但市場經(jīng)過“正常下跌-正常上漲-正常下跌”或“正常上漲-正常下跌-正常上漲”的模態(tài)之后轉(zhuǎn)換為其它模態(tài)的概率較大。
為分析價格網(wǎng)絡(luò)中模態(tài)間的傳導(dǎo)特征,引入k-核的概念。刪除度小于k的節(jié)點(diǎn)以及相應(yīng)的邊,余下子圖的全部頂點(diǎn)定義為高核,高核內(nèi)節(jié)點(diǎn)具有較強(qiáng)的連通性和傳導(dǎo)性,以下對生豬價格波動復(fù)雜網(wǎng)絡(luò)的k-核節(jié)點(diǎn)(k=3)進(jìn)行分析,k核解析后的結(jié)果如下圖所示。
將生豬價格波動復(fù)雜網(wǎng)絡(luò)進(jìn)行k-核解析后,共剩6個高核內(nèi)節(jié)點(diǎn),其中Aaa指向aaA節(jié)點(diǎn)的邊權(quán)重7和aAA指向AAa節(jié)點(diǎn)的邊的權(quán)重5明顯小于其它的邊,進(jìn)而將“Aaa節(jié)點(diǎn)至aaA節(jié)點(diǎn)”和”aAA節(jié)點(diǎn)至AAa節(jié)點(diǎn)”的邊過濾掉,解析后的最終子圖形成一個單向封閉閉環(huán):aaA>aAA>AAA>AAa>Aaa>aaa>aaA,當(dāng)市場生豬價格波動處于此環(huán)路中任意節(jié)點(diǎn)模態(tài)時,有大概率經(jīng)過6個滑動周期的價格波動會回到原該起始節(jié)點(diǎn)。如從節(jié)點(diǎn)aaa代表的模態(tài)作為起始,其本身包括之后6個滑動周期生豬價格波動為:“正常下跌-正常下跌-正常下跌”至“正常下跌-正常下跌-正常上漲”至“正常下跌-正常上漲-正常上漲”至“正常上漲-正常上漲-正常上漲”至“正常上漲-正常上漲-正常下跌”至“正常上漲-正常下跌-正常下跌”再波動回到“正常下跌-正常下跌-正常下跌”從這幾種傳導(dǎo)模態(tài)過程可得,生豬價格具有漸進(jìn)性、持續(xù)性和周期性,持續(xù)的正常漲跌是常態(tài),而異常漲跌和持平在此波動過程中出現(xiàn)概率極小。
圖5 價格復(fù)雜網(wǎng)絡(luò)的k-核解析
對生豬價格網(wǎng)絡(luò)進(jìn)行k-核解析(k=3)后,再將剩余模態(tài)轉(zhuǎn)變時間發(fā)生的頻次歸納統(tǒng)計如表4:
表4 模態(tài)轉(zhuǎn)換發(fā)生時間節(jié)點(diǎn)和頻次
從表4中可以看出,每一種模態(tài)轉(zhuǎn)變發(fā)生的月份頻次有很大不同,aaa模態(tài)至aaA模態(tài)的轉(zhuǎn)變發(fā)生了13次,其中5月份的頻次最高;模態(tài)aaA至aAA模態(tài)轉(zhuǎn)變發(fā)生了14次,其中6月份頻次最高;aAA模態(tài)至AAA模態(tài)轉(zhuǎn)變發(fā)生了11次,其中7月份的頻次最高;AAA模態(tài)至AAa模態(tài)發(fā)生了13次,其中9月份的頻次最高;AAa模態(tài)至AaA模態(tài)則僅在2000年和2001年10月份發(fā)生了2次;AaA模態(tài)至aAA模態(tài)也僅在2000年11月份發(fā)生1次;AAa模態(tài)至Aaa模態(tài)發(fā)生了13次,其中6月份頻次最高;Aaa模態(tài)至aaa模態(tài)發(fā)生了12次,其中2月份的頻次最高。3-核解析后由頻次歸納統(tǒng)計,AAA節(jié)點(diǎn)和aaa節(jié)點(diǎn)具有明顯的自傳導(dǎo)特征,其中轉(zhuǎn)換aaa→aaa共有31次,4月份出現(xiàn)此轉(zhuǎn)換狀態(tài)次數(shù)為6次,出現(xiàn)在3月、5月的次數(shù)均為4次,其中4月份的次數(shù)最大。AAA模態(tài)轉(zhuǎn)換到AAA模態(tài)總共有35次,發(fā)生在8月份的頻次最高,即生豬價格四個月連續(xù)下跌的狀態(tài)發(fā)生過31次,而連續(xù)上漲的次數(shù)發(fā)生過35次,3個月“連續(xù)正常下跌”模態(tài)和3個月“連續(xù)正常上漲”的自傳導(dǎo)現(xiàn)象較多出現(xiàn)在4月和8月,其次“正常下跌-正常下跌-正常上漲”向“正常下跌-正常上漲-正常上漲”模態(tài)轉(zhuǎn)化多數(shù)出現(xiàn)在6月。
文章將生豬價格的原始價格轉(zhuǎn)換為狀態(tài),構(gòu)建有向網(wǎng)絡(luò)圖,分析價格網(wǎng)絡(luò)節(jié)點(diǎn)的加權(quán)出度、核心模態(tài)以及生豬價格波動模態(tài)傳導(dǎo)從而明確生豬價格波動與傳導(dǎo)的動力學(xué)特征,具體結(jié)論如下:
1)AAA,aaa,Aaa,aaA,AAa,aAA六個節(jié)點(diǎn)的加權(quán)出度分布之和高達(dá)79.26%,遠(yuǎn)大于其他節(jié)點(diǎn)。相對符號“B”、“b”和“o”的節(jié)點(diǎn),連續(xù)上漲的模態(tài)AAA和連續(xù)下跌的模態(tài)aaa轉(zhuǎn)出次數(shù)相對較高,且字符“B”和字符“b”在同一節(jié)點(diǎn)中的模態(tài)沒有出現(xiàn),表明2000年1月至2018年9月間生豬價格一般是以正常漲跌進(jìn)行,價格異常上升、異常下降或者持平的狀態(tài)很少。
2)模態(tài)AAA的聚集系數(shù)為0.25,加權(quán)出度為50;aaa的加權(quán)出度為41,聚集系數(shù)為0.2167。一般情況下,當(dāng)生豬價格連續(xù)正常上漲3個月或者連續(xù)下跌3個月后,經(jīng)過市場的調(diào)節(jié)作用,價格的波動狀態(tài)會有所改變。含有“B”、“b”或“o”的模態(tài)加權(quán)出度很低,且加權(quán)集聚系數(shù)都為0,表明價格“異常上漲”和“異常下降”的可能性較小。
3)k-核解析顯示生豬價格波動主要傳導(dǎo)路徑是6個模態(tài)形成的單向閉環(huán)。生豬價格波動有較強(qiáng)的網(wǎng)絡(luò)特性,模態(tài)間的傳導(dǎo)有漸進(jìn)性、持續(xù)性和周期性。
生豬價格網(wǎng)絡(luò)的模態(tài)特征和傳導(dǎo)有一定的規(guī)律性,應(yīng)根據(jù)價格的漲幅發(fā)生在不同月份,加強(qiáng)生豬物聯(lián)網(wǎng)的建設(shè),提前做好價格預(yù)警工作,準(zhǔn)確篩選價格預(yù)警指標(biāo),采用科學(xué)、先進(jìn)的方法確定先行指標(biāo)、同步指標(biāo),優(yōu)化生豬價格預(yù)警體系;盡管研究表明生豬價格仍以正常上漲和正常下降模態(tài)為主,但不能忽視價格的異常波動,近期,疫情等不確定性因素對價格的影響越來越顯著,其對價格的影響要加以關(guān)注;其次,推廣“農(nóng)戶+養(yǎng)殖合作社”模式,加強(qiáng)養(yǎng)殖戶抗風(fēng)險能力。散戶對生豬市場的掌控能力較弱,當(dāng)出現(xiàn)連續(xù)上漲模態(tài)后,盲目跟風(fēng)加大生豬存欄量,導(dǎo)致連續(xù)下跌模態(tài)出現(xiàn),從而導(dǎo)致脫離不了“豬周期”魔咒,建立養(yǎng)殖合作者可以避免資金鏈斷裂,合作共贏。
本文建立了生豬價格轉(zhuǎn)移的復(fù)雜網(wǎng)絡(luò)圖,研究了價格網(wǎng)絡(luò)的模態(tài)特征、核心模態(tài)節(jié)點(diǎn)的加權(quán)集聚系數(shù)、K-核的傳導(dǎo)模態(tài)等特征,雖然有助于從網(wǎng)絡(luò)視角完善價格波動的特征和狀態(tài)間的傳導(dǎo),但由于時間和精力的局限性,對供給、需求和生豬供應(yīng)鏈等因素未納入模型中。同時基于本文得到的模態(tài)對價格進(jìn)行預(yù)測,這將是下一階段的研究方向。