“…Data pre-processing was done using the FMRIB Software Library (FSL, Oxford, UK; Smith et al, 2004; Woolrich et al, 2009), 4dfp, and an in-house pipeline built using Nipype (Gorgolewski et al, 2011), modified based on published functional MRI (fMRI) analysis methods (Fair et al, 2007; Fair et al, 2009; Fair et al, 2012; Iyer et al, 2013; Miranda-Dominguez et al, 2014a), and optimized for macaques (Miranda-Dominguez et al, 2014b), including resting state fMRI (rsfMRI) studies in infant rhesus by our group (Kovacs-Balint et al, 2018; Mavigner et al, 2018). This procedure consisted of (a) quantification and correction of dynamic field map changes, (b) slice-time correction of intensity differences as a result of interleaved slice image acquisition, (c) combined resampling of: within-run rigid-body motion correction, registration of the EPI to the subject's own averaged T1-weighted structural image, and registration of the T1 to age-specific T1-weighted rhesus infant and juvenile brain structural MRI atlases developed by our group (publicly available at: ), using nonlinear registration methods in FSL (FNIRT), (d) BOLD signal normalization to mode of 1000, to scale BOLD values across participants at an acceptable range, (e) BOLD signal detrending, (f) regression of rigid body head motion (six directions), global brain signal, BOLD signal of the ventricles and white matter (derived from manually drawn masks), and first-order derivatives of these signals, and (g) low-pass (f < 0.1 Hz) temporal filter (second order Butterworth filter) (Fair et al, 2007; Fair et al, 2009; Fair et al, 2012; Miranda-Dominguez et al, 2014b).…”