Stunting rate with accurate observation number [message #15090] |
Thu, 31 May 2018 11:46 ![Go to next message Go to previous message](/theme/default/images/down.png) |
anikhpg42@gmail.com
Messages: 38 Registered: December 2017 Location: Bangladesh
|
Member |
![live%3A9d510b086916cb2d](/theme/default/images/skype.png)
|
|
Hi,
In order to find the stunting rate among children under age 5 using BDHS 14 data, I have used the following Stata command.
I used this command both in KR & PR dataset. But the total observation number is not exactly matching with the BDHS-14 report.
Prevalence of Stunting in 2014 (Children under age 5)
My estimation,
Using KR file, 36.24 % (N = 6965)
Using PR file, 36.09 % (N = 7256)
From BDHS'14 Report
36.1 % (N = 7318)
Here, my question is, how can I obtain the same prevalence rate and same observation for calculating child stunting rate?
Where is my mistakes, that I've made in my Stata command?
I have to get the exact observations (N = 7318) and the exact stunting rate (36.10%).
Stata command for stunting (just for the KR file)
*SVY command
gen strata=v023
gen psu=v021
gen sampwt=v005/1000000
svyset psu [pw=sampwt], strata (strata)
//child stunting calculation
codebook hw70
tab hw70 if hw70>9990,m
tab hw70 if hw70>9990, m nolabel
gen HAZ=hw70
replace HAZ=. if HAZ>=9996
*****
gen stunted=.
replace stunted=0 if HAZ~=.
replace stunted=1 if HAZ<-200
svy: tab stunted
Best regards,
Anik
ASIBUL ISLAM ANIK
|
|
|