Ошибка при оценке CI для GLMM с использованием confint () - PullRequest
Новая
       88

Ошибка при оценке CI для GLMM с использованием confint ()

0 голосов
/ 02 ноября 2018

У меня есть набор GLMM с двоичной переменной отклика и набором непрерывных переменных, и я хотел бы получить доверительные интервалы для каждой модели. Я использовал функцию confint(), 95% и метод profile, и она работает без проблем, если она применяется к модели без взаимодействий.

Однако, когда я применяю confint() к модели с взаимодействиями (непрерывно * непрерывно), я получаю эту ошибку: 

> m1CI <- confint(m1, level=0.95, method="profile")
Error in zeta(shiftpar, start = opt[seqpar1][-w]) : 
  profiling detected new, lower deviance

Модель работает без проблем (хотя я применил оптимизатор, потому что у некоторых моделей были проблемы с конвергенцией), и вот окончательная форма одной из них: 

m1 <- glmer(Use~RSr2*W+RSr3*W+RShw*W+RScon*W+RSmix*W+(1|Pack/Year),control=glmerControl(optimizer="bobyqa", optCtrl=list(maxfun=100000)), data = data0516RS, family=binomial(link="logit"))

Кто-нибудь знает, почему это происходит, и как я могу это решить?

Я использую версию R 3.4.3 и lme4 1.1-17

1 Ответ

0 голосов
/ 27 марта 2019

Проблема была решена с помощью следующих инструкций: https://stat.ethz.ch/pipermail/r-sig-mixed-models/2014q3/022394.html

Я использовал confint.merMod из пакета lme4 и увеличил параметр 'devtol', сначала до 1e-8, который не работал для моих моделей, а затем до 1e-7. С этим значением сработало

...