To generate measures of event-related functional connectivity (27,28), we modeled every cue and probe event with a unique delta function, convolved with a canonical hemodynamic response function, using SPM5. All events were modeled, but only cue events for correct trials were included in the analysis, because these represented trial periods in which cognitive control demands were maximized. Subjects who committed no errors had 160 trials (128 A cue; 32 B cue). Separate regressors modeling each event were defined in a general linear model to yield 128 unique A cue and 32 unique B cue beta values for every voxel. Each beta value reflected the magnitude of the hemodynamic response evoked by each event. The beta values were then sorted by condition (i.e., A cue or B cue) and concatenated to generate a beta series for each condition. We then parcellated each participant’s brain into discrete regions of interest to represent network nodes (19,38,39) using an anatomic atlas (40). Regions with less than 25% brain coverage were excluded, yielding 78 anatomic regions (Table S1 in Supplement 1). Pairwise Pearson correlations between
