在 glmer.nb 中指定 theta.ml 参数答案

【问题标题】：Specify theta.ml parameters within glmer.nb在 glmer.nb 中指定 theta.ml 参数
【发布时间】：2021-02-20 16:13:30
【问题描述】：

我正在尝试使用包lme4 中的glmer.nb() 拟合混合效应负二项式模型。

m1 <- glmer.nb(NumberEvents ~ offset(log(days)) + x + (1|ID), data=dfr, nAGQ = 20)

我收到以下警告：

警告信息：在 theta.ml(Y, mu, weights = object@resp$weights, limit = limit, : 达到迭代限制

我想我应该增加theta.ml的迭代次数限制。

theta.ml(limit = 1000)

怎么做？我如何在glmer.nb() 中调用theta.ml()？任何帮助表示赞赏！

【问题讨论】：

我建议你可以将glmer.nb里面的限制设置为initCtrl=list(limit=1000))
嗨@AlvaroMartinez，我遇到了同样的错误消息，所以我按照你的建议在glmer.nb 中添加了initCtrl=list(limit=1000)。但是 R 显示错误消息unused argument (initCtrl = list(limit = 1000))。我想知道你是否知道如何解决它？

【解决方案1】：

虽然很老了，但我还是先回答一下吧。

一旦达到迭代限制，几乎可以肯定您的数据是等分散或分散不足的（即方差 glmer.nb 正在尝试拟合负二项式模型，随着形状参数趋于无穷大，该模型变得等分散（方差 = 均值）。如果您查看估计的形状参数，您通常会发现它非常大。

适当的解决方案通常是恢复到泊松模型。

【讨论】：