of the episode extraction is shown in Fig. 3A and B. Episodes with the same type extracted from trials in the same stimulation condition were concatenated for each subject. A 10 s sliding window with a 1 s step size was employed for PAC computation to improve data utilization. In the comparison of PAC between the three types of episodes, an inner-subject normalization was made by calculating the percentage relative change with respect to NN and scaling to the max value:(1)percentagerelativechangek=PACk−PACNNmax(abs(PACk,NN))×100%
