# upload data
data(us_fiscal_lsuw)
# specify the model and set seed
set.seed(123)
specification = specify_bsvar_mix$new(us_fiscal_lsuw, p = 1, M = 2)
# run the burn-in
burn_in = estimate(specification, 10)
# estimate the model
posterior = estimate(burn_in, 20)
# compute structural shocks' conditional standard deviations
csd = compute_conditional_sd(posterior)
# workflow with the pipe |>
############################################################
set.seed(123)
us_fiscal_lsuw |>
specify_bsvar_mix$new(p = 1, M = 2) |>
estimate(S = 10) |>
estimate(S = 20) |>
compute_conditional_sd() -> csd
Run the code above in your browser using DataLab