Running a national level single-country model

Load your library


Load the data

cleaned_natdata <- get_data(national=TRUE, local=TRUE, mycountry="Nepal")

Get the JAGS model inputs and the cleaned data

pkg_data <- get_modelinputs(startyear=1990, endyear=2025.5, nsegments=12, 
                            raw_data = cleaned_natdata)

Run JAGS model and get posterior point estimates with uncertainty

mod <- run_jags_model(jagsdata = pkg_data, jagsparams = NULL, 
                      n_iter = 50000, n_burnin = 10000, n_thin = 20)

Plot posterior point estimates with uncertainty

plots <- plot_estimates(jagsdata = pkg_data, model_output = mod)