Dear AFNI folks,
I think I've got this figured out, but if not, please let me know.
An all-zero 3d+time dataset (empty+tlrc) was created which had the same grid spacing and number of TRs (229) as a run of data. It was generated by taking a preprocessed 3d+time dataset and using 3dcalc to multiply it by zero. Because the IRF contained 24 time-points, stim_maxlag was set using (n-1)=23.
3dConvolve \
-input empty+tlrc \
-nfirst 0 \
-nlast 228 \
-polort -1 \
-num_stimts 1 \
-stim_maxlag 1 23 \
-stim_file 1 new_impulse_function.1D \
-iresp 1 IRF_file+tlrc \
-output simulated_signal
Correct, no?
Cheers!
David