r - 在 glmer.nb 中指定 theta.ml 参数
问题描述
我正在尝试使用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()
?任何帮助表示赞赏!
解决方案
虽然这很老了,但我会作为记录来回答。
一旦达到迭代限制,几乎可以肯定您的数据是等分散或分散不足的(即方差 <= 均值)。glmer.nb
正在尝试拟合负二项式模型,随着形状参数趋于无穷大,该模型变得等分散(方差 = 均值)。如果您查看估计的形状参数,您通常会发现它非常大。
适当的解决方案通常是恢复到泊松模型。
推荐阅读
- android - 无法在 onDataChange 中获取值
- php - 从 ACF 选项将密件抄送收件人添加到 WooCommerce 电子邮件通知
- flask - python flask base.html 不显示 {{user.username}
- r - 寻找 R 函数来检测 tibble 中的重复值或相同值
- python - How can I generate this plot?
- javascript - RTCPeerConnections(或 mediaStreams?)即使在关闭后也会占用内存和 CPU
- python - 如何让beautifulsoup更快地获取抓取网页中的所有元素?
- c# - NpgsqlConnectionStringBuilder 是否能够验证连接字符串?
- php - 当 PHP array_key_exists 为 FALSE 时,为什么我不能添加键和值?
- python - 学习 Python 并遇到 PIP 安装模块的问题