** DEMOGRAPHIC AND HEALTH SURVEY DATA

/*
		A stata code (do file) to compute skilled birth attendant from DHS data
*/

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

*** YEAR = 2013
** SKILLED BIRTH ATTENDANTS (SBA)

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

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

** SURVEY SET
gen psu =    v021
gen strata = v022

svyset psu [pw = weight], strata(strata) vce(linearized)

keep if v208 !=0 & v208 !=.

*==============================================================================*
** DROP IF NOT WITHIN SAMPLE
qui regr m3a m3b m3c if v208 >0 & v208 !=. [pw=weight]		
drop if e(sample)!=1	

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

// RENAME


rename v190 wealth
rename v025 residence
rename v024 region

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


// GENERATING SKILLED BIRTH ATTENDANT - VARIABLE

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


** SKILLED BIRTH ATTENDANTS

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


// CHECK - [Yes = Percentage delivered by a skilled provider]

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

svy: tab region skilled_birth, count format(%4.0f) miss

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