The DHS Program User Forum: Reproductive Health » Early neonatal mortality variable

Home » Topics » Reproductive Health » Early neonatal mortality variable

Show: Today's Messages :: Show Polls :: Message Navigator

Early neonatal mortality variable [message #13273]

Wed, 11 October 2017 22:00

jahm3675
Messages: 2
Registered: August 2017
Location: Sydney

Member

Hi DHS experts

I am calculating a variable for early neonatal mortality from pkbr61dt data set from DHS suevey 2012-13 of Pakistan.
When I run the following command It tabulates the rates which are matching with the DHS report but the variable that it creates has different rates.
for example:
Stata Command:

* create a child Alive or Died variable using the b5 variable
gen alive=b5
lab def alive 0 "Died" 1 "Alive"
lab val alive alive
lab var alive "Alive or Died by the time of survey"

*** Age at death using variables b6 and b5
gen age_death=.
replace age_death = 0 if b6<=106 & b5==0
replace age_death = 1 if b6>=107 & b6<=130 & b5==0
replace age_death = 2 if (b6> 130 & b6<=999) | b5==1
lab def age_death 0 "Early Neonatal Death" 1 "Late Neonatal Death" 2 "Survived Neonatal Period"
lab val age_death age_death
lab var age_death "Neonatal Mortality Status"

* Set length of calendar to use
gen callen = v018 + 59
* If calendar is aligned right (as in original dataset), use the following:
gen beg = v018
gen end = callen
* If calendar is aligned left (as it is in some datasets), use the following:
*gen beg = 1
*gen end = 60

gen wt=v005/1000000
svyset v021 [pw=wt], strata(v023)

svy: tab age_death if v008-b3<60, per count form(%7.3g)
svy: tab age_death if v008-b3<60, per col form(%7.3g)

----------------------
Neonatal |
Mortality |
Status | count
----------+-----------
Early Ne | 522
Late Neo | 139
Survived | 11316
|
Total | 11977
----------------------

----------------------
Neonatal |
Mortality |
Status | column
----------+-----------
Early Ne | 4.36
Late Neo | 1.16
Survived | 94.5
|
Total | 100
----------------------
Key: column = column percentages

But when I tabulate the new variable it gives following rates:

Neonatal Mortality |
Status | Freq. Percent Cum.
-------------------------+---------------------------------- -
Early Neonatal Death | 2,026 4.03 4.03
Late Neonatal Death | 558 1.11 5.14
Survived Neonatal Period | 47,654 94.86 100.00
-------------------------+---------------------------------- -
Total | 50,238 100.00

Which is for all years. it should give me a variable that already calculates the rates as in the weighted table which is i dont know why a temporary table.

Then I run another command to restrict the rates for six year period:
. keep if b2 >= 2008 & b2 <= 2013

it gives following rates:

Neonatal Mortality |
Status | Freq. Percent Cum.
-------------------------+---------------------------------- -
Early Neonatal Death | 445 3.80 3.80
Late Neonatal Death | 106 0.90 4.70
Survived Neonatal Period | 11,174 95.30 100.00
-------------------------+---------------------------------- -
Total | 11,725 100.00

I am stuck at this stage to get over it but cant figure out a way. the DHS rate is 522 for the 5 years preceding survey.

When I check the data for year of birth and deaths variable:
Neonatal Mortality Year of birth
Status 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013
Early Neonatal Death 92 74 105 105 112 105 133 76 91 99 86 86 80 3
Late Neonatal Death 22 25 12 26 22 12 20 22 23 24 18 21 20 0
Survived Neonatal Per 2,232 1,679 2,360 1,947 2,444 2,431 2,402 2,257 2,284 2,361 2,316 2,073 2,095 45
Total 2,346 1,778 2,477 2,078 2,578 2,548 2,555 2,355 2,398 2,484 2,420 2,180 2,195 48

The rate is still 445. Where are 77 cases going? Why the variable created does not have 522 early neonatal deaths?

Will appreciate your help.Thanks

Report message to a moderator