Statistical analysis
The data were analysed using IBM SPSS Statistics for Windows version 23.0. We followed the modified intention-to-treat analysis; only subjects who were randomized and who received all of the study interventions were included in the final analysis.Data are expressed as the means and standard deviations (SDs) for normally distributed continuous variables and medians and interquartile ranges (IQRs) for nonnormally distributed continuous variables. Categorical variables were reported as frequencies and proportions. Normality testing was performed using the Shapiro‒Wilk test. Student’s t test was used for comparison of normally distributed variables. The Mann‒Whitney U test was used for comparison of nonnormally distributed variables between the two groups. Categorical variables were analysed via the Pearson χ2 test or the Fisher’s exact test as appropriate. All tests were 2-tailed, and 