lmer Получить значения p от anova

Question

lmerTest был разработан как обертка для оценки p-values из анализа смешанной модели lmer, используя оценку Satterthwaite знаменателей степеней свободы (ddf). Но теперь lmerTest теперь сломан. В настоящее время он возвращает сообщение о том, что была внутренняя ошибка вычисления и возвращает только результат lmer (без p-values). Я смог вычислить p-values из функции summary(), используя отличный код Дэн Мирмана для вычисления оценки Kenward-Rogers ddf. Но я не могу найти эквивалентный код для вычисления p-values в вызове anova на модели lmer. Я подозреваю, что нужно просто подать anova() ddf, но я не могу понять, как это сделать.

Заранее благодарим любого, кто может предложить решение этой проблемы.

Ларри Хансикер

Answer 1

lmerTest возвращает выходные данные ANOVA пакета lme4 всякий раз, когда некоторые вычислительные ошибки в получении приближения Satterthwaite (такие, как, например, при расчете асимптотической дисперсии ковариационной матрицы). LmerTest не нарушен, просто можно привести примеры, когда аппроксимация Саттертуэйт не может быть рассчитана. По моему опыту это происходит не часто.