** USE THE BIRTHS RECODE (ETBR70FL) **
** ETHIOPIA DHS 2016

clear all
set matsize 800
set maxvar 9000
set mem 1g
cd "..."
use "ETBR70FL", clear
set more off

********************************************************************************
** WEIGHT VARIABLE
gen weight = v005/1000000

**************************

** SURVEY SET
gen psu =    v021
gen strata = v023
svyset [pweight = weight], psu(psu) strata(strata) vce(linearized)
*svydes

********************************************************************************


// DELIVERY CARE = facility-based deliveries
recode m15 (21/46=1 "Facility-based deliveries") ///
(11/12 96=0 "Non-facility deliverie"), gen(facility_delivery)
label var facility_delivery "Facility-based deliveries"
label val facility_delivery facility_delivery

*******************

cap drop DHS_delivery
recode m15 (21/26=1 "Public sector") ///
(31/36=2 "Private sector") (41/46=3 "NGO") ///
(11/12=4 "Home") (96=5 "Other"), gen(DHS_delivery)
label var DHS_delivery "Place of delivery as in DHS eport"
label val DHS_delivery DHS_delivery

recode m15 (21/26=1 "Public sector") (31/36=2 "Private sector") ///
(41/46=3 "NGO") (11/12 96=4 "Elsewhere"), gen(DHS_deliveries)
label var DHS_deliveries "Replication of DHS table for facility-based deliveries"
label val DHS_deliveries DHS_deliveries

*==============================================================================*

** SKILLED BIRTH ATTENDANT

gen skilled_birth	=	0
label define skilled_birth 0"Unskilled" 1"Skilled"
label var skilled_birth "Birth delivered by skilled birth attendant"
label val skilled_birth skilled_birth

** SKILLED BIRTH ATTENDANTS RECODE

foreach xvar of varlist m3a m3b m3c m3d m3e {
replace skilled_birth=1 if `xvar'==1
}

**

*==============================================================================*

** DROP IF NOT WITHIN SAMPLE
qui regr DHS_delivery [pw=weight]		
drop if e(sample)!=1		

********************************************************************************
** CHECK: TABLE 9.8
tab facility_delivery [iw=weight], m

svy: tab facility_delivery, count format(%4.0f)
svy: tab facility_delivery, percent format(%4.1f) col

********************************************************************************

** CHECK: TABLE 9.9
tab skilled_birth [iw=weight], m

svy: tab skilled_birth, count format(%4.0f)
svy: tab skilled_birth, percent format(%4.1f) col

********************************************************************************
********************************************************************************
********************************************************************************

svy: tab v190 facility_delivery, percent format(%4.1f) miss row
svy: tab v190 skilled_birth, percent format(%4.1f) miss row

svy: tab DHS_deliveries skilled_birth, percent format(%4.1f) miss row

********************************************************************************

exit

** USE THE BIRTHS RECODE (ETBR61FL) **
** ETHIOPIA DHS 2011

clear all
set matsize 800
set maxvar 9000
set mem 1g
cd "..."
use "ETBR61FL", clear
set more off

********************************************************************************
** WEIGHT VARIABLE
gen weight = v005/1000000

**************************

** SURVEY SET
gen psu =    v021
gen strata = v023
svyset [pweight = weight], psu(psu) strata(strata) vce(linearized)
*svydes

********************************************************************************
recode m15 (21/33=1 "Facility-based deliveries") ///
(11/12 96 99=0 "Non-facility deliverie"), gen(facility_delivery)
label var facility_delivery "Facility-based deliveries"
label val facility_delivery facility_delivery

*******************

cap drop DHS_delivery
recode m15 (21/24=1 "Public sector") ///
(31 32=2 "Private sector") (33=3 "NGO") ///
(11/12=4 "Home") (96 99=5 "Other"), gen(DHS_delivery)
label var DHS_delivery "Place of delivery as in DHS eport"
label val DHS_delivery DHS_delivery

recode m15 (21/24=1 "Public sector") (31/32=2 "Private sector") ///
(33=3 "NGO") (11/12 96 99=4 "Elsewhere"), gen(DHS_deliveries)
label var DHS_deliveries "Replication of DHS table for facility-based deliveries"
label val DHS_deliveries DHS_deliveries

*==============================================================================*

** SKILLED BIRTH ATTENDANT

gen skilled_birth	=	0
label define skilled_birth 0"Unskilled" 1"Skilled"
label var skilled_birth "Birth delivered by skilled birth attendant"
label val skilled_birth skilled_birth

** SKILLED BIRTH ATTENDANTS RECODE

foreach xvar of varlist m3a m3b {
replace skilled_birth=1 if `xvar'==1
}

**

*==============================================================================*

** DROP IF NOT WITHIN SAMPLE
qui regr facility_delivery [pw=weight]		
drop if e(sample)!=1		

********************************************************************************
** CHECK: TABLE 9.6
tab facility_delivery [iw=weight], m

svy: tab facility_delivery, count format(%4.0f)
svy: tab facility_delivery, percent format(%4.1f) col

********************************************************************************

** CHECK: TABLE 9.7
tab skilled_birth [iw=weight], m

svy: tab skilled_birth, count format(%4.0f)
svy: tab skilled_birth, percent format(%4.1f) col

********************************************************************************
********************************************************************************
********************************************************************************

svy: tab v190 facility_delivery, percent format(%4.1f) miss row
svy: tab v190 skilled_birth, percent format(%4.1f) miss row

svy: tab DHS_deliveries skilled_birth, percent format(%4.1f) miss row

********************************************************************************

exit