The DHS Program User Forum
Discussions regarding The DHS Program data and results
Home » Topics » Nutrition and Anthropometry » Can't get exact match between one regional calculation and reported stunting value (Stunting of children under 5)
Can't get exact match between one regional calculation and reported stunting value Thu, 17 October 2019 12:38
 dhsLearner Messages: 4Registered: October 2019 Member
I am trying to calculate childhood (under 5) stunting rates by region in Senegal in 2012. For one region, I am getting a different result than what it is reported in the Senegal sumary report. Here is my R code:

# Create a binary variable for children stunted
PR2012\$stuntedUnderFive<-ifelse(PR2012\$hc70< -200, 1, 0)
# Assign NA when z score is >9900 sd (NA) #9998 are don't know responses
PR2012\$stuntedUnderFive<-ifelse(PR2012\$hc70>9900|
PR2012\$hc1> 59, NA, PR2012\$stuntedUnderFive)

# Create weight variable
PR2012\$wt<-PR2012\$hv005/1000000

# Tabluate stunted by region
PR2012%>%
group_by(shzone)%>% #shzone is the four regions (north, east, south west)
summarize(meanStunting = weighted.mean(stuntedUnderFive,wt,na.rm = T))

Results:

1 centre 0.190
2 nord 0.168
3 ouest 0.145
4 sud 0.266

I calculate West (ouest) as 14.5%, but it is reported as 14.7% in the summary (https://dhsprogram.com/pubs/pdf/FR288/FR288eng.pdf).

Can anyone help me figure out what I am doing incorrectly?

 Can't get exact match between one regional calculation and reported stunting value By: dhsLearner on Thu, 17 October 2019 12:38 Re: Can't get exact match between one regional calculation and reported stunting value By: Trevor-DHS on Thu, 17 October 2019 14:16 Re: Can't get exact match between one regional calculation and reported stunting value By: dhsLearner on Thu, 17 October 2019 14:25
 Previous Topic: Anthropometric data for children Turkey 2008 and 2013 Next Topic: 'Don't know' in exclusive breastfeeding
Goto Forum:

Current Time: Fri Apr 12 10:07:34 Coordinated Universal Time 2024