* Calculate measurement date in days
gen mdate = mdy(hc18, hc17, hc19)
* Calculate birth date in days
gen bdate = mdy(hc30, hc16, hc31) if hc16 <= 31
replace bdate = mdy(hc30, 15, hc31) if hc16 > 31
* Calculate age in months with days expressed as decimals.
gen age = (mdate-bdate)/30.4375

//hw3 is height, clean DHS code and converts to cm
gen newh=hc3/10 if hc3!=9999
//hw3 is height in kg to one decimal w/o the decimal; convert to kg w/decimal; clean DHS code and
gen neww=hc2/10 if hc2!=999
* Compute Z-scores (In case of Nepal 2011 DHS, replace hc15=. if hc15==0)
zscore06, a(age) s(hc27) h(newh) w(neww) measure(hc15) male(1) female(2)

*Declaring Survey data
gen wt=hv005/1000000
svyset hv001 [pweight=wt], strata(hv022) vce(linearized) singleunit(missing) || hv002

*Completeness of information 
gen info=0
replace info=1 if hc33~=.

*Stunting
gen stunting=0 if hv103==1               
replace stunting=. if haz06>=6 
replace stunting=1 if haz06<-2 & hv103==1 
*Check with DHS report
svy:tab stunting if info==1

*Underweight
gen underwt=0 if hv103==1              
replace underwt=. if waz06>=6 
replace underwt=1 if waz06<-2 & hv103==1 
*Check with DHS report
svy:tab underwt if info==1

*Wasting
gen wasting=0 if hv103==1 
replace wasting=. if whz06>=6 
replace wasting=1 if whz06<-2 & hv103==1
*Check with DHS report
svy:tab wasting if info==1