and are defined in the frequency domain: delta (2–3.5 Hz), theta (4–7.5 Hz), alpha1 (8–10 Hz), alpha2 (10–12 Hz), beta1 (13–18 Hz), beta2 (18.5–21 Hz), beta3 (21.5–30 Hz) and gamma (30.5–44 Hz). Based on this principle, lagged linear connectivity was calculated. Time-series of current density were extracted for different regions of interests using sLORETA. Power in all 6,239 voxels was normalized to a power of 1 and log transformed at each time point. Region of interest values thus reflect the log-transformed fraction of total power across all voxels, separately for specific frequencies. Regions of interest selected were based on the activity findings of both EEG and fMRI - left and right orbitofrontal cortex, the pregenual anterior cingulate cortex, the dorsal anterior cingulate cortex, the posterior cingulate cortex, the left and right inferior frontal gyrus, the left and right insula, the left parahippocampus.