[問題] AP-GARCH simulation
在用AP-garch(asymmetric power garch)做模擬時遇到了一些問題
請教一下熟悉的版友
我用的模型是apgarch，除了變異數方程不同其餘和一般garch一致
變數代號:
條件標準差σ；殘差u；innovation是e
parameter包含alpha beta gamma omega和δ
變異數方程:
σ(t)^δ = omega + alpha*( |u(t-1)| - gamma*u(t-1))^δ + beta*σ(t-1)^δ
u =σ*e
我自己寫了code估出參數值，然後用來模擬未來十天期的風險值Value at Risk
方法是用樣本最後一天的殘差 u(1) 和條件標準差 σ(1)
反覆帶入求出十天後的條件標準差 σ(11)，再乘上e得出 u(11)
(e的部分假設innovation為t分配，利用估出來的自由度v作出random seed)
以上重複一萬次然後找出最底端的1%分配值
不過作出的結果和文獻的差很多
想請教一下這樣模擬的過程是正確的嗎?
模擬的程式碼如下:
e =trnd(nu,10,10000);
for j=1:10000
sigma(1,j) = sqrt(ht(end));
u(1,j) = data(end);
for i=1:10
sigma(i+1,j) = (omega+alpha*(abs(u(i,j))+gamma*u(i,j))^delta
+beta*sigma(i,j)^delta)^(1/delta);
u(i+1,j) = sigma(i+1,j)*e(i,j);
end
u_simulated(j) = u(11,j);
VaR = quantile(u_simulated,0.01);
end
end
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 111.255.243.215