A general linear model using SPM’s canonical hemodynamic response function (HRF), modeled each condition (rest, 0-, 1-, 2-, and 3-back), and six motion regressors. Linear contrasts compared task load versus 0-back. A second-level one-sample t-test investigated task effect using 2-back versus 0-back contrasts. A second-level linear regression used individual resiliency scores as covariate and the same contrasts as dependent variable. To test our hypotheses, an a priori regions-of-interest mask (frontal-cingulate-parietal-basal ganglia; (Bogacz et al., 2010, O’Reilly and Frank, 2006)) was created using WFU Pickatlas (Maldjian et al., 2003). Regions of significant correlation within this mask were identified using a voxel-wise threshold of p<0.005 uncorrected, combined with cluster size threshold of 61 contiguous voxels. This combined threshold provides protection against type I error (Forman et al., 1995) and was estimated with Monte Carlo simulation using AlphaSim (Douglas Ward, http://afni.nimh.nih.gov/pub/dist/doc/program_help/AlphaSim.html) giving an overall corrected threshold of p<0.05. For identified clusters, activation data were extracted from individual contrast maps for correlation with behavioral measures
