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

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

Платформа Statsmodels LME в настоящее время поддерживает выводы после оценки с помощью тестов Вальда и доверительных интервалов для коэффициентов, анализа правдоподобия профиля, проверки отношения правдоподобия и AIC. [курсив добавлен]

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

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


python statsmodels
person Louis Thibault    schedule 30.05.2015    source источник

Ответы (1)


arrow_upward
6
arrow_downward

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

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

person Kerby Shedden    schedule 04.06.2015
comment
Чтобы быть более точным, llf — это значение критерия REML, если вы используете reml=True для подгонки, и значение логарифмической функции правдоподобия, если вы используете reml=False для подгонки. - person Kerby Shedden; 04.06.2015
comment
Большое спасибо! Как мне изменить оценщик на ML? [Источник] (statsmodels.regression.mixed_linear_model.MixedLM), похоже, не имеет аргументов для указания оценки. Должен ли я вручную установить MixedLM.reml на False? - person Louis Thibault; 05.06.2015
comment
@blz Ты должен это сделать, model.fit(reml=False) - person rightskewed; 05.06.2015