Самый простой пакет R для установки субъекта в качестве случайного эффекта в смешанной модели логита - PullRequest
Новая
       122

Самый простой пакет R для установки субъекта в качестве случайного эффекта в смешанной модели логита

0 голосов
/ 03 марта 2019

У меня есть набор данных, в котором отдельные лица, каждый из которых принадлежит к определенной группе, неоднократно выбирали между несколькими дискретными результатами. 

subID  group   choice
1      Big     A
1      Big     B
2      Small   B
2      Small   B
2      Small   C
3      Big     A
3      Big     B
.       .      .
.       .      .

Я хочу проверить, как членство в группе влияет на выбор, и хочу учитывать не-Независимость наблюдений из-за повторного выбора, сделанного одними и теми же лицами.В свою очередь, я планировал реализовать группу лечения смешанной полиномиальной регрессии в качестве фиксированного эффекта и subID в качестве случайного эффекта.Кажется, что есть несколько вариантов для многочленных логитов в R, и я надеюсь, что некоторые рекомендации, которые могут быть наиболее легко реализованы для этой смешанной модели:

1) multinom -GLM через nnet позволяет использовать функцию multinom.Это выглядит как хороший, понятный, простой вариант ... для фиксированных моделей эффектов.Однако есть ли способ реализовать случайные эффекты с мултином? Предыдущий пост CV предполагает, что multinom способен обрабатывать GLM со смешанными эффектами с распределением Пуассона и лог-ссылкой.Однако я не понимаю (а), почему это так или (б) требуемый синтаксис.Кто-нибудь может уточнить?

2) mlogit - фантастический пакет с невероятно полезными виньетками.Однако документация «смешанный логит» относится к моделям, которые имеют случайные эффекты, связанные с альтернативными специфическими ковариатами (реализовано с помощью аргумента rpar).Моя модель не имеет альтернативных специфических переменных;Я просто хочу учесть случайные перехваты участников.Это возможно с Mlogit?Эта разница автоматически учитывается путем установки subID в качестве id.var при формировании данных в длинную форму с помощью mlogit.data ?РЕДАКТИРОВАТЬ: я только что нашел пример «обманом» mlogit, чтобы обеспечить случайные коэффициенты для переменных, которые различаются между людьми (очень низко здесь ), но я не совсем понимаю синтаксис.

3) MCMCglmm , очевидно, еще один вариант.Однако, как относительный новичок с R и кто-то совершенно незнакомый с байесовской статистикой, мне лично не удобно разбирать пример синтаксиса смешанных логитов с этим пакетом , или, даже следуя синтаксису, угадывать приоры илидругие необходимые аргументы.

Будем весьма благодарны за любые указания относительно наиболее простого подхода и его реализации синтаксиса.Мне также интересно, если случайный эффект subID должен быть вложен в группу (поскольку люди являются членами групп), но это может быть вопросом для CV вместо этого.В любом случае, большое спасибо за любые идеи.

...