конвергенция ошибка lme4

Question

модель, которую я пытался указать было:

M4 <- glmer(CORT_pgmm~AO_Ause+(1|St),data=belcher,family=Gamma(link="log"), control=glmerControl(optimizer="bobyqa")) 

Делая это, я получил следующее предупреждение с помощью glmer с "bobyqa" и "Nelder_Mead" как оптимизаторы

Warning message: 
In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, : 
    Model failed to converge with max|grad| = 0.0146734 (tol = 0.001, component 1) 

.

Я следовал инструкциям, данным г-ном Бен Болкером в предыдущем answer

Как только я c помещайте оптимизаторов в те, что содержатся в пакете optimx, появляется предупреждающее сообщение.

M4 <- glmer(CORT_pgmm~AO_Ause+(1|St),data=belcher,family=Gamma(link="log"), control=glmerControl(optimizer="optimx",optCtrl=list(method="nlminb"))) 

Тем не менее, я не уверен, как узнать, правильно ли это в моем случае. Как я могу быть уверен в этом?

Приношу извинения, если я не написал это правильно, так как это мой первый пост. Заранее благодарю за ваши ответы.

Answer 1

Вы никогда не можете быть уверены (это численная оптимизация случая, о котором мы не можем доказать много в общем случае), но, как правило, я бы сказал, что если у вас есть , те же предполагаемые «оптимальные» оценки параметров с более чем одним разным оптимизатором, вы можете перестать беспокоиться о сбое конвергенции. Насколько близка «приблизительно» зависит от личного вкуса и научных целей. Для моих типичных случаев использования оценки параметров в пределах (скажем) 1% или консервативно 0,1% друг от друга с разными оптимизаторами считались бы «приблизительно равными».