Statistical Analysis
Statistical analyses were performed using SPSS Statistics for Windows version 24.0 (IBM). Analyses adhered to the intent-to-treat principle. Participant characteristics were summarized using descriptive statistics, including means, SDs, frequencies, and proportions. The Shapiro-Wilk test was used to determine whether continuous outcomes were normally distributed. Nonnormally distributed outcomes were log(x) or log(x + 1) transformed for analysis. Linear mixed models with random intercepts were generated for each outcome variable, with fixed effects of group (control and tele-CBT), time (baseline, after intervention, and follow-up), and group-by-time interaction. For reporting, estimated means and SEs from mixed models were back-transformed to their original units as: exp(mean log) ± exp(mean log) × (exp(SE log)-1) or exp(mean log) − 1 ± exp(mean log) × (exp(SE log) − 1) for log(x) or log(x + 1) transformations, respectively (based on the Delta method