# upload data
data(us_fiscal_lsuw)
# specify the model and set seed
specification = specify_bsvar_mix$new(us_fiscal_lsuw, M = 2)
set.seed(123)
# estimate the model
posterior = estimate(specification, 10)
# verify heteroskedasticity
sddr = verify_identification(posterior)
summary(sddr)
# workflow with the pipe |>
############################################################
set.seed(123)
us_fiscal_lsuw |>
specify_bsvar_mix$new(M = 2) |>
estimate(S = 10) |>
verify_identification() |>
summary() -> sddr_summary
Run the code above in your browser using DataLab