R mgcv: автокорреляция bam с пропущенными значениями
У меня есть данные эксперимента. Участники ответили на некоторые вопросы, и я хочу спрогнозировать время ответа на основе некоторых показателей. Время отклика в испытании
t коррелирует со временем отклика в
t-1. Моя -модель (
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не знал бы, что последовательные строки больше не представляют собой последовательные испытания, верно? Итак, как лучше всего не включать отвлекающие факторы для оценки фиксированных и случайных эффектов, а только для оценки автокорреляции?