Priors from pilot
Posted: Tue May 18, 2021 4:13 am
Hi,
I’ve done a pilot study of an unlabeled design on food choice and health labelling, and I’m trying to use the coefficients to get an efficient design with fixed priors. I have a couple of questions I need help with:
Firstly, I’ve realized I made a mistake when designing the pilot… I added an alternative-specific constant to my unlabeled design. I think this has meant that Ngene checked for repeated alternatives and row repetition, but NOT for dominance. Therefore, there were some dominated alternatives in my design (4 out of 24 rows, 2 per block, so I hope not disastrous). I’d checked, but there were always some dominated alternatives, so I thought that was just the best that could be done with my syntax!
The coefficients are of the expected sign and mostly significant…
Would you still recommend using them for the pilot?
Many thanks for the help!
I add below the original syntax, and the coefficients
Coefficients are shown as effects coded, alt3 is an opt-out
; alts = alt1*, alt2*, alt3
; rows = 24
; block = 2
; eff = (mnl, d)
; alg = swap(stop=total(10 mins))
; model:
U(alt1) = b1 + b2[-0.0001]*price[60,80,100,120]
+ b3.effects[0.0001]*low_salt[1,0]
+ b4.effects[0.0001]*safety_certified[1,0]
+ b5.effects[0.0001]*no_antibiotics[1,0]
+ b6.effects[0.0001]*originSA[1,0]
/
U(alt2) = b2*price
+ b3*low_salt
+ b4*safety_certified
+ b5*no_antibiotics
+ b6*originSA
$
VARIABLES Priors
Price -0.00170 (0.00264)
SaltContentD 0.307*** (0.0649)
SafetycertificationD 0.377*** (0.118)
RoutineantibioticuseD 0.0141 (0.0534)
OriginD 0.182*** (0.0566)
Standard errors in parentheses
*** p<0.01, ** p<0.05, * p<0.1
I’ve done a pilot study of an unlabeled design on food choice and health labelling, and I’m trying to use the coefficients to get an efficient design with fixed priors. I have a couple of questions I need help with:
Firstly, I’ve realized I made a mistake when designing the pilot… I added an alternative-specific constant to my unlabeled design. I think this has meant that Ngene checked for repeated alternatives and row repetition, but NOT for dominance. Therefore, there were some dominated alternatives in my design (4 out of 24 rows, 2 per block, so I hope not disastrous). I’d checked, but there were always some dominated alternatives, so I thought that was just the best that could be done with my syntax!
The coefficients are of the expected sign and mostly significant…
Would you still recommend using them for the pilot?
Many thanks for the help!
I add below the original syntax, and the coefficients
Coefficients are shown as effects coded, alt3 is an opt-out
; alts = alt1*, alt2*, alt3
; rows = 24
; block = 2
; eff = (mnl, d)
; alg = swap(stop=total(10 mins))
; model:
U(alt1) = b1 + b2[-0.0001]*price[60,80,100,120]
+ b3.effects[0.0001]*low_salt[1,0]
+ b4.effects[0.0001]*safety_certified[1,0]
+ b5.effects[0.0001]*no_antibiotics[1,0]
+ b6.effects[0.0001]*originSA[1,0]
/
U(alt2) = b2*price
+ b3*low_salt
+ b4*safety_certified
+ b5*no_antibiotics
+ b6*originSA
$
VARIABLES Priors
Price -0.00170 (0.00264)
SaltContentD 0.307*** (0.0649)
SafetycertificationD 0.377*** (0.118)
RoutineantibioticuseD 0.0141 (0.0534)
OriginD 0.182*** (0.0566)
Standard errors in parentheses
*** p<0.01, ** p<0.05, * p<0.1