The first 10 volumes of each functional time series were discarded from the analysis to allow for magnetization equilibrium and for the acclimatization of the subjects to the scanning environment. The remaining 190 volumes were corrected for intravolume acquisition time differences between slices and were also corrected for intervolume geometrical displacement due to head movement. The functional images were then normalized to the standard space of the Montreal Neurological Institute (MNI) and resampled to a voxel size of 3 mm × 3 mm × 3 mm, and spatially smoothed with a 4-mm full-width at half-maximum (FWHM) Gaussian kernel. After that, the linear trend of the time courses was removed and temporal band pass filtering (0.01–0.08 Hz) was carried out to reduce the effects of low frequency drifts and physiological high-frequency noise. Participants included in this study were restricted to head motion of less than 2.0 mm in any direction and 2.0° of angular motion during the resting state scan. In addition, the volume-level mean framewise displacement (FD) (Van Dijk et al., 2012) was calculated and used as an extra global
