However, when I run the code below (which has 2 models, and reject conditions ), it seems to impose the restrictions appropriately. Has there been an update? Should I trust the design that's being generated? I am using version 1.4.0
thanks
Michael
- Code: Select all
Design
;alts(m1)=alt1, alt2,alt3
;alts(m2)=alt1, alt2,alt3
;rows=1000
;block=200
;eff=m1(mnl,d)+m2(mnl,d)
;require:
alt2.m<=alt1.m,
alt2.cost<alt3.cost,
alt2.time<=alt1.time,
alt2.risk<=alt1.risk
;reject:
alt1.risk >= 2 and alt2.risk <1,
alt1.time > 0.360 and alt2.time <0.180
;alg = mfederov
;model(m1):
U(alt1)=
mb[-0.1]*m[5,4,3,2,1]+
t[-0.1]*time[0.007,0.014,0.03,0.045,0.09,0.182,0.365,0.73,1.825,3.65,7.3]+
r[-0.1]*risk[0.01,0.02,0.04,0.08,0.15,0.26,1,2,3, 5,10,15,26]
/
U(alt2)=
mb[-0.1]*m[5,4,3,2,1]+
t[-0.1]*time[0.007,0.014,0.03,0.045,0.09,0.182,0.365,0.73,1.825,3.65,7.3]+
r[-0.1]*risk[0.01,0.02,0.04,0.08,0.15,0.26,1,2,3, 5,10,15,26]+
cost[-0.0002]*cost[5,10,15,25, 50, 100, 250,500,1000,2000,5000,10000]
/
U(alt3)=cost*cost
;model(m2):
U(alt1)=
mb[-0.2]*m[5,4,3,2,1]+
t[-0.3]*time[0.007,0.014,0.03,0.045,0.09,0.182,0.365,0.73,1.825,3.65,7.3]+
r[-0.4]*risk[0.01,0.02,0.04,0.08,0.15,0.26,1,2,3, 5,10,15,26]
/
U(alt2)=
mb[-0.2]*m[5,4,3,2,1]+
t[-0.3]*time[0.007,0.014,0.03,0.045,0.09,0.182,0.365,0.73,1.825,3.65,7.3]+
r[-0.4]*risk[0.01,0.02,0.04,0.08,0.15,0.26,1,2,3, 5,10,15,26]+
cost[-0.0002]*cost[5,10,15,25, 50, 100, 250,500,1000,2000,5000,10000]
/
U(alt3)=cost*cost
$