Я работаю над набором данных панели и пытаюсь запустить регрессию логита с фиксированными эффектами.
Я обнаружил, что модели glmer из пакета lme4 и пакета bife подходят для такой работы.
Однако, когда я запускаю регрессию с каждой моделью, у меня не совпадают результаты (оценки, стандартные ошибки и т. Д.)
Вот код и результаты для модели блеска с перехватом:
glmer_1 <- glmer(CVC_dummy~at_log + (1|year), data=own, family=binomial(link="logit"))
summary(glmer_1)
Estimate Std. Error zvalue Pr(>|z|)
(Intercept) -6.43327 0.09635 -66.77 <2e-16 ***
at_log 0.46335 0.01101 42.09 <2e-16 ***
Без перехвата:
glmer_2 <- glmer(CVC_dummy~at_log + (1|year)-1, data=own, family=binomial(link="logit"))
summary(glmer_2)
Estimate Std.Error z value Pr(>|z|)
at_log 0.46554 0.01099 42.36 <2e-16 ***
И с пакетом bife:
bife_1 <- bife(CVC_dummy~at_log | year, data=own, model="logit")
summary(bife_1)
Estimate Std. error t-value Pr(> t)
at_log 0.4679 0.0110 42.54 <2e-16 ***
Почему оценочные коэффициенты at_log отличаются между двумя пакетами?
Какой пакет я должен использовать?