Hi Josh,
Note that if you are using GAM, BLOCK, SPMG*, or similar, then you are actually telling 3dREMLfit what the HRF is supposed to look like. The cases where the program is figuring that out are with multiple basis functions, an particularly in the case of TENT.
And in order for the program to get a useful estimate of what the HRF looks like, it is very important to model drift and motion effects, or they will corrupt the very estimates that you seek. It does not seem reasonable to remove those model parameters.
Maybe there is an aspect I am missing. Please feel free to provide more details and possibly even a model (e.g. 3dDeconvolve command).
- rick