# upload data
data(us_fiscal_lsuw)
# specify the model and set seed
set.seed(123)
specification = specify_bsvar_msh$new(us_fiscal_lsuw)
# run the burn-in
burn_in = estimate(specification, 10)
# estimate the model
posterior = estimate(burn_in, 20)
# compute regime probabilities
rp = compute_regime_probabilities(posterior)
rp_summary = summary(rp)
# workflow with the pipe |>
############################################################
set.seed(123)
us_fiscal_lsuw |>
specify_bsvar_msh$new() |>
estimate(S = 10) |>
estimate(S = 20) |>
compute_regime_probabilities() |>
summary() -> rp_summary
Run the code above in your browser using DataLab