Statistical analysis
All statistical analyses were performed using SPSS 21.0 (SPSS Inc., USA). The method of exploratory analysis (normality test and distribution curve) was used to test the distribution state of the data. The chi-square test (categorical variables), independent-samples t-test (normally distributed variables), and Wilcoxon Rank Sum Test (not normally distributed variables) were used to examine the difference between groups on baseline demographic and outcome variables. A repeated-measures ANOVA was performed for 3-time points to evaluate the effect of group factors, time factors, and the interaction between time × group on the outcomes variables. A simple effect test was used to examine the difference among groups within each time point and between the 3-time points within each group when the interactions of time×group were significant. Effect size, one of the critical indicators to explain the potential clinical significance, was calculated using Cohen’s d [