Analysis
Quantitative data were analyzed with SPSS-Version 22.0. Data met the assumptions for planned statistical tests, including normality and absence of multicollinearity. We
used descriptive statistics to identify sample characteristics and the intervention implementation profile (e.g. hours,duration). Repeated measures analysis of variance
(RM-ANOVA) was used to examine: a) multivariate change in all outcomes from baseline to 6 months to 12 months, and b) univariate change in each outcome from baseline to 6 months and 12 months. When the assumption of sphericity was violated for an outcome, the Greenhouse-Geisser Epsilon F statistic was interpreted. Post-hoc pairwise
comparisons were conducted using the Bonferroni test.