est_df <- boot_bw(
x = indicatorsHH, w = villageData, statistic = bootClassic,
params = "anc1", strata = "region", replicates = 9, parallel = TRUE
) |>
boot_bw_estimate()
## Add population ----
pop_df <- somalia_population |>
subset(select = c(region, total))
names(pop_df) <- c("strata", "pop")
estimate_total(est_df, pop_df, strata = "region")
Run the code above in your browser using DataLab