Using these four approaches (MVL-MI-Canolty; MVL-MI-Özkurt; KL-MI-Tort; PLV-MI-Cohen) we calculated PAC between phases 7–13 Hz (in 1 Hz steps) and amplitudes 34–100 Hz (in 2 Hz steps), for the time-period 300–1,500 ms following grating presentation and a 1,200 ms baseline period. PAC-values were calculated separately for each trial and then averaged to obtain a single MI-value per amplitude and phase. This was repeated using surrogate data, created by shuffling trial and phase-carrying information (200 surrogates), to normalize MI-values. On a PC with 32 GB of RAM, and Intel(R) Core™ i7-4790 processor, the computation time for these procedures was 4.5 h (script: 5_visual_PAC_four_methods.m).