*load in dataset
	use "$jordhsraw\1990\JOKR21DT\JOKR21FL", clear

*create an age var
	gen age_ch=v008-b3

*drop all vars with bf in it
	capture drop *bf*

*exclusive breastfeeding
	*generate an othfood var to indicate if consumed other foods.
	capture drop othfood
	egen othfood=rowtotal(v409 v409a v410 v410a v411 v411a v412 v413a v413 v414), missing

*create ebf based on age range, last-born, living with mother, current bf status, and other food consumed
	*age_ch<6 smaller than 6 month of age 
	*b9==0 lives with mother 
	*m4==95 currently breastfeeding for the specific child 
	*othfood==0 did not consume any other foods 
	gen ebf=.
	replace ebf=0 if age_ch<6 & b9==0 & midx==1
	replace ebf=1 if age_ch<6 & b9==0 & m4==95 & midx==1 & othfood==0
	bigtab ebf age_ch othfood 
	
	*validate the estimate against DHS report 
	*23.95% Sample size 724.7 Not validated. 
	tab ebf [iweight=v005/1000000]