如何解决R mgcv:bam自相关缺失值
我有一个实验的数据。参与者对一些项目做出了回应,我想根据一些措施来预测响应时间。试用 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
不会知道连续行不再代表连续试验,对吗?那么在估计固定和随机效应时不包括干扰因素,而只包括自相关估计的最佳做法是什么?
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。