Statistical analysis
The normality of all demographic, clinical and fMRI variables was tested by using the Kolmogorov–Smirnov method and histogram inspection. For baseline data, nonparametric tests for data with a skewed distribution, Pearson’s chi-squared tests or Fisher’s exact tests for categorical data, and independent sample t-tests for data with a normal distribution were applied to assess group differences. The presented analyses were performed in the intention-to-treat (ITT) population, and missing data were mitigated by using the last observation carried forward method.Functional MRI time series were modeled using a general linear model with SPM12 software (The effect of the home-based TR approach on motor function was evaluated by using linear mixed-effects models of restricted maximum likelihood estimation within individuals. Time was analyzed as a repeated variable. Group effect, time effect and group-by-time interaction effect were inputted as fixed effects, and age, sex, educational level of patients, the number of treatment sessions, OT/PT and ETNS durations were entered as random effects. The changes in motor function scores from baseline to the end of the 12-week rehabilitative intervention and from baseline to the end of the follow-up period were analyzed. Partial correlation analyses were conducted in each group between the mean dALFF/dReHo/dDC value of each cluster, showing significant differences in the comparisons of the main effect of group or the interaction effect of group × time and motor function scores. Age, gender, educational level, the number of treatment sessions, OT/PT and ETNS durations were also considered as covariates, and the significance levels for correlation analyses were set at uncorrected 