Gang,
Thanks for pointing out 3dLME. I think this is what I was looking for.
I've included my attempt at a model.txt file below. The goal was to create one that would test whether the estimated HRFs differ between two groups, with a covariate in the model as well. The lines I'm least confident in are the "Model:" line and the "CorStr:" line. Am I on the right track?
------------------------------
Data:Volume
Output:OutputFileName
MASK:../data/Mask+tlrc.BRIK
Model:-1+MyCov+Time:Group
COV:MyCov
RanEff:1
VarStr:0
CorStr:1~order|Subj
SS:sequential
Clusters:4
GroupOne-GroupTwo
0*0*0*One-0*0*0*Two
Subj MyCov Time order group InputFile
S1 23 Coeff0 0 One S1Coeff0+tlrc
S1 23 Coeff1 1 One S1Coeff1+tlrc
S1 23 Coeff2 2 One S1Coeff2+tlrc
S1 23 Coeff3 3 One S1Coeff3+tlrc
S1 23 Coeff4 4 One S1Coeff4+tlrc
S1 23 Coeff5 5 One S1Coeff5+tlrc
S1 23 Coeff6 3 One S1Coeff6+tlrc
S1 23 Coeff7 4 One S1Coeff7+tlrc
S1 23 Coeff8 5 One S1Coeff8+tlrc
... <-- Additional group One subjects entered here...
S12 58 Coeff0 0 Two S12Coeff0+tlrc
S12 58 Coeff1 1 Two S12Coeff1+tlrc
S12 58 Coeff2 2 Two S12Coeff2+tlrc
...
----------------------------
Thanks again for your assistance,
-John