李阿勇

[摘 要] 本文在貝葉斯理論框架下對常見的兩類模型：GARCH模型和SV模型進(jìn)行了比較研究，同時基于上證綜指和深證成指數(shù)據(jù)對兩類模型進(jìn)行了實證研究。無論在理論上還是實證分析，SV模型對資產(chǎn)收益波動性的刻畫能力要強于GARCH模型。

[關(guān)鍵詞] 貝葉斯； 隨機波動； GARCH

doi ： 10 . 3969 / j . issn . 1673 - 0194 . 2014 . 06. 044

[中圖分類號] F832.5 [文獻(xiàn)標(biāo)識碼] A [文章編號] 1673 - 0194（2014）06- 0073- 03

0 引 言

波動性模型通常用來描述資產(chǎn)收益率的條件方差，目前對波動性建模可以分為兩類： ① ARCH模型及其擴(kuò)展形式，這類模型采用確定性函數(shù)描述資產(chǎn)收益率的條件方差； ② SV模型及其擴(kuò)展形式，這類模型采用隨機方程來描述資產(chǎn)收益率的條件方差。Engle在對英國的經(jīng)濟(jì)通貨膨脹進(jìn)行研究的過程中提出了ARCH模型，Engle假設(shè)資產(chǎn)收益率的條件方差是時變的，同時可以用誤差滯后項的平方進(jìn)行解釋。Bollerslev在Engle研究的基礎(chǔ)上，將資產(chǎn)收益率的條件方差同時可以由其滯后項作為解釋變量納入模型中，提出了GARCH模型。GARCH模型是對ARCH模型的重要擴(kuò)展，在研究金融和資本市場時間序列的特征中得到了廣泛運用。GARCH模型的擴(kuò)展形式有很多，Bollerslev、Chou、Kroner等學(xué)者都曾對GARCH模型進(jìn)行過詳細(xì)的論述，Li等學(xué)者也對GARCH模型族進(jìn)行過綜述。SV模型最早由Taylor等學(xué)者提出，Taylor假設(shè)收益率的擾動項是不可觀測的，可以用一個隨機過程進(jìn)行刻畫，這個假設(shè)決定了隨機波動模型是一個具有動態(tài)波動特性的模型。Clark、Harvey、Polson、Rossi等學(xué)者在早期均對隨機波動模型進(jìn)行過研究，Ghysels等學(xué)者對SV模型族進(jìn)行了詳細(xì)的論述。本文選取應(yīng)用最廣泛的GARCH - N，SV - N兩個模型進(jìn)行建模，統(tǒng)一采用基于Gibbs抽樣技術(shù)的MCMC方法對這兩個模型的參數(shù)進(jìn)行估計，并利用中國股市數(shù)據(jù)進(jìn)行實證研究，分析這些模型對我國股票市場的刻畫能力。

1 GARCH模型及其貝葉斯統(tǒng)計推斷

ARCH形式比較簡單，但其在實際使用過程中往往需要用很高的階，在應(yīng)用時受到較大的限制。Bollerslev提出了廣義ARCH模型（GARCH模型），解決了ARCH模型參數(shù)過多的缺點。GARCH模型同時把均值修正后資產(chǎn)收益率的過去值和條件方差的過去值作為條件方差的解釋標(biāo)量，是對ARCH模型的重要擴(kuò)展，在研究金融和資本市場時間序列的特征中得到了廣泛運用。GARCH模型的具體形式如下：

at = σtεt，εt ～ i.i.dN（0，1） （1）

σ2t = α0 + ■αia2t - i + ■βj σ2t - j （2）

式中，p ≥ 0，q ≥ 0，α0 > 0，αi ≥ 0，βj ≥ 0。

at為時間t的均值修正收益率，σ2t為條件方差，εt服從標(biāo)準(zhǔn)正態(tài)分布，p ≥ 0，q ≥ 0，α0 > 0，αi ≥ 0，βj ≥ 0。這個模型被稱為GARCH - N模型。

2 SV模型及其貝葉斯統(tǒng)計推斷

SV模型最早由Taylor等學(xué)者提出，Taylor假設(shè)收益率的擾動項是不可觀測的，可以用一個隨機過程進(jìn)行刻畫，這個假設(shè)決定了隨機波動模型是一個具有動態(tài)波動特性的模型。SV模型具體形式如下：

at = σtεt，εt ～ i.i.dN（0，1） （3）

lnσ2t = α0 + βlnσ2t - 1 + ηt，ηt ～ i.i.dN（0，σ2η） （4）

式中，at表示消去均值后第t期的收益。

3 實證研究

為了分析GARCH模型和SV模型對中國股票市場波動率的刻畫能力，本文選取上證綜指和深圳成指來對兩個模型進(jìn)行實證研究。我國股票市場與1996年12月16日起實行了漲跌停板限價交易制度，為了保證數(shù)據(jù)的一致性，本文樣本數(shù)據(jù)選擇1997年1月2日至2013年3月1日，共3 908天的日線數(shù)據(jù)，同時考察日收益率，以每日收盤價為基準(zhǔn)。本文采用連續(xù)復(fù)利計算收益率，即rt = ln（St） - ln（St - 1），St為第t日的收盤價，rt為第t日收益率。

表1為上證綜指收益率與深證成指收益率的描述性統(tǒng)計結(jié)果。從表中可以看出上證綜指收益率還是深證成指收益率的峰度分別為7.329 8、6.439 7，遠(yuǎn)大于正態(tài)分布的峰度3，說明收益率序列分布呈現(xiàn)明顯尖峰肥尾特征。表中J - B為Jarque - Bera檢驗統(tǒng)計量，用以檢驗序列是否服從正態(tài)分布，H = 1則拒絕序列服從正態(tài)分布的假設(shè)，J - B統(tǒng)計量也表明收益率分布并非正態(tài)分布。LM（q）為Engle提出的拉格朗日乘子（LM）q階檢驗數(shù)，H = 1則拒絕波動率自相關(guān)系數(shù)全為0的假設(shè)，說明序列有ARCH效應(yīng)。小括號內(nèi)的數(shù)值為檢驗p值。LM檢驗表明無論之后階數(shù)取10或20，LM統(tǒng)計量都較大，波動率自相關(guān)系數(shù)全為0的假設(shè)被拒絕，說明模型具有明顯的ARCH效應(yīng)，即明顯的波動聚集性。

為了進(jìn)一步對收益率的情況進(jìn)行分析，本文對上證綜指收益率序列和深證成指收益率序列的平穩(wěn)性進(jìn)行研究，分別對兩組收益率序列作ADF單位根檢驗，滯后階數(shù)為5。在顯著性水平1%下，上證綜指收益率的ADF統(tǒng)計量為 -25.353 6，遠(yuǎn)小于臨界值 -2.863 6，接受上證綜指收益率序列為隨機游走的概率為0.001，因此拒絕上證綜指收益率為隨機游走的假設(shè)，即該收益率序列是平穩(wěn)的。深證成指收益率的ADF統(tǒng)計量為 -25.115 9，也遠(yuǎn)小于臨界值，說明深證成指收益率也是平穩(wěn)的。

目前對SV模型最有效的估計方法是基于MCMC技術(shù)的貝葉斯統(tǒng)計推斷方法，MCMC方法最重要的軟件包是BUGS和WinBUGS。BUGS是Bayesian inference using gibbs sampling的縮寫，最初由英國劍橋大學(xué)生物統(tǒng)計研究所開發(fā)，是目前進(jìn)行MCMC計算最方便的軟件。WinBUGS是BUGS的Windows版本，可以免費使用，本文運用WinBUGS軟件完成GARCH模型和SV模型的參數(shù)估計工作。對兩個模型模擬30 000次，首先對模型的收斂性進(jìn)行判別。

圖1為樣本迭代軌跡歷史，由圖1可以看出迭代歷史基本趨于穩(wěn)定，說明迭代過程是收斂的，從自相關(guān)函數(shù)圖也可對收斂性進(jìn)行判斷。由自相關(guān)函數(shù)圖可以看出，模型中各參數(shù)的自相關(guān)函數(shù)很快接近于0，說明迭代過程已經(jīng)收斂。

得到模型的各參數(shù)估計結(jié)果見表2。

以2013年3月1日至2013年4月1日共21個交易日的上證綜指和深證成指日收益數(shù)據(jù)為樣本外數(shù)據(jù)對兩類模型的樣本外預(yù)測能力進(jìn)行評價。

評價指標(biāo)選取均方根誤差統(tǒng)計量（RMSE）和絕對誤差統(tǒng)計量（MAE），兩個統(tǒng)計指標(biāo)如下：

RMSE = ■ （5）

MAE = ■| ri - ■i | （6）

式中，ri為樣本點i的實際收益率，■i為預(yù)測收益率。

由RMSE和MAE的表達(dá)式可以看出兩個統(tǒng)計指標(biāo)值越小說明預(yù)測結(jié)果越精確。通過前文所述的模型進(jìn)行樣本外數(shù)據(jù)預(yù)測效果。

由表可以看出無論是RMSE指標(biāo)還是MAE指標(biāo)，SV模型的預(yù)測效果均好于GARCH模型。

4 結(jié) 論

本章對常見的兩類模型：GARCH模型和SV模型進(jìn)行了比較研究，同時基于上證綜指和深證成指數(shù)據(jù)對兩類模型進(jìn)行了實證研究。無論在理論上還是實證分析都可以看出SV模型對資產(chǎn)收益波動性的刻畫能力要強于GARCH模型。上海股市收益率與深圳股市收益率存在明顯的波動聚集性和尖峰厚尾性。GARCH模型和SV模型對這兩種性質(zhì)均具有較強的刻畫能力，而且SV模型對兩個股票市場波動性的描述更精確。

主要參考文獻(xiàn)

[1] R F Engle. Autoregressive Conditional Heteroscedasticity with Estimates of the Variance of United Kingdom Inflations [J]. Econometrics，1982（50）：987-1007.

[2] 王鵬，王建瓊，魏宇. 自回歸條件方差—偏度—峰度：一個新的模型[J]. 管理科學(xué)學(xué)報，2009（5）：121-129.

[3] C W J Granger，J Hallman. Long Memory Series with Attractors [J]. Oxford Bulletin of Econometrics and Statistics，1991，53（1）：11-26.

[4] Y Hammao，R Masulis. Correlations in Price Changes and Volatility Across International Stock Markets[J]. Review of Financial Studies，1990（3）：281-307.

[5] P Christoffersen，K Jacobs. Which GARCH Model for Option Valuation？[J]. Management Science，2004，50（9）：1204-1221.

[6] T Bollerslev. Generalized Autoregressive Conditional Heteroskedasticity [J]. Journal of Econometrics，1986（31）：307-327.