Fatal编程技术网
  • r/
  • R Bam无法生成稀疏矩阵？

R Bam无法生成稀疏矩阵？

r

R Bam无法生成稀疏矩阵？,r,gam,mgcv,R,Gam,Mgcv,我正在尝试制作一个包含大约630万行数据的二项GAM（因此是bam）。我有两个随机效应（40级和7级），使用s（…，“re”），3个分类变量（5+5+4级）和一个循环平滑（s（Hour，by=GroupAB，bs=“cc”））。GroupAB有25个级别，是前面提到的两个变量之间的连接，我在平滑外部包含了这些变量之间的额外交互。类似于（y~s（…，，re）+s（…，，re）+s（Hour，by=GroupAB，“cc”）+A*B+C，family=“binomial”） 运行模型时，我收到以下错

我正在尝试制作一个包含大约630万行数据的二项GAM（因此是
bam
）。我有两个随机效应（40级和7级），使用
s（…，“re”）
，3个分类变量（5+5+4级）和一个循环平滑（
s（Hour，by=GroupAB，bs=“cc”）
）。GroupAB有25个级别，是前面提到的两个变量之间的连接，我在平滑外部包含了这些变量之间的额外交互。类似于
（y~s（…，，re）+s（…，，re）+s（Hour，by=GroupAB，“cc”）+A*B+C，family=“binomial”）

运行模型时，我收到以下错误消息：

Error in predict.gam(G, newdata = mf[ind, ], type = "lpmatrix", newdata.guaranteed = TRUE,  : 
  number of items to replace is not a multiple of replacement length
我追踪到bam试图生成稀疏矩阵（
sparse.matrix.model（）
in）。看到它使用了
bam（）
的
chunk.size=
参数，我开始玩弄它。将chunck尺寸从10k增加到75K后，我使模型运行并给出合理的结果

但首先为什么会出现错误呢？是否有技巧避免此错误消息？还是与数据有关