Pre-processing of the data was performed as previously described [22] and included three-dimensional motion correction, temporal Gaussian smoothing of 4 s, spatial Gaussian smoothing of 4 mm, linear detrending and temporal high pass filtering using 0.01 Hz.