lme4 glmm模型收敛性问题

Question

我正在尝试使用lme4包进行glmm，并得到收敛代码0和一条声明: Model to converge with max|grad| = 0.00791467 (tol = 0.001，component 1)。我对使用lme4包很感兴趣，因为我希望在添加其他协变量时使用AIC值来确定适当的模型。

两周前，当我尝试相同的方法时，我收到了一条警告消息，指出由于max|grad|问题，模型无法收敛，但这次没有收到警告消息，只收到了汇总输出末尾的语句。

这是否意味着模型不收敛？我还使用了glmmPQL方法。两种模型之间的系数参数估计是相似的。

下面是glmer (lme4)模型代码。我增加了最大乐趣来处理我上次运行模型时遇到的其他问题。

l1<-glmer(Meat_Weight~logsh+SAMS_region_2015+(1|StationID),
        family="Gamma"(link="log"),data=datad,control=glmerControl(optCtrl=list(maxfun=100000)))

这是glmmPQL代码。

m1<-glmmPQL(fixed=Meat_Weight~logsh+SAMS_region_2015,random=~1|StationID,
        family=Gamma(link="log"),data=datad)

我确信这不是诊断问题的信息，但如果有人有建议，我可以提供更多的数据。

谢谢

Answer 1

尝试更改优化器

    l1<-glmer(Meat_Weight~logsh+SAMS_region_2015+(1|StationID),
    family="Gamma"(link="log"),data=datad, control = glmerControl(optimizer="bobyqa"))