[請益] GARCH-GPD Monte Carlo simulation

[naturalsmen (日日夜夜)]

各位大大好

小弟最近在看McNeil and Frey (2000)出的論文

Estimation of tail-related risk measures for heteroscedastic

financial time series: an extreme value approach

http://www.sciencedirect.com/science/article/pii/S0927539800000128


對最後他估計multiple day returns的算法有些疑問

我照著他的做法做了一遍 用他網頁提供的DAX index歷史資料

http://www.macs.hw.ac.uk/~mcneil/ftp/DAX.txt

但是10-day VaR的結果在0.99 quantile的地方差距很大

他的violations只有48筆 期望值是51筆

我做出來的violations卻有109筆@@

0.95 quantile跟作者的差距不大 但感覺是錯的

然後又去看了國內的人做的

http://140.123.5.6/deptfin/course_data/Data/Sup7.pdf

最後的GARCH-GPD的violations也跟期望值差距不大



我照McNeil(2000)上面的算法:

取1000天的窗格

1. 先用AR(1)-GARCH(1,1)預測未來10天DAX index negative log return ->

   得到1000筆殘差值然後標準化 跟 未來10天估計的條件平均數條件變異數

2. 取0.9 quantile和0.1 quantile作為左右尾的門檻值 ->

   左右尾分別fit一個GPD model -> 得到兩組(xi, beta)

3. 隨機從1.得到的1000筆標準化殘差取1筆 取後放回

4. 該殘差若>右尾門檻值 隨機從右尾fit得的GPD抽一個值 ->

   得右尾門檻值+右尾GPD抽到的值

5. 4.的判斷改成左尾 -> 得左尾門檻值-左尾GPD抽到的值

6. 4. 5.的條件不合則該殘差不變

7. 重複3.~6.

我共重複5000次 得到5000筆新的殘差值的混合分配

長這樣：http://imgur.com/veey8XC

上面是GARCH model得的標準化殘差分配

下面是新殘差的混合分配

...為啥沒有什麼變啊啊啊 我以為會出現厚尾之類的

不確定是不是那邊出錯

8. 從新殘差隨機抽10筆 模擬未來10天對數報酬 <-

   條件平均數 + 條件變異數*這10筆殘差

   然後加總

9. 8.重複1000次得到1000筆新的報酬率

10. 再fit一次GPD 門檻設0.9 quantile -> 得到t+1+...+t+10的VaR

然後重複上面的步驟到樣本資料結束

...就出事了@@

不知道版上有沒有大大曾經做過相關的東西

希望有大大能夠協助解惑

小弟願意以一頓飯報答

感謝!!!


--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 132.208.73.221
※ 文章網址: https://www.ptt.cc/bbs/Economics/M.1461284855.A.2FE.html
※ 編輯: naturalsmen (132.208.73.221), 04/22/2016 08:28:58
→ naturalsmen: 附帶一提 小弟寄信問過作者 但他說我敘述的步驟沒錯
→ naturalsmen: 但英文信我不太確定對方是不是真的懂我的意思
→ naturalsmen: 小弟是用R跑的 code:http://pastie.org/10807554
推 calvinhobbes: 在這理問這個很難得到解答吧... 但祝你好運
推 bearching: 要不要去統計版或R板問問看?
推 djching: 還真的不知道XD
→ naturalsmen: 那我轉去統計版看看 感謝~
※ naturalsmen:轉錄至看板 Statistics                               04/25 02:07
推 bearching: 個人直覺上是軟體差異的問題，例如黃裕烈就有比較VAR在
→ bearching: JMulti,Eviews,R的輸出結果的比較~
→ DIDIMIN: 各軟體給予的起始值差異使得估計結果不同