Analysis
T1–T3
The dataset included all participants who provided data at T1. Baseline characteristics are summarized as means and standard deviations for continuous variables and frequencies and proportions for categorical variables. For the outcome measures and process evaluation measures, we first compared the means of the three time points. Because this program was implemented on a group basis, data is nested by group. Therefore, we fitted mixed models for repeated measures (MMRM), with time points (T1, T2, and T3) as fixed effects and program sites as random effects. If the outcome measures showed significant changes in MMRM over time, T1–T2 and T1–T3 multiple comparison tests (Bonferroni correction) and effect sizes (Cohen’s d) were performed to examine the differences between time points. A The interviews were recorded and converted into textual data. The textual data were summarized for each statement made by the participant. The participants were classified according to their study IDs. The number of participants who reported about what impressed them regarding the program was counted based on the summarized statements for each program’s structural elements and content themes. The data on changes after program participation were categorized based on the similarity of the summarized statements, and the corresponding number of participants was counted. A principal researcher established the program’s structural elements, content themes, and categories. Counts were conducted separately by two researchers (MK and MN) for purposes of reliability (rate of agreement of 94.4%), and areas where there were differences in classification were reviewed by a third researcher (RO) to decide on the category.