sort age
egen sequence=seq(),by(caseid)
tab sequence
drop if sequence==1
drop sequence

*Generate weight
iw=v005/1000000

*DHS-7 and later surveys will have the b19 variable, age in months 
*To create age in months in older surveys, run the following lines
               capture confirm variable b19
               if _rc { // b19 does not exist, so create equivalent for old calculation method
                              gen b19 = v008 - b3
                              label variable b19 "Age of child in months or months since birth"
               }



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

*PNC - 2008- DHS surveys
               
               *pnc timing         
               recode m71 ( 100/202 =2 "PNC within 2 days") (203/305= 1 "PNC 3-41 days") (900/1000 =9 "Don't know, missing") (else=0 "No PNC or after 41 days") if bidx==1, gen(pnctimebaby) 
                recode pnctimebaby (1 0 9= 0 "No PNC, pnc 2+days or dk/miss") (2=1 "PNC within 2 days") if bidx==1, gen(pnc2dbaby) 
                lab var pnctimebaby "Timing of PNC for child"
               lab var pnc2dbaby "Postnatal care within 2 days"
               gen pnc2dbaby = 0
               ta pnc2dbaby if bidx==1 & b19<24 [iw=wt], m
               
               *provider of pnc 
                recode m72 (11/29 = 1 "Skilled provider") (96 =0 "Non-skilled provider") ( else = 9 "Don't know, missing") if bidx==1, gen (pncprovbaby) 
                gen pncprovbaby = 0 
               
               *pnc within 2 days
               *PNC calculated if it was given by provider                           
               gen skillpnc2dbaby = 0 
               replace skillpnc2dbaby = 1 if pnc2dbaby==1 & pncprovbaby==1 
               lab var skillpnc2dbaby "Postnatal care within 2 days by skilled provider"
               lab val skillpnc yesno
               ta skillpnc2dbaby if bidx==1 & b19<24 [iw=wt],m

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

*For DHS7 surveys:

*!RUN ABOVE CODE AND THE FOLLOWING

               *generate general variable for if they had a check at all or not
               gen check =0 
               replace check =1 if (m70==1 |m74==1)
               
               *for those newer surveys
               replace skillpnc2dbaby = 0 if check ==0 
               
               *Now look at babies born in a facility, if checked before discharge             
               *provider of pnc check - health facility-born babies
               recode m76 (11/29 = 1 "Skilled provider") (96 =0 "Non-skilled provider") ( else = 9 "Don't know, missing")  if check==1 , gen (pncprovbaby_hf)

               *code timing of check if health facility delivery
               recode m75 ( 100/202 =2 "PNC within 2 days") (203/305= 1 "PNC 3-41 days") (900/1000 =9 "Don't know") (else=0 "No PNC or after 41 days")  if check==1, gen(pnctimebaby_hf)
               ta pnctimebaby_hf  [iw=wt]
               
               recode pnctimebaby_hf (1 0 9= 0 "No PNC, pnc 2+days or dk/miss") (2=1 "PNC within 2 days") if check==1 , gen(pnc2dbaby_hf)
               lab var pnc2dbaby_hf "Postnatal care within 2 days"
               ta pnc2dbaby_hf  bidx if b19<24 [iw=wt], m

               *PNC calculated if it was given by medically trained provider                       
               gen skillpnc2dbaby_hf = 0 if bidx==1 
               replace skillpnc2dbaby_hf = 1 if pnc2dbaby_hf==1 & pncprovbaby_hf==1 &check==1
               lab var skillpnc2dbaby_hf "Postnatal care within 2 days by skilled provider"
               ta skillpnc2dbaby_hf bidx if b19<24 [iw=wt], m

               ************
               *combine two timing variables
               gen newpnctimebaby =  skillpnc2dbaby_hf
               *if they delivered at home or only got checked after discharge
               replace newpnctimebaby = 1 if skillpnc2dbaby==1 & (facdel==0 | facdel == 9 | (facdel ==1 & m74==0 ))
               label var newpnctimebaby "Postnatal care within 2 days by skilled provider"
               tab newpnctimebaby bidx if b19<24  [iw=wt], col