by nastay » Fri Mar 04, 2022 7:39 pm
Dear Michiel,
Thank you very much for your replay and suggestion. I have revised the Ngene codes, would you please see if all looks fine with the codes? I really appreciate your help.
I am working with following attributes and levels:
1- Age (10, 50, 75 years)
2- Life years without treatment- LY-UNTREAT (3, 5, 10)
3- Quality of life- QoL (50,75,90)
4- Effect of tretament as life years after treatment-LY-TREAT (0.5,1,3)
5- Distribution of effect- DISTRIBUTION (all get the same effect in average, most are under average while 5% get 3xaverage effect (1.5, 3,9), half get minimal effect while the other half get double (1,2,6))
6- Risk of known and unknown side effect- RISK (No, 10% risk of known and no risk of unknown, 10% risk of known and there is risk of unknown as well).
I applied a nonlinear transformation on age and QoL. Both life years attributes are continuous and I assume there is linear relationship in attribute levels. DISTRIBUTION and RISK are categorical and I used effects coding.
Design
;alts = altA, altB
;rows = 24
;block = 3
;eff = (mnl, d)
;model:
U(altA) = b1[0]
+ b2 [0] * age[2.30,3.91,4.32]
+ b3 [0] * LY-UNTREAT[3,5,10]
+ b4 [0] * QoL[3.91,4.32,4.50]
+ b5[0] * LY-TREAT[0.5,1,3]
+ b6.effects[0|0] * DISTRIBUTION[1,2,3]
+ b7.effects[0|0] * RISK[1,2,3] /
U(altB) = b2 * age
+ b3 * LY-UNTREAT
+ b4 * QoL
+ b5 * LY-TREAT
+ b6.effects * DISTRIBUTION
+ b7.effects * RISK
$