Skip to contents

Running a subnational level single-country model

subnational_local_models.Rmd

Load your library

Load the data 

cleaned_data <- get_data(national=FALSE, 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_data)

Run JAGS model and get posterior point estimates with uncertainty 

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

Plot posterior point estimates with uncertainty 

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