inst/doc/subnational_singlecountry_customdata_models.R

In mcmsupply: Estimating Public and Private Sector Contraceptive Market Supply Shares

## ----setup, include=FALSE-----------------------------------------------------
knitr::opts_chunk$set(echo = TRUE)

## ---- include=TRUE, message=FALSE, eval=FALSE---------------------------------
#  library(mcmsupply)
#  library(dplyr)
#  set.seed(1209)

## ---- include=TRUE, message=FALSE, warning=FALSE, error=TRUE------------------

cleaned_data <- get_data(national=FALSE, local=TRUE, 
                         surveydata_filepath = "~/Documents/R/mcmsupply/inst/data-raw/my_custom_subnational_data_bad.xlsx",
                         mycountry="Ethiopia")

## ---- include=TRUE, message=FALSE, warning=FALSE, eval=FALSE------------------
#  
#  cleaned_data <- get_data(national=FALSE, local=TRUE,
#                           surveydata_filepath = "~/Documents/R/mcmsupply/inst/data-raw/my_custom_subnational_data_good.xlsx",
#                           mycountry="Ethiopia")

## ---- include=TRUE, message=FALSE, warning=FALSE, eval=FALSE------------------
#  pkg_data <- get_modelinputs(startyear=1990, endyear=2025.5,
#                              nsegments=12, raw_data = cleaned_data)
#  

## ---- include=TRUE, message=FALSE, eval=FALSE---------------------------------
#  mod <- run_jags_model(jagsdata = pkg_data, jagsparams = NULL,
#                        n_iter = 40000, n_burnin = 10000, n_thin = 15)

## ---- include=TRUE, message=FALSE, eval=FALSE---------------------------------
#  plots <- plot_estimates(jagsdata = pkg_data, model_output = mod)

## ---- include=TRUE, message=FALSE, eval=FALSE---------------------------------
#  estimates_2018 <- pull_estimates(model_output = mod, country = 'Ethiopia', year=2018)
#  
#  head(estimates_2018)