Смешанная биноминальная смешанная модель с посредничеством - PullRequest
Новая
       77

Смешанная биноминальная смешанная модель с посредничеством

0 голосов
/ 24 марта 2020

Я использую пакет lme4 в R, чтобы подогнать отрицательную биномиальную модель, которую я хотел бы использовать в посредничестве (подходит лучше, чем Пуассон).

Я получаю ошибку: "семейство glmer для модель результата не поддерживается "

Для простоты я просто переключил класс на lmer; например, от 

> class(model.Y)
[1] "glmerMod"
attr(,"package")
[1] "lme4"

до 

> class(model.Ynew)
[1] "lmerMod"
attr(,"package")
[1] "lme4"

и все работает. Очевидно, это может создать проблемы, но я проверил невязки, подгонял и предсказал значения, и никаких изменений не произошло. Результаты также согласуются с результатами модели Пуассона. 

> summary(residuals(model.Y)-residuals(model.Ynew))
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
      0       0       0       0       0       0 
> summary(fitted(model.Y)-fitted(model.Ynew))
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
      0       0       0       0       0       0 
> summary(predict(model.Y)-predict(model.Ynew))
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
      0       0       0       0       0       0

Есть ли у кого-нибудь понимание последствий этого или возможного пути продвижения вперед, который не настолько ad ho c?

...