Изменяет ли ChoiceModelR порядок дел?
Я использую ChoiceModelR, чтобы проанализировать совместный дизайн. Каждый участник должен был ответить на 12 наборов выбора, каждый из которых состоял из 3 вариантов выбора плюс отсутствие выбора. 6 переменных описали варианты выбора.
Мои imput-данные для 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.
Мой вопрос:
То, что меня беспокоит, это порядок моих участников. Я ожидаю, что ничего не изменилось, но я не уверен. Я ожидал бы первые 36 рядов dataChoice_train
(принадлежащий участнику 12628880) будет представлен в верхнем слое hb.post.baseline$betadraw
(так в hb.post.baseline$betadraw[1,,]
).
Таким образом, я мог бы использовать unique(dataChoice_train)
назначить идентификатор участника моему файлу бета-розыгрышей.
Может ли кто-нибудь это подтвердить? Есть ли лучший / более прямой метод назначения рассчитанных бета-версий участникам?
Заранее спасибо!