Как я могу выполнить тест отношения правдоподобия на линейной модели смешанного эффекта?

Документация для линейных моделей смешанных эффектов Statsmodels утверждает, что

Структура StatModels LME в настоящее время поддерживает вывод после оценки через тесты Wald и доверительные интервалы на коэффициенты, анализ правдоподобия, тестирование отношения правдоподобия и AIC. [выделено мной]

Я заметил метод MixedLM.loglike, но я не могу найти функцию/метод для запуска теста отношения правдоподобия.

Может кто-нибудь любезно указать мне в правильном направлении?

ОТВЕТЫ

Ответ 1

Я запускаю ветвь разработки, поэтому все может измениться, но класс результатов, возвращаемый MixedLM.fit(), должен иметь атрибут, называемый 'llf'. Это значение функции логарифмического правдоподобия при оцененных параметрах. Если у вас есть две вложенные модели и в 2 раза больше разницы в их значениях llf, то при нулевой гипотезе, где более простая модель истинна, это будет chi ^ 2 случайная величина со степенями свободы, равная разнице в степенях свободы для двух моделей.

Обратите внимание, что многие люди считают, что вы должны переключить оценку на ML (не по умолчанию REML) при использовании тестов LR.