| Re: Discrepancy with estimates from Tables 11.1 and 11.2 in Niger 2012 DHS [message #22089 is a reply to message #22017] |
Fri, 29 January 2021 09:52   |
Shireen-DHS
Messages: 140
Registered: August 2020
Location: USA
|
Senior Member |
|
|
Dear Deepali,
The code below will match the tables. This code is available on on code share library hosted on the GitHub website (https://github.com/DHSProgram). Please read the readme file first. The indicators you are interested in are found here ( https://github.com/DHSProgram/DHS-Indicators-Stata/tree/mast er/Chap11_NT). Before running any code please read the notes in the main file or the specific do file you are using.
For Table 11.1:
Actually the denominator for the early initiation of breastfeeding is 5,143 (see the 4th column in the table).
The code below matches the table.
*open KR file
**** child's age ****
gen age = v008 - b3
*code copied from the NT_BF_INIT.do file
//Start breastfeeding within 1 hr
gen nt_bf_start_1hr= (m4!=94 & m4!=99) & (inrange(m34,0,100)) if (midx==1 & age<24)
label values nt_bf_start_1hr yesno
label var nt_bf_start_1hr "Started breastfeeding within one hour of birth - last-born in the past 2 years"
gen wt=v005/1000000
tab nt_bf_start_1hr [iw=wt]
*I obtain 52.9% and 5143 cases as in the table.
For Table 11.2:
For this indicator you need to select for the youngest child under 24 months living with the mother. As shown in the NTmain.do file on line 97 this is done as follows:
* Selecting for youngest child under 24 months and living with mother
keep if age < 24 & b9 == 0
* if caseid is the same as the prior case, then not the last born
keep if _n == 1 | caseid != caseid[_n-1]
*note this will drop cases so it should only be done if you are analyzing IYCF indicators.
*the following code was copied from the NT_IYCF.do file from line 49 to 86 to get the breastfeeding status indicator.
//currently breastfed
gen nt_bf_curr= m4==95
label values nt_bf_curr yesno
label var nt_bf_curr "Currently breastfeeding - last-born under 2 years"
//breastfeeding status
gen water=0
gen liquids=0
gen milk=0
gen solids=0
*Child is given water
replace water=1 if (v409>=1 & v409<=7)
*Child given liquids
foreach xvar of varlist v409a v410 v410a v412c v413*{
replace liquids=1 if `xvar'>=1 & `xvar'<=7
}
*Given powder/tinned milk, formula, or fresh milk
foreach xvar of varlist v411 v411a {
replace milk=1 if `xvar'>=1 & `xvar'<=7
}
*Given any solid food
foreach xvar of varlist v414* {
replace solids=1 if `xvar'>=1 & `xvar'<=7
}
replace solids=1 if v412a==1 | v412b==1 | m39a==1
gen nt_bf_status=1
replace nt_bf_status=2 if water==1
replace nt_bf_status=3 if liquids==1
replace nt_bf_status=4 if milk==1
replace nt_bf_status=5 if solids==1
replace nt_bf_status=0 if nt_bf_curr==0
label define bf_status 0"not bf" 1"exclusively bf" 2"bf & plain water" 3"bf & non-milk liquids" 4"bf & other milk" 5"bf & complemenatry foods"
label values nt_bf_status bf_status
label var nt_bf_status "Breastfeeding status for last-born child under 2 years"
//Age categories
cap recode age (0/1=1 " 0-1") (2/3=2 " 2-3") (4/5=3 " 4-5") (6/8=4 " 6-8") (9/11=5 " 9-11") (12/17=6 " 12-17") (18/23=7 " 18-23") , gen(agecats)
tab agecats nt_bf_status [iw=wt], row
* This matches the final report (Table 11.2) by age categories. For instance I obtain for age 0-1 I obtain 2.7% for not breastfed and 415 cases.
Thank you.
Shireen Assaf
The DHS Program
|
|
|
|
The DHS Program User Forum
Members
Search
Help
Register
Login
Home
