# Access a subset of the usmacro_growth dataset
data <- usmacro_growth[,c("GDPC1", "CPIAUCSL", "FEDFUNDS")]
# Estimate a model
mod <- bvar(data, sv_keep = "all", quiet = TRUE)
# Extract posterior draws of VAR coefficients
bvar_coefs <- coef(mod)
# Compute summary statistics
summary_stats <- summary(bvar_coefs)
# Compute summary statistics of VAR coefficients without using coef()
summary_stats <- summary(mod$PHI)
# Test which list elements of 'mod' are of class 'bayesianVARs_draws'.
names(mod)[sapply(names(mod), function(x) inherits(mod[[x]], "bayesianVARs_draws"))]
Run the code above in your browser using DataLab