Since this was a randomized, blind clinical trial, the statistical analyses for treatment effect evaluation of the primary and secondary outcomes are relatively straightforward. Baseline-to-end-point score changes in the three global SCL-90-R indices and rates of response in the GSI were computed as the primary endpoints. For defining rate of response, patients with a reduction of at least 30% from the baseline GSI score were classified as "much improved"; at least 50%, as "very much improved." Subanalyses of the baseline-toend-point score changes of the nine SCL factors and sleep quality score (the average of the scores of the three SCl-90-R items on sleep quality) were secondary endpoints. Statistical analysis on both primary and secondary outcomes was done using intention-to-treat analysis (ITT) with statistical software SPSS. Missing values in the SCL-90-R questionnaire for the patients who withdrew from the study before the eighth week were imputed using the last-observation-carried-forward method. For primary outcomes, effect sizes (for three global indices) and number needed to treat (NNT, for rate of response in the GSI), as well as the P values from two sample t-tests and chi-square tests, are reported in the treatment effect assessment. The same analytic approaches were applied to the secondary outcomes. Additionally, Fisher's exact test was used to compare the difference in dropout rate and AEs between the two treatment groups.