To create functional connectivity (FC) matrices, time series of voxels within each of the 20 ROIs were averaged and correlated with the averaged time series of all other ROIs resulting in 190 correlation coefficients which were then transformed using Fisher's z calculation. Functional connectivity matrices were extracted using the conn toolbox.