The DHS Program User Forum
Discussions regarding The DHS Program data and results
Home » Countries » Bangladesh » Bangladesh infant mortality replication
Bangladesh infant mortality replication [message #5570] Tue, 09 June 2015 13:30 Go to previous message
busbyj2 is currently offline  busbyj2
Messages: 3
Registered: June 2015
Location: Austin Texas
Member
We're trying to replicate the infant mortality data from 2011 DHS survey for Bangladesh for the Divisions but we're getting different numbers than what's in the report.

Here is a piece of our STATA code. Would welcome any troubleshooting.

gen prev5yr = v008 - b3 < 60
gen prev3yr = v008 - b3 < 36

gen infdeath = b6 <= 211
replace infdeath = . if b6 == 997 | b6 == 998 | b6 == 999

gen chldeath = b6 <= 305
replace chldeath = . if b6 == 997 | b6 == 998 | b6 == 999

gen attbirth = m3a == 1
replace attbirth = 1 if m3b == 1
replace attbirth = 1 if m3c == 1
replace attbirth = 1 if m3d == 1
replace attbirth = 1 if m3e == 1
replace attbirth = 1 if m3f == 1

gen delhfac = m15 != 11
replace delhfac = 1 if m3g == 1

************** 1 - INFANT MORTALITY (PREV 5 YEARS) **********************

gen indic=1

svy: mean infdeath, subpop(prev5yr)
* show size of subpopulations
estat size
* store Part_Results in various variables
matrix m=e(b)
matrix s=e(V)
matrix su=e(_N)
estat effects, deft
matrix d=e(deft)

gen natl=m[1,1]
gen senatl=s[1,1]
gen sizenatl=su[1,1]
gen deftnatl=d[1,1]

svy: mean infdeath, over(region) subpop(prev5yr)
estat size
matrix mq=e(b)
matrix sq=e(V)
matrix suq=e(_N)
estat effects, deft
matrix dq=e(deft)

forvalues i = 1/`numreg'{
gen reg`i'=mq[1,`i']
}
forvalues i = 1/`numreg'{
gen sereg`i'=sq[`i',`i']
}
forvalues i = 1/`numreg'{
gen sizereg`i'=suq[1,`i']
}
forvalues i = 1/`numreg'{
gen deftreg`i'=dq[1,`i']
}
savesome indic-deftreg`numreg' if _n==1 using birth_result1, replace
matrix drop m s su mq sq suq d dq
drop natl senatl sizenatl deftnatl indic

forvalues i = 1/`numreg'{
drop reg`i'
}
forvalues i = 1/`numreg'{
drop sereg`i'
}
forvalues i = 1/`numreg'{
drop sizereg`i'
}
forvalues i = 1/`numreg'{
drop deftreg`i'
}
 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: Deal with missing data in BMI
Next Topic: Number of households
Goto Forum:
  


Current Time: Fri Mar 29 07:26:02 Coordinated Universal Time 2024