Leveraging machine learning to estimate individualized treatment effects in cluster-randomized trials

TL;DR

This paper develops a unified mixed-effects machine learning framework to estimate individualized treatment effects in cluster-randomized trials, accounting for both individual and cluster-level covariates.

Contribution

It introduces a comprehensive approach combining various machine learning methods with mixed-effects models for causal inference in CRTs.

Findings

Methods evaluated across diverse simulations.

Application demonstrated in Ghana hypertension trial.

Provided practical guidance and reproducible code.

Abstract

Cluster-randomized trials (CRTs) are widely used to evaluate interventions delivered at the clinic, practice, or community level. Although standard analyses typically target average treatment effects, such summaries mask potentially meaningful variation in treatment response across individuals and clusters. This work addresses the estimation of conditional average treatment effects (CATEs) for continuous outcomes in two-arm parallel CRTs by defining causal estimands that incorporate both individual- and cluster-level baseline covariates while marginalizing over unobserved cluster heterogeneity. To estimate these quantities, we develop a unified framework based on mixed-effects machine learning, integrating and extending a range of existing approaches, including Bayesian additive regression trees with random effects, multilevel Bayesian causal forests, mixed-effects random forests, several mixed-effects gradient boosting procedures, and generalized additive mixed models, while incorporating cluster-specific random intercepts to account for within-cluster dependence. We evaluate these methods across diverse simulation scenarios and demonstrate their use in the Task Shifting and Blood Pressure Control in Ghana CRT, which investigates strategies for improving hypertension management. Drawing on these investigations, we provide practical guidance for applying mixed-effects machine learning to quantify treatment-effect heterogeneity in CRTs, together with reproducible code that enables investigators to implement all methods within a coherent workflow.