我使用ChoiceModelR分析联合设计。 每个参与者必须回答12个选择集,每个选择集包括3个选择选项和无选择项。 6个变量描述了选择选项。
我对ChoiceModelR的估算数据如下:
> head(dataChoice_train)
participantID ques alt V_1 V_2 V_3 V_4 V_5 V_6 choice.cbc
1 12628880 1 1 3 6 3 2 1 1 2
2 12628880 1 2 1 5 3 1 5 2 0
3 12628880 1 3 2 3 4 2 5 1 0
4 12628880 2 1 4 2 1 2 1 1 2
5 12628880 2 2 6 7 3 3 3 2 0
6 12628880 2 3 1 5 4 1 5 2 0
参与者12628880拥有前12 * 3 = 36行,接下来的36属于参与者12628881,依此类推。
我跑步
hb.post.baseline <- choicemodelr(data=dataChoice_train, xcoding=rep(0, 6),
mcmc=list(R=20000, use=10000),
options=list(save=TRUE,none=TRUE))
> dim(hb.post.baseline$betadraw)
[1] 846 23 1000
我有我所有的846名参与者和我的23个变量实现(不包括每个变量的参考实现)。 我得到了1000个估算值,它与标准保持值10相符。
我担心的是参与者的订购。
我希望一切都没有改变,但我不确定。
我希望dataChoice_train
的前36行(属于参与者12628880)将在hb.post.baseline$betadraw
的顶层(因此在hb.post.baseline$betadraw[1,,]
中表示)。
这样,我可以使用unique(dataChoice_train)
将参与者ID分配给我的betadraws文件。
有人可以确认吗? 是否有更好/更直接的方法将计算出的beta分配给参与者?
谢谢!