gives the amplitude of the higher frequency band oscillation at each phase of the lower frequency band rhythm. Next, the phases of the lower band were binned into 18 frequency bins, and the mean of the amplitude of the higher frequency band over each phase bin was calculated. We then normalize the mean amplitudes by dividing each bin value by the sum of the mean amplitudes over all the bins. Lastly, we calculated the associated entropy measure for the normalized mean amplitudes and used them to compute the MI (see the Supplementary Material for the MATLAB code). First, we calculated the gamma/theta CFC with gamma between 30 and 40 Hz and theta between 4 and 7 Hz. We did not include 40–50 Hz gamma due to contamination with line noise. For the exploratory analysis, we performed CFC for gamma–beta, gamma–alpha, gamma–delta, beta–alpha, beta–theta, beta–delta, alpha–theta, alpha–delta, and theta–delta. The MATLAB script used for calculating CFC in the current study can be found in the Supplementary Material.
