Statistical analyses
The baseline characteristics of the participants and the PCCs were summarised using the means for continuous data and percentages for categorical data. We tested for baseline differences to ensure that cluster randomisation had led to comparable treatment groups. We used the R software version 4.1.1 [The intention-to-treat analyses included all participants, regardless of the number of sessions received. The use of mixed models with time, treatment, interaction between time and treatment, and baseline covariate as parameters obviates separate multiple imputation methods for missing data [