plot function

Graphical Summary of Output Created with bgvar

Plotting function for fitted values, residuals, predictions, impulse responses and forecast error variance decompositions created with the BGVAR package.

## S3 method for class 'bgvar'
plot(x, ..., resp = NULL, global = TRUE)

## S3 method for class 'bgvar.resid'
plot(x, ..., resp = NULL, global = TRUE)

## S3 method for class 'bgvar.pred'
plot(x, ..., resp = NULL, cut = 40, quantiles = c(0.1, 0.16, 0.5, 0.84, 0.9))

## S3 method for class 'bgvar.irf'
plot(
  x,
  ...,
  resp = NULL,
  shock = 1,
  quantiles = c(0.1, 0.16, 0.5, 0.84, 0.9),
  cumulative = FALSE
)

## S3 method for class 'bgvar.fevd'
plot(x, ..., resp, k.max = 10)

Arguments

• x: Either an object of class bgvar, bgvar.res, bgvar.irf, bgvar.predict or bgvar.fevd.
• ...: Additional arguments; set graphical parameters.
• resp: Specify either a specific variable, a specific country or a specific variable in a specific country which should be plotted. If set to NULL all countries is plotted.
• global: If global=TRUE global residuals are plotted, otherwise country residuals.
• cut: Length of series to be plotted before prediction begins.
• quantiles: Numeric vector with posterior quantiles. Default is set to plot median along with 68%/80% confidence intervals.
• shock: Specify the shock which should be plotted.
• cumulative: Default is set to FALSE. If cumulative=TRUE cumulative impulse response functions are plotted.
• k.max: plots the k series with the highest for the decomposition of resp.

Returns

No return value.

Examples

library(BGVAR)
data(testdata)
model.ssvs <- bgvar(Data=testdata,W=W.test,plag=1,draws=100,burnin=100,
                    prior="SSVS",eigen=TRUE)

# example for class 'bgvar'
plot(model.ssvs, resp=c("EA.y","US.Dp"))

# example for class 'bgvar.resid'
res <- residuals(model.ssvs)
plot(res, resp="EA.y")

# example for class 'bgvar.pred'
fcast <- predict(model.ssvs,n.ahead=8)
plot(fcast, resp="y", cut=20)

# example for class 'bgvar.irf'
shockinfo <- get_shockinfo("chol")
shockinfo$shock <- "US.stir"; shockinfo$scale <- +1
irf.chol<-irf(model.ssvs, n.ahead=24, shockinfo=shockinfo)
plot(irf.chol, resp="US")

# example for class 'bgvar.fevd'
fevd.us=fevd(irf.chol,var.slct=c("US.stir"))
plot(fevd.us, resp="US.stir", k.max=10)

Author(s)

Maximilian Boeck, Martin Feldkircher