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