*Early initiation of breast feeding
gen early_init=0 if b19<24
replace early_init=1 if m34<101 & b19<24
la val early_init A
label var early_init "Early initiation of BF"
tab v013 early_init

*Exclusive breast feeding
gen water=0
gen liquids=0
gen milk=0
gen solids=0
gen breast=0
gen bottle=0

*TO DETERMINE IF CHILD IS GIVEN WATER, SUGAR WATER, JUICE, TEA OR OTHER.
replace water=1 if (v409>=1 & v409<=7)
                                                                   
* IF GIVEN OTHER LIQUIDS
foreach xvar of varlist v409a v410 v410*  v413* {
replace liquids=1 if `xvar'>=1 & `xvar'<=7
}

cap replace liquids=1 if v412c>=1 & v412c<=7
                                                                                                
* IF GIVEN POWDER/TINNED milk, FORMULA OR FRESH milk
foreach xvar of varlist v411 v411a v412 v414p {
replace milk=1 if `xvar'>=1 & `xvar'<=7
}

* IF STILL BREASTFEEDING
replace breast=1 if m4==95 

* IF WAS EVER BOTTLE FED
replace bottle=1 if m38==1 

*IF GIVEN ANY SOLID FOOD
foreach xvar of varlist v414* {
replace solids=1 if `xvar'>=1 & `xvar'<=7
}
replace solids=1 if v412a==1 | v412b==1 
                
gen diet=7
replace diet=0 if water==0 & liquids==0 & milk==0 & solids==0
replace diet=1 if water==1 & liquids==0 & milk==0 & solids==0
replace diet=2 if            liquids==1 & milk==0 & solids==0 
replace diet=3 if                         milk==1 & solids==0 
replace diet=4 if                         milk==0 & solids==1 
replace diet=5 if                         milk==1 & solids==1 
replace diet=6 if breast==0 

*diet=1: given only water (full bf)
*diet=2: given only liquids (bf & liquids)
*diet=3: given only milks (bf & milk)
*diet=4: given only solids (bf & solids)
*diet=5: given only milk and solids (bf & milk & solids)
*diet=6: not still breastfeeding (weaned)
*diet=7: not now being breastfed (m4~=95)

gen ebf=0 if b19<6
replace ebf=1 if diet==0 & b19<6 
la val ebf A
lab var ebf "Exclusive BF under 6 months"

gen water_bf=0 if b19<6
replace water_bf=1 if diet==1 & b19<6 
la val water_bf A
lab var water_bf "Water associated BF under 6 months"

*Survey setting
gen wt=v005/1000000
svyset psu [pweight=wt], strata(v023) vce(linearized) singleunit(centered) 

*Tabulate 
svy:tab survey ebf, row percent
svy:tab survey ebf, count 

svy:tab survey water_bf, row percent
svy:tab survey water_bf, count

svy:tab survey early_init, row percent
svy:tab survey early_init, count 

		My output				From DHS Reports		
	EBF	Water associated BF	Count		EBF	Water associated BF	Count
2010-11	38.60	39.03	1200		39.0	39.2	1204
2012-13	37.57	45.36	639		37.5	46.4	640
2014	32.18	40.89	587		33.0	41.2	588
2015	33.03	37.68	603		33.3	38.4	604
2016	36.21	38.56	559		36.4	38.3	560
2017	41.66	37.26	1050		42.1	37.6	1052

	Early initiation of BF				
	My output			From DHS Reports	
	%	Count		%	Count
2015	29.69	2384		30.5	2383
2016	28.61	2296		29.4	2281
2017	33.05	4373		33.6	4368