******************************************************************
*STATA syntax code for the computation of index*
*values with the Australian value set (TTO)*
******************************************************************;
gen disut_mo = .
replace disut_mo = 0 if missing(disut_mo) & mobility == 1
replace disut_mo = 0.068 if missing(disut_mo) & mobility == 2
replace disut_mo = 0.374 if missing(disut_mo) & mobility == 3

gen disut_sc = .
replace disut_sc = 0 if missing(disut_sc) & selfcare == 1
replace disut_sc = 0.087 if missing(disut_sc) & selfcare == 2
replace disut_sc = 0.267 if missing(disut_sc) & selfcare == 3

gen disut_ua = .
replace disut_ua = 0 if missing(disut_ua) & activity == 1
replace disut_ua = 0.053 if missing(disut_ua) & activity == 2
replace disut_ua = 0.139 if missing(disut_ua) & activity == 3

gen disut_pd = .
replace disut_pd = 0 if missing(disut_pd) & pain == 1
replace disut_pd = 0.068 if missing(disut_pd) & pain == 2
replace disut_pd = 0.449 if missing(disut_pd) & pain == 3

gen disut_ad = .
replace disut_ad = 0 if missing(disut_ad) & anxiety == 1
replace disut_ad = 0.097 if missing(disut_ad) & anxiety == 2
replace disut_ad = 0.397 if missing(disut_ad) & anxiety == 3

gen disut_mo3sc3 = 0
replace disut_mo3sc3 = -0.064 if mobility == 3 & selfcare ==3
gen disut_mo3ua3 = 0
replace disut_mo3ua3 = 0.025 if mobility == 3 & activity ==3
gen disut_mo3pd3 = 0
replace disut_mo3pd3 = -0.092 if mobility == 3 & pain ==3
gen disut_mo3ad3 = 0
replace disut_mo3ad3 = -0.013 if mobility == 3 & anxiety ==3
gen disut_sc3ua3 = 0
replace disut_sc3ua3 = 0.055 if selfcare == 3 & activity == 3
gen disut_sc3pd3 = 0
replace disut_sc3pd3 = -0.09 if selfcare == 3 & pain == 3
gen disut_sc3ad3 = 0
replace disut_sc3ad3 = -0.105 if selfcare == 3 & anxiety ==3
gen disut_ua3pd3 = 0
replace disut_ua3pd3 = -0.025 if activity ==3 & pain ==3
gen disut_ua3ad3 = 0
replace disut_ua3ad3 = -0.043 if activity ==3 & anxiety ==3
gen disut_pd3ad3 = 0
replace disut_pd3ad3 = -0.185 if pain == 3 & anxiety ==3

gen disut_total = disut_mo + disut_sc + disut_ua + disut_pd + disut_ad + disut_mo3sc3 + disut_mo3ua3 + disut_mo3pd3 + disut_mo3ad3 + disut_sc3ua3 + disut_sc3pd3 + disut_sc3ad3 + disut_ua3pd3 + disut_ua3ad3 + disut_pd3ad3 

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