Horseshoe Forests for High-Dimensional Causal Survival Analysis

Abstract

We develop a Bayesian tree ensemble model to estimate heterogeneous treatment effects in censored survival data with high-dimensional covariates. Instead of imposing sparsity through the tree structure, we place a horseshoe prior directly on the step heights to achieve adaptive global-local shrinkage. This strategy allows flexible regularisation and reduces noise. We develop a reversible jump Gibbs sampler to accommodate the non-conjugate horseshoe prior within the tree ensemble framework. We show through extensive simulations that the method accurately estimates treatment effects in high-dimensional covariate spaces, at various sparsity levels, and under non-linear treatment effect functions. We further illustrate the practical utility of the proposed approach by a re-analysis of pancreatic ductal adenocarcinoma (PDAC) survival data from The Cancer Genome Atlas.

References in corpus (9)

• Recursive Partitioning for Heterogeneous Causal Effects
• BART: Bayesian additive regression trees
• Sparsity information and regularization in the horseshoe and other shrinkage priors
• Robust Inference on Average Treatment Effects with Possibly More Covariates than Observations
• Discussion on "Bayesian Regression Tree Models for Causal Inference: Regularization, Confounding, and Heterogeneous Effects" by Hahn, Murray and Carvalho
• A simple sampler for the horseshoe estimator
• Contributed Discussion to Uncertainty Quantification for the Horseshoe by Stéphanie van der Pas, Botond Szabó and Aad van der Vaart
• Estimating heterogeneous survival treatment effect in observational data using machine learning
• Semiparametric Bayesian causal inference