****************************************************************** *STATA syntax code for the computation of index* *values with the Argentinian value set (VAS)* ******************************************************************; gen m1 =0 gen m2 =0 gen s1 =0 gen s2 =0 gen u1 =0 gen u2 =0 gen p1 =0 gen p2 =0 gen a1 =0 gen a2 =0 replace m1 = 1 if mobility == 2 replace m2 = 1 if mobility == 3 replace s1 = 1 if selfcare == 2 replace s2 = 1 if selfcare == 3 replace u1 = 1 if activity == 2 replace u2 = 1 if activity == 3 replace p1 = 1 if pain == 2 replace p2 = 1 if pain == 3 replace a1 = 1 if anxiety == 2 replace a2 = 1 if anxiety == 3 gen o12 = 1 if (m1 + s1 + u1 + p1 + a1 > 0) & (m2 + s2 + u2 + p2 + a2 == 0) replace o12 = 0 if o12 == . gen o13 = 1 if (m2 + s2 + u2 + p2 + a2 > 0) & (m1 + s1 + u1 + p1 + a1 == 0) replace o13 = 0 if o13 == . gen a23 = 1 if (m1 + s1 + u1 + p1 + a1 > 0) & (m2 + s2 + u2 + p2 + a2 > 0) replace a23 = 0 if a23 == . gen n2g3 = (m1 + s1 + u1 + p1 + a1) if (m2 + s2 + u2 + p2 + a2 > 0) replace n2g3 = 0 if n2g3 == . gen n2sq = (m1 + s1 + u1 + p1 + a1)^2 gen n3sq = (m2 + s2 + u2 + p2 + a2)^2 gen disut_total = .248*m1 + .247*m2 + .184*s1 + .178*s2 + .209*u1 + .148*u2 + .185*p1 + .157*p2 + .150*a1 + .116*a2 + 0.028*o12 + 0.388*o13 + 0.232*a23 + (-0.086*n2g3) + (-0.02*n2sq) + (-0.008*n3sq) gen EQ_index = 1 - disut_total replace EQ_index = round(EQ_index, .001)