Functional data were preprocessed to correct for slice scan time differences (using sinc interpolation), 3D motion artifacts (Trileaner sinc interpolation), linear drifts, and low-frequency non-linear drifts (high pass filter less than 2 cycles/time course).