Nothing
"irf.svecest" <-
function(x, impulse=NULL, response=NULL, n.ahead=10, ortho=TRUE, cumulative=FALSE, boot=TRUE, ci=0.95, runs=100, seed=NULL, ...){
if(!is(x, "svecest")){
stop("\nPlease provide an object of class 'svecest', generated by 'SVEC()'.\n")
}
y.names <- colnames(x$var@x)
if(is.null(impulse)){
impulse <- y.names
} else {
impulse <- as.vector(as.character(impulse))
if(any(!(impulse %in% y.names))) {
stop("\nPlease provide variables names in impulse\nthat are in the set of endogenous variables.\n")
}
impulse <- subset(y.names, subset = y.names %in% impulse)
}
if(is.null(response)){
response <- y.names
} else {
response <- as.vector(as.character(response))
if(any(!(response %in% y.names))){
stop("\nPlease provide variables names in response\nthat are in the set of endogenous variables.\n")
}
response <- subset(y.names, subset = y.names %in% response)
}
## Getting the irf
irs <- .irf(x = x, impulse = impulse, response = response, y.names = y.names, n.ahead = n.ahead, ortho = ortho, cumulative = cumulative)
## Bootstrapping
Lower <- NULL
Upper <- NULL
if(boot){
ci <- as.numeric(ci)
if((ci <= 0)|(ci >= 1)){
stop("\nPlease provide a number between 0 and 1 for the confidence interval.\n")
}
ci <- 1 - ci
if(!(is.null(seed))) set.seed(abs(as.integer(seed)))
BOOT <- .bootirfsvec(x = x, n.ahead = n.ahead, runs = runs, ortho = ortho, cumulative = cumulative, impulse = impulse, response = response, ci = ci, seed = seed, y.names = y.names)
Lower <- BOOT$Lower
Upper <- BOOT$Upper
names(Lower) <- impulse
names(Upper) <- impulse
}
result <- list(irf=irs, Lower=Lower, Upper=Upper, response=response, impulse=impulse, ortho=ortho, cumulative=cumulative, runs=runs, ci=ci, boot=boot, model = class(x))
class(result) <- "varirf"
return(result)
}
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.