孫廣彪 鄭軍

【摘 要】 金融風(fēng)險(xiǎn)管理者最關(guān)注的是VaR技術(shù)的精度，目的就是通過構(gòu)造VaR的置信區(qū)間和期望損失預(yù)測(cè)來評(píng)估常見的動(dòng)態(tài)模型的精度，并量化誤差的大小。要構(gòu)造合適的置信區(qū)間，關(guān)鍵問題在于解決投資收益條件方差的動(dòng)態(tài)行為。文章通過構(gòu)造VaR的預(yù)測(cè)區(qū)間和期望損失來評(píng)估常見的動(dòng)態(tài)模型的精度，并量化誤差的大小，通過Monte Carlo方法來提供定量依據(jù)。它不僅可以作為金融機(jī)構(gòu)評(píng)估和管理個(gè)別資產(chǎn)或資產(chǎn)組合市場(chǎng)風(fēng)險(xiǎn)的工具，而且可以作為金融監(jiān)管部門監(jiān)管金融機(jī)構(gòu)和評(píng)估市場(chǎng)風(fēng)險(xiǎn)的重要手段。

【關(guān)鍵詞】 金融風(fēng)險(xiǎn)； 管理； VaR； 評(píng)估

一、關(guān)于VaR的介紹

風(fēng)險(xiǎn)價(jià)值模型（Value—at—Risk，VaR）是近年發(fā)展起來的用于測(cè)量和控制金融風(fēng)險(xiǎn)的量化模型。VaR技術(shù)越來越廣泛地用于投資組合風(fēng)險(xiǎn)計(jì)量、風(fēng)險(xiǎn)資本配置和績效評(píng)價(jià)。金融風(fēng)險(xiǎn)管理者當(dāng)然最關(guān)注VaR技術(shù)的精確度。

VaR從統(tǒng)計(jì)的意義上講，本身是個(gè)數(shù)字，是指面臨“正常”的市場(chǎng)波動(dòng)時(shí)“處于風(fēng)險(xiǎn)狀態(tài)的價(jià)值”。即在給定的置信水平（通常是1%或5%）和一定的持有期限內(nèi)（通常是一天或一周），預(yù)期的最大損失量（可以是絕對(duì)值，也可以是相對(duì)值）。期望損失（Expected Shortfall，ES）指位于超出VaR損失的條件期望。VaR技術(shù)在度量尾部風(fēng)險(xiǎn)時(shí)是無用的，誤差較大。

本文寫作目的有兩方面：一是評(píng)估計(jì)算VaR和ES中產(chǎn)生的潛在損失；二是通過VaR和ES的置信區(qū)間來量化誤差的嚴(yán)重性。

Jorion和Pritsker考慮過VaR風(fēng)險(xiǎn)值的估計(jì)。但構(gòu)造合適的VaR和ES的置信區(qū)間，關(guān)鍵問題在于解決投資收益條件方差的動(dòng)態(tài)行為，本文將運(yùn)用著名的GARCH模型來量化這些動(dòng)態(tài)行為。GARCH模型適用于波動(dòng)性的分析和預(yù)測(cè)，已經(jīng)成為金融風(fēng)險(xiǎn)管理中的主力，在GARCH—VaR模型和ES預(yù)測(cè)中很少產(chǎn)生參數(shù)估計(jì)錯(cuò)誤。

Pascual，Romo和Ruiz利用GARCH產(chǎn)生的時(shí)間序列獲得預(yù)測(cè)密度，提出了一個(gè)新的Bootstrap重采樣技術(shù)。本文發(fā)展了該重采樣技術(shù)，提出的重采樣技術(shù)相對(duì)來說更容易實(shí)現(xiàn)，并可擴(kuò)展到多元風(fēng)險(xiǎn)模型。

二、模型的構(gòu)建和風(fēng)險(xiǎn)措施

本文對(duì)一個(gè)給定的金融資產(chǎn)或投資組合建立每日損失（負(fù)回報(bào)）的動(dòng)態(tài)模型：

Lt=σtεt，t=1，…，T （1）

其中，εt是獨(dú)立同分布的，均值為0，方差為1，分布函數(shù)為G。這里考慮G為標(biāo)準(zhǔn)的Students t分布，自由度為d?！靓舤～t（d）為模擬波動(dòng)動(dòng)態(tài)，使用對(duì)稱的GARCH（1，1）模型，σ■■=ω+αL■■+βσ■■，其中α+β<1。需要說明的是，本文方法適用于更復(fù)雜的σ■■模型和其他εt分布。

本文關(guān)注損失分布的尾部情形，為此考慮兩個(gè)主流的風(fēng)險(xiǎn)措施：VaR技術(shù)和期望損失ES。前者簡單說就是損失分布的條件分位數(shù)，后者是超過VaR的那部分損失的期望。

已知T時(shí)期信息情況下，VaR以覆蓋率p度量T+1時(shí)期，用VaR■■表示這個(gè)正值：

Pr（LT+1>VaR■■■F■）=p （2）

這里FT表示在時(shí)間T可用的信息；p通常是一個(gè)小數(shù)字，如p=0.01或p=0.05。

類似，已知T時(shí)期信息情況下，ES以覆蓋率p度量T+1時(shí)期，用ES■■表示這個(gè)正值：

ES■■=E（LT+1■LT+1>VaR■■，F(xiàn)T） （3）

給定模型（1），可以得到VaR■■和ES■■的簡化的表達(dá)式：

VaR■■=σT+1G■■≡σT+1c1，p （4）

其中，G■■表示G的（1—p）分位數(shù)、標(biāo)準(zhǔn)損失分布εt=Lt/σt，σT+1是T+1時(shí)期的條件波動(dòng)。例如，若G是標(biāo)準(zhǔn)正態(tài)分布Φ，p=0.05，則G■■=Φ■■=1.645，從而有VaR■■=1.645σT+1。在一般情況下，當(dāng)ε～G，方程（4）表明，可以將VaR■■表示為σT+1和常數(shù)c1，p=G■■的乘積，其值取決于G和p。

類似地，給定模型（1）：

ES■■=σT+1E（ε■ε>G■■）

≡σT+1c2，p（5）

其中，ε是獨(dú)立同分布的隨機(jī)變量，均值為0，方差為1，分布函數(shù)為G。若ε～N（0，1），則對(duì)任意常數(shù)a，有E（ε■ε>a）=■，其中φ和Φ表示標(biāo)準(zhǔn)正態(tài)隨機(jī)變量的密度函數(shù)和分布函數(shù)。此時(shí)有ES■■=σT+1■和c2，p

≡■。若ε服從標(biāo)準(zhǔn)的Students t分布，自由度為d，則c2，p由不同公式給出。為描述這個(gè)公式，令td服從標(biāo)準(zhǔn)的Students t分布，自由度為d，Andreev和Kanto給出，對(duì)任何常數(shù)a，有E（td■td>a）=（1+■）■■，其中f和F表示td的概率密度和累積密度函數(shù)。于是有：

c2，p≡E（ε■ε>G■■）=1+（■G■■）2/d■■■

其中，G■■是ε分布的（1—p）分位數(shù)。特別地，G■■=■t■■，其中t■■是td分布的（1—p）分位數(shù)。

實(shí)際上，無法計(jì)算VaR■■和ES■■的真實(shí)值，因?yàn)樗鼈円蕾囉跀?shù)據(jù)生成的過程（也就是說，它們依賴于G和條件方差模型σ■■）。因此，需要估計(jì)它們，從而估計(jì)風(fēng)險(xiǎn)。本文的最終目標(biāo)就是要通過建立置信區(qū)間（或預(yù)測(cè)區(qū)間）來量化風(fēng)險(xiǎn)估計(jì)。

三、蒙特卡羅的結(jié)果

對(duì)GARCH方差預(yù)測(cè)，雖然可以考慮利用逼近思想例如三角方法來計(jì)算預(yù)測(cè)區(qū)間，但是在筆者所考慮的非參數(shù)模型中并不能直接計(jì)算VaR和ES的預(yù)測(cè)區(qū)間，甚至在參數(shù)情況下，這種三角逼近的方法也有可能比重采樣方法差。下面通過構(gòu)造VaR的預(yù)測(cè)區(qū)間和期望損失ES來評(píng)估常見的動(dòng)態(tài)模型的精度，并量化誤差的大小。通過蒙特卡羅（Monte Carlo）方法來提供定量依據(jù)，重點(diǎn)是GARCH模型下的時(shí)間變異風(fēng)險(xiǎn)的實(shí)際情況。

本文考慮三種不同的估計(jì)方法：

一是歷史模擬法（Historical Simulation）：下面稱為HS法，它利用損失的經(jīng)驗(yàn)分布來計(jì)算VaR和ES。

二是正態(tài)條件分布（Normal Conditional Distribution）：下面稱為Normal法。

三是極值理論中的Hill估計(jì)方法（The Hill's estimates）：下面稱為Hill法。

（一）GARCH損失情況

現(xiàn)在考慮GARCH—t（d）數(shù)據(jù)生成過程（DGP）的四種情況。設(shè)α=0.10，ω=（202/252）*（1—α—β）。因此，無條件波動(dòng)為每年20%。選擇如下四個(gè)參數(shù)：

1.Benchmark參數(shù)：β=0.80，d=8；

2.High Persistence參數(shù)：β=0.89，d=8；

3.Low Persistence參數(shù)：β=0.40，d=8；

4.Normal Distribution參數(shù)：β=0.80，d=500。

在Hill估計(jì)極值分布前，需要選擇一個(gè)截止點(diǎn)Tu，它定義了要估計(jì)的尾部指數(shù)參數(shù)極值的子樣本。為了挑選這個(gè)重要參數(shù)，從上述四個(gè)DGP參數(shù)來進(jìn)行Monte Carlo模擬實(shí)驗(yàn)，估計(jì)尾部指數(shù)的臨界值，并最終計(jì)算出為期一天的VaR和ES預(yù)測(cè)的結(jié)果偏差和均方根誤差（RMSE）。圖1和圖2顯示500和1 000的樣本估計(jì)總體的情況。

在這兩種情況下，選擇截?cái)帱c(diǎn)要符合極值子樣本的0.5%至10%的最大損失。圖1和圖2 的橫軸表示極值觀察量（超出500和1 000的數(shù)量），縱軸表示偏差和均方根誤差（RMSE）。從減少均方根誤差以實(shí)現(xiàn)偏差接近于零的角度來看，對(duì)這四個(gè)DGP，不管是VaR還是ES，2%的截止是合理的。

表1—表4給出了對(duì)應(yīng)上述四個(gè)DGP的蒙特卡羅結(jié)果。表1—表4的上半部分計(jì)算VaR，下半部分計(jì)算ES。表的左側(cè)給出了相應(yīng)風(fēng)險(xiǎn)的點(diǎn)估計(jì)，右側(cè)給出了基于Bootstrap法的VaR區(qū)間估計(jì)。對(duì)VaR和ES預(yù)測(cè)，考慮兩個(gè)估計(jì)樣本大小，T = {500，1 000}。在這些實(shí)驗(yàn)中，點(diǎn)估計(jì)進(jìn)行100 000次Monte Carlo計(jì)算。對(duì)Bootstrap置信區(qū)間，只進(jìn)行5 000次Monte Carlo計(jì)算，每次帶有999次Bootstrap復(fù)制。

（二）VaR與ES的點(diǎn)預(yù)測(cè)

本文主要目的是給出VaR和ES的有限樣本的置信區(qū)間，首先來考慮各種模型準(zhǔn)確預(yù)測(cè)風(fēng)險(xiǎn)的能力。VaR和ES的點(diǎn)預(yù)測(cè)結(jié)果見表1—表4左側(cè)的偏差和均方根誤差。

1.Benchmark參數(shù)情況

表1頂部給出樣本大小為T=500時(shí)，Benchmark DGP的VaR的結(jié)果。

首先考慮VaR估計(jì)的偏差，主要注意的是HS的向上偏差和Normal的向下偏差。后者是可以預(yù)料的，正態(tài)分布的尾部相對(duì)于1%的覆蓋率的確太小了，其他顯示的偏差都很小。對(duì)VaR估計(jì)的均方根誤差，HS最高且遠(yuǎn)遠(yuǎn)高于其他，而Hill要低得多。Normal的均方根誤差也比較低，如前所述，它有相當(dāng)大的偏差。當(dāng)增加樣本大小為1 000，表1顯示一般偏差變?。籋S仍向上偏差，Normal仍向下偏差。在均方根誤差方面，Hill表現(xiàn)非常好。

下面來看ES點(diǎn)預(yù)測(cè)。可以發(fā)現(xiàn)Normal有非常大的向下偏差。與VaR結(jié)果進(jìn)行比較，各種ES估計(jì)模型都有更大的均方根誤差。均方根誤差的增加部分緣于偏差的增加。

2.High Persistence參數(shù)情況

表2的上半部給出High Persistence的DGP的VaR結(jié)果?？梢钥吹?，結(jié)果類似于表1，但HS卻不同，HS現(xiàn)在偏差更大，均方根誤差已超過VaR真實(shí)值的50%，近似為2.71。Hill再次表現(xiàn)非常好。

表2的下半部給出High Persistence的DGP的ES結(jié)果?？梢钥吹?，結(jié)果類似于表1，但HS不同。HS的偏差和均方根誤差都很大。

3.Low Persistence參數(shù)情況

表3的上半部給出Low Persistence的DGP的VaR結(jié)果。HS現(xiàn)在表現(xiàn)不錯(cuò)，Hill再次表現(xiàn)非常好。

表3的下半部給出Low Persistence的DGP的ES結(jié)果。HS現(xiàn)在表現(xiàn)不錯(cuò)。

4.Conditional Normal情況

表4的上半部給出Normal的DGP的VaR結(jié)果。和表1中相比，偏差和均方根誤差都相當(dāng)小。HS仍舊表現(xiàn)較差，Normal表現(xiàn)非常好，Hill一直表現(xiàn)較好。

表4的下半部給出Normal的DGP的ES結(jié)果。和表1中相比，偏差和均方根誤差都相當(dāng)小。ES模型的偏差和均方根誤差均與表4中VaR結(jié)果的偏差和均方根誤差相近。

（三）VaR和ES的Bootstrap預(yù)測(cè)區(qū)間

上面討論VaR和ES點(diǎn)預(yù)測(cè)的精度?，F(xiàn)在討論VaR和ES的Bootstrap預(yù)測(cè)區(qū)間。也就是說，下面要評(píng)估各種方法的Bootstrap預(yù)測(cè)區(qū)間的能力。預(yù)測(cè)區(qū)間的結(jié)果在每個(gè)表的右側(cè)，給出了真實(shí)覆蓋率、置信區(qū)間的平均限值、置信區(qū)間平均寬度相對(duì)于VaR點(diǎn)預(yù)測(cè)的百分比。

1.Benchmark參數(shù)情況

再看表1頂部，相對(duì)于90%的覆蓋率而言，HS的覆蓋率非常低。此外，平均置信區(qū)間很寬。HS方法忽略了DGP在覆蓋率和寬度方面的動(dòng)態(tài)。

在VaR方面，Normal與HS一樣不好，但Normal有更小的寬度。不幸的是Normal覆蓋率又過小。當(dāng)增加樣本大小為1 000，HS覆蓋率變差。Normal覆蓋率變差，寬度較好。Hill有較好的覆蓋率和寬度。

下面來看ES的Bootstrap預(yù)測(cè)區(qū)間?？梢园l(fā)現(xiàn)HS有較低的覆蓋率，置信區(qū)間很寬。Hill有最好的覆蓋率但是又很寬。

從表1整體來看，Hill覆蓋率較好，而HS、Normal的覆蓋率差。Hill比較可取，一般來說它ES的覆蓋率比VaR的覆蓋率要差點(diǎn)。

2.High Persistence參數(shù)情況

先看表2上部的VaR結(jié)果。相較于表2中的Benchmark情況，HS覆蓋率變差（低），寬度變差（寬）。Normal覆蓋率更好。

表2底部的ES的結(jié)果。比較表3中的VaR和ES的結(jié)果，可以看到ES比VaR的覆蓋率通常是更糟的。與表1中Benchmark情況的ES比較，HS覆蓋率更糟糕且寬度更差。Normal仍然覆蓋率很差。Hill有更好的覆蓋率但區(qū)間更寬。

3.Low Persistence參數(shù)情況

表3上部右側(cè)為VaR結(jié)果。HS模型具有更好的覆蓋率和稍好的寬度。Normal覆蓋率變差但寬度更好。Hill覆蓋率相近且寬度比以前更好。

表3下部為ES結(jié)果。HS表現(xiàn)更好。

4.Conditional Normal情況

比較表4和表1中的VaR，HS覆蓋率變差但寬度比以前低。Normal覆蓋率和寬度均更好。Hill寬度比以前更好。Hill模型比較可行。

表4下部為ES結(jié)果。此時(shí)覆蓋率一般都變好。HS和Normal比其他表現(xiàn)要差。Normal覆蓋率很好。

（四）結(jié)果概述

根據(jù)表1—表4的結(jié)果，可以得出結(jié)論，對(duì)VaR和ES估計(jì)，HS不僅在點(diǎn)估計(jì)方面差，而且有很差的置信區(qū)間，即使在波動(dòng)的持久程度相對(duì)溫和時(shí)也是如此；Normal當(dāng)正態(tài)假設(shè)接近真實(shí)時(shí)當(dāng)然表現(xiàn)很好，否則就不行了。Hill的表現(xiàn)相當(dāng)不錯(cuò)，即使在有條件的正態(tài)分布。

一般情況下，計(jì)算ES的均方根誤差相比于VaR的均方根誤差要高得多（相對(duì)真實(shí)值）。因此，雖然在理論上ES會(huì)傳達(dá)更多的損失分布的尾部信息，但它難以準(zhǔn)確估計(jì)。當(dāng)討論兩個(gè)風(fēng)險(xiǎn)措施的優(yōu)點(diǎn)時(shí)要考慮這一點(diǎn)。

不幸的是，相較于VaR而言，很難在事前評(píng)估ES的準(zhǔn)確度。

當(dāng)Hill為VaR點(diǎn)估計(jì)90%的置信區(qū)間提供相當(dāng)可靠的覆蓋率，但ES相應(yīng)的覆蓋率通常遠(yuǎn)低于90%，因此是不可靠的。從一個(gè)保守風(fēng)險(xiǎn)管理者的角度來說，過度覆蓋不如低覆蓋。

四、結(jié)論

風(fēng)險(xiǎn)管理者和投資組合管理者躑躅于VaR估計(jì)的精度。量化點(diǎn)估計(jì)的精度可以使風(fēng)險(xiǎn)管理者做出更明智的決策。需要指出的是VaR 技術(shù)是一種非常復(fù)雜的方法，VaR計(jì)算方法的選用以及歷史數(shù)據(jù)的處理方式等都將顯著地影響到?jīng)Q策過程的復(fù)雜程度和最終的精度。本文采用VaR技術(shù)可以對(duì)選定的財(cái)務(wù)風(fēng)險(xiǎn)措施進(jìn)行預(yù)測(cè)與分析，在今后的研究與實(shí)踐中，還需進(jìn)一步將其他風(fēng)險(xiǎn)估計(jì)法的優(yōu)點(diǎn)綜合于VaR估計(jì)技術(shù)中?！?/p>

【參考文獻(xiàn)】

[1] A Review of Financial Market Events in Autumn 1998：BIS Committee on the Global Financial System.BIS，Basel，1999c，October.

[2] Jorion，P.Risk2：Measuring the Risk in Value—At—Risk[J].Financial Analysts Journal，1996，52：47—56.

[3] Pritsker，M.Towards Assessing the Magnitude of Value—at—Risk Errors Due to Erros in the Correlation Matrix[J].Financial Engineering News，1997（10/11）.

[4] Pascual，L.，Romo，J.，and Ruiz，E.Forecasting Returns and Volatilities in GARCH Processes Using the Bootstrap，manuscript，Departamento de Estadistica y Econometria，Universidad Carlos III de Madrid，2001.

[5] Andreev，A.，and Kanto，A.Conditional Value—at—Risk estimation using non—integer degrees of freedom of Students t—distribution，forthcoming[J].Journal of Risk，2004.