|Sampling Weigths for children under 5 [message #19318]
||Fri, 29 May 2020 11:16
Registered: May 2020
I am new to survey data and I am having trouble trying to find my error in my R-Code. I chose Philippines 2017 and I am trying to calculate the percentage of vitamin A supplementation for this country for children under 5.
This is my R-Code
phi17 <- read.dta("PHKR71FL.DTA")
# Using labels as column names
var.labels <- attr(phi17, "var.labels")
names(phi17) <- var.labels
For speficying the model I use the following code:
# Simplifying column names
phi17$weight <- as.numeric(phi17$`women's individual sample weight (6 decimals)`/1000000)
phi17$psu <- phi17$`primary sampling unit`
phi17$strata <- phi17$`stratification used in sample design`
phi17$vi_a = phi17$`vitamin a in last 6 months`
phi_17_design <- svydesign(~ psu,
strata = ~strata,
weights = ~ weight,
data = phi_17)
So, I used the weights for the women but when I calculate the percentage of children who received a vitamin A supplement in the last 6 months in Philippines 2017, i get the following results:
svymean(~ vi_a, phi_17_design, na.rm = TRUE) --> yes: 69.47%
HOWEVER, when I read the final report, the result should actually be 76%.
Does anyone recognize the mistake? Did I specified the weigths correctly??
Thank you in advance!