Я выполнил логистическую оценку со взвешенными данными опроса и кластерными SE, используя svyglm из пакета Survey.Поскольку я хочу использовать тест отношения правдоподобия для сравнения двух моделей, мне нужна вероятность журнала, как мне извлечь его из объекта svyglm?
logLik(mod1)
возвращает значение 15426.78, которое кажется отключенным, и дополнительно предупреждает, что svyglmне соответствует максимальному правдоподобию.
Когда я использую stargazer для экспорта своего вывода, он обеспечивает логарифмическое правдоподобие и AIC, которые равны -7,350.223 и 14,724.450 соответственно.Как Stargazer получает эти значения?Тем более, что
mod1$aic
возвращает 14 780,21
Любая помощь приветствуется.