金融時(shí)間序列分析第部分時(shí)間序列分析基礎(chǔ)波動(dòng)率模型

1、金融時(shí)間序列分析金融時(shí)間序列分析 陸貴斌2012年年10月月波動(dòng)率模型波動(dòng)率模型 （一）（一）問(wèn)題的提出問(wèn)題的提出 xt = xt-1 + ut 其中其中, , ut 為白噪聲過(guò)程。為白噪聲過(guò)程。計(jì)量經(jīng)濟(jì)學(xué)模型中計(jì)量經(jīng)濟(jì)學(xué)模型中的異方差的異方差通常通常屬于遞增型異方差屬于遞增型異方差,但利率，匯率，股票收益等時(shí)間序列中存在的異方差但利率，匯率，股票收益等時(shí)間序列中存在的異方差卻不卻不是是遞增型異方差。遞增型異方差。例如，例如，匯率，匯率，股票價(jià)格常常用隨機(jī)游走描述股票價(jià)格常常用隨機(jī)游走描述： 1995-2000年日元兌美元匯率時(shí)間序列及差分序列見(jiàn)下圖年日元兌美元匯率時(shí)間序列及差分序列見(jiàn)下圖80

2、100120140160200400600800100012001400JPY (1995-2000)-8-6-4-20246200400600800100012001400D(JPY) (1995-2000)日元兌美元匯率序列日元兌美元匯率序列JPY(1995-2000) 日元兌美元匯率差分序列（收益）日元兌美元匯率差分序列（收益）D(JPY) 02468200400600800100012001400Volatility of returns0102030405060200400600800100012001400DJPY2收益絕對(duì)值序列收益絕對(duì)值序列 (1995-2000) D(JPY)

3、的平方的平方 (1995-2000) 這種序列的特征是：這種序列的特征是：（1）過(guò)程的方差不僅隨時(shí)間變化，而且有時(shí)變化得很激烈。）過(guò)程的方差不僅隨時(shí)間變化，而且有時(shí)變化得很激烈。 （2）按時(shí)間觀察，表現(xiàn)出）按時(shí)間觀察，表現(xiàn)出 “波動(dòng)集群波動(dòng)集群” （volatility clustering）特征，）特征， 方差在一定時(shí)段中比較小，而在另一時(shí)方差在一定時(shí)段中比較小，而在另一時(shí)段中比較大。段中比較大。（3）從取值的分布看表現(xiàn)的則是）從取值的分布看表現(xiàn)的則是 “高峰厚尾高峰厚尾” 特征，特征，（leptokurtosis and fat-tail）即均值附近與尾區(qū)的概率值比）即均值附近與尾區(qū)的概率

4、值比正態(tài)分布大，而其余區(qū)域的概率比正態(tài)分布小。正態(tài)分布大，而其余區(qū)域的概率比正態(tài)分布小。 高高峰厚尾峰厚尾分布曲線分布曲線 正態(tài)分布正態(tài)分布曲線曲線 高高峰厚尾峰厚尾分布分布特征特征示意圖示意圖 顯然現(xiàn)期方差與前期的顯然現(xiàn)期方差與前期的“波動(dòng)波動(dòng)”有關(guān)系。有關(guān)系。自回歸條件異方差模型自回歸條件異方差模型（Engle 1982）通常有兩類(lèi)：通常有兩類(lèi)：1）用確定的函數(shù)來(lái)刻畫(huà)異方差的演變，）用確定的函數(shù)來(lái)刻畫(huà)異方差的演變，GARCH模型模型2）用隨機(jī)方程來(lái)描述異方差。）用隨機(jī)方程來(lái)描述異方差。隨機(jī)波動(dòng)率模型隨機(jī)波動(dòng)率模型 （二）自回歸條件異方差模型（二）自回歸條件異方差模型（ARCH模型）模型）(

5、Autoregressive Conditional Heteroskedasticity Model)一、模型的提出一、模型的提出考察一個(gè)考察一個(gè)AR(p)過(guò)程過(guò)程1122tttptptycyyy其中：其中： 是白噪聲：是白噪聲：t0,tE2,()0,tsstEelse 的條件期望是：的條件期望是：ty121122(,)tttttptpE y yycyyy 的無(wú)條件期望是：的無(wú)條件期望是：ty12()(1)tpE ycA. 的無(wú)條件期望是常數(shù)，但的無(wú)條件期望是常數(shù)，但 的條件期望卻是隨時(shí)的條件期望卻是隨時(shí)間而變化的。間而變化的。B. 的無(wú)條件方差是一常數(shù)的無(wú)條件方差是一常數(shù) ，但，但 的條件

6、方差卻的條件方差卻可能隨時(shí)間而變化?？赡茈S時(shí)間而變化。一種方法是一種方法是將將 視為服從視為服從 AR(m) 過(guò)程過(guò)程。tytyt2t2t22221122tttmt mtcu 其中：其中： 是一新的白噪聲：是一新的白噪聲：tu0,tEu 2,()0,tsstE u uelse此時(shí)，此時(shí)，222212122221122(,)(,) ttttttttmt mEEc 即：即：212122221122(,)(,) ttttttttmt mVarEc 考慮考慮 k 變量回歸模型（或變量回歸模型（或AR(p)過(guò)程）過(guò)程）01 1ttkkttyxx滿足滿足121122201122()0,()()tttttt

7、tttqt qEVarEh 二、模型表達(dá)形式二、模型表達(dá)形式01212(1),0(2),1qq 其中：其中：則稱(chēng)則稱(chēng) 服從服從 q 階自回歸條件異方差過(guò)程。記為階自回歸條件異方差過(guò)程。記為t( )tARCH q保證條件方差為正數(shù)保證條件方差為正數(shù)平穩(wěn)序列平穩(wěn)序列101 1101 1()()tttkkttttkktE yExxxx均值方程均值方程 條件均值模型條件均值模型注意注意：1t表示表示 t 時(shí)刻之前可獲得的信息集。時(shí)刻之前可獲得的信息集。例例 ARCH(1)模型為模型為21011()tttthVar ARCH(q)過(guò)程的另一種表達(dá)方式過(guò)程的另一種表達(dá)方式01 1ttkkttyxxttth

8、 v22201122tttqt qh . .(0,1)i i dtvN2111()0,()()1ttttttE vE vVar v三、三、ARCH模型的性質(zhì)模型的性質(zhì)1、 的無(wú)條件均值和方差分別為的無(wú)條件均值和方差分別為t201( )0( )( ) 1tttEVarE21011()tttthVar ttth v主要討論主要討論ARCH(1)的性質(zhì)的性質(zhì)22122011011011( )()()()()()ttttttttVarEE EEhEEVar 2t因?yàn)橐驗(yàn)?是平穩(wěn)的，所以是平穩(wěn)的，所以221101( )()()( )1ttttVarVarEVar顯然要求顯然要求1012、令、令221()

9、tttttthVar則則 是白噪聲過(guò)程。是白噪聲過(guò)程。t21221()()()0tttttttEEE VarEE E222122()()()() () ()0tstssstssttsststtsststsEE EEEEEhEEE EEEEst 這樣這樣22011ttttth 過(guò)程過(guò)程 的形式類(lèi)似于的形式類(lèi)似于AR(1)。過(guò)程過(guò)程 前后不相關(guān)，但過(guò)程前后不相關(guān)，但過(guò)程 卻是前后相關(guān)的卻是前后相關(guān)的（因而也（因而也是不獨(dú)立的）。是不獨(dú)立的）。2tt2t3、ARCH(1)模型的尾部特征模型的尾部特征42222211011()3 ()33()ttttttEEh 44221011222400111122

10、4001111() ()3 ()3 (2)32()()ttttttttEE EEEVarE 殘差平方服從一個(gè)殘差平方服從一個(gè)異方差的異方差的AR(1)過(guò)程過(guò)程隨機(jī)變量之間不相關(guān)，只能說(shuō)明它們之間沒(méi)有線性關(guān)系，不能說(shuō)隨機(jī)變量之間不相關(guān)，只能說(shuō)明它們之間沒(méi)有線性關(guān)系，不能說(shuō)明它們之間沒(méi)有非線性關(guān)系。明它們之間沒(méi)有非線性關(guān)系。422400012112241011240142111()32()13(12)3()13(1)(),(1)(13)(130)ttttEEEmE 4411()(),()()ttttEEVarVar4220112221102121()3(1)(1)()(1)(13)13313ttE

11、Var四、四、ARCH模型的建立模型的建立ARCH模型的建立的基本步驟模型的建立的基本步驟（1）建立時(shí)間序列的計(jì)量模型，去掉任何線性關(guān)系，）建立時(shí)間序列的計(jì)量模型，去掉任何線性關(guān)系， 對(duì)估計(jì)的殘差檢驗(yàn)對(duì)估計(jì)的殘差檢驗(yàn)ARCH效果；效果；（2）識(shí)別）識(shí)別ARCH模型的階數(shù)，估計(jì)模型；模型的階數(shù)，估計(jì)模型；（3）檢驗(yàn)）檢驗(yàn)ARCH模型，根據(jù)情況修改模型。模型，根據(jù)情況修改模型。1、建立收益率序列的計(jì)量模型，檢驗(yàn)、建立收益率序列的計(jì)量模型，檢驗(yàn)ARCH效果效果（1）選擇條件均值模型并估計(jì)出參數(shù)）選擇條件均值模型并估計(jì)出參數(shù)可以是回歸模型可以是回歸模型 或或ARMA模型模型tttyX( )( )ttB

12、 yB （2）計(jì)算出殘差值）計(jì)算出殘差值( )( )tttttByXyB或檢驗(yàn)殘差的平方檢驗(yàn)殘差的平方 是否存在自相關(guān)。是否存在自相關(guān)。2t檢驗(yàn)序列相關(guān)常用的方法有檢驗(yàn)序列相關(guān)常用的方法有：(a) D.W統(tǒng)計(jì)量檢驗(yàn)統(tǒng)計(jì)量檢驗(yàn)(b) QLB統(tǒng)計(jì)量檢驗(yàn)統(tǒng)計(jì)量檢驗(yàn)(c) 拉格朗日乘數(shù)檢驗(yàn)拉格朗日乘數(shù)檢驗(yàn)(LM檢驗(yàn)檢驗(yàn))注意：注意： D.W統(tǒng)計(jì)量檢驗(yàn)的不足統(tǒng)計(jì)量檢驗(yàn)的不足1、僅僅檢驗(yàn)殘差序列是否存在一階序列相關(guān)；、僅僅檢驗(yàn)殘差序列是否存在一階序列相關(guān)；2、回歸方程右邊如果存在滯后因變量，、回歸方程右邊如果存在滯后因變量， D.W檢驗(yàn)不再有效；檢驗(yàn)不再有效；3、 D.W統(tǒng)計(jì)量的擾動(dòng)項(xiàng)在原假設(shè)下依賴(lài)于系數(shù)

13、矩陣。統(tǒng)計(jì)量的擾動(dòng)項(xiàng)在原假設(shè)下依賴(lài)于系數(shù)矩陣。 因此，因此，檢驗(yàn)序列相關(guān)常用檢驗(yàn)序列相關(guān)常用 QLB統(tǒng)計(jì)量檢驗(yàn)或拉格朗日統(tǒng)計(jì)量檢驗(yàn)或拉格朗日乘數(shù)檢驗(yàn)乘數(shù)檢驗(yàn)(LM檢驗(yàn)檢驗(yàn))。QLB統(tǒng)計(jì)量統(tǒng)計(jì)量(Ljung-Box Q統(tǒng)計(jì)量統(tǒng)計(jì)量) 檢驗(yàn)法檢驗(yàn)法1、計(jì)算出殘差值、計(jì)算出殘差值( )( )tttttByXyB或2、計(jì)算出殘差平方的均值、計(jì)算出殘差平方的均值2211TttT3、計(jì)算出殘差平方的樣本自相關(guān)系數(shù)、計(jì)算出殘差平方的樣本自相關(guān)系數(shù)222212221()()()Ttt it iiTtt 4、計(jì)算、計(jì)算QLB統(tǒng)計(jì)量統(tǒng)計(jì)量21(2)piLBiQT TTiT為樣本容量，為樣本容量，p為設(shè)定的滯后階數(shù)

14、。為設(shè)定的滯后階數(shù)。原假設(shè)原假設(shè)H0：序列不存在：序列不存在p階自相關(guān)；階自相關(guān)；備選假設(shè)：備選假設(shè)：H1為序列存在為序列存在p階自相關(guān)。階自相關(guān)。221(2)( )piLBiQT TpTi近似當(dāng)當(dāng)12( )LBQp時(shí)時(shí)拒絕原假設(shè)拒絕原假設(shè)H0拉格朗日乘數(shù)檢驗(yàn)拉格朗日乘數(shù)檢驗(yàn)(LM檢驗(yàn)檢驗(yàn))法法1、使用最小二乘法估計(jì)最適當(dāng)?shù)摹⑹褂米钚《朔ü烙?jì)最適當(dāng)?shù)腁R(n)模型模型01122tttnt ntyyyy2、計(jì)算出殘差值、計(jì)算出殘差值 ，進(jìn)行下面的回歸分析，進(jìn)行下面的回歸分析t222201122tttptptu 3、零假設(shè)、零假設(shè)012:0pH備擇假設(shè)備擇假設(shè)112:,pH 不全為零。不全為零。

15、（不存在（不存在ARCH 效果）效果）4、檢驗(yàn)統(tǒng)計(jì)量：、檢驗(yàn)統(tǒng)計(jì)量：LM=T R2其中，其中，T為回歸式的樣本容量，為回歸式的樣本容量，R2為擬合優(yōu)度為擬合優(yōu)度21ESSRSSRTSSTSS 6、判別法則：、判別法則：若若LM 2 (p)，接受，接受H1（說(shuō)明殘差是（說(shuō)明殘差是ARCH過(guò)程）。過(guò)程）。5、統(tǒng)計(jì)量、統(tǒng)計(jì)量LM的分布的分布在零假設(shè)成立的條件下，統(tǒng)計(jì)量在零假設(shè)成立的條件下，統(tǒng)計(jì)量LM近似服從近似服從2( )p2、識(shí)別、識(shí)別ARCH模型的階數(shù)，估計(jì)模型模型的階數(shù)，估計(jì)模型階數(shù)識(shí)別：階數(shù)識(shí)別：可以分析時(shí)間序列可以分析時(shí)間序列 的自相關(guān)函數(shù)和的自相關(guān)函數(shù)和偏自相關(guān)函數(shù)。也可以使用偏自相關(guān)函

16、數(shù)。也可以使用AIC或或BIC準(zhǔn)則來(lái)確定序準(zhǔn)則來(lái)確定序列適當(dāng)?shù)臏箝L(zhǎng)度。列適當(dāng)?shù)臏箝L(zhǎng)度。 2t參數(shù)估計(jì)：參數(shù)估計(jì)：一般使用極大似然法估計(jì)：一般使用極大似然法估計(jì)：21111( ,)exp22TtTptpttfhh 為了簡(jiǎn)單，可以使用條件似然函數(shù)：為了簡(jiǎn)單，可以使用條件似然函數(shù)：21111( ,)exp( ,)22TtTptpttffhh 3、檢驗(yàn)、檢驗(yàn)ARCH模型模型對(duì)于對(duì)于ARCH模型，把殘差標(biāo)準(zhǔn)化模型，把殘差標(biāo)準(zhǔn)化tttvh則則 應(yīng)該是正態(tài)白噪聲序列。所以應(yīng)該是正態(tài)白噪聲序列。所以（1） Ljung-Box Q 統(tǒng)計(jì)量檢驗(yàn)序列是否存在自相關(guān)，統(tǒng)計(jì)量檢驗(yàn)序列是否存在自相關(guān)，（2）利用）利用

17、 JB檢驗(yàn)檢驗(yàn) 或或 Q-Q 圖檢驗(yàn)序列的正態(tài)性圖檢驗(yàn)序列的正態(tài)性 tv 此時(shí)，此時(shí)， 各期之間是獨(dú)立的，因而不可能各期之間是獨(dú)立的，因而不可能ARCH存在效果。存在效果。 也可以直接對(duì)也可以直接對(duì)ARCH模型的殘差進(jìn)行模型的殘差進(jìn)行ARCH檢驗(yàn)。檢驗(yàn)。 tv（三）（三）GARCH模型模型一、一、ARCH(q)模型的缺點(diǎn)模型的缺點(diǎn)1、當(dāng)、當(dāng) q 較大時(shí)，參數(shù)估計(jì)很難做到精確較大時(shí)，參數(shù)估計(jì)很難做到精確2、為了保證條件方差為正，參數(shù)要求為正。、為了保證條件方差為正，參數(shù)要求為正。 當(dāng)參數(shù)過(guò)多時(shí)，用實(shí)際數(shù)據(jù)估計(jì)出的模型往往不能滿足當(dāng)參數(shù)過(guò)多時(shí)，用實(shí)際數(shù)據(jù)估計(jì)出的模型往往不能滿足 這一要求，從而，模

18、型不具實(shí)用性。這一要求，從而，模型不具實(shí)用性。3、ARCH模型會(huì)高估波動(dòng)率。模型會(huì)高估波動(dòng)率。二、二、GARCH模型的定義模型的定義01 1ttkkttyxx考慮考慮 k 變量回歸模型變量回歸模型滿足滿足12( ,)ttttyf t yy稱(chēng)稱(chēng)( , )tGARCH q p121011()0()ttqptttit ijtjijEVarhh 注：注：GARCH模型的優(yōu)點(diǎn)在于可用低階模型的優(yōu)點(diǎn)在于可用低階GARCH模型來(lái)模型來(lái)代表高階的代表高階的ARCH模型，從而使模型的識(shí)別和估計(jì)都變模型，從而使模型的識(shí)別和估計(jì)都變得比較容易。得比較容易。例如，用例如，用GARCH（1，1）擬合）擬合 的變化特征與

19、用的變化特征與用ARCH（20）擬合相比較，兩者的效果十分相近。）擬合相比較，兩者的效果十分相近。tARCH項(xiàng)系數(shù)項(xiàng)系數(shù)GARCH項(xiàng)系項(xiàng)系數(shù)數(shù)三、三、GARCH模型的另一種定義模型的另一種定義01 1ttkkttyxx若若ttth v2011qptit ijtjijhh . .(0,1)i i dtvN則稱(chēng)則稱(chēng)( , )tGARCH q p四、四、GARCH模型的性質(zhì)模型的性質(zhì)1、當(dāng)、當(dāng)p=0時(shí)，時(shí)，GARCH過(guò)程變?yōu)檫^(guò)程變?yōu)锳RCH過(guò)程。過(guò)程。2、 GARCH過(guò)程的含義是過(guò)程的含義是 條件方差條件方差 ht 是是 ht-1, ,ht-p 和和1,tt q 的函數(shù)。的函數(shù)。3、參數(shù)、參數(shù) 是保

20、證條件方差為正的充分條件。是保證條件方差為正的充分條件。4、若、若則則5、 的平穩(wěn)性條件為的平穩(wěn)性條件為 0(1,2, ),0(1,2, )iiiqip( , ),tGARCH q p2( , ),max,tARMA p r rp q 2t111,qp這時(shí)這時(shí) 也是寬平穩(wěn)的。也是寬平穩(wěn)的。t五、五、GARCH(1,1)模型模型2、GARCH(1,1) 的條件方差和無(wú)條件方差的條件方差和無(wú)條件方差條件方差是條件方差是 ht ，通過(guò)對(duì)上式兩邊取期望可得無(wú)條件方差通過(guò)對(duì)上式兩邊取期望可得無(wú)條件方差ttth v2011111111, (0,0,1)ttthh 10()ttttEVarh20111ttD

21、E1、大的、大的 會(huì)緊跟著另一個(gè)大的會(huì)緊跟著另一個(gè)大的 ，這樣就會(huì)產(chǎn)生，這樣就會(huì)產(chǎn)生在金融時(shí)間序列中有名的波動(dòng)率聚類(lèi)現(xiàn)象。在金融時(shí)間序列中有名的波動(dòng)率聚類(lèi)現(xiàn)象。21t2t3、GARCH(1,1)模型與模型與ARCH模型的關(guān)系模型的關(guān)系 把把GARCH(1,1)過(guò)程遞推展開(kāi)，等價(jià)于一個(gè)過(guò)程遞推展開(kāi)，等價(jià)于一個(gè)ARCH()。( , ),tGARCH q p4、若、若 則則2( , ),max,tARMA p rrp q例如例如(1,2),tGARCH20112211tttthhh 即即上式兩邊加上上式兩邊加上 得，得，2t2220112211tttttthhh 222011221122222201

22、1122211221122222011122211122()()()() ()ttttttttttttttttttttthhhhhhh 所以所以222011122222211122()()()ttttttttthhh 令令2,ttthw則則 是白噪聲，且是白噪聲，且 2220112211122()ttttttwww tw即即2(2,2),tARMA5、GARCH(1,1) 模型的尾部特征模型的尾部特征可以證明：當(dāng)可以證明：當(dāng) 時(shí)，時(shí)， 221111 2()04211222111()1()33()1()2ttEVar因此，因此， GARCH(1,1) 模型的尾部比正態(tài)分布的尾部厚。模型的尾部比正

23、態(tài)分布的尾部厚。特別，當(dāng)特別，當(dāng) 時(shí)，時(shí)，1042()3()ttEVar說(shuō)明系數(shù)說(shuō)明系數(shù) 在決定在決定 的尾部行為時(shí)起關(guān)鍵作用。若的尾部行為時(shí)起關(guān)鍵作用。若1t10則則 的尾部行為與正態(tài)噪聲相似，若的尾部行為與正態(tài)噪聲相似，若 ，則，則 有厚尾性。有厚尾性。10tt向前多步預(yù)測(cè)：向前多步預(yù)測(cè)：011( )()(1),2TTh kh kk事實(shí)上，因?yàn)槭聦?shí)上，因?yàn)?2ttthv21011201120111()(1)ttttttttthhhvhhh v 注意到，注意到，21(1)0ttE v所以所以T時(shí)向前時(shí)向前k步預(yù)測(cè)為：步預(yù)測(cè)為：011( )()(1)TTh kh k21011(1)TTTThh

24、h 向前一步預(yù)測(cè)：向前一步預(yù)測(cè)：6、 GARCH(1,1) 模型的預(yù)測(cè)模型的預(yù)測(cè)011( )( ),()1Tth kVark 反復(fù)迭代可得反復(fù)迭代可得11011111 ()( )()(1)1kkTTh kh當(dāng)當(dāng) 時(shí)有時(shí)有111非對(duì)稱(chēng)模型非對(duì)稱(chēng)模型 資產(chǎn)的收益率和波動(dòng)性存在一個(gè)特點(diǎn)資產(chǎn)的收益率和波動(dòng)性存在一個(gè)特點(diǎn)不對(duì)稱(chēng)性，收益不對(duì)稱(chēng)性，收益率上升時(shí)，波動(dòng)性較小，收益率下降時(shí)，波動(dòng)性較大。率上升時(shí)，波動(dòng)性較小，收益率下降時(shí)，波動(dòng)性較大。一、一、TARCH模型模型(門(mén)限門(mén)限ARCH模型模型)ttth v222011111ttptptqt qthhhS 其中其中 是一個(gè)虛擬變量，是一個(gè)虛擬變量，S11

25、1,00,0ttS上式表明：好消息上式表明：好消息 和壞消息和壞消息 對(duì)對(duì)條件方差的有不同的影響。條件方差的有不同的影響。1(0)t1(0)t例：例： 估計(jì)估計(jì)1995.1.2-2001.12.31滬市每日股票價(jià)格收盤(pán)指數(shù)滬市每日股票價(jià)格收盤(pán)指數(shù)的的TARCH模型。模型。均值方程：均值方程：1ln()1.000024ln()tttspsp方差方程：方差方程：6221112.57 100.070.0630.909tttthShARCH(1)GARCH(1)RESID*ARCH(1)二、二、EGARCH模型模型 (指數(shù)指數(shù)GARCH模型模型)ttth v01110111ln( )ln()ln() ln()ln() ttptpqit it iit iitptpqt it it iiiit it it ihhhvE vg vhhEghhh在條件方差在條件方差 ht 中引入了參數(shù)中引入了參數(shù) gi ，使得，使得 ht 在沖擊在沖擊 vt 取正、負(fù)值時(shí)取正、負(fù)值時(shí)有不同程度的變化，從而可以描述沖擊對(duì)價(jià)格波動(dòng)的非對(duì)稱(chēng)影響。有不同程度的變化，從而可以描述沖擊對(duì)價(jià)格波動(dòng)的非對(duì)稱(chēng)影響。在在Eviews中，條件方差方程為：中，條件方差方程為：11111ln( )ln()tttttthchhh例：例： 估計(jì)估計(jì)1995.1.2-2001.12.31滬市每日股