基于貝葉斯Wishart波動模型的原油市場與股市動態(tài)相依性研究

朱慧明，彭 成，游萬海，鄧 超

（湖南大學(xué)工商管理學(xué)院，湖南長沙410082）

基于貝葉斯Wishart波動模型的原油市場與股市動態(tài)相依性研究

朱慧明，彭 成，游萬海，鄧 超

（湖南大學(xué)工商管理學(xué)院，湖南長沙410082）

針對時變相關(guān)系數(shù)矩陣在多變量隨機波動模型的估計問題，構(gòu)建了貝葉斯動態(tài)相關(guān)Wishart波動模型。在CC-MSV模型的基礎(chǔ)上，設(shè)置精度矩陣服從Wishart分布，使得模型的相關(guān)系數(shù)矩陣具有時變特征。通過模型的統(tǒng)計結(jié)構(gòu)分析，選擇參數(shù)先驗分布，設(shè)計相應(yīng)的Gibbs-MTM-ARMS混合算法，據(jù)此估計模型參數(shù)；并利用上證綜合指數(shù)、標(biāo)普500指數(shù)與原油期貨價格數(shù)據(jù)進行實證分析。研究結(jié)果表明：模型能夠有效地刻畫原油市場與股票市場的動態(tài)相依性；金融危機期間，股票市場與原油市場的相關(guān)性較強，并且難以判斷正負(fù)方向；金融危機后，中國股票市場與原油市場呈現(xiàn)極微弱的相關(guān)性，而美國股票市場與原油市場的正相關(guān)性較為明顯。

動態(tài)相依性；隨機波動；貝葉斯分析；Wishart分布；Gibbs-MTM-ARMS混合算法

1 引言

原油作為基礎(chǔ)的能源和化工原料，是最主要的生產(chǎn)要素之一，在經(jīng)濟與金融發(fā)展中發(fā)揮著重要的作用。股票市場是經(jīng)濟發(fā)展水平的重要標(biāo)志，研究國際原油市場與股票市場的關(guān)系無論是從國家能源安全，還是從風(fēng)險防范、資源配置角度都具有重要的理論意義與現(xiàn)實意義。隨著世界經(jīng)濟對原油依賴性越來越強，原油市場與股票市場的關(guān)系也日益成為學(xué)者們關(guān)注的焦點問題［1-7］。金融資產(chǎn)收益率與收益率的波動是金融市場兩個基本的概念。目前，大量文獻主要從收益的視角研究原油市場與股票市場之間的相關(guān)性，如Park和Ratti［8］利用多元向量自回歸模型發(fā)現(xiàn)油價對股票市場具有顯著的影響；Miller和Ratti［9］運用協(xié)整向量誤差相關(guān)模型發(fā)現(xiàn)原油價格波動與OECD國家的股票市場收益是負(fù)相關(guān)的。同樣的，收益率的波動性作為體現(xiàn)金融市場質(zhì)量與效率的有效指標(biāo)，直接影響著市場的不確定性與風(fēng)險，描述原油市場與股票市場的波動率特征，測度波動率的大小及刻畫它們的動態(tài)特征也具有極為重要的理論與實際意義。因此，從波動模型的角度研究原油市場與股票市場的聯(lián)系也越來越受到青睞，如姬強和范英［10］利用動態(tài)條件相關(guān)的多元GARCH模型，分析美國次信貸危機發(fā)生前后國際原油市場和中、美股票市場間的協(xié)動性；Vo［11］利用常數(shù)相關(guān)和動態(tài)相關(guān)多變量隨機波動模型研究美國股票市場與原油市場之間的相關(guān)性。多元GARCH模型與多變量隨機波動模型在考慮波動方程時均假設(shè)變量的波動僅受到前期同一變量影響，在一定程度上限制了模型的解釋能力。Philipov和Glickman［12］、Asai［13］等人則通過設(shè)置精度矩陣服從Wishart分布，解決了這一問題，使得模型的適應(yīng)性更加廣泛。

估計和預(yù)測多變量金融時間序列的相關(guān)系數(shù)是進行最優(yōu)投資組合選擇，管理風(fēng)險，預(yù)測VaR等金融市場活動中一個極為重要的環(huán)節(jié)。為了估計和預(yù)測時變性的相關(guān)系數(shù)，學(xué)者們主要構(gòu)建了三種波動模型：多變量GARCH模型［14］、多變量隨機波動（Multivariate Stochastic Volatility，MSV）模型［15］和多變量Wishart波動模型［16-18］。在實際經(jīng)濟問題分析中，變量的協(xié)方差陣是不可直接觀察的，GARCH類模型利用前期已知信息進行測度，多變量隨機波動（Multivariate Stochastic Volatility，MSV）模型則將其視為內(nèi)在的潛狀態(tài)變量。Bau-wens和Laurent［19］指出GARCH類模型面臨維數(shù)詛咒難題；Asai等人［15］構(gòu)建多變量隨機波動模型的貝葉斯MCMC算法，減少了參數(shù)個數(shù)。然而，它們的波動方程僅僅考慮變量各自的時序影響，而忽略了變量之間的相互影響。變量之間的波動溢出效應(yīng)以均值方程為媒介間接體現(xiàn)的，而未直接反映在波動方程。Wishart波動模型［16-18］設(shè)置收益率序列的協(xié)方差陣服從Wishart過程，使得模型波動方程直接受到自身與其他變量的前期波動的影響，從而體現(xiàn)變量之間的波動溢出效應(yīng)，進而研究多個變量之間的波動與風(fēng)險特征。

為了估計原油市場與股票市場的動態(tài)相關(guān)系數(shù)，構(gòu)建了動態(tài)相關(guān)Wishart多變量隨機波動模型，設(shè)計Gibbs-MTM-ARMS混合抽樣的MCMC數(shù)值計算方法對模型參數(shù)進行估計；并且，利用國際原油期貨價格（WTI）、上證綜合指數(shù)（SSECI）和標(biāo)普500指數(shù)（S&P500）進行實證分析，刻畫原油市場與中、美股票市場各自的波動特征以及它們之間動態(tài)相依性。

2 貝葉斯動態(tài)多變量Wishart隨機波動模型的構(gòu)建

2.1 動態(tài)多變量Wishart隨機波動模型

多變量隨機波動模型可以看作是一種對多變量ARCH/GARCH類模型的擴展。它通過對資產(chǎn)收益波動加入一個隨機干擾項，使波動與收益的變化同時具有隨機不確定性。考慮常數(shù)相關(guān)多變量隨機波動（Constant Correlation MSV，CC-MSV）模型：

此處，yt=（y1t，y2t，…，ynt）′表示n種資產(chǎn)在t時刻的收益，相關(guān)系數(shù)矩陣為Pε=（ρij）n×n是常數(shù)矩陣（對角線元素為1，即ρii=1）。潛在的狀態(tài)變量ht=（h1t，h2t，…，hnt）′和μ=（μ1，μ2，…，μn）′均表示n×1維的向量，并且指定h0=μ。第一個等式描述收益的變化特征，稱為收益方程；第三個等式稱為波動方程，描述收益的波動特征，ηt～N（0， diag為隨機擾動項。CC-MSV模型能夠有效的估計協(xié)方差矩陣Var（yt）=Σt=（σij，t），刻畫潛在的波動特征，但是不能體現(xiàn)金融收益序列時變的相關(guān)系數(shù)特征，而在一定程度上限制了它的應(yīng)用。

為了得到收益序列的動態(tài)相關(guān)系數(shù)矩陣，Engle［20］提出了動態(tài)條件相關(guān)（Dynamic Conditional Correlation，DCC）模型，分析多變量GARCH模型中相關(guān)系數(shù)矩陣的動態(tài)行為。假設(shè)Dt= diag為收益向量在t時刻的協(xié)方差矩陣，則動態(tài)相關(guān)系數(shù)矩陣Pt=（ρt，ij）n×n為：

顯然，t時刻收益向量的協(xié)方差陣為Σt= DtPtDt。若協(xié)方差矩Σt是正定的，存在精度矩陣，且假設(shè)Φt服從自由度為v，尺度矩陣為St-1的Wishart分布，那么，構(gòu)建DC-W-MSV模型：

顯然，參數(shù)d和Θ={v，l0：d，M0：d｝決定了協(xié)方差矩陣的動態(tài)行為。Mi衡量了協(xié)方差矩陣的跨期敏感度，反映協(xié)方差受到前期的影響。v的大小決定Φt估計值的離散程度。滯后階數(shù)d表明Σt的持續(xù)效應(yīng)。li揭示了協(xié)方差跨期關(guān)系的強度，稱為持久參數(shù)。li越接近0，當(dāng)期波動對未來資產(chǎn)回報的波動影響越小；li越接近1，影響越大；li＞1表明波動結(jié)構(gòu)是非平穩(wěn)的。

2.2 模型參數(shù)的貝葉斯后驗推斷

多變量隨機波動模型中因為包含不可觀測的隱波動變量，導(dǎo)致難以表達(dá)精確的似然函數(shù)形式，使得估計潛在狀態(tài)變量和參數(shù)都極為困難。MCMC算法因為處理高維積分問題方面的獨特優(yōu)勢，成為時變波動模型中應(yīng)用最廣泛的方法。

按照前文的分析，分為兩個階段估計模型中的相關(guān)系數(shù)矩陣{Pt｝：

（Ⅱ）令γt=（γ1t，γ2t，…，γnt）′～N（0，Pt），同時結(jié)合模型（3），從而估計參數(shù)Θ和動態(tài)相關(guān)系數(shù)矩陣{Pt｝，t=1，2，…，T。

CC-MSV模型的MCMC算法設(shè)計參考文獻［15］，下文主要討論如何設(shè)計步驟（Ⅱ）的MCMC算法。記Ar：m={Ar，Ar+1，…，Am｝，給定Pt的條件下，γt服從均值向量為0和協(xié)方差為Pt的多元正態(tài)分布；而給定參數(shù)Θ和Φ1：t-1的條件下，Φt服從自由度為v，尺度矩陣為St-1的Wishart分布。那么，γ1：T和Φ1：T的聯(lián)合密度函數(shù)為邊緣似然函數(shù)與潛變量Φ1：T的密度函數(shù)之積，即：

那么，模型的似然函數(shù)為一個高維積分過程：

它的具體形式難以獲得，使得極大似然法難以實現(xiàn)。在貝葉斯理論的框架下，參數(shù)Θ與潛變量Φ1：T的后驗聯(lián)合分布π（Φ1：T，Θ｜γ1：T）可以利用MCMC抽樣方法進行模擬，因此可以將參數(shù)空間擴展為包括潛變量在內(nèi)的新的參數(shù)空間，從而估計模型參數(shù)和潛變量，因此問題轉(zhuǎn)化為從完全條件后驗分布π（Φ1：T｜Θ，γ1：T）與π（Θ｜Φ1：T，γ1：T）中抽取隨機樣本。從式（5）中可以看出，模型參數(shù)Θ只與π（Φ1：T｜Θ）有關(guān)，因此，可以獨立地推斷它的完全條件后驗密度函數(shù)；而推導(dǎo){Φt｝的完全條件后驗密度需要考慮相關(guān)系數(shù)矩陣Pt。

定理1 參數(shù)Θ和潛變量Φ1：T的后驗聯(lián)合分布正比于聯(lián)合密度函數(shù)π（γ1：T，Φ1：T｜Θ）與參數(shù)先驗分布π（Θ）之積，且

證明：參數(shù)Θ和潛變量Φ1：T的后驗聯(lián)合分布：

顯然，∫π（γ1：T，Φ1：T，Θ）dΦ1：TdΘ是與Θ和Φ1：T無關(guān)的常數(shù)因子，因此：

π（Φ1：T，Θ｜γ1：T）∝π（γ1：T，Φ1：T｜Θ）∝π（γ1：T，Φ1：T｜Θ）π（Θ）=π（γ1：T｜Φ1：T）π（Φ1：T｜Θ）π（Θ）

證畢。

為了進行貝葉斯分析，需要設(shè)置模型參數(shù)Θ的先驗分布。根據(jù)Philipov和Glickman［12］，Asai和Mc Aleer［13］等人的觀點，模型參數(shù)先驗分布的設(shè)置如下：

此處，Gamma表示伽瑪分布。因為Wishart分布中，自由度v必須大于尺度矩陣的維數(shù)n，所以設(shè)置v-n的先驗分布為伽瑪分布，從而保證v＞n。同樣，為了滿足矩陣參數(shù)Mi的對稱正定性，這里選擇逆Wishart分布作為先驗分布，當(dāng)d=0時，逆Wishart分布是共軛先驗。li設(shè)置為擴散先驗，使得對li可能的取值具有相等概率。

下面給出參數(shù)Θ和潛在狀態(tài)變量Φt的后驗完全條件密度函數(shù)的具體形式，進而設(shè)計MCMC抽樣算法，估計模型參數(shù)。

定理2 潛在狀態(tài)變量Φt的完全條件后驗密度函數(shù)為：

此處：

證明：根據(jù)定理1，可知Φt的完全條件后驗密度函數(shù)：

將上式化簡即得定理2的結(jié)論。證畢。

根據(jù)Φt的條件后驗密度函數(shù)進行抽樣之前，需要指定（Φ0，Φ-1，…，Φ-d）。采用MTM（Multiple-Try Metropolis）算法，估計參數(shù)Φt，假設(shè)當(dāng)前值為

（1）獨立從分布ξ～Wishartn（1/ΔΦt，ΔΦtI）生成（ξ1，ξ2，…，ξk），令：布，令：

為對稱函數(shù)，則：

（2）隨機地從 {Φt（1），Φt（2），…，Φt（k）｝按照比例w（Φt（i），抽取候 選值，并 令 W1=）。

參數(shù)集Θ={v，l0：d，M0：d｝含有d+2個實參數(shù)和d+1個矩陣參數(shù)，下面記M=M0：d。

推論1 參數(shù)空間Θ={v，l0：d，M｝的聯(lián)合后驗密度為：

證明：略（可參考定理2的證明過程）。

值得注意的是，表達(dá)式（10）中參數(shù){v，l0：d，M｝的先驗分布設(shè)置是相互獨立的。各參數(shù)的完全條件后驗分布密度函數(shù)如下：

定理3 參數(shù)M的完全條件后驗密度為：

證明：略（可參考定理2的證明過程）。

從表達(dá)式（11）可知，參數(shù)難以直接進行抽樣，為此，同樣采用MTM對M進行抽樣分析，詳細(xì)步驟參考參數(shù)Φt的抽樣方法。當(dāng)d=0時，M的完全條件后驗分布是自由度為vT+fM0，尺度矩陣為的逆Wishart分布，可直接利用Gibbs抽樣。

定理4 自由度v的完全條件后驗密度為：

π（v｜Φ1：T，l0：d，M）∝exp{（vΔ（v）+nv T log（v）-T log（Γn（v/2）+（αv-1）log（v-n））/2｝ （12）

此處：

證明：略（可參考定理2的證明過程）。

定理5 參數(shù){li｝，i=0，1，2，…，d的完全條件聯(lián)合后驗密度為：

證明：略（可參考定理2的證明過程）。

因為v和{li｝，i=0，1，2，…，d的條件后驗分布的復(fù)雜性，應(yīng)用ARMS（Adaptive Rejection Metropolis Sampling）抽樣方法，從v和{li｝，i=0，1，2，…，d的條件后驗分布進行抽樣。

2.3 Gibbs-MTM-ARMS混合抽樣算法設(shè)計

根據(jù)DCC-W-MSV模型中參數(shù)的條件后驗分布，設(shè)計Gibbs-MTM-ARMS混合抽樣的數(shù)值計算方法，以此估計模型參數(shù)及置信區(qū)間。記，為參數(shù)的初始值；假設(shè)為第r次抽樣的估計值，則第（r+1）次抽樣的迭代步驟如下：

（1）利用MCMC算法，估計CC-MSV模型中

的參數(shù)，并得到DC-W-MSV的樣本數(shù)據(jù)γt；

（6）重復(fù)步驟（2）～（5），直至參數(shù)的后驗條件分布為平穩(wěn)分布。

3 實證分析

3.1 指標(biāo)選取及統(tǒng)計特征分析

為了分析原油市場與股票市場的聯(lián)系，選取上證綜合指數(shù)（SSECI）、標(biāo)準(zhǔn)普爾500指數(shù)（S&P500）和紐約商業(yè)交易所推出的輕質(zhì)低硫原油（WTI）期貨價格這三個指標(biāo)，在2006年1月3日至2012年12月31日的日收盤價作為樣本。舍去節(jié)假日差異等因素造成的缺失數(shù)據(jù)，共有1704個樣本，其中，上證指數(shù)和S&P500指數(shù)來源于雅虎金融數(shù)據(jù)庫；原油期貨價格來源于EIA數(shù)據(jù)庫。數(shù)據(jù)的選取的跨度包含了美國“次貸危機”引發(fā)的全球金融危機和隨后發(fā)生“歐債危機”；在此期間，原油價格波動也尤為劇烈。收益序列：Rit=100（ln（pi，t）-ln（pi，t-1）），Rit表示收益率；pit表示指標(biāo)i在時刻t的收盤價。各指標(biāo)的時間趨勢圖和對數(shù)收益率時間圖如下：

圖1 各指數(shù)與對數(shù)收益率時間圖

圖2 TAR（1）模型的殘差時間圖

從圖1可知，各指標(biāo)的收益率體現(xiàn)波動率聚集的性質(zhì)，而且呈現(xiàn)明顯的自相關(guān)和異方差性質(zhì)。為了減弱趨勢項及自相關(guān)性的影響，根據(jù)Asai等人［13］的觀點，采用門限自回歸模型對收益序列進行處理，即：

那么，殘差υ∧it=Rit-E（Rit｜Ii，t-1）=Rit-就是CC-MSV模型中均值方程的收益序列，Ii，t-1表示過去的信息集。殘差估計如下：

從圖2可知，殘差出現(xiàn)大幅波動的時期是不同于對數(shù)收益率的，而且SSECI的殘差波動比S&P500和WTI更加頻繁。

3.2 CC-MSV模型的MCMC估計

下面分別利用SSECI與WIT、S&P500與WTI經(jīng)過處理的數(shù)據(jù)，構(gòu)建二元CC-MSV模型。迭代次數(shù)設(shè)置為30000次，舍去初始10000次不平穩(wěn)數(shù)據(jù)，同時為了減少鏈條自相關(guān)性，在剩余鏈條，每2個隨機數(shù)只保留一個，實際用于分析的數(shù)據(jù)為10000個。模型估計結(jié)果如表1所示：

表1 CC-MSV模型參數(shù)的MCMC估計

此處，Q表示分位點。從表1中可以看出，各參數(shù)MC誤差遠(yuǎn)小于標(biāo)準(zhǔn)差，且Geweke統(tǒng)計量（Z值）的絕對值小于1.96，在95%的置信水平下，可判斷迭代初的樣本均值與迭代末的樣本均值不存在顯著性的差異，抽樣獲得的Markov鏈?zhǔn)鞘諗康摹?/p>

圖3a-3f為CC-MSV模型中參數(shù)φ11、φ22和ρ12的后驗分布密度函數(shù)圖。從圖中可以看出，各參數(shù)的邊緣后驗分布核密度估計的曲線平滑，有明顯的單峰對稱特征，說明參數(shù)貝葉斯估計值的誤差非常小。CC-MSV模型中，波動率估計如圖4所示：

從波動率估計圖4可以看出，兩個模型對WTI波動率的估計無顯著差異；SSECI比S&P500和WTI的波動頻率更大，WTI的波動水平是最大的，S&P500與WTI的波動具有同期性，而SSECI表現(xiàn)不明顯，這基本與門限自回歸的殘差圖2吻合。在CC-MSV模型中，=0.1233＞0，= 0.3605＞0且，說明（1）S&P500與WTI的相關(guān)性要強于SSECI和WTI的相關(guān)性，美國對石油的依賴大于中國導(dǎo)致原油價格波動對美國市場的影響要高于中國。（2）WTI和SSECI、S&P500是正相關(guān)的，但是SSECI表現(xiàn)不明顯。從圖1和圖4中，可以看出S&P500與WTI的波動具有協(xié)同性，但是SSECI與WTI的波動則毫無規(guī)律。因此，從靜態(tài)角度研究SSECI與WTI之間的相關(guān)性是不全面的。

3.3 DC-W-MSV模型的MCMC估計

圖3 CC-MSV模型與DC-W-MSV模型中各參數(shù)的后驗密度圖

圖4 CC-MSV模型的波動率估計圖

利用2.2中的步驟進行濾波，獲得標(biāo)準(zhǔn)化處理后的數(shù)據(jù)γit，構(gòu)建DC-W-MSV模型（d=1），估計動態(tài)相關(guān)系數(shù)。模型參數(shù)的MCMC估計結(jié)果如表2所示：

表2中，各參數(shù)MC誤差遠(yuǎn)小于標(biāo)準(zhǔn)差，Geweke統(tǒng)計量的絕對值也小于1.96，因此，可判斷抽樣得到的鏈條是平穩(wěn)收斂的。各參數(shù)的后驗密度如圖3g-3p所示，從圖中可知，除了參數(shù)l（3）、v（4）和l（4），其余參數(shù)的邊緣后驗分布核密度估計的曲線平滑，有明顯的單峰對稱特征，說明參數(shù)貝葉斯估計值的誤差非常小。而l（3）、v（4）和l（4）的后驗密度圖出現(xiàn)雙峰，甚至三峰的情況，可能的是因為數(shù)據(jù)具有多峰的特征。

DC-W-MSV模型中各指標(biāo)的波動率（標(biāo)準(zhǔn)差）與動態(tài)相關(guān)系數(shù)圖如下：

表2 貝葉斯DC-W-MSV模型參數(shù)的MCMC估計

圖5 DC-W-MSV模型的波動率及動態(tài)相關(guān)系數(shù)圖

從圖5中可以看出，DC-W-MSV模型的波動率估計圖能夠很好地擬合的對數(shù)收益率時間圖1。（1）SSECI從2005年6月起開始上漲，并于2007年10月16日達(dá)到最高點6124，隨后一路下跌至2008年10月28日的最低1664點，之后觸底反彈。（2）2007年下半年美國次貸危機爆發(fā)，投資者開始對按揭證券的價值失去信心，引發(fā)流動性危機，直到2008年9月，這場金融危機開始失控，向全球蔓延，并導(dǎo)致多間相當(dāng)大型的金融機構(gòu)倒閉或被政府接管，從而引發(fā)了2007年至2009年的環(huán)球金融危機S&P500在此期間波動尤為劇烈。（3）WTI具有波動幅度大，起伏時間快的特征：國際原油價格自2002年開始進入單邊爆發(fā)式上漲階段，2006年超過60美元/桶，2008年7月更是達(dá)到145美元/桶的歷史最高點；此后，受到美國“次貸危機”沖擊，出現(xiàn)“斷崖式”急跌，12月份為30.28美元/桶，半年時間跌去80%；緊接著，各國政府聯(lián)合救市，全球經(jīng)濟緩慢復(fù)蘇，原油價格又節(jié)節(jié)攀升，2011年2月再次突破百元大關(guān)。

從動態(tài)相關(guān)系數(shù)圖5c中看出，在世界金融危機期間，無論是SSECI，還是S&P500，與WTI的相關(guān)系數(shù)都比較大的（絕對值），說明在此期間原油市場與股票市場的相依性是比較強的，但是正負(fù)難以判斷。金融危機后，SSECI與WTI的相關(guān)系數(shù)較小，而S&P500與WTI在呈現(xiàn)明顯的正相關(guān)性。從指數(shù)（價格）的時間圖1看出，在全球經(jīng)濟逐漸回暖，國內(nèi)經(jīng)濟增長減緩的環(huán)境下，SSECI一直表現(xiàn)低迷，而S&P500與WTI總體上表現(xiàn)上漲的趨勢。

表3分析了SSECI與WTI、S&P500與WTI的動態(tài)相關(guān)系數(shù)的統(tǒng)計特征。從表中可知，SSECI與WTI的相關(guān)性較小，而S&P500與WTI呈現(xiàn)正相關(guān)性。模型的動態(tài)相關(guān)系數(shù)分布圖如下：

從圖8可知，動態(tài)相關(guān)系數(shù)具有正態(tài)分布的性質(zhì)。因此，可以利用正態(tài)分布的性質(zhì)估計相關(guān)系數(shù)置信水平，進行風(fēng)險管理。

4 結(jié)語

針對多變量隨機波動模型中時變相關(guān)系數(shù)的估計問題，在CC-MSV模型的基礎(chǔ)上，設(shè)置精度矩陣服從Wishart分布，構(gòu)建了貝葉斯動態(tài)相關(guān)Wishart多變量隨機波動（DC-W-MSV）模型，設(shè)計了Gibbs-MTM-ARMS混合抽樣算法估計模型參數(shù)。利用國際原油價格（WTI）、上證綜合指數(shù)（SSECI）和標(biāo)普500指數(shù)（S&P500）進行原油市場與股票市場的動態(tài)相依性的實證分析。實證結(jié)果表明：模型能夠有效刻畫原油市場與股票市場的動態(tài)相依性。股票市場與原油市場的相關(guān)性在金融危機期間表現(xiàn)較強，并且難以判斷正負(fù)方向；金融危機后，中國股票市場與原油市場相關(guān)性極為微弱，而美國股票市場與原油市場則表現(xiàn)較強的正相關(guān)性。

表3 DC-W-MSV模型的動態(tài)相關(guān)系數(shù)分析

圖8 DC-W-MSV模型中動態(tài)相關(guān)系數(shù)分布圖

［1］Kilian L，Park C.The impact of oil price shocks on the US stock market［J］.International Economic Review，2009，50（4）：1267-1287.

［2］Wen Xiaoqian，Wei Yu，Huang Dengshi.Measuring contagion between energy market and stock market during financial crisis：A copula approach［J］.Energy Economics，2012，34（5）：1435-1446.

［3］Masih R，Peters S，De Mello L.Oil price volatility and stock price fluctuations in an emerging market：Evidence from South Korea［J］.Energy Economics，2011，33（5）：975-986.

［4］Arouri M E H，Rault C.Oil prices and stock markets in GCC countries：Empirical evidence from panel analysis［J］.International Journal of Finance&Economics，2012，17（3）：242-253.

［5］Nguyen C C，Bhatti M I.Copula model dependency between oil prices and stock markets：Evidence from China and Vietnam［J］.Journal of International Financial Markets，Institutions and Money，2012，22（4）：758-773.

［6］Kilian L.Exogenous oil supply shocks：How big are they and how much do they matter for the US economy？［J］.The Review of Economics and Statistics，2008，90（2）：216-240.

［7］劉紅忠，何文忠，李治平.A股市場上得“中石油魔咒”現(xiàn)象及其解釋［J］.財經(jīng)研究，2012，38（8）：109-121.

［8］Park J，Ratti R A.Oil price shocks and stock markets in the USand 13 European countries［J］.Energy Economics，2008，30（5）：2587-2608.

［9］Miller J I，Ratti R A.Crude oil and stock markets：Stability，instability，and bubbles［J］.Energy Economics，2009，31（4）：559-568.

［10］姬強，范英.次貸危機前后國際原油市場與中美股票市場間的協(xié)動性研究［J］.中國管理科學(xué)，2010，18（6）：42-50.

［11］Vo M.Oil and stock market volatility：A multivariate stochastic volatility perspective［J］.Energy Economics，2011，33（5）：956-965.

［12］Philipov A，Glickman M E.Multivariate stochastic volatility via Wishart processes［J］.Journal of Business&Economic Statistics，2006，24（3）：313-328.

［13］Asai M，Mc Aleer M.The structure of dynamic correlations in multivariate stochastic volatility models［J］. Journal of Econometrics，2009，150（2）：182-192.

［14］Bauwens L，Laurent S，Rombouts J V K.Multivariate GARCH models：A survey［J］.Journal of applied econometrics，2006，21（1）：79-109.

［15］Asai M，Mc Aleer M，Yu Jun.Multivariate stochastic volatility：A review［J］.Econometric Reviews，2006，25（2-3）：145-175.

［16］Gouriéroux C，Jasiak J，Sufana R.The Wishart autoregressive process of multivariate stochastic volatility［J］. Journal of Econometrics，2009，150（2）：167-181.

［17］Golosnoy V，Gribisch B，Liesenfeld R.The conditional autoregressive Wishart model for multivariate stock market volatility［J］.Journal of Econometrics，2012，167（1）：211-223.

［18］Jin Xin，Maheu J M.Modeling realized covariances and returns［J］.Journal of Financial Econometrics，2013，11（2）：335-369.

［19］Bauwens L，Laurent S，Rombouts J V K.Multivariate GARCH models：A survey［J］.Journal of applied econometrics，2006，21（1）：79-109.

［20］Engle R.Dynamic conditional correlation：A simple class of multivariate generalized autoregressive conditional heteroskedasticity models［J］.Journal of Business&Economic Statistics，2002，20（3）：339-350.

The Dynamic Dependence Between Crude Oil Market and Stock Market
Based on Bayesian Wishart Multivariate Stochastic Volatility model

ZHU Hui-ming，PENG Cheng，YOU Wan-hai，DENG Chao
（College of Business Administration，Hunan University，Changsha 410082，China）

Bayesian dynamic correlation Wishart volatility model is established in this paper to address estimation problem of time-varying coefficient matrix in multivariate stochastic volatility.In order to make the correlation coefficient matrix incorporated time-varying characteristics，the precision matrices in CC-MSV models are set to following the Wishart distribution.Based on the analysis of statistic structure of model and the selection of parameters prior，the Gibss-MTM-ARMS sampling algorithm method is utilized to estimate model parameters.The empirical research applies the data of Shanghai Composite Index，S&P500 and crude oil future price.The research results show that the correlation between crude oil market and stock market is strong during the financial crisis，but it is difficult to identify the direction.Further，after the financial crisis，the U.S.stock market is obviously positive correlated with crude oil market，while the correlation between China stock market and crude oil market is very weak.All those indicate that models used in this paper can effectively depict the dynamic dependency between crude oil market and stock market.

dynamic dependence；stochastic volatility；Bayesian analysis；Wishart distribution；Gibbs-MTM-ARMS algorithm

F224

A

1003-207（2014）07-0001-09

2013-07-16；

2014-02-19

國家自然科學(xué)基金項目（71221001，71031004，7171075）；教育部博士點基金項目（20110161110025）；湖南省自然科學(xué)基金項目（11JJ3090）

朱慧明（1966-），男（漢族），湖南湘潭人，湖南工商大學(xué)管理學(xué)院，教授，博士生導(dǎo)師；研究方向：貝葉斯計量經(jīng)濟模型.