Hello AFNI experts,

I have been doing research on using 3dDeconvolve vs. 3dREMLfit, and it seems like 3dREMLfit is the ideal option because fMRI time series are temporally autocorrelated, and 3dREMLfit will attempt to account for this autocorrelation. It also sounds like 3dREMLfit will soon become the new default in first level analysis. Therefore, I would like to use 3dREMLfit stats when I run 3dMVM and 3dttest++, but I am running into the issue of there not being an -iresp command for 3dREMLfit. I use the -iresp command in 3dDeconvolve to model my HDR, which I use to compare brain activation in certain ROIs across two groups. I'm wondering if there is a similar function in 3dREMLfit that will allow me to model the HDR as I do with the -iresp command in 3dDeconvolve?