Реализация субъекта как случайного эффекта в модели смешанного логита в R
Извиняюсь, если это повторная публикация, но я не мог найти, что это покрыто непосредственно и на языке мирян в другом месте. Кроме того, поскольку это вопрос реализации, а не статистики, я подумал, что лучше отправлять сообщения в SO, а не в CV, но, пожалуйста, дайте мне знать, если я должен поступить иначе.
У меня есть набор данных, в котором индивидуумы, каждый из которых принадлежит к определенной группе, неоднократно выбирали между несколькими дискретными результатами. Что-то похожее на:
subID group choice
1 Red A
1 Red B
2 Blue B
2 Blue B
2 Blue C
3 Red A
3 Red B
Меня интересует, как членство в группах предсказывает выбор, и я хочу объяснить несамостоятельность наблюдений из-за неоднократных выборов, сделанных одними и теми же людьми. В свою очередь, я планировал реализовать смешанную полиномиальную регрессию, рассматривая subID как случайный эффект.
Казалось бы, mlogit - это общий маршрут для смешанных логитов. Однако мне немного неясно, как установить желаемый случайный эффект. Я знаком с установкой предметов в качестве случайных факторов для линейных моделей в lmer() (например, Y ~ group + (1|subID), data=base), но при рассмотрении обсуждений, которые были здесь и здесь, я немного скинул по схеме млогит. Могу ли я действительно учитывать subID через аргумент rpar?
Любое руководство будет тщательно оценено. Благодарю.