Случайные эффекты с моделями count

Я пытаюсь создать модель препятствий со случайными эффектами либо в r, либо в stata. Я просмотрел пакет glmmADMB, но у меня возникли проблемы с его загрузкой в ​​​​R, и я не могу найти документацию по пакету в Cran. Этот пакет еще доступен? Кто-нибудь успешно использовал его для оценки барьерной модели со случайными эффектами?

В качестве альтернативы, есть ли способ оценить это в stata? Есть ли способ оценить случайные эффекты с любым типом данных подсчета в stata?

Мы будем очень признательны за любые советы.

Дженнифер


person user2418436    schedule 24.05.2013    source источник
comment
Не знаю насчет Stata, но пакет cmp может на это способен: stata-journal.com/article.html?article=st0224 glmmADMB нет в CRAN по техническим причинам, он есть на R-forge. См. glmmadmb.r-forge.r-project.org... MCMCglmm Пакет также может создавать модели с препятствиями со случайными эффектами.   -  person Ben Bolker    schedule 24.05.2013
comment
Кстати, вы предпочитаете R или Stata?   -  person Ben Bolker    schedule 25.05.2013


Ответы (1)


В Stata xtnbreg и xtpoisson по умолчанию используется оценка случайных эффектов. Вы всегда можете оценить две части отдельно вручную. См. главу о подсчетах в книге Cameron and Trivedi's Stata для перекрестных примеров.

У вас также есть написанные пользователем hplogit и hnlogit для моделей подсчета препятствий. Они используют логит/пробит для первого этапа и усеченный нулем пуассоновский/отрицательный бином для второго этапа. Кроме того, хорошим подходом может быть конечная смешанная модель (см. написанный пользователем fmm). Есть еще ztpnm. Все это модели в разрезе.

person Dimitriy V. Masterov    schedule 24.05.2013