[英]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 邮件试图解决同样的问题:
看完邮件,我把关于集群的所有代码都删掉了,比如bam()
function中的参数cluster
。 然后错误消息消失。
我不知道细节,但我希望这个技巧能帮助你。
一种可能的原因
qr.qty(qrx, f) 中的错误:右侧应该有 14195 行而不是 14196 行
内存不足。 这可以解释为什么您看到了某些数据集的错误,而没有看到其他数据集。 这在使用大型集群时尤其常见。
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.