| Sampling Weigths for children under 5 [message #19318] | 
			Fri, 29 May 2020 11:16   | 
		 
		
			
				
				
				
					
						  
						rafaaqp1990
						 Messages: 2 Registered: May 2020 
						
					 | 
					Member  | 
					 | 
		 
		 
	 | 
 
	
		Dear all, 
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! 
 
Best regards 
 
rafaaqp1990 
 
 
 
 
		
		
		
 |  
	| 
		
	 | 
 
 
 | 
	| 
		
 | 
	| 
		
 |