作者letibe (remember the fate)
站内Statistics
标题[问题] AP-GARCH simulation
时间Tue Dec 3 16:09:53 2013
在用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