# \donttest{
# Access a subset of the fred_qd dataset
data <- fred_qd[, c("CPIAUCSL", "UNRATE", "FEDFUNDS")]
# Transform it to be stationary
data <- fred_transform(data, codes = c(5, 5, 1), lag = 4)
# Estimate a BVAR using one lag, default settings and very few draws
x <- bvar(data, lags = 1, n_draw = 1000L, n_burn = 200L, verbose = FALSE)
# Get coefficent values at the 10%, 50% and 90% quantiles
coef(x, conf_bands = 0.10)
# Only get the median of the variance-covariance matrix
vcov(x, conf_bands = 0.5)
# }
Run the code above in your browser using DataLab