Analysis of variance revealed no significant differences between the intervention group and the waiting control group in terms of demographic characteristics or pre-intervention SPIN, LASK, and CORE (see Table 1).
Table 1
Demographic/clinical characteristics of participants per condition at pre-intervention
Demographic/clinical variable | Intervention group (N = 25) | Control group (N = 29) |
Gender % Female | 22 (88.00%) | 25 (86.21%) |
Mean Age(SD) | 22.32 (3.01) | 23.03 (2.87) |
Education, n(%) | | |
College | 12 (48.00%) | 11 (37.93%) |
Undergraduate | 11 (44.00%) | 15 (51.72%) |
Postgraduate | 2 (8.00%) | 3 (10.34%) |
CORE-OM(SD) | 46.64 (14.39) | 47.41 (14.77) |
SPIN(SD) | 40.08 (10.62) | 41.48 (9.25) |
LSAS(SD) | 70.96 (19.54) | 75.52 (22.55) |
PANAS-P(SD) | 27.12 (4.48) | 26.52 (5.57) |
PANAS-N(SD) | 27.84 (6.67) | 28.93 (7.16) |
Note: CORE-OM = The Clinical Outcomes in Routine Evaluation - Outcome Measure; SPIN = The Social Phobia Inventory; LSAS = The Liebowitz Social Anxiety Scale; PANAS-P = The Positive and Negative Affect Schedule-Positive; PANAS-N = The Positive and Negative Affect Schedule-Negative. |
INSERT Table ABOUT HERE
Analyses of SPIN and LSAS scores revealed that there were statistically significant main effects for GROUP (SPIN: F(1, 52) = 6.528, p = 0.014, η2 = 0.112; LSAS: F(1, 52) = 4.133, p = 0.047, η2 = 0.074), TIME (SPIN: F(1, 52) = 32.970, p < 0.001, η2 = 0.388; LSAS: F(1, 52) = 21.26, p < 0.001, η2 = 0.29) and interaction between GROUP and TIME (SPIN: F(1, 52) = 18.182, p < 0.001, η2 = 0.259; LSAS: F(1, 52) = 9.678, p < 0.03, η2 = 0.157).
The simple effect analysis found that the intervention group had a significant decrease on the SPIN (t(1, 52) = 46.608, p < 0.001,η2 = 0.473, MSPINpre = 40.08, SD = 10.62; MSPINpost = 28.64, SD = 10.21)and LSAS scores (t(1, 52) = 27.757, p < 0.001, η2 = 0.348, MLSASpre = 70.96 SD = 19.54; MLSASpost= 55.16, SD = 18.39); however, the control group had no significant improvement on SPIN (t(52) = 1.179, p = 0.282, η2 = 0.022; MSPINpre = 41.48, SD = 9.25; MSPINpost = 39.79, SD = 9.74) and LSAS scores (t(52) = 1.215, p = 0.275, η2 = 0.023; MLSASpre = 75.52, SD = 22.55; MLSASpost = 72.45, SD = 22.87).
Analysis of CORE-OM scores revealed a significant TIME effect, F(1, 52) = 11.123, p = 0.002, η2 = 0.176. The pre-score of CORE-OM (Mcorepre = 46.64, SD = 14.39) was significantly higher than the post-score of CORE-OM (Mcorepost = 39.04, SD = 12.49); however, the main effect of GROUP and the interaction between GROUP and TIME did not reached significance (all ps < 0.05).
Analysis of PANAS scores revealed a significant TIME effect (positive emotions: F(1, 52) = 9.896, p = 0.003, η2 = 0.160; negative emotions: F(1, 52) = 27.877, p < 0.001, η2 = 0.349) and a non-significant GROUP effect (p > 0.05). In details, the pre-score of PANAS-P (MPANAS−Ppre = 27.12, SD = 4.48) was significantly lower than the post-score of PANAS-P (MPANAS−Ppost = 29.40, SD = 4.84). Meanwhile, the pre-score of PSNAS-N (MPANAS−Npre = 27.84, SD = 6.67) was significantly higher than the post-score of PSNAS-N (MPANAS−Npost = 21.12, SD = 5.78).
We also found a significant interaction effect of TIME and GROUP on PANAS-N, F(1, 52) = 9.545, p = 0.003, η2 = 0.155. The simple effect analysis found that the intervention group had a significant decrease on the PANAS-N scores (MPANAS−Npre= 27.84, SD = 6.67; MPANAS−Npost = 21.12, SD = 5.78); however, the control group had no significant improvement (MPANAS−Npre = 28.93, SD = 7.16; MPANAS−Npost = 27.17, SD = 5.28). The interaction effect of TIME and GROUP on PANAS-P had not reached significance (p = 0.337).
The analyses of heart rate variability and low blood pressure variability found no variance homogeneity.The analysis of high blood pressure variability revealed a significant effect of TIME, F(1, 52) = 4.926, p = 0.031, η2 = 0.087. In details, the variability of high blood pressure in the pre-test (MPre = 4.07%, SD = 7.06%) was significantly higher than the variability of high blood pressure in the post-test (MPre = 0.60%, SD = 6.55%). However, we found no significant effect of GROUP and the interaction of GROUP and TIME on high blood pressure variability (see Fig. 3). All the comparation data and ANOVA results were listed in supplementary materials (see Tables 2 and 3). Additionally, no participants reported physical discomfort experienced after viewing virtual reality imagery.
INSERT FIGURE 3 ABOUT HERE
The one-way ANOVA analysis (n = 21) showed a significant main effect of TIME on LSAS scores (F(1, 28) = 18.224, p < 0.001, η2 = 0.477) and SPIN scores (F(1, 28) = 26.519, p < 0.001, η2 = 0.570). The post-hoc analyses revealed that the post-intervention scores (MSPINpost = 29.25, SD = 9.26; MLSASpost = 58.48, SD = 15.28) were significantly decreased than pre-intervention scores (MSPINpre = 40.00, SD = 11.41; MLSASpre = 69.90, SD = 20.53), as well as the follow-up scores (MSPINfollow−up = 26.62, SD = 8.90; MLSASfollow−up = 50.19, SD = 17.66) were significantly decreased than pre-intervention scores (all ps < 0.05).
We also found a significant main effect of TIME on CORE-OM scores (F(1, 23) = 8.459, p < 0.01, η2 = 0.297), PANAS-P scores, and PANAS-N scores. The post-hoc analyses revealed that the post-intervention scores (MCORE−OMpost = 40.86, SD = 10.97; MPANAS−Ppost = 29.05, SD = 5.05; MPANAS−Npost = 21.86, SD = 5.43) are significantly improved than pre-intervention scores (MCORE−OMpre = 47.62, SD = 14.88; MPANAS−Ppre = 26.90, SD = 4.55; MPANAS−Npre = 28.05, SD = 6.88), as well as the follow-up scores (MCORE−OMfollow−up = 36.95, SD = 10.39; MPANAS−Pfollow−up = 29.24, SD = 6.00; MPANAS−Nfollow−up = 20.38, SD = 5.47) are significantly improved than pre-intervention scores (all ps < 0.05).
The main effects of TIME on high and low blood pressure variability reached significance (high blood pressure variability: F(1, 20) = 8.451, p < 0.01, η2 = 0.297; low blood pressure variability: F(1, 20) = 18.742, p < 0.001, η2 = 0.484). The post-hoc analyses revealed that the post-intervention scores (MHPVpost = 1.05%, SD = 6.59%; MLPVpost = -3.63%, SD = 7.31%) are significantly improved than pre-intervention scores (MHPVpre = 5.20%, SD = 6.90%; MLPVpre = 2.31%, SD = 6.49%) (see Fig. 4). All the comparation data and ANOVA results were listed in supplementary materials (see Tables 4 and 5).
INSERT FIGURE 4 ABOUT HERE