******************************************************************
*STATA syntax code for the computation of index*
*values with the Iranian value set (TTO)*
******************************************************************;

gen disut_mo = .
replace disut_mo = 0 if missing(disut_mo) & mobility == 1
replace disut_mo = 0.093 if missing(disut_mo) & mobility == 2
replace disut_mo = 0.22 if missing(disut_mo) & mobility == 3

gen disut_sc = .
replace disut_sc = 0 if missing(disut_sc) & selfcare == 1
replace disut_sc = 0.103 if missing(disut_sc) & selfcare == 2
replace disut_sc = 0.235 if missing(disut_sc) & selfcare == 3

gen disut_ua = .
replace disut_ua = 0 if missing(disut_ua) & activity == 1
replace disut_ua = 0.085 if missing(disut_ua) & activity == 2
replace disut_ua = 0.127 if missing(disut_ua) & activity == 3

gen disut_pd = .
replace disut_pd = 0 if missing(disut_pd) & pain == 1
replace disut_pd = 0.075 if missing(disut_pd) & pain == 2
replace disut_pd = 0.149 if missing(disut_pd) & pain == 3

gen disut_ad = .
replace disut_ad = 0 if missing(disut_ad) & anxiety == 1
replace disut_ad = 0.098 if missing(disut_ad) & anxiety == 2
replace disut_ad = 0.205 if missing(disut_ad) & anxiety == 3

gen i3 = (mobility == 3) + (selfcare == 3) + (activity == 3) + (pain == 3) + (anxiety == 3)
replace i3 = i3 - 1
replace i3 = 0 if i3 <=0

gen disut_i3 = 0
replace disut_i3 = i3 *0.024

gen disut_total = disut_mo + disut_sc + disut_ua + disut_pd + disut_ad + disut_i3

gen EQ_index = . 
replace EQ_index = 1 - 0.081 - disut_total
replace EQ_index = 1 if mobility == 1 & selfcare == 1 & activity == 1 & pain == 1 & anxiety == 1 
replace EQ_index = round(EQ_index, .001)