繁体   English   中英

拟合随机效果Z20F35E630DAF44DBFA4C3F68F53999999999999999999999999999999999999999DAFAM()而不是GAM()ZC1C425268E68E6855174C174F174140278E608ENENENENENENENENENENENENENENED时,错误

[英]error when fitting random effects model using bam() rather than gam() function in mgcv package, R

I am fitting a model with many random effects using the bam() function within the mgcv package for R. 我的基本 model 结构如下所示:

fit <- bam(y ~ s(x1) + s(x2) + s(xn) + s(plot, bs = 're'), data = dat)

这个 function 适用于我的数据的 4 个子集,但不适用于第五个,这令人惊讶。 相反,它会引发此错误:

Error in qr.qty(qrx, f) : 
  right-hand side should have 14195 not 14196 rows

如果我改用gam()而不是bam() function,这个错误就会消失。 如果我从 model 中删除随机效果,它也会消失。 我真的不确定是什么导致了这个错误,或者如何处理它。 不幸的是,生成一个可重现的示例需要传递一个非常大的数据集,因为它不清楚为什么在这个特定的数据集上抛出这个错误,而其他 4 个数据集拟合完全相同的 model。

任何想法为什么会抛出这个错误,以及如何克服它,将不胜感激。

我有同样的问题,我发现这个 r-help 邮件试图解决同样的问题:

[R] bam (mgcv) 未使用指定数量的核心

看完邮件,我把关于集群的所有代码都删掉了,比如bam() function中的参数cluster 然后错误消息消失。

我不知道细节,但我希望这个技巧能帮助你。

一种可能的原因

qr.qty(qrx, f) 中的错误:右侧应该有 14195 行而不是 14196 行

内存不足。 这可以解释为什么您看到了某些数据集的错误,而没有看到其他数据集。 这在使用大型集群时尤其常见。

暂无
暂无

声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.

 
粤ICP备18138465号  © 2020-2024 STACKOOM.COM