Abstract:We recently developed a new method random-intercept accelerated failure time model with Bayesian additive regression trees (riAFT-BART) to draw causal inferences about population treatment effect on patient survival from clustered and censored survival data while accounting for the multilevel data structure. The practical utility of this method goes beyond the estimation of population average treatment effect. In this work, we exposit how riAFT-BART can be used to solve two important statistical questions with clustered survival data: estimating the treatment effect heterogeneity and variable selection. Leveraging the likelihood-based machine learning, we describe a way in which we can draw posterior samples of the individual survival treatment effect from riAFT-BART model runs, and use the drawn posterior samples to perform an exploratory treatment effect heterogeneity analysis to identify subpopulations who may experience differential treatment effects than population average effects. There is sparse literature on methods for variable selection among clustered and censored survival data, particularly ones using flexible modeling techniques. We propose a permutation-based approach using the predictor's variable inclusion proportion supplied by the riAFT-BART model for variable selection. To address the missing data issue frequently encountered in health databases, we propose a strategy to combine bootstrap imputation and riAFT-BART for variable selection among incomplete clustered survival data. We conduct an expansive simulation study to examine the practical operating characteristics of our proposed methods, and provide empirical evidence that our proposed methods perform better than several existing methods across a wide range of data scenarios. Finally, we demonstrate the methods via a case study of predictors for in-hospital mortality among severe COVID-19 patients and estimating the heterogeneous treatment effects of three COVID-specific medications. The methods developed in this work are readily available in the R ${\textsf {R}}$ package riAFTBART $\textsf {riAFTBART}$ .

Incorporating External Data into the Analysis of Clinical Trials via Bayesian Additive Regression Trees

Estimation of Causal Effects of Multiple Treatments in Observational Studies with a Binary Outcome

Nonparametric competing risks analysis using Bayesian Additive Regression Trees (BART)

Bayesian Additive Regression Trees using Bayesian model averaging

Co-data Learning for Bayesian Additive Regression Trees

Seemingly unrelated Bayesian additive regression trees for cost-effectiveness analyses in healthcare

PAM-HC: A Bayesian Nonparametric Construction of Hybrid Control for Randomized Clinical Trials Using External Data

BASIC: A Bayesian Adaptive Synthetic-Control Design for Phase II Clinical Trials.

Adaptively leverage multiple real-world data sources for treatment effect estimation based on similarity

Incorporating External Data for Analyzing Randomized Clinical Trials: A Transfer Learning Approach

Dynamic Treatment Regimes Using Bayesian Additive Regression Trees for Censored Outcomes

Bayesian Nonparametric Modeling for Causal Inference

A new method for clustered survival data: Estimation of treatment effect heterogeneity and variable selection

Digital twins and Bayesian dynamic borrowing: Two recent approaches for incorporating historical control data

Estimation of causal effects of multiple treatments in healthcare database studies with rare outcomes

Leveraging machine learning: Covariate-adjusted Bayesian adaptive randomization and subgroup discovery in multi-arm survival trials

BEATS: Bayesian hybrid design with flexible sample size adaptation for time-to-event endpoints

A propensity score‐integrated approach for leveraging external data in a randomized controlled trial with time‐to‐event endpoints

A Bayesian machine learning approach for estimating heterogeneous survivor causal effects: Applications to a critical care trial

Bayesian Nonparametric Common Atoms Regression for Generating Synthetic Controls in Clinical Trials

Identification of subpopulations with distinct treatment benefit rate using the Bayesian tree.