Сравните glmer с glm после вменения мышей - PullRequest
1 голос
/ 09 июня 2019

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

fm1 <- с (dti.mice1, glmer (Лечение ~ (1 | Больница) + Возраст)) </p>

fm2 <- с (dti.mice1, glm (Лечение ~ Возраст))</p>

pool.compare (fm1, fm2)

Ошибка: модель «fit1» не больше, чем «fit0»

1 Ответ

0 голосов
/ 12 июня 2019

Похоже, проблема в том, что fm1 и fm2 не распознаются как вложенные модели (при этом fm1 является более крупной моделью).

Из документации мыши :

fit1
Объект класса 'mira', созданный with.mids ().

fit0
Объект класса'mira', созданный with.mids ().Модель в fit0 представляет собой вложенную форму fit0 для fit1.

Таким образом, полная / более крупная модель должна соответствовать fit1.Который имеет место для вашего примера.Возможно, вы не можете использовать glm и glmer вместе.Что произойдет, если вы используете glm для обеих моделей?Это работает?

...