Abstract:In randomized controlled trials, forest plots are frequently used to investigate the homogeneity of treatment effect estimates in subgroups. However, the interpretation of subgroup-specific treatment effect estimates requires great care due to the smaller sample size of subgroups and the large number of investigated subgroups. Bayesian shrinkage methods have been proposed to address these issues, but they often focus on disjoint subgroups while subgroups displayed in forest plots are overlapping, i.e., each subject appears in multiple subgroups. In our approach, we first build a flexible Cox model based on all available observations, including categorical covariates that identify the subgroups of interest and their interactions with the treatment group variable. We explore both penalized partial likelihood estimation with a lasso or ridge penalty for treatment-by-covariate interaction terms, and Bayesian estimation with a regularized horseshoe prior. One advantage of the Bayesian approach is the ability to derive credible intervals for shrunken subgroup-specific estimates. In a second step, the Cox model is marginalized to obtain treatment effect estimates for all subgroups. We illustrate these methods using data from a randomized clinical trial in follicular lymphoma and evaluate their properties in a simulation study. In all simulation scenarios, the overall mean-squared error is substantially smaller for penalized and shrinkage estimators compared to the standard subgroup-specific treatment effect estimator but leads to some bias for heterogeneous subgroups. We recommend that subgroup-specific estimators, which are typically displayed in forest plots, are more routinely complemented by treatment effect estimators based on shrinkage methods. The proposed methods are implemented in the R package bonsaiforest.

Comparing Approaches to Treatment Effect Estimation for Subgroups in Clinical Trials

Subgroup identification in clinical trials: an overview of available methods and their implementations with R

Estimating the treatment effect in a subgroup defined by an early post-baseline biomarker measurement in randomized clinical trials with time-to-event endpoint

Using shrinkage methods to estimate treatment effects in overlapping subgroups in randomized clinical trials with a time-to-event endpoint

Detecting critical treatment effect bias in small subgroups

A Utility-Based Design for Randomized Comparative Trials with Ordinal Outcomes and Prognostic Subgroups

Tutorial in biostatistics: data‐driven subgroup identification and analysis in clinical trials

Identifying treatment response subgroups in observational time-to-event data

A Machine-Learning Approach for Estimating Subgroup- and Individual-Level Treatment Effects: An Illustration Using the 65 Trial

Testing for a treatment effect in a selected subgroup

Exploring Heterogeneity in Clinical Trials with Latent Class Analysis

Evaluating the impact of treating the optimal subgroup

Clinical trials with nested subgroups: Analysis, sample size determination and internal pilot studies

SUBGROUP IDENTIFICATION AND VARIABLE SELECTION FOR TREATMENT DECISION MAKING

Propensity Score Weighting for Causal Subgroup Analysis

Comparison of methods that combine multiple randomized trials to estimate heterogeneous treatment effects

A comparison of methods for model selection when estimating individual treatment effects

Online Testing of Subgroup Treatment Effects Based on Value Difference

Efficient Targeted Learning of Heterogeneous Treatment Effects for Multiple Subgroups

Selecting optimal subgroups for treatment using many covariates

Predicting subgroup treatment effects for a new study: Motivations, results and learnings from running a data challenge in a pharmaceutical corporation