基于隨機波動率假設(shè)的權(quán)證定價理論評述

摘要：隨機波動率（SV）模型在衍生品定價和風(fēng)險管理中的應(yīng)用開始發(fā)揮越來越重要的作用，但是，由于很難得到似然函數(shù)的閉型表達式，SV模型的參數(shù)估計問題嚴重限制了它在金融實踐領(lǐng)域的普及應(yīng)用。不過，近年來，學(xué)者們提出了許多旨在解決SV模型參數(shù)估計問題的有效且可行的新方法，大大推進了SV模型的應(yīng)用化進程。本文將在權(quán)證定價分析的框架內(nèi)，重點評述SV模型的參數(shù)估計方法，并從理論和實證的角度對它們的優(yōu)點和不足進行簡要評介和比較。

關(guān)鍵詞：SV模型；權(quán)證定價；估計方法

Abstract：Stochastic Volatility（SV）models are increasingly important in practical derivatives pricing applications， yet the empirical application of SV models in financial has been limited due to the difficulties involved in the evaluation of the likelihood function. However， recently there has been fundamental progress in this area due to the proposal of several new estimation methods that try to overcome this problem， being at the same time， empirically feasible. In this paper， we analysis the main estimators of the parameters and the volatility of univariate SV models under the framework of warrant pricing. We describe the main advantages and limitations of each of the methods both from the theoretical and empirical point of view.

Key Words：SV models，warrant pricing，estimating methods

中圖分類號：F830文獻標(biāo)識碼：A文章編號：1674-2265（2009）05-0026-05

一、引言

Black-Scholes期權(quán)定價理論在金融領(lǐng)域的里程碑地位早已在理論界和實踐行業(yè)得到廣泛認可。在實踐領(lǐng)域，正是因為Black-Scholes 期權(quán)定價公式的問世，期權(quán)類金融衍生產(chǎn)品市場在過去的四十五年中得到了迅猛發(fā)展，并已成為了整個金融市場的重要組成部分。與此同時，針對Black-Scholes定價理論的局限性而展開的學(xué)術(shù)研究層出不窮。和任何其他理論一樣，該定價理論的局限性主要表現(xiàn)在偏離現(xiàn)實的種種假設(shè)條件上。

由于波動率是金融產(chǎn)品定價和風(fēng)險管理中最重要的考慮因素之一，所以針對標(biāo)的資產(chǎn)價格常數(shù)波動率假設(shè)的研究最為普遍和深刻。標(biāo)的資產(chǎn)價格波動率是常數(shù)的假設(shè)和現(xiàn)實之間的沖突是非常明顯的，通常資產(chǎn)價格波動率表現(xiàn)出典型的時變特征，比如有厚尾、波動聚集性、杠桿效應(yīng)、長記憶性和持續(xù)性等。目前，對金融波動的以上特征進行刻畫的模型可以分兩大類：第一類是ARCH/GARCH類波動率模型，Engel（1982）、Bollerslev（1986）將當(dāng)期波動率表示成過去波動率和平方誤差項的確定性函數(shù)；第二類是由Taylor（1986）、Ghyselst（1996）和Shephard（1996）提出的隨機波動率模型（SV），它把波動率假設(shè)成一個服從某種隨機過程的潛在變量。

相比較而言，SV模型的理論基礎(chǔ)更加貼近現(xiàn)實，能夠包含更多的波動率時變特征，如果可以把它和權(quán)證定價結(jié)合起來，在理論上可以大大提高權(quán)證定價的精確合理程度。但是，隨機波動率模型的應(yīng)用卻因為其估計方法上的困難而沒有ARCH/GARCH類波動率模型那樣普遍。不過隨著計算機模擬技術(shù)和計量經(jīng)濟學(xué)的發(fā)展，有越來越多的隨機波動率模型的估計方法被提出來，大大促進了SV模型在實踐領(lǐng)域的應(yīng)用，進而也使人們能更精確合理地對衍生產(chǎn)品進行定價。本文主要就隨機波動率模型的參數(shù)估計方法進行評述，并在此基礎(chǔ)上對權(quán)證定價進行理論上的分析與探討。

二、隨機波動率模型的類型和特征

早期研究這一領(lǐng)域的有Clark(1973)、Eppis和Eppis(1976)、Tauchen和Pitts(1983)，采用混合分布假設(shè)來近似資產(chǎn)價格行為，為收益波動率的結(jié)構(gòu)化建模提供了一個理論框架，但是并沒有建立起能夠刻畫具體特征的波動率模型。直到1982年才由Taylor提出反映波動率群集的離散型SV模型，這一模型后來大多被稱為是基本SV模型。后來，SV模型經(jīng)歷了從基本模型到能反映更多波動率特征的擴展型模型、從離散時間模型到連續(xù)時間模型、從一元模型到多元模型的發(fā)展路徑。

（一）離散時間的SV模型

Taylor（1986）提出標(biāo)準(zhǔn)離散時間SV模型來刻畫金融市場中的有關(guān)波動問題，來解釋金融收益序列波動的波動率群集效應(yīng)。

基本的離散SV模型如下：

式中，是第t期的去掉均值后的收益， 是對數(shù)波動率； 是一個常數(shù)，表示平均的波動水平； 是結(jié)構(gòu)參數(shù)，反映了波動率的群集效應(yīng)；，

誤差過程 與互不相關(guān)，二者都是不可觀測的，

；參數(shù)是常數(shù)，度量了波動擾動的標(biāo)準(zhǔn)誤差。

基本的SV 模型中的許多假定在實踐中都不盡合理，修改這些假設(shè)，可得到各種擴展的SV模型。同時，單元的SV模型可以擴展到多元，更利于考察各收益序列之間的相關(guān)性。比如，Ruiz（1994），Harvey、Ruiz和Shepard（1994）通過假設(shè)收益率誤差項服從t分布或者GED分布（廣義誤差分布）而非正態(tài)分布，提出了肥尾類SV模型；Harvey和Shephard（1996）并不事先假設(shè)擾動項 和 兩者之間的相關(guān)系數(shù)為0，而是一個未知參數(shù)，從而提出了杠桿SV模型；Breidt、Crato與DeLima（1998），Granger（1999），Diebold 與Inoue（1999）和Granger 與Hyung（1999）把長記憶引入SV模型，提出了長記憶SV模型。

（二）連續(xù)時間隨機波動率模型

雖然在現(xiàn)實中，觀測值總是離散的，這給連續(xù)時間SV模型的估計帶來了困難，但是連續(xù)時間SV模型還是由于在理論上的明顯優(yōu)勢而幾乎和離散時間SV模型同時引起學(xué)者們的關(guān)注。隨著計算機技術(shù)的發(fā)展，獲取高頻數(shù)據(jù)開始變得越來越方便，計量經(jīng)濟學(xué)的發(fā)展為連續(xù)時間SV模型的估計提供了更多方法選擇的余地。著名的有Hull -White 模型、Stein-Stein 模型和Heston模型。

三、隨機波動率模型的主要統(tǒng)計推斷方法分析

由于似然函數(shù)精確表達式導(dǎo)致的參數(shù)估計難問題，SV模型在很長一段時間里同ARCH相比，并沒有在實踐領(lǐng)域表現(xiàn)出應(yīng)有的吸引力。然而，最近幾年，總體上關(guān)于非線性潛變量模型的估計已經(jīng)有了顯著的進步，特別是對SV模型的估計，學(xué)者們基于計量經(jīng)濟學(xué)和仿真模擬等技術(shù)，提出了許多可行的估計方法。根據(jù)各種估計方法的特點，可以把這些方法分成如下四大類：

（一）矩估計方法——廣義矩估計法（GMM）

Andersen和S?觟rensen（1993）以及Jacquier、Polson和Rossi（1994）針對該對數(shù)正態(tài)SV模型提出了GMM估計。主要思想是運用SV模型的平穩(wěn)性和遍歷性性質(zhì)，這些性質(zhì)使得樣本矩收斂于它們的無條件期望。對一個容量為T的樣本，令表示每一個樣本矩和其理論矩之間差距的一個向量，m是需要計算的矩的個數(shù)。廣義矩（GMM）估計量是通過最小準(zhǔn)則函數(shù)來構(gòu)造的：

其中 是一個的權(quán)重矩陣，它反映了賦予相應(yīng)的每一個矩的重要性。

GMM方法簡單易行，但是并不有效，因為它是通過犧牲有效性來換取相對簡單性的。它的主要問題是如何合理地確定權(quán)重矩陣。Greet Dphaene 和 Olivier Vergote (2004) 得到了針對標(biāo)準(zhǔn)離散SV模型的GMM估計中最優(yōu)權(quán)重矩陣以及最優(yōu)GMM估計量漸近協(xié)方差矩陣的閉型表達式。利用這兩個閉型表達式可以對GMM估計量和其他估計量的有效性進行比較，并且可以從大范圍的有效矩條件中選擇最優(yōu)的小范圍矩條件。另外，當(dāng)接近于1時，收斂于無條件矩會相當(dāng)緩慢，這時只有大樣本能夠解決這種問題。

Hansen和Scheinkman（1995）利用無窮小生成元來構(gòu)建隱含于穩(wěn)態(tài)馬爾可夫過程中的矩條件，進而通過這些矩條件來構(gòu)造GMM估計量。George Chacko 和 Luis M. Viciera（1999）基于隨機過程的條件特征函數(shù)，提出了改進連續(xù)時間過程估計的譜廣義矩估計方法。以上兩種改進后的GMM方法都為連續(xù)時間隨機波動率模型的參數(shù)估計提供了有效參考。最近，Whilla Chen 和Rohit S. Deo（2005）利用GMM對長記憶隨機波動率模型進行了估計，并對于樣本矩的收斂速度進行了討論。

在國內(nèi)，李傳樂、王美今（2006）采用模擬廣義矩估計方法（GMM），以上證綜合指數(shù)為樣本，考察了漲跌停板制度對滬市股票收益波動的影響，并將SV模型的實證結(jié)果與GARCH模型進行了比較，發(fā)現(xiàn)SV模型的估計更符合實際。

（二）基于極大似然函數(shù)的估計方法——準(zhǔn)最大似然估計法（QML）

由于SV模型是非線性的，這使基于線性狀態(tài)空間模型的估計是無效的。不過，Harvey、Ruiz和Shephard（1994）以及Nelson（1988）通過對標(biāo)準(zhǔn)SV模型中的式（1）的觀測值進行平方對數(shù)變換，把非線性（高斯）SV模型轉(zhuǎn)化為線性非高斯?fàn)顟B(tài)空間模型，并據(jù)此提出了準(zhǔn)最大似然估計法（QML）。他們假設(shè) 是正態(tài)分布，利用Kalman濾子得到的準(zhǔn)極大似然函數(shù)，并使之最大化來得到準(zhǔn)極大似然估計值。

傳統(tǒng)的極大似然估計方法假設(shè)不存在模型設(shè)定誤差，視觀測變量的密度函數(shù)為其真實的似然函數(shù)，QML估計量則允許模型設(shè)定誤差的存在，并不事先決定似然函數(shù)的具體形式。Ruiz（1994）指出，QML的估計量是一致的，并且是漸近正態(tài)的。但是，由于沒有基于 的真實似然函數(shù)，QML方法并不有效，而且近似效果會隨著波動率的波動率的減小而變得更差。Jacquier Polson和Rossi（1994）指出，較高的持續(xù)性 和較小的 會使QML估計量產(chǎn)生更大的MSE，所以直接把QML應(yīng)用于SV模型中會出現(xiàn)重大問題，因為金融數(shù)據(jù)通常都會有上述兩個特征。另外，當(dāng)觀測值為零時，QML方法還會遇到模型內(nèi)樣本點問題（inliers problem），因為此時，隨機波動率模型的線性空間模型轉(zhuǎn)化（對數(shù)轉(zhuǎn)化）便無法實現(xiàn)。Breidt和Carriquiry（1996）、Fuller （1996）以及 Bollerslev 與 Wright（2001） 通過線性轉(zhuǎn)化的方式給觀測值附加一個向上的平移項，改進了傳統(tǒng)QML方法，克服了模型內(nèi)樣本點問題。

盡管QML方法有許多局限性，但是它的靈活性是不容置疑的，可以方便地對肥尾類SV模型和多元SV模型進行參數(shù)估計。Andrew C. Harvey和Neil Shephard (1996）通過在構(gòu)成準(zhǔn)似然函數(shù)的條件密度函數(shù)中引入觀測變量的符號變量，把QML估計方法延伸到杠桿隨機波動率模型的參數(shù)估計。

在國內(nèi)，蘇衛(wèi)東、張世英（2002）通過引入禁忌遺傳算法對準(zhǔn)極大似然估計方法(QML)進行了改進，很好地克服了傳統(tǒng)QML方法在最大化準(zhǔn)似然函數(shù)時遇到的初點選擇問題，并且通過對滬深股市的實證分析表明，改進后的方法較傳統(tǒng)QML方法和GMM方法都有效。

（三）輔助模型估計法

1. 間接推斷法。Smith（1993），Gourieroux、 Monfort 和Renault（1993）通過引入一個向量（比如說）參數(shù)化的輔助模型，提出了關(guān)于復(fù)雜結(jié)構(gòu)模型基于模擬技術(shù)的間接推斷方法。這種方法論允許利用錯誤設(shè)定的輔助模型，因為基于正確設(shè)定的結(jié)構(gòu)模型的模擬過程，以及對觀測值的模擬路徑通過同一個輔助模型所發(fā)揮的校正作用，可以對模型設(shè)定錯誤所產(chǎn)生的偏差進行有效的調(diào)整。

在確定輔助模型后，首先利用樣本數(shù)據(jù)得到一個輔助模型的參數(shù)估計量 。接著基于模擬技術(shù)，從給定結(jié)構(gòu)參數(shù) 的結(jié)構(gòu)模型中獲得H次模擬樣本

，并借此產(chǎn)生 的新估計量。最后基于H個重復(fù)模擬和T個觀測值的樣本，運用加權(quán)矩陣 來選擇的一個間接估計量，使得二次距離最小化，即：

上面提到的結(jié)構(gòu)參數(shù) 是由鏈結(jié)函數(shù)給定的。正是因為需要對在數(shù)值最優(yōu)化算法中出現(xiàn)的各個計算鏈結(jié)函數(shù)，使得間接推斷這種方法的計算負擔(dān)較重，降低了實現(xiàn)效率。

2. 有效矩估計法（EMM）。有效矩估計的主要思想和間接推斷估計方法相一致，都是通過引入一個輔助模型來對隨機模型進行相關(guān)分析和估計，只是兩者的處理方法不同。Gourieroux、Monfort 和Renault（1993）研究表明，間接推斷估計量和EMM估計量是近似等同的，都是一致的且服從漸近正態(tài)分布。由于有效矩估計法避免了對每個模擬樣本都用得分函數(shù)來求得估計量 ，并且不用最小化類似于（3）式的二次距離，所以和間接推斷法相比，大大提高了計算效率。但是和間接推斷法一樣，EMM也會出現(xiàn)在模型設(shè)定檢驗中的過度識別問題。此外，EMM方法無法解決波動率的過濾和平滑問題，所以，還需要用附加的方法來對波動率進行估計和預(yù)測。

Andersen、Chung和S?觟renson(1999)利用蒙特卡羅方法對SV模型的EMM估計量的有限樣本特性進行了分析，發(fā)現(xiàn)EMM的處理效率低于MCMC方法，但要優(yōu)于GMM方法。在概念上，EMM方法可以被看成是介于GMM方法和直接對精確似然函數(shù)進行推斷的方法兩者之間。

Engle（1994）、Ghysels和Jasiak（1994）已經(jīng)利用EMM方法對連續(xù)時間的SV模型進行了參數(shù)估計。Gallant 和Long（1997）、Jiang 和 van der Sluis （2000）對一些離散SV模型進行了估計。在后來的發(fā)展中，EMM主要是被用來對連續(xù)時間SV模型進行參數(shù)估計，比如，Andersen和Lund（1997）、 Chernov和Ghysels（2000）以及 Dai 和 Singleton （2001）。

（四）貝葉斯MCMC方法（Bayesian analysis based on the Markov Chain Monte Carlo method）

之所以把MCMC方法列入完全模擬方法中，是因為這種方法直接把參數(shù)的后驗期望作為其估計值。在一般情況下，結(jié)構(gòu)參數(shù) 的后驗密度函數(shù)

以及其后驗期望函數(shù)的閉型表達式是無法得到的。在這種復(fù)雜的情況下，模擬技術(shù)就成了一種必要的方法，其中貝葉斯MCMC技術(shù)是非常重要的一類。MCMC方法是一種動態(tài)的抽樣方法，克服了傳統(tǒng)Monte Carlo 模擬的高維、靜態(tài)的缺陷， 提高了估算精度。

Jacquier、Polson和Rossi(1994)首次將貝葉斯方法用于對標(biāo)準(zhǔn)隨機波動率模型的分析。他們在貝葉斯框架中，利用基于數(shù)據(jù)增擴原理（Tanner和Wong，1987）的馬爾可夫鏈蒙特卡羅技術(shù)，獲取潛在變量（收益波動率）和所有參數(shù)之間的聯(lián)合后驗分布，進而從該聯(lián)合后驗分布中提取樣本，計算后驗矩估計值。這種方法的一個重要特點是把SV模型看成是由先驗密度函數(shù) 和兩個邊際分布

構(gòu)成的分層結(jié)構(gòu)模型，然后根據(jù)貝葉斯原理得到

的聯(lián)合后驗分布。從聯(lián)合后驗分布中抽取樣本的方法主要有Gibbs抽樣方法和Metropolis-Hasting抽樣方法。事實上，Gibbs抽樣方法是特殊的Metropolis-Ha-

sting抽樣方法，也可以看成是EM算法（Expection-

Maximum algorithm，期望值最大算法）的隨機版本。

Jacquier、Polson和Rossi（1994）提出的這種方法采用的是單步蒙特卡羅方法，在每個時點上，潛在變量從僅抽樣一次。這種方法，雖然在概念上很簡單，但是從模擬的角度來看，它并不是特別有效。Kim、Shephard和Chib（1998）提出了基于多步蒙特卡羅模擬的替代方法。這種方法可以對所有潛在波動率進行一次性抽樣，所以大大提高了處理效率。

借助于能夠?qū)崿F(xiàn)這種方法的Winbugs和OXmetrics等軟件包的發(fā)展，MCMC方法在對SV模型的實證分析中得到了廣泛應(yīng)用。Chib、Nardari和Shephard(2002) 提出一種有效的貝葉斯MCMC算法，對由包含肥尾t分布、模型局外變量和跳躍成份的擴展SV模型進行了參數(shù)估計。Abanto-Valle、Bandyop-

adhyay、Lachos和Enriquez（2008）基于一系列對稱正態(tài)分布的混合尺度，利用魯棒貝葉斯MCMC方法對肥尾SV進行參數(shù)估計。Jacquier、Polson和Rossi（2004）把基于MCMC的貝葉斯方法擴展到帶有肥尾特征和杠桿效應(yīng)的SV模型中。Jun Yu（2005）對Jacquier、Polson和Rossi（2004）提到的擴展SV模型做了適當(dāng)?shù)男拚⒗肕CMC方法對其進行了重新分析，改進的效果非常明顯。

貝葉斯MCMC方法不僅可以應(yīng)用于離散SV模型，而且也可以廣泛適用于連續(xù)SV模型。Roberts、Papaspiliopoulos和Dellaportas（2004），Griffin 和Steel（2003），Griffin和Steel（2008）以及Gander和Stephens（2007） 將貝葉斯MCMC方法應(yīng)用于波動率服從O-U過程的連續(xù)SV模型的參數(shù)估計和特征分析。Fruhwirth-Schnatter和Sogner（2008）將此方法用于Heston連續(xù)SV模型的相關(guān)分析。

在國內(nèi)，孟利鋒、張世英、何信（2004）采用貝葉斯框架下的MCMC方法對具有杠桿效應(yīng)的SV模型進行了參數(shù)估計， 并借助Wingugs對上海和深圳股市的指數(shù)收益時間序列進行實證分析，捕捉到顯著的杠桿效應(yīng)。周彥、張世英（2007）利用MCMC方法對連續(xù)時間SV模型進行了估計。

四、隨機波動率假設(shè)下的權(quán)證定價分析與展望

（一）權(quán)證定價分析

統(tǒng)計推斷估計方法的不斷發(fā)展大大促進了基于SV模型的衍生品定價研究的發(fā)展，進而也為在隨機波動率假設(shè)下的權(quán)證定價研究奠定了基礎(chǔ)。

第一，在某些SV模型下，比如Heston模型，是可以得到歐式期權(quán)價格的閉型解的。在這種情況下，只要通過有效可行的估計方法得到SV模型參數(shù)和波動率的估計值就可以求出歐式期權(quán)的價格，在做適當(dāng)修正后，便可以直接對權(quán)證進行定價。

第二，雖然對于大部分SV模型而言，歐式期權(quán)價格的閉型表達式是無法獲得的，但是，在得到SV模型的參數(shù)估計值之后，通過不太復(fù)雜的模擬方法，比如蒙特卡羅模擬方法，我們就可以得到性能良好的歐式期權(quán)價格的近似值。對于離散時間SV模型，可以直接進行模擬，而對于連續(xù)時間SV模型，則要先通過離散化來構(gòu)建模擬路徑，然后再進行模擬定價。目前的研究中廣泛采用的離散化方法主要是歐拉離散化方法和Milstein離散方法，雖然這些方法簡單易行，但會由于忽視掉波動率不能為負的直觀假設(shè)而產(chǎn)生較大的離散誤差。因此，如果能提出一種能確保波動率為非負且操作性強的離散化方法，必定能更合理地對SV模型進行仿真模擬和參數(shù)估計。

目前，SV模型和資產(chǎn)定價模型或風(fēng)險管理模型的結(jié)合，只是簡單地分為兩步，即先通過SV模型估計出波動率，然后將其納入到資產(chǎn)定價模型和風(fēng)險管理模型。但是，這樣會忽略掉模型與模型之間可能存在的可融合問題，比如，假設(shè)條件或者模型設(shè)定沖突等，所以有必要對SV模型和其他模型的嵌套問題進行深入研究。

第三，每種SV模型都有各自在模型設(shè)定上的特點，比如假設(shè)條件和模型結(jié)構(gòu)的不同，而以不同資產(chǎn)為標(biāo)的物的歐式期權(quán)又有各自不同的市場微觀結(jié)構(gòu)，且不同的市場會有不同的市場機制，所以是否存在一種對任何歐式期權(quán)定價都是最優(yōu)選擇的SV模型，在理論上是值得懷疑的。如果不存在，那么在不同的歐式期權(quán)或權(quán)證和SV模型之間是否存在一一對應(yīng)的最優(yōu)選擇關(guān)系呢？鑒于我國權(quán)證市場上權(quán)證種類稀少的實際情況，我們有條件用不同的SV模型對每只權(quán)證進行擬合分析，并探索其中可能存在的內(nèi)在規(guī)律。如果在我國權(quán)證市場不斷發(fā)展之后，權(quán)證的種類數(shù)量超出了逐一分析的限度，那么在理論邏輯上，可以依據(jù)一定的標(biāo)準(zhǔn)，比如標(biāo)的股票的規(guī)模或者是市場交易的結(jié)果和理論推測的結(jié)果偏差，對權(quán)證進行分類研究。

（二）SV模型估計方法研究的展望

通過對隨機波動率模型的參數(shù)估計方法的評述，我們發(fā)現(xiàn)該領(lǐng)域今后的研究需要在以下兩個方面有所突破：

第一，雖然估計方法之間關(guān)于實現(xiàn)難易程度、精確度和處理速度等方面的比較研究取得了一些成果，但是都還只是局部比較。而且即使是局部比較，所得到的結(jié)果也是不一致的，有的甚至是截然相反的。所以有必要對不同SV模型估計方法的優(yōu)點和缺點進行全面的比較研究。

第二，進一步簡化實現(xiàn)手段，開發(fā)出更具操作性的SV模型估計軟件。雖然近年來SV模型估計的軟件化取得不小成果，但是，相比于ARCH/GARCH模型的估計還是顯得過于復(fù)雜。這也是SV模型未能在實踐領(lǐng)域得到進一步普及的主要原因，所以在提高估計方法實現(xiàn)效率這方面還有很大的發(fā)展空間，其中相關(guān)開放性軟件的開發(fā)是最主要的一塊。

參考文獻：

[1]孟利鋒，張世英，何信.SV模型參數(shù)估計的經(jīng)驗特征函數(shù)方法[J].系統(tǒng)工程，2004，22（12）.

[2]周彥，張世英.基于MCMC方法的連續(xù)時間SV模型建模研究[J].工業(yè)工程，2007，10（1）.

[3]毛明來，陳通，徐正國.SV類模型體系探討[J].西北農(nóng)林科技大學(xué)學(xué)報（社會科學(xué)版），2006，6（4）.

[4]蔣祥林，王春峰.基于貝葉斯原理的隨機波動率模型分析及其應(yīng)用[J]. 系統(tǒng)工程，2005，23（10）.

[5]Chib.Marginal Likelihood from the Gibbs Output[J]. Journal of the American Statistical Association，1995，Vol.90，No.432，1313-1321.

[6]Chib，Ivan Jeliazkov.Marginal likelihood from the Metropolis-Hastings output[J]. Journal of the American Statistical Association，2001，96，453，270-281.

[7]Chib，S.， Greenberg，E. Markov chain Monte Carlo simulation methods in Econometrics[J]. Econometric Theory，1996，12，409-431.

[8]Davide Raggu.Adaptive MCMC methods for inference on affine stochastic volatility models with jumps[J]. Econometrics Journal，2004，volume 8，235-250.

[9]Hansen，L. P. and J. A. Scheinkman. Back to the future：Generating moment implications for continuous-time Markov processes[J]. Econometrics，1995，63，767-804.

[10]Harvey， A. C，E. Ruiz and E. Shepherd. Multivar-

-iate stochastic variance models. Rev. Econom. Stud. 1994，61， 247-264.

[11]Heston，S. L. A closed-form solution for options with stochastic volatility with applications to bond and currency options. Rev. Financ. Stud. 1993，6， 327-343.

（特約編輯 齊稚平）