基于Copula函數(shù)的滬深300股指期貨相依風(fēng)險(xiǎn)研究

摘 要：本文對(duì)滬深300股指和股指期貨仿真交易收益率極端風(fēng)險(xiǎn)和相依關(guān)系進(jìn)行了研究，用DCC-GARCH模型描述了股指期貨和現(xiàn)貨之間動(dòng)態(tài)的條件相關(guān)系數(shù)，并以極值分布為邊際分布對(duì)四種常用的Copula函數(shù)進(jìn)行了擬合，發(fā)現(xiàn)Frank Copula的擬合效果最好，其次為Clayton Copula。在此基礎(chǔ)之上，對(duì)不同組合的VaR和CVaR進(jìn)行測(cè)度，發(fā)現(xiàn)投資組合比例與風(fēng)險(xiǎn)之間呈現(xiàn)“U”型特征，這也為股指期貨套期保值提供了一種新的研究方式。

關(guān)鍵詞：極值分布；Copula函數(shù)；套期保值

Abstract：This paper study extreme risk and the relationship between the China 300 stock index and stock trading simulation yield. DCC-GARCH model used to describe the dynamic conditions of the correlation coefficient between the stock index futures and spot goods，and extreme value distribution for the marginal distribution of the four commonly used Copula function fitting.We found Frank Copula is the best fit，followed by Clayton Copula. On this basis，with different combinations of VaR and CVaR for measure，we find“U”-type features between the proportion of portfolio and risk，which provides a new approach for hedge research.

Key Words：extreme value distribution，Copula function，hedge

中圖分類號(hào)：F830 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1674-2265（2012）12-0003-09

一、引言

金融資產(chǎn)的相關(guān)性塑造了這樣一種風(fēng)險(xiǎn)機(jī)制：它為風(fēng)險(xiǎn)的傳播提供了便捷的渠道，使之在與其相關(guān)聯(lián)的金融產(chǎn)品間蔓延，不斷地積累能量并觸發(fā)新的風(fēng)險(xiǎn)以致于金融風(fēng)波，如美國(guó)的次貸危機(jī)所引發(fā)的全球經(jīng)濟(jì)動(dòng)蕩。因此，如何度量金融資產(chǎn)的相關(guān)性歷來為學(xué)者所關(guān)注。

股指期貨不僅可以作為風(fēng)險(xiǎn)管理工具，有效分散和轉(zhuǎn)移企業(yè)的金融風(fēng)險(xiǎn)，而且還可以發(fā)揮套期保值的作用，為投資者規(guī)避風(fēng)險(xiǎn)。但股指現(xiàn)貨和期貨之間特殊的關(guān)聯(lián)性使得杠桿值進(jìn)一步增加，這又成為金融市場(chǎng)上新的風(fēng)險(xiǎn)來源。本文考察我國(guó)股指期貨正式推出之前期貨和現(xiàn)貨的動(dòng)態(tài)相關(guān)關(guān)系，分析其風(fēng)險(xiǎn)特征并進(jìn)行套期保值的研究。

目前，對(duì)金融資產(chǎn)收益時(shí)間序列的風(fēng)險(xiǎn)研究眾多，估計(jì)方法主要包括歷史模擬法、參數(shù)方法和非參數(shù)方法。歷史模擬是利用收益序列的經(jīng)驗(yàn)分布來近似真實(shí)分布，該方法雖然簡(jiǎn)單，但不能對(duì)過去觀察不到的數(shù)據(jù)進(jìn)行外推，在運(yùn)用中受到限制。參數(shù)方法假定資產(chǎn)收益服從某一特定過程，基于某一特定分布，得出的風(fēng)險(xiǎn)值多是對(duì)金融資產(chǎn)收益的總體風(fēng)險(xiǎn)的度量，并未很好地考慮到極端風(fēng)險(xiǎn)。由于金融市場(chǎng)上的收益率存在尖峰厚尾的特征，極端事件的發(fā)生雖然稀少，損失卻很巨大，人們最為關(guān)注的風(fēng)險(xiǎn)就是這種極端風(fēng)險(xiǎn)。極值分布作為一種非參數(shù)方法，不需設(shè)定模型，而是讓數(shù)據(jù)去選擇，相對(duì)于一般的橢圓分布，它更能捕捉到市場(chǎng)的極端風(fēng)險(xiǎn)；而極端風(fēng)險(xiǎn)間的相關(guān)是一種非線性相關(guān)，由于金融收益率具有的“波動(dòng)叢集性”的特征，使得一般的線性相關(guān)無法準(zhǔn)確描述金融資產(chǎn)間的關(guān)聯(lián)關(guān)系。Copula作為一種數(shù)學(xué)函數(shù)可以用來度量金融市場(chǎng)上的非線性相關(guān)，正確設(shè)定研究對(duì)象的邊緣分布是構(gòu)造Copula函數(shù)的關(guān)鍵。Copula這一概念最早由斯克拉（SKlar，1959）提出，最近幾年才發(fā)現(xiàn)它能應(yīng)用在金融風(fēng)險(xiǎn)管理中，對(duì)這一方法比較系統(tǒng)的理論介紹可以參見喬（Joe，1997）、內(nèi)爾森（Nelsen，1999）和恩布雷希茨等（Embrechts等，1999）的研究。隆吉（Longi，2000）應(yīng)用Copula方法對(duì)信用風(fēng)險(xiǎn)及信用衍生產(chǎn)品定價(jià)進(jìn)行研究。弗雷和麥克尼爾（Frey和McNeil，2003）使用Copula方法來分析尾部相關(guān)性，采用比較一般化的相關(guān)系數(shù)進(jìn)行定義。近年來使用Copula方法對(duì)組合風(fēng)險(xiǎn)進(jìn)行度量的研究還有波耶、戈塞爾和薩蒙（Bouye、Gaussel和Salmon，2002），隆吉和索林克（Longi和Solink，2001），格拉澤曼等（Glasserman等，2002），恩布雷希茨等（2003），羅森堡和舒曼（Rosenberg和Schuermann，2005）等。恩布雷希茨等（2003）和凱路比尼等（Cherubini等，2004）對(duì)Copula在金融中的相關(guān)應(yīng)用給出了比較全面的綜述。

本文在綜合現(xiàn)有研究的基礎(chǔ)上，以極值分布作為構(gòu)造Copula函數(shù)的邊緣分布，分析滬深300股指現(xiàn)貨和期貨收益率的極端風(fēng)險(xiǎn)，找出能夠精確度量極端風(fēng)險(xiǎn)的相依結(jié)構(gòu)，并對(duì)相依關(guān)系做出合理解釋。

二、模型介紹

（一）VaR及CVaR模型

1. VaR模型。VaR也稱在險(xiǎn)價(jià)值，是一定時(shí)期內(nèi)，在一定的置信度下，投資組合可能出現(xiàn)的最大損失，是國(guó)際上新近發(fā)展起來的一種卓有成效的風(fēng)險(xiǎn)量化技術(shù)。假設(shè) 代表某一金融資產(chǎn)的損失，其密度函數(shù)為 ，則VaR可以表示為：

（1）

當(dāng)密度函數(shù) 為連續(xù)函數(shù)時(shí)，也可以表示為： ，其中 為損失分布 的反函數(shù)。

2. CVaR的計(jì)算模型。CVaR是指損失超過VaR 的條件均值，也稱為期望短缺 CVaR（Expected Shortfall VaR）或尾部VaR（Tail VaR），其數(shù)學(xué)表達(dá)式為： ，X 表示損失（即負(fù)的收益）。CVaR彌補(bǔ)了 VaR 的缺陷。對(duì)于一個(gè)連續(xù)性的隨機(jī)變量，CVaR是超過VaR的損失的期望值。此時(shí)有：

（2）

因?yàn)閂aR往往會(huì)低估風(fēng)險(xiǎn)，大多數(shù)金融情景下的風(fēng)險(xiǎn)衡量，會(huì)帶來災(zāi)難性的后果。而這個(gè)問題可以被CVaR所解決。CVaR 滿足平移不變性、正齊次性、次可加性和單調(diào)性，因而是一致性的風(fēng)險(xiǎn)度量。CVaR 代表了超額損失的平均水平，反映了損失超過 VaR 閾值時(shí)可能遭受的平均潛在損失的大小，較之 VaR 更能體現(xiàn)潛在的風(fēng)險(xiǎn)價(jià)值。

（二）DCC-GARCH模型

DCC-GARCH模型由恩格爾（Engle，2002）提出，是在波勒斯勒夫（Bollerslev，1990）提出的常相關(guān)模型的基礎(chǔ)上發(fā)展出來的。設(shè) 為具有零均值的收益序列，則DCC-GARCH的具體模型為：

（3）

其中 是一個(gè)對(duì)角矩陣，對(duì)角線的元素即為各變量的條件方差， 為標(biāo)準(zhǔn)化殘差 的條件協(xié)方差矩陣，S為標(biāo)準(zhǔn)化殘差 的無條件協(xié)方差矩陣， 為條件相關(guān)系數(shù)矩陣，符號(hào)“ ”代表Hadamard乘積，即兩矩陣對(duì)應(yīng)元素相乘， ， 為所要估計(jì)的參數(shù)。

收益服從多元正態(tài)分布是使用極大似然估計(jì)的前提，若此前提假設(shè)不成立，只能使用擬極大似然估計(jì)（Quasi-Maximum Likelihood）。估計(jì)過程分為兩階段：（1）使用一元GARCH模型對(duì)各變量進(jìn)行估計(jì)。（2）使用前一步驟所得的標(biāo)準(zhǔn)化殘差來估計(jì)條件相關(guān)系數(shù)（條件協(xié)方差）。

（三）POT模型及極端風(fēng)險(xiǎn)測(cè)度

極值理論作為度量極端風(fēng)險(xiǎn)的一種方法，具有超越樣本數(shù)據(jù)的估計(jì)能力，并可以準(zhǔn)確地描述分布尾部的分位數(shù)。它主要包括兩類模型：BMM模型和POT模型。其中BMM模型是一種傳統(tǒng)的極值分析方法，主要用于處理具有明顯季節(jié)性數(shù)據(jù)的極值問題，POT模型是一種新型的模型，對(duì)數(shù)據(jù)要求的數(shù)量比較少，是目前經(jīng)常使用的一類極值模型。本文將采用POT模型進(jìn)行估計(jì)。

1. POT模型的理論基礎(chǔ)。假設(shè)序列 的分布函數(shù)為 ，定義 為隨機(jī)變量 超過閾值 的條件分布函數(shù)，它可以表示為：

（4）

根據(jù)條件概率公式我們可以得到：

（5）

定理（皮金斯，1975）：對(duì)于一大類分布 （幾乎包括所有的常用分布），存在一個(gè) ， 使得條件超限分布函數(shù) ：

（6）

當(dāng) 時(shí)， ；當(dāng) 時(shí)， 。函數(shù) 為廣義帕累托分布。

2. POT模型的建立和求解。本文采用最大概似估計(jì)法（Maximum Likelihood CVaRtimation）估計(jì)參數(shù) 與 ，求極值分布的方法如下：

第一步， 找出GPD概率密度函數(shù)：

（7）

第二步，求其對(duì)數(shù)似然函數(shù)：

（8）

第三步，確定閾值 ：

迪普伊（Dupuis，1998）給出了對(duì)閾值 的估計(jì)方法，一般有兩種：根據(jù)Hill圖或根據(jù)樣本的超限期望圖。本文采用樣本的超限期望圖確定閾值 ，令 ，樣本的超限期望函數(shù)定義為：

（9）

超限期望圖為點(diǎn) 構(gòu)成的曲線，選取充分大的 作為閾值，使得當(dāng) 時(shí)， 為近似線性函數(shù)：

（10）

第四步，估計(jì)結(jié)果：

當(dāng) 確定以后，利用 的值，根據(jù)公式（6）進(jìn)行最大似然估計(jì)得到 和 。同時(shí)，我們得到 的值中比閾值 大的個(gè)數(shù)，記為 ，根據(jù)公式（4），用頻率代替 的值，可以得到 的表達(dá)式：

（11）

3. 基于極值分布的VaR估計(jì)。由上式取反函數(shù)，可以得出VaR的估計(jì)式，對(duì)于給定某個(gè)置信水平 ，可以由 的分布函數(shù)公式（8）得到：

（12）

根據(jù)GPD的條件分布函數(shù)公式（5）可以得到：

（13）

（四）Copula函數(shù)及相依風(fēng)險(xiǎn)測(cè)度

1. 四種常用的Copula函數(shù)。尤里（Juri，2002）提出的Copula函數(shù)包含了變量尾部相關(guān)的全部信息，可以使分析者更加全面地了解變量之間的尾部相關(guān)關(guān)系。為了對(duì)變量之間的尾部相依關(guān)系進(jìn)行刻劃，我們?cè)趯?duì)滬深300股指現(xiàn)貨和期貨相依關(guān)系的研究中使用了四種Copula函數(shù)，分別是Gumbel、Guass、Clayton和Frank函數(shù)。

（1）Gauss Copula函數(shù)。

（14）

為一元標(biāo)準(zhǔn)正態(tài)分布， 。 事實(shí)上是邊際分布為標(biāo)準(zhǔn)正態(tài)分布的函數(shù)。諸多研究表明，在數(shù)據(jù)的中間部分，模擬 與實(shí)際數(shù)據(jù)吻合較好，但在尾部差別相當(dāng)明顯，當(dāng)要度量尾部風(fēng)險(xiǎn)相依性時(shí)，Gauss Copula是無能為力的。

（2）Gumbel Copula函數(shù)。Gumbel Copula 函數(shù)的分布函數(shù)和密度函數(shù)分別為：

（15）

（16）

其中 。當(dāng) 時(shí)， ，即隨機(jī)變量 、 獨(dú)立；當(dāng) 時(shí)，隨機(jī)變量 、 趨向于完全相關(guān)，且 ，即當(dāng) 時(shí)，Gumbel Copula函數(shù)趨向于Fréchet上界 。

Gumbel Copula的密度函數(shù)具有非對(duì)稱性，上尾高下尾低，其密度分布呈“J”字型。Gumbel Copula 函數(shù)上尾的相關(guān)性較強(qiáng)，可用于描述在上尾處具有較強(qiáng)相關(guān)關(guān)系的現(xiàn)象。收益率呈左偏分布時(shí)，多用Gumbel Copula去描述相關(guān)關(guān)系。

（3）Clayton Copula函數(shù)。Clayton Copula函數(shù)的分布函數(shù)與密度函數(shù)分別為：

（17）

（18）

其中 。當(dāng) 時(shí)， ，即隨機(jī)變量 、 趨向于獨(dú)立；當(dāng) 時(shí)，隨機(jī)變量 、 趨向于完全相關(guān)，且 ，即當(dāng) 時(shí)，Clayton Copula函數(shù)趨向于Fréchet上界 。

Clayton Copula的密度函數(shù)同樣具有非對(duì)稱性，上尾低下尾高，其密度分布呈“L”型。Clayton Copula 函數(shù)下尾的相關(guān)性較強(qiáng)，可用于描述在下尾處具有較強(qiáng)相關(guān)關(guān)系的現(xiàn)象。收益率呈右偏分布時(shí)，多用Clayton Copula去描述相關(guān)關(guān)系。

（4）Frank Copula函數(shù)。Frank Copula函數(shù)的分布函數(shù)和密度函數(shù)分別為：

（19）

（20）

其中 。 時(shí)，隨機(jī)變量 、 正相關(guān)， 時(shí)隨機(jī)變量 、 趨向于獨(dú)立， 時(shí)隨機(jī)變量 、 負(fù)相關(guān)。Frank Copula的密度函數(shù)具有對(duì)稱性，其密度分布呈“U”型。

2. Copula模型的估計(jì)。Copula模型的參數(shù)估計(jì)采用兩階段的極大似然估計(jì)方法。極大似然方程為：

（21）

為聯(lián)合密度函數(shù)， 是邊際分布 的參數(shù)， 是Copula分布密度函數(shù)的參數(shù)，所有的參數(shù)集為：

（22）

第一步，估計(jì)邊緣分布函數(shù)參數(shù)：

第二步，在已知第一步參數(shù)的情況下，對(duì) 進(jìn)行估計(jì)：

（23）

以上方法又稱邊際推斷函數(shù)法（IFM），參數(shù)的估計(jì)稱為邊際推斷估計(jì)（IFME），IFME和MLE的效率是非常接近的，具有一致性和漸近正態(tài)性。

3. 基于Copula的VaR算法。結(jié)合前面的極值分布，測(cè)定VaR的算法主要步驟如下：

第一步，估計(jì)滬深300股指現(xiàn)貨和期貨收益率的極值分布參數(shù)。

第二步，以GPD為邊際分布，估計(jì)四種Copula函數(shù)的參數(shù)。

第三步，根據(jù)一定準(zhǔn)則選取擬合最優(yōu)的Copula函數(shù)。

第四步，使用恰當(dāng)?shù)倪\(yùn)算法則從最優(yōu)的Copula函數(shù)產(chǎn)生pseudo隨機(jī)數(shù)。

第五步，用蒙特卡羅模擬，產(chǎn)生10000組模擬的資產(chǎn)組合估計(jì)VaR。

三、實(shí)證分析

（一）樣本選取及說明

滬深300指數(shù)期貨合約于2006年10月30日開始在中國(guó)金融期貨交易所進(jìn)行仿真交易，從2007年起進(jìn)行考察。本文選取2007年1月4日到2008年9月26日間期現(xiàn)指數(shù)收盤價(jià)作為分析對(duì)象。數(shù)據(jù)來源于中國(guó)金融期貨交易所（www.cffex.com.cn）。滬深300指數(shù)期貨推出的是當(dāng)月、下月和隨后兩個(gè)季月的合約，每個(gè)期貨合約都有到期日，為了克服期貨價(jià)格的不連續(xù)性，我們把每一天離到期日最近的合約品種的收盤價(jià)格連接起來。收益率定義為： ， 表示第t日指數(shù)收盤價(jià)，為了凸顯研究對(duì)象的數(shù)字特征，通過單調(diào)變換把收益率定義為 。滬深300股指現(xiàn)貨和期貨收益率分別用s和f表示。

（二）樣本描述統(tǒng)計(jì)

首先對(duì)滬深股指300現(xiàn)貨和期貨收益率進(jìn)行了描述統(tǒng)計(jì)，結(jié)果表明，滬深股指300現(xiàn)貨收益率的標(biāo)準(zhǔn)差為2.7080，明顯低于仿真股指期貨收益率的標(biāo)準(zhǔn)差，期貨收益率的峰度為7.5163，現(xiàn)貨收益率的峰度為3.9259，這說明期貨收益率的整體波動(dòng)要大一些。從J-B正態(tài)檢驗(yàn)來看，無論是期貨還是現(xiàn)貨的收益率都顯著拒絕原正態(tài)分布的假設(shè)。雖然現(xiàn)貨收益率的均值和中位數(shù)都大于期貨收益率，但現(xiàn)貨收益率的偏度為-0.3590，期貨的為1.1641。

采用DCC-GARCH模型對(duì)滬深300股指現(xiàn)貨和期貨收益率序列進(jìn)行擬合，圖1顯示兩個(gè)收益率序列動(dòng)態(tài)相關(guān)的特征示例，其中縱坐標(biāo)為動(dòng)態(tài)時(shí)變相關(guān)系數(shù)。兩指數(shù)收益序列動(dòng)態(tài)相關(guān)性非常強(qiáng)，相關(guān)系數(shù)最小為0.0293，最高為0.8909，平均為0.7077，標(biāo)準(zhǔn)差為0.0751。

（三）極值分布建模及極端風(fēng)險(xiǎn)測(cè)度

1. 極值分布的閾值選取。基于極值理論中的POT模型，我們需要確定恰當(dāng)?shù)拈撝?對(duì)超限分布進(jìn)行GPD擬合。迪穆謝爾（DuMouchel，1983）的研究表明，在 允許的情況下選取10%左右的數(shù)據(jù)作為極值數(shù)據(jù)組是比較合適的，否則可能會(huì)出現(xiàn)樣本內(nèi)過度擬合，樣本外不適用。圖2和圖3分別描繪了滬深300股指現(xiàn)貨和期貨收益率下尾和上尾以及各自的閾值變動(dòng)范圍。圖2為超限期望圖及下尾閾值范圍，左上為現(xiàn)貨收益率下尾，右上為期貨收益率下尾，左下為現(xiàn)貨收益率下尾閾值范圍，右下為期貨收益率下尾閾值范圍。圖3為滬深300股指現(xiàn)貨和期貨收益率上尾閾值范圍，左上為現(xiàn)貨收益率上尾，右上為期貨收益率上尾，左下為現(xiàn)貨收益率上尾閾值范圍，右下為期貨收益率上尾閾值范圍。以現(xiàn)貨收益率下尾分布為例：樣本的平均超限函數(shù)圖在-3.3附近近似直線，具有明顯的帕累托分布特征，并且我們結(jié)合閾值平緩變動(dòng)的范圍，經(jīng)過多次試驗(yàn)，最終確定當(dāng) 、 時(shí)能夠較好地刻劃下尾特征。

2. GPD參數(shù)估計(jì)結(jié)果。根據(jù)選取的閾值，進(jìn)行極大似然估計(jì)，通過尾部的擬合得到參數(shù)估計(jì)結(jié)果（見表1）。由圖4、圖5可以看出，下尾尾部的大部分點(diǎn)基本都落在圖中的曲線上。由于數(shù)據(jù)分布的原因，上尾擬合不可避免地存在一些偏差，總體來看GPD的擬合效果是比較理想的。

3. 基于極值分布的VaR和CVaR估計(jì)結(jié)果。為了凸顯極值分布對(duì)極端風(fēng)險(xiǎn)的刻劃能力，表2同時(shí)給出了基于正態(tài)分布和極值分布的VaR和CVaR估計(jì)值，由表2不難看出：期貨收益率的風(fēng)險(xiǎn)要明顯大于現(xiàn)貨收益率的風(fēng)險(xiǎn)，極值分布下的期貨風(fēng)險(xiǎn)值大概是現(xiàn)貨的兩倍；在0.975和0.99的置信水平下，基于極值分布的風(fēng)險(xiǎn)估計(jì)值要大于正態(tài)分布下的風(fēng)險(xiǎn)值，說明越是在高分位數(shù)，極值分布越能捕捉到尾部風(fēng)險(xiǎn)。

（四）Copula函數(shù)建模和相依風(fēng)險(xiǎn)測(cè)度

1. 經(jīng)驗(yàn)Copula函數(shù)。經(jīng)驗(yàn)Copula函數(shù)是基于非參數(shù)建模的方法對(duì)金融時(shí)間序列的真實(shí)分布做出較為精確的擬合，這種估計(jì)能夠比其他方法更接近數(shù)據(jù)的實(shí)際分布，因此，在將滬深300股指現(xiàn)貨和期貨收益率的邊際分布（極值分布）轉(zhuǎn)化為均勻分布進(jìn)行不同的Copula函數(shù)估計(jì)之前，有必要考察經(jīng)驗(yàn)Copula函數(shù)的一些性狀。

由于Gumbel Copula在上尾有著較高的相關(guān)性，即主要反映市場(chǎng)交易活躍時(shí)的相關(guān)，不妨稱之為“牛式相關(guān)”；Clayton Copula在下尾有著較高的相關(guān)性，即主要反映市場(chǎng)交易不活躍時(shí)的相關(guān)，稱之為“熊式相關(guān)”；Frank Copula是對(duì)稱分布，屬于“對(duì)稱相關(guān)”。由經(jīng)驗(yàn)Copula概率密度的等高線可以看出，該分布兩端概率的密集程度大體相當(dāng)，或下尾比上尾處的概率稍微密集了些，因此有理由猜想利用對(duì)稱分布的Frank Copula或Clayton Copula可以較為真實(shí)地反映這一相關(guān)結(jié)構(gòu)。

圖6：經(jīng)驗(yàn)Copula密度函數(shù)和概率密度的等高線

由表3可以看出，F(xiàn)rank Copula的極大似然值最大，并且AIC和BIC是最小的，因此選取Frank Copula是比較恰當(dāng)?shù)摹８鶕?jù)以上估計(jì)結(jié)果對(duì)函數(shù)進(jìn)行擬合，圖7給出經(jīng)驗(yàn)Copula和四種Copula的擬合圖的等高線，左上為Gumbel Copula，右上為Gauss Copula，左下為Frank Copula，右下為Clayton Copula。不難看出，F(xiàn)rank Copula的擬合效果最好，其次為Clayton Copula。為了進(jìn)一步驗(yàn)證這一結(jié)論的可靠性，利用Frank Copula函數(shù)對(duì)滬深300股指現(xiàn)貨和期貨收益率序列進(jìn)行模擬，對(duì)照真實(shí)收益率，可以推知基于Frank Copula對(duì)收益率分布的擬合非常理想（見圖8）。

3. 基于Frank Copula的VaR和CVaR估計(jì)結(jié)果。

現(xiàn)在考慮極端事件的相依風(fēng)險(xiǎn)的測(cè)度問題，s和f分別表示滬深300股指現(xiàn)貨和期貨的收益率，設(shè) 和 分別為股指現(xiàn)貨和期貨的投資比例，這兩種資產(chǎn)的組合投資的對(duì)數(shù)收益率可以表示為：

選取極值分布作為邊際分布，F(xiàn)rank Copula為連接函數(shù)，模擬出聯(lián)合分布F的隨機(jī)數(shù)，對(duì)于給定的組合投資收益率R，在給定的置信水平q下，用蒙特卡羅模擬，產(chǎn)生10000組模擬的資產(chǎn)組合估計(jì)VaR。

以上給出了五種組合下的VaR和CVaR值，不難發(fā)現(xiàn)，隨著現(xiàn)貨投資比例的增加，風(fēng)險(xiǎn)值不斷降低，然后又隨之增加，大致呈現(xiàn)“U”型特征。以上面數(shù)據(jù)為例，在0.975和0.99的置信水平下，現(xiàn)貨的投資比例為0.7、期貨的投資比例為0.3時(shí)，風(fēng)險(xiǎn)值最低，不斷地縮短步長(zhǎng)，總可以找到風(fēng)險(xiǎn)最低的投資組合，該思路可以運(yùn)用到套利上面去，從而為研究股指期貨的最優(yōu)套期保值提供了新的范式，這也是本文下一步的研究重點(diǎn)。

四、結(jié)論

第一，通過DCC-GARCH模型對(duì)股指仿真期貨和現(xiàn)貨收益率的條件相關(guān)的動(dòng)態(tài)描述可以看出，滬深300股指期貨和現(xiàn)貨收益率的變動(dòng)趨勢(shì)大致相同，具有較強(qiáng)的動(dòng)態(tài)相關(guān)關(guān)系，投資者可以利用現(xiàn)貨和期貨收益波動(dòng)的相關(guān)性進(jìn)行有效的套期保值，規(guī)避系統(tǒng)性風(fēng)險(xiǎn)；而監(jiān)管者需要更加密切地關(guān)注證券市場(chǎng)的動(dòng)態(tài)變化，防止股指期貨推出后市場(chǎng)的惡意操縱。

第二，以極值分布為邊際分布構(gòu)造的Copula函數(shù)極大限度地捕捉到了極端風(fēng)險(xiǎn)的相依關(guān)系。進(jìn)一步的實(shí)證表明：期貨風(fēng)險(xiǎn)要明顯大于現(xiàn)貨風(fēng)險(xiǎn)，極值分布下的期貨風(fēng)險(xiǎn)值大概是現(xiàn)貨的兩倍。滬深300股指現(xiàn)貨和期貨收益率的尾部相關(guān)結(jié)構(gòu)可以用Frank函數(shù)進(jìn)行描述，即無論市場(chǎng)瘋狂上漲還是急劇下跌過程中，股指現(xiàn)貨和期貨收益率之間存在著同時(shí)相互影響、相互加強(qiáng)的雙邊風(fēng)險(xiǎn)關(guān)聯(lián)關(guān)系，同時(shí)，尾部相關(guān)結(jié)構(gòu)也可以用Clayton函數(shù)進(jìn)行刻劃，這一結(jié)果意味著當(dāng)資本市場(chǎng)出現(xiàn)持續(xù)下跌、劇烈波動(dòng)等極端風(fēng)險(xiǎn)事件時(shí)，股指期貨的風(fēng)險(xiǎn)關(guān)聯(lián)性相互加強(qiáng)。因此，應(yīng)設(shè)定恰當(dāng)?shù)姆€(wěn)定機(jī)制，盡可能不干擾市場(chǎng)正常波動(dòng)和交易情況，特別是對(duì)于預(yù)防股災(zāi)等極端股市危機(jī)是極其必要的。

第三，在利用Frank Copula進(jìn)行極端風(fēng)險(xiǎn)測(cè)度的同時(shí)，我們發(fā)現(xiàn)，在投資組合比例的步長(zhǎng)為0.2時(shí)，隨著現(xiàn)貨投資比例的增加，風(fēng)險(xiǎn)值先是不斷降低后又逐漸增加，組合投資風(fēng)險(xiǎn)大致呈“U”型分布，縮短步長(zhǎng)，可以得到風(fēng)險(xiǎn)最低時(shí)的投資組合，至少可以得到對(duì)應(yīng)風(fēng)險(xiǎn)局部最小值的投資組合，其實(shí)這是一種風(fēng)險(xiǎn)對(duì)沖的思想，從而為最優(yōu)套期保值比率的研究提供了一種新的思路。

參考文獻(xiàn)：

[1]王志誠(chéng)，周春生.金融風(fēng)險(xiǎn)管理研究進(jìn)展：國(guó)際文獻(xiàn)綜述[J].管理世界，2006，（4）.

[2]馬林.金融風(fēng)險(xiǎn)度量及其有效前沿[D].華中科技大學(xué)碩士學(xué)位論文，2006年1月.

[3]Danielsson J.，CG.de Vries.1997.Value at Risk and Extreme Returns，London School of Economics[M]，F(xiàn)inancial Markets Group Discussion Paper 1997（272-273）.

[4]韋艷華，張世英，郭焱.金融市場(chǎng)相關(guān)程度與相關(guān)模式的研究[J].系統(tǒng)工程報(bào).2004，19（4）.

[5]Sklar，A.，1959，“Fonctions de Repartition a n Dimensionset Leurs Marges”，Publications de l'Institut de Statistiquede l' Universite de Paris，Vol.8， pp.229～231.

[6]Joe，H.，1997，“Multivariate Models and Dependence Concepts”，Chapman Hall.

[7]Nelsen，R.，1999，“An Introduction to Copulas”，Springer，New York.

[8]Embrechts，P.，Lindskog，F(xiàn).and McNeil，A.，2003，“Modeling Dependence with Copulas and Applications to Risk Management”，In：Rachev，S.（Ed.），Handbook of Heavy Tailed Distributions in Finance，Elsevier，Rotterdam， pp.329～384.

[9]Longin，F(xiàn).，2000，“From Value at Risk to Stress Testing：the Extreme Value Approach”， Journal of Banking and Finance，Vol.24，pp.1097～1130.

[10]Frey，R.and McNeil，A，.2003，“Modeling Dependent Defaults”，Journal of Risk，Vol.6，No.1，pp.59～92.

[11]Bouye，E.，Gaussel，N. and Salmon，M.，2002，“Investigating Dynamic Dependence Using Copula”， Manuscript， Financial Econometrics Research Center.

[12]Longin，F(xiàn). and Solink，B.，2001，“Extreme Correlation of International Equity Markets”，Journal of Finance， Vol.56，pp.649～676.

[13]Glasserman，P.，Heidelberger，P.and Shahabuddin， P.，2002，“Portfolio Value-at-risk with Heavy-tailed Risk Factors”， Mathematical Finance，Vol.12，pp.239～270.

[14]Embrechts，P.，Lindskog，F(xiàn).and McNeil，A.，2003，“Modeling Dependence with Copulas and Applications to Risk Management”，In：Rachev，S.（Ed.），Handbook of Heavy Tailed Distributions in Finance，Elsevier，Rotterdam， pp.329～384.

[15]Rosenberg，J.and Schuermann，T.，2005，“A General Approach to Integrated Risk Management with Skewed， Fat-tailed Risk”，F(xiàn)RB of New York Staff Report No.185.

[16]Embrechts，P.，Lindskog，F(xiàn).and McNeil，A.，2003，“Modeling Dependence with Copulas and Applications to Risk Management”，In：Rachev，S.（Ed.），Handbook of Heavy Tailed Distributions in Finance，Elsevier，Rotterdam， pp.329～384.

[17]鄭振龍，王保合.基于極值理論的風(fēng)險(xiǎn)價(jià)值度量[J].金融學(xué)季刊.2005，1（1）

[18]McNeil.A．J，F(xiàn)rey.R.2000.Estimation of Tail-Rlated Risk Measures for Heteroscedastic Financial Time Series：An Extreme Value Approach[J].Journal of Empirical Finance 2000，7：27（1-300）．

[19]DuMouchel W.M.1983.Estimating the Stable Index _ in Order to Measure Tail Thickness： A Critique[M]. Annals of Statistics 11.1983.（1019-1031）.

[20]李竹渝，魯萬(wàn)波，龔金國(guó).經(jīng)濟(jì)金融計(jì)量學(xué)中的非參數(shù)估計(jì)技術(shù)[M].科學(xué)出版社.2007.

[21]歐陽(yáng)資生.信用風(fēng)險(xiǎn)相依模型及其應(yīng)用研究[M].知識(shí)產(chǎn)權(quán)出版社.2008年2月.

（特約編輯 齊稚平；校對(duì) GX）