| Bangladesh infant mortality replication [message #5570] |
Tue, 09 June 2015 13:30  |
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'
}
|
|
|
|
|
|
|
|
|
|
|
|
|
|
The DHS Program User Forum
Members
Search
Help
Register
Login
Home
