use "...IAIR7EFL.DTA" , clear

gen age=b19_01
gen period=60

* tetanus toxoid injections
recode m1_1 (0 1 8 9 . = 0 "No") (2/7 = 1 "Yes"), gen(rh_anc_toxinj)
replace rh_anc_toxinj = . if age>=period
label var rh_anc_toxinj "Received 2+ tetanus injections during last pregnancy"
       
* protection against neonatal tetanus
       
* temporary vars needed to compute the indicator
gen totet = 0 
gen ttprotect = 0                          
replace totet = m1_1 if (m1_1>0 & m1_1<8)
replace totet = m1a_1 + totet if (m1a_1 > 0 & m1a_1 < 8)
                              
* generating variable for date of last injection - will be 0 for women with at least 1 injection at last pregnancy
gen lastinj = 9999
replace lastinj = 0 if (m1_1>0 & m1_1 <8)
gen int ageyr = (age)/12 

* years ago of last shot - (age at of child), yields some negatives
replace lastinj = (m1d_1 - ageyr) if m1d_1 <20 & (m1_1==0 | (m1_1>7 & m1_1<9996)) 

* generate summary variable for protection against neonatal tetanus 
replace ttprotect = 1 if tet2lastp ==1 

* at least 2 shots in last 3 years
replace ttprotect = 1 if totet>=2 &  lastinj<=2 

* at least 3 shots in last 5 years
replace ttprotect = 1 if totet>=3 &  lastinj<=5

* at least 4 shots in last 10 years
replace ttprotect = 1 if totet>=4 &  lastinj<=9 

* at least 2 shots in lifetime
replace ttprotect = 1 if totet>=5  
label variable ttprotect "Full neonatal tetanus Protection"
                              
gen rh_anc_neotet = ttprotect
replace rh_anc_neotet = . if  bidx_01!=1 | age>=period 
label variable rh_anc_neotet "Protected against neonatal tetanus"
       
tab rh_anc_toxinj [iweight=v005/1000000]
tab rh_anc_neotet [iweight=v005/1000000]