You can do it by extracting the cubic column(s) from an xmat.1D file created with 3dDeconvolve -polort 3 -x1D_stop and adding that as a -ortvec to a 3dDeconvolve -polort 1, though there's probably an easier way. I'm not convinced it's a good idea, though, since it effectively amounts to assuming that the inflection point of your scanner's drift occurred halfway through each run.
My intuition on this is probably not as valuable as that of other people on this board, though.
ijs
Edited 1 time(s). Last edit at 06/12/2015 05:22PM by Isaac Schwabacher.