張克榮,劉武藝
(阜陽師范學(xué)院a.經(jīng)濟與商業(yè)學(xué)院;b.科研處,安徽阜陽236041)
價格波動的研究方法及其模型
張克榮a,劉武藝b
(阜陽師范學(xué)院a.經(jīng)濟與商業(yè)學(xué)院;b.科研處,安徽阜陽236041)
在價格理論與實踐中,經(jīng)濟學(xué)家們試圖建立嚴(yán)密的數(shù)學(xué)模型來研究和分析時間序列經(jīng)濟周期和誤差修正等價格波動問題。文章參考有關(guān)時間序列分析方面的研究成果,系統(tǒng)闡釋了應(yīng)用于價格波動相關(guān)的計算方法及其模型。
價格;波動;研究方法;模型
其基本原理即算術(shù)平均,包括簡單移動平均、加權(quán)移動平均、項和項移動平均,對稱的亨德森移動平均、PA(階段平均)等方法。該方法直接采用時間序列的移動平均值來代表經(jīng)濟序列的長期趨勢,優(yōu)點是計算簡便、方法客觀,適用于長期趨勢較為復(fù)雜且隨機波動很大的時間序列數(shù)據(jù)的處理;同時也便于不同時間序列波動幅度大小變化的比較研究。
簡單移動平均法采用的方法是取一定數(shù)量時期的數(shù)據(jù)平均,按時間順序逐次推進,每推進一次,就舍去前一個數(shù)據(jù),同時增加一個后續(xù)相鄰的數(shù)據(jù),再進行平均,依次類推,最后形成一個新的序列。若原時間序列沒有明顯的不穩(wěn)定變動的話,則可用最近一次移動平均數(shù)作為下一個時期預(yù)測值。此方法的特點是只能用于近期預(yù)測,即只能對于后續(xù)相鄰的那一項預(yù)測,而且也僅適用于時間序列變化比較平穩(wěn)的近期預(yù)測。
季節(jié)模型是反映具有季節(jié)變動規(guī)律的時間序列模型。季節(jié)變動指以一年為一個周期的變化,引起季節(jié)變動的首要因素是四季更迭,適合季節(jié)分析的數(shù)據(jù)是季度數(shù)據(jù)、月度數(shù)據(jù)、周數(shù)據(jù)、日數(shù)據(jù)等等。季節(jié)調(diào)整方法是從時間序列觀測值中消除季節(jié)變動的一種方法。隨時間變動的經(jīng)濟變量的影響因素通常歸為四類:第一類是長期變動趨勢(T),它是由各個時期普遍的、持續(xù)的、決定性的基本因素的作用,使指標(biāo)在長時期內(nèi)沿著逐漸增加、逐漸減少或平穩(wěn)的趨勢演變。第二類是季節(jié)變動(S),指在年內(nèi)的周期變動,并且這種變動每年重復(fù)出現(xiàn)。第三類是循環(huán)變動(C),又稱為周期波動。指時間序列在為期較長的時間內(nèi)起伏的變動。反映的是呈漲落交替的波浪式起伏變動。第四類是不規(guī)則變動(I),又稱隨機變動。
四種變動與原序列的關(guān)系可以被概括為兩種模型:
乘法模型:
加法模型:
其中,乘法模型適于T、S、C相關(guān)的情形,比如季節(jié)變動的幅度隨趨勢上升而增大。加法模型適于T、S、C相互獨立的情形。模型中的S部分也被稱為季節(jié)因子(seasonal factor),反映序列隨時間變化過程中,受季節(jié)因素影響的程度。其表現(xiàn)形式因模型類型不同而異。如果從乘法模型出發(fā),季節(jié)因子以季節(jié)指數(shù)(seasonal index)形式出現(xiàn);如果從加法模型出發(fā),季節(jié)因子以季節(jié)變差形式出現(xiàn)。但這兩種模型只是形式上的不同,因為對乘數(shù)模型取對數(shù),就成為加法模型,即:
對原序列進行季節(jié)調(diào)整,就是將季節(jié)變動從序列中去除,以消除序列里存在季節(jié)波動帶來的干擾?;舅悸肥牵?/p>
乘法模型:
加法模型:
在社會經(jīng)濟統(tǒng)計中,主要采用乘法模型。
在季節(jié)調(diào)整的方法中,目前主要有四種:美國商務(wù)部人口普查局研究開發(fā)的X11方法,在X11上改進的CensusX12方法,移動平均法和Tramo/Seats法。
采用計量經(jīng)濟模型,先作散點圖,然后選擇能擬合原序列的曲線方程?;貧w分析是計量經(jīng)濟分析中使用最多的方法,主要用來分析兩個及以上的變量之間因果關(guān)系。研究經(jīng)濟變量的變化影響因素及作用程度。定義隨機變量:
式中y是因變量,亦稱被解釋變量;x1,x2,…,xp是自變量,亦稱解釋變量;f(x1,x2,…,xp)是回歸函數(shù);ε是隨機誤差,表示受隨機因素影響而未能觀察到的偶然因素。y由自變量和隨機誤差共同決定。根據(jù)模型表達式是否是線性,可將回歸模型分為線性回歸模型和非線性回歸模型;根據(jù)模型中自變量的個數(shù),還可以分為一元回歸模型和多元回歸模型。
經(jīng)典線性回歸模型有三個基本假設(shè):①自變量是確定性變量且彼此不相關(guān),即cov(xi,xj)=0(i≠j);②隨機誤差項服從相互獨立且期望為零、標(biāo)準(zhǔn)差為σ的正態(tài)分布;③樣本容量個數(shù)多于參數(shù)個數(shù)即n>p+1。
回歸函數(shù)為線性函數(shù)的模型是線性回歸模型,其一般形式為:
如(7)式中只含有一個自變量,即為一元線性回歸模型。現(xiàn)代線性回歸分析是由一組探求變量之間關(guān)系的技術(shù)組成,其中心主題是建立模型、評價擬合精度及可靠性,并做出結(jié)論。
消除趨勢法更多地借用濾波技術(shù),主要包括HP濾波法,是Hodrick,Prescott在1980年和Kydland,Prescott在1992年提出的一種分離經(jīng)濟波動、分析經(jīng)濟波動特征的方法,該方法將經(jīng)濟活動分解為趨勢因素和周期因素,剔除趨勢因素,以周期因素作為研究經(jīng)濟波動特性的依據(jù)。此外還有Band-Pass濾波法,是Baxter和King在1995年提出,也被稱為BK濾波法。主要是利用了時間序列時域分析和頻域分析,其基本做法是利用移動平均濾波將經(jīng)濟數(shù)據(jù)中的周期成分分離出來。另外還有狀態(tài)空間模型下的Kalman濾波法。
目前比較常用的是HP濾波法(Hodrick和Prescott,1990)。
HP濾波法假定Yt是一個包含趨勢和周期波動的時間序列,是趨勢波動因素,是周期波動因素,則:
計算HP濾波是要從Yt中分離出,最小化:
HP濾波因為建立在對經(jīng)濟變量較為合理的基礎(chǔ)上而得到廣泛應(yīng)用。從統(tǒng)計學(xué)的觀點來看,λ的選取必須是隨意的,因為任何一個非平衡的時間序列都可以分解為無數(shù)個非平衡趨勢因素與平穩(wěn)周期因素的組合。不同的λ值決定了不同的周期方式和平滑度,λ=0時,滿足最小化問題的趨勢等于序列yt,λ增加時,估計趨勢中的變化總數(shù)相對于序列中的變化減少,即λ越大,估計趨勢越光滑;λ趨于無窮大時,估計趨勢將接近線性函數(shù)。目前對于月度數(shù)據(jù),一般而言λ是使用14400(Hodrick和Prescott,1998);對于季度數(shù)據(jù),學(xué)者們一般認(rèn)同λ=1600;但如果使用年度數(shù)據(jù)分析,Backus和Kehoe(1992)認(rèn)為λ應(yīng)為100,而OECD(經(jīng)合組織)則認(rèn)為25比較合適。
食品價格的形成不是獨立的,而是在復(fù)雜的社會經(jīng)濟環(huán)境中各因素互相作用下形成并波動的,食品價格與很多因素有著相互儲存的因果關(guān)系和相關(guān)性。多變量的價格因果模型是利用各因素之間的相關(guān)性,通過對相關(guān)因素的分析,找出受這些變量影響的食品價格發(fā)展的規(guī)律性,從而進行分析。這類模型主要有幾何分布滯后模型、多項式分布滯后模型、自回歸分布滯后模型等動態(tài)時間序列模型,結(jié)構(gòu)化與非結(jié)構(gòu)化的向量自回歸模型和誤差修正模型,條件異方差模型等。在這類模型中都首先涉及單位根檢驗以及協(xié)整檢驗等,以便進行模型識別和確定。
由于很多時間序列都是非平穩(wěn)的,而時間序列模型主要是針對平穩(wěn)時間序列的建模方法,因此,先要對時間序列進行平穩(wěn)性檢驗,這就是檢驗時間序列是否具有單位根。
如果一個時間序列的均值或自協(xié)方差函數(shù)隨時間而改變,那么這個序列就是非平穩(wěn)的時間序列。隨機過程{yt,t=1,2,…},若
式中,ρ=1,εt為一個穩(wěn)定過程,且E(εt)=0,cov(εt,εt-s)=μt<∞,這里s=0,1,2,…,則稱該過程為單位根過程。單位根過程經(jīng)過一階差分成為平穩(wěn)過程,即
這樣的時間序列yt稱為一階單整序列,記作I(1)。一般地,如果非平穩(wěn)時間序列xt經(jīng)過d次差分達到平穩(wěn),則稱其為d階單整序列,記作I(d)。其中,d表示單整階數(shù),是序列包含的單位根個數(shù)。
ADF檢驗?zāi)P陀腥N設(shè)定模式,選擇正確的設(shè)定模式是十分重要的。ADF檢驗法是對DF檢驗法的改進,在DF檢驗中,考慮一個AR(1)過程
其中,εt是白噪聲。若參數(shù)-1<ρ<1,則序列yt是平穩(wěn)的。而當(dāng)ρ>1或ρ<1,序列是爆炸性的,沒有實際意義。所以,只需檢驗ρ是否在(-1,1)區(qū)間內(nèi)。在實際檢驗時,將(12)式寫為:△yt=γyt-1+εt,其中,γ=ρ-1。檢驗假設(shè)為:
H0∶γ=0,h1∶γ<0
在DF檢驗中,常常因為序列存在高階滯后相關(guān)而破壞了隨機擾動項是白噪聲的假設(shè),ADF檢驗對此做了改進。它假定序列yt服從AR(p)過程。檢驗方程為:
協(xié)整檢驗是Johansen在1988年及在1990年與Juselius一起提出的一種以VAR模型為基礎(chǔ)的檢驗回歸系數(shù)的方法。Johansen協(xié)整檢驗的基本思想和具體研究方法如下。
首先建立一個VAR(p)模型
其中y1t,y2t,…,ykt都是非平穩(wěn)的I(1)變量;xt是一個確定的d維的外生向量,代表趨勢項、常數(shù)項等確定性項;εt是k維擾動向量。將式(14)經(jīng)過差分變換后,得到式(15)。
由于I(1)過程經(jīng)過差分變換將變成I(0)過程,即式(13)中△yt,△yt-j(j=1,2,…,p)都是I(0)變量構(gòu)成的向量,那么只要ΠYt-1是I(0)的向量,即y1,t-1,y2,t-1,…,yk,t-1之間具有協(xié)整關(guān)系,就能保證△yt是平穩(wěn)過程。變量y1,t-1,y2,t-1,…,yk,t-1之間是否具有協(xié)整關(guān)系主要依賴于矩陣Π的秩。將yt的協(xié)整檢驗變成對矩陣Π的分析問題,是Johansen協(xié)整檢驗的基本原理。因為矩陣Π的秩等于它的非零特征根的個數(shù),因此可以通過對非零特征根個數(shù)的檢驗來檢驗協(xié)整關(guān)系和協(xié)整向量的秩。
Johansen協(xié)整檢驗主要包括兩種檢驗方法,一種是特征根跡檢驗,另一種是最大特征值檢驗。
在特征根跡檢驗中,由r個最大特征根可得到r個協(xié)整向量,而對于其余k-r個非協(xié)整組合來說,λr+1,…,λk應(yīng)該為0,于是可得到原假設(shè)、備擇假設(shè)為:
Hr0∶λr>0,λr+1=0
Hr1∶λr+1>0,r=0,1,…,k-1
相應(yīng)的檢驗統(tǒng)計量為:
ηx為特征根跡統(tǒng)計量。并依次檢驗這一系列統(tǒng)計量的顯著性。在最大特征值檢驗中,原假設(shè)、備選假設(shè)為:
Hr0∶λr+1=0,Hr1∶λr+1>0
檢驗統(tǒng)計量是基于最大特征值的,其形式為:
其中ξr為最大特征根統(tǒng)計量。并依次檢驗這一系列統(tǒng)計量的顯著性。就兩個時間序列的聯(lián)系而言,應(yīng)當(dāng)對yt=(p1t,p2t)、k=2進行驗證。當(dāng)ξ0<臨界值,接受H00(r=0),表明有0個協(xié)整向量,表明不存在協(xié)整關(guān)系。當(dāng)ξ0>臨界值,拒絕H00(r=0),表明至少有1個協(xié)整向量,表明存在協(xié)整關(guān)系。
格蘭杰因果檢驗實質(zhì)上是檢驗一個變量的滯后變量是否可以引入到其他變量方程中。一個變量如果受到其他變量的滯后影響,則稱它們具有格蘭杰因果關(guān)系。
在一個二元p階的VAR模型中:
其統(tǒng)計量為:
S1服從F分布。如果S1大于F的臨界值,則拒絕原假設(shè);否則接受原假設(shè):x不是y的格蘭杰原因。
從20世紀(jì)70年代末起,以計量經(jīng)濟學(xué)家Hendry為代表,一些經(jīng)濟學(xué)家將理論和數(shù)據(jù)信息有效結(jié)合,提出了動態(tài)計量經(jīng)濟學(xué)模型的理論與方法,為時間序列模型帶來了重要的發(fā)展。Hendry認(rèn)為,模型的建立應(yīng)該是從一個能夠代表數(shù)據(jù)生成過程的自回歸分布滯后模型(Autoregressive Distributed Lag,ADL)逐漸簡化,最后得到包含變量之間長期穩(wěn)定關(guān)系的簡單模型。包括幾何分布滯后模型、多項式分布滯后模型、自回歸分布滯后模型等。
幾何分布滯后模型反映了變量的影響程度隨滯后期的延長而按幾何級數(shù)遞減,模型為:
幾何分布滯后模型中的隨機擾動項通常存在一階負(fù)相關(guān)關(guān)系,參數(shù)估計比較復(fù)雜。為解決幾何分布滯后模型存在的問題,Almon提出了多項式分布滯后模型(Polynomial Distributed Lag,PDL)。對于方程右邊只含有外生變量及其滯后項的分布滯后模型:
由于解釋變量之間明顯存在多重共線性,不能采用OLS估計。為此,將βi分解為
式中,i=0,1,2,…,p,并且
即將每個參數(shù)用一個多項式表示,式(21)可改寫為多項式分布滯后模型:
一個PDL模型由三個因素確定,滯后期p,多項式的次數(shù)q,約束條件。
與DL和PDL兩種分布滯后模型比較,Jorgenson(1966)提出的自回歸分布滯后(Auto-regressive Distributed Lag,ADL)模型應(yīng)用更加廣泛。(p,q)階自回歸分布滯后模型的基本表達式為:
式中,xt-1是滯后i期的外生變量向量(維數(shù)與變量個數(shù)相同),且每個外生變量的最大滯后階數(shù)為τi;βi是參數(shù)向量。
誤差修正模型(ECM)的基本形式是由Davidson,Hendry,Srba和Yeo于1978年提出的,稱為DHSY模型。模型形式為:
式中,ecm是誤差修正項。模型的意義可以通過ADL(1,1)模型加以解釋。ADL(1,1)模型為:yt=β0+β1x1+β2yt-1+β3xt-1+εt,移項后整理可得:
[1]高鐵梅.計量經(jīng)濟分析方法與建模:EVIEWS應(yīng)用及實例[M].北京:清華大學(xué)出版社,2006.
[2]張曉峒.EViews使用指南與案例[M].北京:機械工業(yè)出版社,2007.
[3]易丹輝.數(shù)據(jù)分析與EVIES應(yīng)用[M].北京:中國統(tǒng)計出版社,2008.
F304.2
A
1002-6487(2011)07-0044-04
1 時間序列的趨勢及波動周期測定方法
安徽省教育廳人文社科項目(2007sk262);安徽省高校青年教師資助計劃項目(2008jqw108)
張克榮(1980-),女,安徽懷遠人,碩士,講師,研究方向:農(nóng)業(yè)經(jīng)濟管理。劉武藝(1979-),男,安徽安慶人,博士,副教授,研究方向:農(nóng)業(yè)經(jīng)濟管理。
(責(zé)任編輯/亦民)
在價格波動的研究方法上,經(jīng)濟學(xué)家們不斷建立嚴(yán)密的數(shù)學(xué)模型來研究時間序列問題,如ARMA模型、VAR模型、多元統(tǒng)計分析方法、狀態(tài)空間模型、HP濾波、BP濾波等各種濾波方法都被廣泛用來分析時間序列和經(jīng)濟周期等問題。本文主要參考高鐵梅[1]、張曉峒[2]和易丹輝[3]等有關(guān)時間序列分析方面的研究,闡釋應(yīng)用于價格波動的一些研究方法及其模型。目前,有關(guān)價格波動的研究方法,主要有季節(jié)調(diào)整方法、時間序列分析法、協(xié)整分析、向量自回歸模型、投入產(chǎn)出模型、誤差修正模型和脈沖響應(yīng)模型等。
對于月度、季度等時間序列的波動特征分析,其測度的基本思想是:首先從時間序列中消除季節(jié)因素和不規(guī)則因素;再利用趨勢分解方法把長期趨勢和循環(huán)要素分離開,研究經(jīng)濟變量的長期趨勢波動和景氣循環(huán)變動。由于食品消費價格指數(shù)的月度時間序列總是呈現(xiàn)一定的季節(jié)波動特征,這種季節(jié)變化常常會掩蓋波動的真實特征。因此,在進行數(shù)據(jù)分析前需要對時間序列進行季節(jié)調(diào)整,從而顯示出序列潛在的趨勢循環(huán)因素,趨勢循環(huán)因素能夠真實地反映經(jīng)濟時間序列運動的規(guī)律。目前測定長期趨勢的方法主要有四類。