In FPRgroup/FPEMcountry: fpemlocal

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)
devtools::load_all()
library(ggplot2)
library(grid)
library(gridExtra)

Set output directory and divs for the set of results

# outputdir <- "output_upto2013"
# filenamecreated <- "results_upto2013.csv"
divs <- read.csv("C:/Users/desktop-g/Downloads/Track20Database033020.csv") 
divs <- divs %>% dplyr::select(division_numeric_code) %>% unique()
popcountdiv <- population_counts
divs <- divs %>% dplyr::filter(division_numeric_code %in% popcountdiv$division_numeric_code)
divnames <- divs %>% 
  dplyr::left_join(divisions) %>%
  dplyr::select(name_country) %>%
  dplyr::pull()
divs <- divs %>% dplyr::pull()
# outputdir <- "output"
# filenamecreated <- results_unpd2020.csv"
# # vector of countries where we have data
# global_divs <- fpemdata:::index_m$index_area_df$division_numeric_code
# divs <- dplyr::tibble(division_numeric_code = as.numeric(as.character(global_divs[global_divs %in% fpemdata::divisions$division_numeric_code])))
# divs <- divs %>% dplyr::filter(division_numeric_code %in% fpemdata::population_counts$division_numeric_code)
#  # dplyr::pull()
# divnames <- divs %>% 
#   dplyr::left_join(fpemdata::divisions) %>%
#   dplyr::select(name_country) %>%
#   dplyr::pull()
# divs <- divs %>% dplyr::pull()

for (div in divs) {
  runname <- paste0("avenir_", div)
  fit <- fit_fp_c_autosave(
    runname = runname,
    surveydata_filepath = "C:/Users/desktop-g/Downloads/Track20Database033020.csv",
    is_in_union = "ALL",
    division_numeric_code = div,
    first_year = 1970,
    last_year = 2030,
    diagnostic = FALSE
)
  calc_fp_c_autosave(runname = runname)
}

unions <- c("Y", "N", "ALL")

allreslist <- list()
for (j in 1:length(divs)) {
  divname <- divnames[j]
  div <- divs[j]
  runname <- paste0("avenir_", div)
  reslist <- readRDS(file.path("output/results", paste0(runname, ".rds")))

  newreslist <- list()
  for (i in 1:3) {
    is_in_union <- unions[i]
    fit <- fitlist[[is_in_union]]
    res <- reslist[[is_in_union]]
    newreslist[[i]] <- data.frame(
              name_country = divname,
              division_numeric_code = div,
              is_in_union = is_in_union,
              year = res$contraceptive_use_any$year ,
              percentile = res$contraceptive_use_any$percentile,
              contraceptive_use_any = res$contraceptive_use_any$value,
              contraceptive_use_modern = res$contraceptive_use_modern$value,
              contraceptive_use_traditional = res$contraceptive_use_traditional$value,
              unmet_need_any = res$unmet_need_any$value,
              unmet_need_modern = res$unmet_need_modern$value,
              demand_satisfied_modern = res$demand_satisfied_modern$value
        )
  }
  # rbind across the list
  allreslist[[div]] <- do.call(rbind, newreslist)
}
df <- do.call(rbind, allreslist)
write.csv(df, filenamecreated)