use BDBR51FL, clear

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

** WEIGHT VARIABLE
gen weight = v005/1000000

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

** SURVEY SET
gen psu =    v021
gen strata = v023
svyset psu [pw = weight], strata(strata)

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

// RENAME

rename v013 age
rename v106 education
rename v190 wealth
rename v025 residence
rename v024 region
//rename sdist district

///////////////////////////////////////////////////////////////////////////

// GENERATING DEPENDENT VARIABLES


** Doctor
recode m3a (1=1) (else=0), gen(assisted_doctor)
label define assisted_doctor 0"No" 1"Yes"
label var assisted_doctor "Birth delivered by a doctor"
label val assisted_doctor assisted_doctor
svy: tab wealth assisted_doctor, percent format(%4.1f) row

** Nurse (I AM NOT SURE WHICH VARIABLE IS MEANT TO BE A PARAMEDIC)
cap drop nurse
gen nurse = .
replace nurse = 0 if m3b !=1
replace nurse = 0 if m3d !=1
replace nurse = 0 if m3e !=1
replace nurse = 1 if m3b ==1
replace nurse = 1 if m3d ==1
replace nurse = 1 if m3e ==1
svy: tab residence nurse, percent format(%4.1f) row

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

** 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 m3d m3e {
replace skilled_birth=1 if `xvar'==1
}


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

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

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

** CHECK - DOES NOT MATCH WHAT'S IN THE REPORT
tab skilled_birth
svy: tab skilled_birth, count format(%4.0f)
svy: tab skilled_birth, percent format(%4.1f)

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

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

exit