Nothing
"Psi.varest" <-
function(x, nstep=10, ...){
if(!is(x, "varest")){
stop("\nPlease provide an object of class 'varest', generated by 'VAR()'.\n")
}
nstep <- abs(as.integer(nstep))
Phi <- Phi(x, nstep = nstep)
Psi <- array(0, dim=dim(Phi))
params <- ncol(x$datamat[, -c(1:x$K)])
sigma.u <- crossprod(resid(x)) / (x$obs - params)
P <- t(chol(sigma.u))
dim3 <- dim(Phi)[3]
for(i in 1:dim3){
Psi[, , i] <- Phi[, , i] %*% P
}
return(Psi)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.