我正在尝试实现有关广义预测误差方差分解的R代码
> fevd_generalised <- function(model, n.ahead=10, normalize=TRUE) {
if (class(model) != "varest") {
return("The model class is not varest!")
}
A <- Phi(model, n.ahead)
epsilon <- residuals(model)
Sigma <- t(epsilon)%*%epsilon / (model$obs - model$p*model$K)
gi <- array(0, dim(A))
sigmas <- sqrt(diag(Sigma))
for (j in 1:dim(A)[3]) {
gi[,,j] <- t( t( A[,,j]%*%Sigma ) / sqrt(sigmas) )
}
d <- array(0, dim(A)[c(2,3)])
for (j in 1:dim(d)[2]) {
d[,j] <- diag(A[,,j]%*%Sigma%*%t(A[,,j]))
}
num <- apply(gi^2,1:2,sum)
den <- c(apply(d,1,sum))
fevd <- num/den
if (normalize) {
return(fevd/apply(fevd, 1, sum))
} else {
return(fevd)
}
}
> Blockquote
然而，我想得到广义预测误差方差分解的整个矩阵超过10年，知道我在R软件中首次亮相！请您帮助我，并提前致谢