Hz with 1 Hz bandwidth. The time series of the signal i, i.e., xi(t), was set to the time series of the 5 Hz oscillator while the time series of the signal j was set to a time-delayed copy of the time series of signal i, i.e., xj(t) = xj(t − τ), with a time delay τ = 75 ms. Uncorrelated white Gaussian noise was finally added to both the generated signals. For each signal realization, the circular phase difference between signals i and j was calculated as the phase of their complex-valued coherency (see section Univariate methods for mathematical details). For scenario 2 (Figure 3B) the two signals i and j were independently generated as two 5-Hz oscillators. A narrow peak in the circular phase difference distribution can be observed in the case of phase-coupled signals (A), as opposed to the case of independent signals (B).