Statistical analysis
All data were analyzed using SPSS version 24.0 (IBM Corp., Armonk, NY, USA). The data are expressed as mean ± standard deviation (SD). To compare the general demographic data, one-way ANOVA was used for comparing means between the four groups, and nonparametric data were compared using the Chi-square test.Repeated-measure of covariance was used for most of the statistical analysis. The scores of THINC-it at W0, W2, and W4 were evaluated as the outcome variables. We conducted four types of comparisons to analyze the effects of the rTMS intervention mode (active/sham) and the stimulating mode (ACC/DLPFC) on cognitive function. The first comparison included all participants with the targeting of the ACC; the between-participant factor was set as the group (active/sham), and the within-participant factor was set as the time point (W0/W2) and (W2/W4), respectively, to examine the interaction between time and intervention mode. The second comparison included all participants with the DLPFC target; the between-participant and within-participant factors were the same as those used for the ACC target. The third comparison included all participants who received active-sham stimulation; the between-participant factor was the group (ACC/DLPFC), and the within-participant factor was the time point (W0/W2) and (W2/W4). For the last comparison, we included all participants who received sham-active stimulation, using the same between-participant and within-participant factors as the active-sham stimulation. Significance was assessed by the interaction between the time point and the mode of stimulation, indicating whether or not the change in outcome variables over time differed between the mode of stimulation.The mood-symptom data were analyzed using repeated-measures ANOVA, with the treatment group (A1, DLPFC active-sham rTMS Group; A2, DLPFC sham-active rTMS Group; B1, ACC active-sham rTMS Group; B2, ACC sham-active rTMS Group) as the between-participant factor and time point when cognitive and clinical measurements were conducted (W0, W2, W4) as the within-participant factor. The threshold for significance level was defined as The FC between the voxels in the whole brain and two ROIs, DLPFC and ACC, was determined separately by Pearson’s correlation coefficients in DPABI (