Question

我正在尝试将lme4包用于glmm并且收集代码为0并且声明：模型未能与max | grad汇合= 0.00791467（tol = 0.001，组分1）。我有兴趣使用lme4包，因为我希望有AIC值来确定适当的模型，因为我添加了其他协变量。

两周前，当我尝试相同的方法时，我收到一条警告信息，表明由于max | grad |，模型无法收敛。问题，但这次没有收到警告消息，只是摘要输出结束时的语句。

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

这是glmer（lme4）型号代码。我增加了maxfun来处理上次运行模型时遇到的其他问题。

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"))

lme4 glmm模型收敛问题

1 个答案: