Preprocessing of the imaging data was done in the same way as for the resting state scans. For the analysis we used SPM5 (http://www.fil.ion.ucl.ac.uk/spm/software/spm5/). Trial-related activity was assessed by convolving a vector of the onset times of the stimuli with a synthetic hemodynamic response function (HRF). The general linear model (GLM), as implemented in SPM5, was used to model the effects of interest as well as other confounding effects (scanner drift and motion). Statistical Parametrical Maps were identified for each participant by applying linear contrasts to the parameter estimates (beta weight) applying to the events of interest, resulting in a t-statistic for every voxel. Random effects analyses were employed to calculate group effects.