Run JAGS model and get posterior point estimates with
uncertainty
mod <- run_jags_model(jagsdata = pkg_data, jagsparams = NULL,
n_iter = 80000, n_burnin = 10000, n_thin = 35)
Plot posterior point estimates with uncertainty
plots <- plot_estimates(jagsdata = pkg_data, model_output = mod)