Statistical Analysis
The change in outcomes between the T1 and T2 surveys was calculated and compared by group using a Data from T1 to T3 in the intervention group were analyzed using a single-group, pre–post design. Due to missing values (loss to follow-up) and intraindividual correlations in the outcomes, a mixed model for repeated-measures was used in this analysis. Each time point was entered as a dummy variable (T1 vs T2 and T1 vs T3), and its fixed effects were estimated using the restricted maximum likelihood method. Degrees of freedom were adjusted using the Satterthwaite method. All analyses were performed using Stata 16 software (StataCorp LP, College Station, TX).