by johnr » Tue Jan 10, 2012 11:11 am
The good news is that we working on setting up Ngene to support this type of problem. We hope to have the next version of Ngene (V1.1.1) released sometime soon. For the present, the syntax proposed for this type of problem will look something like this:
Design
;alts = alt1, alt2
;rows = 12
;eff = (mnl,d)
;model:
U(alt1) = b1[0.2] * pbEa[0.2,0.4] * Early[10,12,14] + b2[0.5]* pbOn[0.5,0.3] * Ontime[20,22,24]
+ b3[-0.4] * pbLa[fcn(1 - alt1.pbEa - alt1.pbOn)] * Late[25,27,29] /
U(alt2) = b1 * pbEa * Early + b2 * pbOn * Ontime + b3 * pbLa[f(1-alt2.pbEa -alt2.pbOn)] * Late $
Note that this is somewhat more complex than at first glance. This is because your W (or in the above my pbLa) must sum to 1 (or 100) – which is a new feature in Ngene. This will be handled via the fcn command above which makes the level of an attribute a function of the levels of another attribute. Secondly, it also requires the estimation of an interaction term without a main effect being estimated (b2 * A*W where A and W do not have priors - in the above my pbLa and Late for example are interactions but not main effects). This also is a new feature.
As you can hopefully appreciate, this requires substantive testing in that we need to test thoroughly any new feature to ensure that it is compatible with all other functionality – e.g., what will happen if the design is orthogonal and you want to do this – this attribute cannot be orthogonal hence we need to ensure that the program doesn’t crash when users do things that are not possible. We also have to test that it is compatible with other model types, with other features such as pivoting, etc. This is a non-trivial task
Thus stay tuned. It is coming and hopefully, subject to there being no substantive compatibility issues, sometime very soon.
John