1

我正在尝试使用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()?任何帮助表示赞赏!

4

1 回答 1

1

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

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

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

于 2022-01-06T01:13:07.127 回答