# NOT RUN {
library(BGVAR)
data(eerData)
model.ssvs <- bgvar(Data=eerData,W=W.trade0012,plag=1,draws=100,burnin=100,
prior="SSVS")
# example for class 'bgvar'
plot(model.ssvs, resp="EA")
# }
# NOT RUN {
# example for class 'bgvar.resid'
res <- residuals(model.ssvs)
plot(res, resp="EA")
# }
# NOT RUN {
# example for class 'bgvar.pred'
fcast <- predict(model.ssvs,n.ahead=8,save.store=TRUE)
plot(fcast, resp="US.Dp", Cut=20)
# }
# NOT RUN {
# example for class 'bgvar.irf'
shocks<-list();shocks$var="stir";shocks$cN<-"US";shocks$ident="chol";shocks$scal=-100
irf.chol.us.mp<-irf(model.ssvs,shock=shocks,n.ahead=24)
plot(irf.chol.us.mp,resp="US.y")
# }
# NOT RUN {
# example for class 'bgvar.fevd'
fevd.us.mp=fevd(irf.chol.us.mp,var.slct=c("US.Dp","EA.y"))
plot(fevd.us.mp, resp="US.Dp", k.max=10)
# }
Run the code above in your browser using DataLab