Сравните glmer с glm после импутации мышей.

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

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

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

pool.compare (FM1, FM2)

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


person Eveline Wiegers    schedule 09.06.2019    source источник


Ответы (1)


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

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

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

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

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

person Steffen Moritz    schedule 11.06.2019