time window was rather short, which might lead to inaccurate PAC estimation, the z-score of the MVL (zMVL) was calculated using the mean and standard deviation of 200 surrogates, created by recombining the instantaneous phase and randomized shuffled amplitudes. zMVL values not larger than 1.96 (equivalent to 95% confidence interval) were assigned a value 0. Therefore, for each 3-s epoch in a component, we had 55 zMVL values with 50 ms time resolution.
