Chapter 5 Malawi and Zambia

zambia_sample_200 <- oversample_wrapper_pair(name = "Zambia",
                            bins = zambia_sampling_bin,
                            raster_1 = zambia_population,
                            raster_2 = zambia_population_1k,
                            raster_3 = zambia_population_5k,
                            random_number = 680,
                            original_min_unit3 = 1,
                            original_sample_unit3 = 100,
                            by_factor=1.1,
                            verbose=2)

malawi_sample_200 <- oversample_wrapper_pair(name = "Malawi",
                            bins = maw_sampling_bin,
                            raster_1 = malawi_population,
                            raster_2 = malawi_population_1k,
                            raster_3 = malawi_population_5k,
                            random_number = 680,
                            original_min_unit3 = 1,
                            original_sample_unit3 = 100,
                            by_factor=1.1,
                            verbose=2)

respondents_per_bin <- function(sp_obj,name,onlydf=FALSE) {
  if (any("Unit_3"==sp_obj$type)) {

    tab1 <- sp_obj %>% 
    st_as_sf() %>% 
    filter(pick=="Sample"&type=="Unit_3") %>% 
    mutate(Bin = sapply(Name,function(x) {
      strsplit(x,"@") %>% 
        lapply(function(y) y[2]) %>% 
        unlist()
    })) %>% 
    as.data.frame() %>% 
    group_by(Bin) %>% 
    summarise(`Sampling Units`=n()) %>% 
    mutate(`Number of Households` = `Sampling Units` * 50) %>% 
    mutate(Bin = as.numeric(Bin)) %>% 
    arrange(Bin)
} else {
  tab1 <- sp_obj %>% 
    st_as_sf() %>% 
    filter(pick=="Sample"&type!="Bin") %>% 
    mutate(Bin = sapply(Name,function(x) {
      strsplit(x,"@") %>% 
        lapply(function(y) y[2]) %>% 
        unlist()
    })) %>% 
    as.data.frame() %>% 
    group_by(Bin) %>% 
    summarise(`Sampling Units`=n()) %>% 
    mutate(`Number of Households` = `Sampling Units` * 25) %>% 
    mutate(Bin = as.numeric(Bin)) %>% 
    arrange(Bin)
}
  if(onlydf) return(tab1)
  
  kable(tab1,caption = paste0("Number of Respondents and sampling units in ",
                              name,". (n = ",
                              tab1 %$% 
                                sum(`Number of Households`),")")) %>% 
  kable_styling(full_width = FALSE,position = "left")
  
}

nairobi_sample_150 %>% 
  respondents_per_bin("Nairobi") 

lusaka_sample_150 %>% 
  respondents_per_bin("Lusaka") 

lilongwe_sample_150 %>% 
  respondents_per_bin("Lilongwe") 

zambia_sample_200 %>% 
  respondents_per_bin("Zambia")

malawi_sample_200 %>% 
  respondents_per_bin("Malawi")