Я подбираю модель со множеством случайных эффектов, используя функцию bam() в пакете mgcv для R. Моя основная структура модели c выглядит так:
fit <- bam(y ~ s(x1) + s(x2) + s(xn) + s(plot, bs = 're'), data = dat)
Эта функция работает для 4 подмножества моих данных, но не пятое, что удивительно. Вместо этого возникает следующая ошибка:
Error in qr.qty(qrx, f) :
right-hand side should have 14195 not 14196 rows
Эта ошибка исчезнет, если я переключусь на использование функции gam(), а не bam(). Он также уходит, если я отбрасываю случайный эффект от модели. Я действительно не уверен, что вызывает эту ошибку, и что с ней делать. К сожалению, создание воспроизводимого примера потребует передачи очень большого набора данных, так как непонятно, почему эта ошибка возникает в этом конкретном наборе данных по сравнению с 4 другими наборами данных, соответствующими той же модели.
Есть идеи, почему это возникает ошибка, и мы будем очень признательны за то, как ее преодолеть.