r - 在 glmer.nb 中指定 theta.ml 参数

Question

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

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)

怎么做？我怎么能theta.ml()在里面打电话glmer.nb()？任何帮助表示赞赏！

Answer 1

虽然这很老了，但我会作为记录来回答。

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

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