******************************************************************
*STATA syntax code for the computation of index*
*values with the US (Pickard et al) value set (TTO)*
******************************************************************;

gen disut_mo = .
replace disut_mo = 0 if missing(disut_mo) & mobility == 1
replace disut_mo = 0.096 if missing(disut_mo) & mobility == 2
replace disut_mo = 0.122 if missing(disut_mo) & mobility == 3
replace disut_mo = 0.237 if missing(disut_mo) & mobility == 4
replace disut_mo = 0.322 if missing(disut_mo) & mobility == 5

gen disut_sc = .
replace disut_sc = 0 if missing(disut_sc) & selfcare == 1
replace disut_sc = 0.089 if missing(disut_sc) & selfcare == 2
replace disut_sc = 0.107 if missing(disut_sc) & selfcare == 3
replace disut_sc = 0.22 if missing(disut_sc) & selfcare == 4
replace disut_sc = 0.261 if missing(disut_sc) & selfcare == 5

gen disut_ua = .
replace disut_ua = 0 if missing(disut_ua) & activity == 1
replace disut_ua = 0.068 if missing(disut_ua) & activity == 2
replace disut_ua = 0.101 if missing(disut_ua) & activity == 3
replace disut_ua = 0.255 if missing(disut_ua) & activity == 4
replace disut_ua = 0.255 if missing(disut_ua) & activity == 5

gen disut_pd = .
replace disut_pd = 0 if missing(disut_pd) & pain == 1
replace disut_pd = 0.06 if missing(disut_pd) & pain == 2
replace disut_pd = 0.098 if missing(disut_pd) & pain == 3
replace disut_pd = 0.318 if missing(disut_pd) & pain == 4
replace disut_pd = 0.414 if missing(disut_pd) & pain == 5

gen disut_ad = .
replace disut_ad = 0 if missing(disut_ad) & anxiety == 1
replace disut_ad = 0.057 if missing(disut_ad) & anxiety == 2
replace disut_ad = 0.123 if missing(disut_ad) & anxiety == 3
replace disut_ad = 0.299 if missing(disut_ad) & anxiety == 4
replace disut_ad = 0.321 if missing(disut_ad) & anxiety == 5

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

gen EQ_index = . 
replace EQ_index = 1 - disut_total
replace EQ_index = round(EQ_index, .001)