Statistical analysis
We used the PASS software (PASS 15. NCSS, USA) to calculate the total sample size with α = 0.05 (two-sided) and β = 0.1. Based on our pre-experimental results and previous literature [Statistical analysis adopted SPSS19.0 statistical software (IBM, Armonk, USA). The independent enumeration data were expressed as number of cases and percentage (%), and comparison between the two groups were analyzed with the Chi square or Chi square correction or Fisher’s exact test according to the condition of N and T. The independent approximately normal distribution measurement data (assessed with the Shapiro-Wilk test) was expressed as mean ± standard deviation and then analyzed with the independent-samples t test between the two groups. Extremely non-normally distributed data were expressed as median (25th to 75th percentile) and analyzed with the Mann-Whitney U test. The repeated normal distribution measurement data within the groups were analyzed by one-way repeated-measures ANOVA, with Bonferroni correction test was used to assess downstream time points (T