Data were analyzed using SPM8 (Wellcome Department of Imaging Neuroscience, University College, London) and PASW Statistics v. 20.0 (SPSS, Inc.). Functional volumes were corrected for differences in slice acquisition timing and rigid-body realigned to the initial volume of the first functional imaging scan to account for residual movement after prospective motion correction (Thesen et al., 2000). Each subject’s high-resolution anatomic image was co-registered to the reference functional volume and segmented into tissue classes. Nonlinear spatial transformation parameters from this segmentation were utilized to convert BOLD volumes to the Montreal Neurological Institute (MNI) space, with the resulting volumes re-sampled to isotropic (2 mm per side) voxels and smoothed by a 6 mm full-width at half-maximum isotropic Gaussian Kernel.
