indicators <- c(
  "demo", "anthro", "food", "hunger", "adl", "disability",
  "mental", "dementia", "health", "oedema", "screening", "income",
  "wash", "visual", "misc"
)
classicIndicators <- indicators[indicators != "anthro"]
## Bootstrap classic
classicEstimates <- estimate_classic(
  x = indicators.ALL, w = testPSU, 
  indicators = classicIndicators, replicates = 9
)
probitEstimates <- estimate_probit(
  x = indicators.ALL, w = testPSU, replicates = 9
)
merge_op(x = classicEstimates, y = probitEstimates)
Run the code above in your browser using DataLab