//** ETHIOPIA DHS 2016 **// clear all set matsize 800 set mem 1g set maxvar 9000 cd "C:\Users\User88\Desktop\ETHIOPIA DHS 2016\ETKR70DT" use "ETKR70FL", clear set more off ******************************************************************************** **WEIGHT VARIABLE gen weight = v005/1000000 ******************************************************************************** ******************************************************************************** ** SURVEY SET gen psu = v021 gen strata = v023 svyset psu [pw = weight], strata(strata) vce(linearized) ******************************************************************************** ******************************************************************************** rename v013 age rename v106 education rename v190 wealth rename v025 residence rename v024 region ******************************************************************************** **PLACE OF DELIVERY cap drop DHS_delivery recode m15 (21/27=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 cap drop DHS_delivery2 recode m15 (21/27=1 "Public sector")(31/36=2 "Private sector") /// (11/12=4 "Home") (41/46 96=5 "NGO/Other"), gen(DHS_delivery2) label var DHS_delivery2 "Place of delivery" label val DHS_delivery2 DHS_delivery2 ******************************************************************************** **DELIVERY IN HEALTH FACILITY cap drop facility_delivery 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 ******************************************************************************** ** ASSISTED BY MEDICALLY TRAINED PROVIDER DURING DELIVERY 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 foreach xvar of varlist m3a m3b m3c m3d m3e { replace skilled_birth=1 if `xvar'==1 } * cap drop skilled_provider egen skilled_provider = rowmax(m3a m3b m3c m3d m3e) label define skilled_provider 0"Unskilled" 1"Skilled" label var skilled_provider "Percentage delivered by a skilled provider" label val skilled_provider skilled_provider *tab skilled_provider [iw=weight], m ******************************************************************************** ** DROP MISSING CASES keep if DHS_delivery!=. ******************************************************************************** ** CHECK svy: tab wealth DHS_delivery, percent format(%4.1f) miss row svy: tab wealth facility_delivery, percent format(%4.1f) miss row svy: tab wealth skilled_birth, percent format(%4.1f) miss row *svy: tab wealth skilled_provider, percent format(%4.1f) miss row ******************************************************************************** ** RECODE VARIABLES FOR ANALYSIS ** * Mother’s age at birth cap drop agebirth gen agebirth=(b3-v011)/12 cap drop age_at_birth recode agebirth (min/19.91667=1 "<20") (20/34.91667=2 "20-34") /// (35/max=3 "35-49"), gen(age_at_birth) label var age_at_birth "Mother's age at birth" label val age_at_birth age_at_birth * Birth order gen birth_order1 = bord replace birth_order1 = bord-1 if b0 == 2 replace birth_order1 = bord-2 if b0 == 3 recode birth_order1 (1=1 "1") (2/3=2 "2-3") (4/5=3 "4-5") /// (6/20=4 "6+"), gen(birth_order) label var birth_order "Birth order" label values birth_order birth_order * Antenatal care visits cap drop antenatal_dhs recode m14 (0=0 "None") (1/3=1 "1-3") (4/20=4 "4+") /// (98=88 "DNK") (else =99 "Missing"), gen(antenatal_dhs) label var antenatal_dhs "Antenatal visits - as in DHS report" label val antenatal_dhs antenatal_dhs ** ========================================================================== ** ** ========================================================================== ** ** REPLICATING TABLE 9.8 (page 149) svy: tab age_at_birth DHS_delivery, percent format(%4.1f) row svy: tab birth_order DHS_delivery, percent format(%4.1f) row svy: tab antenatal_dhs DHS_delivery, percent format(%4.1f) row // report excludes DNK & Missing svy: tab residence DHS_delivery, percent format(%4.1f) row svy: tab region DHS_delivery, percent format(%4.1f) row svy: tab education DHS_delivery, percent format(%4.1f) row svy: tab wealth DHS_delivery, percent format(%4.1f) row **************** /* ** COUNTS svy: tab age_at_birth DHS_delivery, count format(%4.0f) miss svy: tab birth_order DHS_delivery, count format(%4.0f) miss svy: tab antenatal_dhs DHS_delivery, count format(%4.0f) miss svy: tab residence DHS_delivery, count format(%4.0f) miss svy: tab region DHS_delivery, count format(%4.0f) miss svy: tab education DHS_delivery, count format(%4.0f) miss svy: tab wealth DHS_delivery, count format(%4.0f) miss */ ** ========================================================================== ** ** REPLICATING TABLE 9.8 (page 149) - Facility Births svy: tab age_at_birth facility_delivery, percent format(%4.1f) row svy: tab birth_order facility_delivery, percent format(%4.1f) row svy: tab antenatal_dhs facility_delivery, percent format(%4.1f) row svy: tab residence facility_delivery, percent format(%4.1f) row svy: tab region facility_delivery, percent format(%4.1f) row svy: tab education facility_delivery, percent format(%4.1f) row svy: tab wealth facility_delivery, percent format(%4.1f) row **************** /* svy: tab age_at_birth facility_delivery, count format(%4.0f) miss svy: tab birth_order facility_delivery, count format(%4.0f) miss svy: tab antenatal_dhs facility_delivery, count format(%4.0f) miss svy: tab residence facility_delivery, count format(%4.0f) miss svy: tab region facility_delivery, count format(%4.0f) miss svy: tab education facility_delivery, count format(%4.0f) miss svy: tab wealth facility_delivery, count format(%4.0f) miss */ ** ========================================================================== ** ** REPLICATING TABLE 9.9 (page 150) - Percentage delivered by a skilled provider svy: tab age_at_birth skilled_provider, percent format(%4.1f) row svy: tab birth_order skilled_provider, percent format(%4.1f) row svy: tab antenatal_dhs skilled_provider, percent format(%4.1f) row svy: tab residence skilled_provider, percent format(%4.1f) row svy: tab region skilled_provider, percent format(%4.1f) row svy: tab education skilled_provider, percent format(%4.1f) row svy: tab wealth skilled_provider, percent format(%4.1f) row **************** /* svy: tab age_at_birth skilled_birth, count format(%4.0f) miss svy: tab birth_order skilled_birth, count format(%4.0f) miss svy: tab antenatal_dhs skilled_birth, count format(%4.0f) miss svy: tab residence skilled_birth, count format(%4.0f) miss svy: tab region skilled_birth, count format(%4.0f) miss svy: tab education skilled_birth, count format(%4.0f) miss svy: tab wealth skilled_birth, count format(%4.0f) miss */ ** ========================================================================== ** ** ================================= GRAPHS ================================= ** ** ========================================================================== ** ssc install hplot, replace //needs internet connection ssc install catplot, replace //needs internet connection catplot DHS_delivery2 [aw=weight], by(wealth, compact note("") col(1)) /// bar(1, blcolor(gs8) bfcolor(brown*.4)) blabel(bar, format(%9.1f) /// pos(base)) percent(wealth) subtitle(, pos(9) ring(1) bcolor(none) nobexpand place(e)) /// ytitle(Place of delivery (%)) var1opts(gap(*0.1) axis(noline)) /// var2opts(gap(*.2)) ysize(5) yla(none) ysc(noline) /// plotregion(lcolor(none)) catplot DHS_delivery2 wealth [aw=weight], percent(wealth) stack asyvars /// bar(1, bcolor(ltblue)) bar(2, bcolor(pink*.2)) bar(3, bcolor(brown*.4)) /// bar(4, bcolor(bluishgray)) bar(5, bcolor(gs10)) bar(6, bcolor(red*.4)) /// bar(7, bcolor(sand)) blabel(bar, format(%9.1f) pos(center)) /// blabel(bar, format(%9.1f) pos(center) size(tiny)orientation(vertical) color(black)) /// ytitle(%) title(Place of delivery by HH wealth) subtitle(Ethiopia DHS 2016) catplot facility_delivery wealth [aw=weight], percent(wealth) stack asyvars /// bar(1, bcolor(gold)) bar(2, bcolor(ltblue)) bar(3, bcolor(gs0)) /// blabel(bar, format(%9.1f) pos(center)) bar(4, bcolor(g)) ytitle(%) /// title(Percentage delivered in a health facility) subtitle(Ethiopia DHS 2016) catplot skilled_provider wealth [aw=weight], percent(wealth) stack asyvars /// bar(1, bcolor(gold)) bar(2, bcolor(ltblue)) bar(3, bcolor(gs0)) /// blabel(bar, format(%9.1f) pos(center)) bar(4, bcolor(g)) ytitle(%) /// title(Prevalence of skilled birth attendance) subtitle(Ethiopia DHS 2016) ** ========================================================================== ** exit