Я получил домашнее задание на тест LRT, но не могу понять, как получить тест LRT в проблеме VaR.
Мне очень жаль.
Вопрос в том, что «отклонить или нет VaR с помощью LRT»
VaR
vol.fit <- vol.garch
resid <- resid.garch
x <- sort(resid[-1])
library(fBasics)
plot(density(x), lwd = 3, main = "GARCH residuals")
lines(x, dnorm(x), col = "gray", lwd = 2) # normal fit
est.hyp <- hypFit(x, trace = F, doplot = F)@fit$estimate # hyperbolic
lines(x, dhyp(x, alpha = est.hyp[1], beta = est.hyp[2], delta = est.hyp[3], mu = est.hyp[4]), col = 2, lwd = 2)
est.nig <- nigFit(x, trace = F, doplot = F)@fit$estimate # normal inverse gaussian
lines(x, dnig(x, alpha = est.nig[1], beta = est.nig[2], delta = est.nig[3], mu = est.nig[4]), col = 4, lwd = 3, lty = 2)
p <- c(0.005, 0.01, 0.05)
q.n <- qnorm(p)
q.hyp <- qhyp(p, alpha = est.hyp[1], beta = est.hyp[2], delta = est.hyp[3], mu = est.hyp[4])
q.nig <- qnig(p, alpha = est.nig[1], beta = est.nig[2], delta = est.nig[3], mu = est.nig[4])
# 95% VaR
VaR.n <- vol.fit*q.n[3]
VaR.hyp <- vol.fit*q.hyp[3]
VaR.nig <- vol.fit*q.nig[3]
mean(r[-1] < VaR.n[-1])
Ожидаемый результат - VaR с Normal будет отклонен.