Customized square waveforms for each participant were generated for the specific counterbalanced order of experimental conditions and then convolved with a double-gamma hemodynamic response function. For each participant, hemodynamic parameters for the different explanatory variables (e.g., high arousal and low arousal) were estimated, and statistical maps of interest were generated for the contrasts across these explanatory variables. The contrast maps from individual subjects were warped into common stereotaxic space before mixed-effects group analyses were performed using the FMRIB’s Local Analysis of Mixed Effects module to obtain the group mean of brain activation.
