我有一个实验的数据。参与者对一些项目做出了回应，我想根据一些措施来预测响应时间。试用t中的响应时间与 中的响应时间相关t-1。我的bam-model ( mgcv-package) 如下所示：
bam(RT ~ s(length) + s(color) + s(length, subject, bs='re'),
data=dat, discrete=TRUE, AR.start=dat$AR.start, rho=0.8)
在实验中，提出了目标和干扰物。我只想分析目标。在lmer中，我会删除带有干扰项的行dat <- dat[dat$type == 'target',]。在目前的情况下，恐怕这不合适，因为bam不知道连续的行不再代表连续的试验，对吧？那么，在估计固定和随机效应时不包括干扰项，而仅用于估计自相关的最佳做法是什么？