Statistical analyses
Comorbidity
Descriptive statistics included mean and standard deviation for numeric variables and proportion for categorical variables. Mann–Whitney-U tests and Chi-square were used for unadjusted comparison of (continuous and categorical, respectively) baseline characteristics between control and intervention arms.This study used generalised linear mixed-effects models (GLMMs) to compare outcome variables between intervention and control groups at baseline and endpoint. Logistic GLMM was used for binary outcome variables (e.g., taking at least one PIM: yes/no), Poisson was used for discrete count variables (e.g., number of regular medications), and gamma distributions were used for continuous positively skewed variables (e.g., daily dose of psychotropic medicines). The models adjusted for potential confounders (age, gender, Charlson Comorbidity Index (CCI), presence of other educational intervention