Hi all Ngeners
I am working on a design for a choice experiment with the following properties...
- Three alternatives (two plus one status quo (SQ))
- 4 attributes, with 3, 4, 3 and 4 levels respectively
- The SQ is defined by the same attributes but the levels are fixed - the SQ gives a fixed utility.
- The levels in the SQ are fixed so that there can be both better and worse scenarios. The SQ is defined by the second worse levels on all attributes except for the cost, which is zero.
- It is not relevant with an alternative without improvement in any attribute and a zero or positive cost.
I have some issues that I really need some advice on...
- I have tried to use the ;require command to secure that my SQ has the correct levels. The problem is that Ngene mostly hangs up and I have to restart the software. If I reduce the number of candidates quite a lot (e.g.=10) it sometimes works but it is not really stable. It does not feel okay to reduce the number of candidates that much! Any other suggestions?
- I can of course just use a constant for the SQ alternative (set the constant equal to the sum of priors and attribute levels), but is that in any way problematic? I have used the ;reject command to get rid of "dominating" alternatives but I don't know if this is a good approach.
- If I use a constant for the SQ and the Bayesian approach, should I calculate the constant utility by adding the "means"? Without the Bayesian approach it is straightforward to just "add the priors".
Please, any advice and/or suggestions are very welcome!
Below you find my syntax
Design
;alts=A*,B*,SQ
;rows=8
;alg=mfederov(candidates=100)
;eff=(mnl,d,mean)
;con
;reject:
A.clear<=1 and A.fish<=1 and A.bottom<=1,
B.clear<=1 and B.fish<=1 and B.bottom<=1,
A.clear>1 and A.fish=1 and A.bottom=1 and A.cost=0,
A.clear=1 and A.fish>1 and A.bottom=1 and A.cost=0,
A.clear=1 and A.fish=1 and A.bottom=2 and A.cost=0,
B.clear>1 and B.fish=1 and B.bottom=1 and B.cost=0,
B.clear=1 and B.fish>1 and B.bottom=1 and B.cost=0,
B.clear=1 and B.fish=1 and B.bottom=2 and B.cost=0
;require:
SQ.clear=1,
SQ.fish=1,
SQ.bottom=1,
SQ.cost=0
;model:
U(A)=b2.dummy[(u,0.2,0.4)|(u,0.1,0.3)]*clear[2,1,0]
+b3.dummy[(u,0.2,0.4)|(u,0.1,0.3)|(u,0,0.1)]*fish[3,2,1,0]
+b4.dummy[(u,0.1,0.3)|(u,0,0.1)]*bottom[2,1,0]
+b5[(n,-0.004,0.002)]*cost[0,30,50,70]
/
U(B)=b2*clear+b3*fish+b4*bottom+b5*cost
/
U(SQ)=b1[0]
+b2*clear+b3*fish+b4*bottom+b5*cost
$