Statistical analyses
RPE, affective valence, MP
Data are presented as the mean ± standard deviation (M ± SD). The normal distribution of each data set was evaluated by the Shapiro-Wilk normality test. Two-way repeated measures ANOVA (3 × 3 factorial design; 3 stimulation conditions and 3-time points) was used to analyze PP, MP, FI, EMG, and RPE at each time point. Two-way repeated measures ANOVA (3 × 2 factorial design; 3 stimulation conditions and 2-time points) was also used to analyze CWST and CRT at pre and post-time. Bonferroni post hoc test was used for the pairwise comparisons. In case of a violation in the assumption of sphericity, the Greenhouse-Geisser epsilon correction was applied. When the assumption of normality was not met (affective valence and arousal), the Friedman test was adopted and if significant results were obtained, Bonferroni correction was used for pairwise comparisons. Partial eta squared (η