Phase frequency by amplitude frequency plots of PAC estimates, or comodulograms, were calculated using EEG recordings from midfrontal channel FZ for every 30 seconds of EEG data to capture the dynamics of the signal over time. PAC estimates were calculated for phase frequencies between 0.1-13 Hz, and amplitude frequencies between 4-50 Hz across a 224 x 224 matrix using 1 Hz window width at 0.1 Hz step to generate high resolution comodulograms. To remove bias in estimates of PAC strength, a null distribution was generated by swapping amplitude time blocks in each comodulogram (200 permutations). TensorPAC package was used to generate PAC comodulograms (Combrisson et al., 2020).