GARCH模型kupiec回测检验在eviews里怎么做
答案:1 悬赏:80 手机版
解决时间 2021-01-19 04:52
- 提问者网友:战皆罪
- 2021-01-18 17:53
GARCH模型kupiec回测检验在eviews里怎么做
最佳答案
- 五星知识达人网友:何以畏孤独
- 2021-01-18 18:18
关于Kupiec的检验方法,可以使用R进行了重新编程作图检验,
参考代码如下:
n <- 984;
m <- seq(from=1, to=150, length=150);
p <- .05; # 例外发生的概率
conf <- .95; # chi^2 分布的置信水平
LR=-2*log(((1-p)^(n-m) )* (p^m))+2*log(((1-(m/n))^(n-m))*((m/n)^m));
critical_value <- qchisq(conf, 1);
plot(m,LR, type='l', col='blue', lwd=2, main='Kupiec Test (1995)', xlab='出现风险的次数m')
abline(h=critical_value, col='red', lwd=3)
LR_C <- which(LR<=critical_value)
abline(v=c(LR_C[1],LR_C[length(LR_C)]), lty='dashed')
LR_C
其中的n就是回顾测试的历史数据总数;m是这些数据中超过VaR的天数,这里的critical_value就是chi^2的临界值。
最后的LR_C输出的是满足VaR回测要求的m的所有可能取值。
参考代码如下:
n <- 984;
m <- seq(from=1, to=150, length=150);
p <- .05; # 例外发生的概率
conf <- .95; # chi^2 分布的置信水平
LR=-2*log(((1-p)^(n-m) )* (p^m))+2*log(((1-(m/n))^(n-m))*((m/n)^m));
critical_value <- qchisq(conf, 1);
plot(m,LR, type='l', col='blue', lwd=2, main='Kupiec Test (1995)', xlab='出现风险的次数m')
abline(h=critical_value, col='red', lwd=3)
LR_C <- which(LR<=critical_value)
abline(v=c(LR_C[1],LR_C[length(LR_C)]), lty='dashed')
LR_C
其中的n就是回顾测试的历史数据总数;m是这些数据中超过VaR的天数,这里的critical_value就是chi^2的临界值。
最后的LR_C输出的是满足VaR回测要求的m的所有可能取值。
我要举报
如以上问答信息为低俗、色情、不良、暴力、侵权、涉及违法等信息,可以点下面链接进行举报!
大家都在看
推荐资讯