r - 使用混合效应模型时的收敛警告
问题描述
我进行了实地观察,计算了一个地块中的鸟类数量,我在不同的月份重复了 4 次观察。这样做的目的是了解土地利用是否会影响观察到的鸟类数量。
我的公式是bird.count ~ lu + (1|plot)
bird.count 一个计数数据,有很多零观测值,计数平均值 = 2,方差 = 5.9,lu 是一个有六个级别的因子,而 plot 是一个有 36 个级别的因子,我在其中进行了重复观察。我总共有 144 个观察结果。
基于所有这些,我进行了负二项式混合效应回归,但收到了关于未能收敛的警告,并且系数的 Std.Error 非常高并且都是一样的。
然后,我尝试平均每个地块的 bird.count,然后总共进行了 36 次观察。我这样做是为了不必使用混合效果模型。我假设一个高斯分布运行了一个 glm,但是系数的 Std.Error 非常高并且都是一样的。
我尝试过的事情,
fit.x10 <- glmer.nb (bird.count ~ 1 + (1|plot) , data= birds.data)
There were 13 warnings (use warnings() to see them)
fit.x11 <- glmer.nb (bird.count ~ lu + (1|plot), data= birds.data)
There were 32 warnings (use warnings() to see them)
这没有发出警告
fit.y10 <- glmmTMB(
bird.count ~ (1|plot), data= birds.data,
family = nbinom2, ziformula = ~0, se = TRUE,
verbose = FALSE, doFit = TRUE
)
这没有给出警告,但截距和系数的 Std.Error 很高并且值相同。对于随机效应图,var= 0.1946 和 std= 0.4412
fit.y11 <- glmmTMB(
bird.count ~ lu+ (1|plot), data= birds.data,
family = nbinom2, ziformula = ~0, se = TRUE,
verbose = FALSE, doFit = TRUE
)
这没有发出警告,但截距和系数的 Std.Error 很高并且值相同
fit.z11 <- glm (bird.count ~ lu , data= birds.data.avg, family = gaussian)
关于如何诊断问题或我做错了什么的任何建议?
解决方案
推荐阅读
- python - Pip 升级失败 (Windows 10)
- php - PHP遍历csv的每一行并返回匹配元素的行
- java - 通过 Selenium 和 Java 使用 GeckoDriver 和 Firefox 无法点击链接
- db2 - 根据 DB2 中一个月的第一条和最后一条记录获取字段的值
- c# - HttpRequestMessage 扩展中的 C# TimeSpan 转换
- c - glUseProgram 是否有任何原因可能会使我的程序崩溃?(C Win32)
- apache-kafka - 从 Kafka 到 Vertica 的数据导入(vkconfig Kafka 连接器)
- powerbi - 时间范围内每个用户 ID 的订单数
- haskell - 列表是列还是行?
- python-3.x - 如何在 Windows 上安装 ELKI?