A Causal Transformation Model for Time-to-Event Data Affected by Unobserved Confounding

Motivated by studies investigating causal effects in survival analysis, we propose a structural bivariate transformation model to quantify the impact of a binary treatment on a time-to-event outcome. The model equations are connected through a bivariate Gaussian distribution, with the dependence parameter capturing unobserved confounding, and are specified as functions of additive predictors to flexibly account for the impacts of observed confounders. Moreover, the baseline survival function is estimated using monotonic P-splines, the effects of binary or factor instruments can be regularized through a ridge penalty, and interactions between treatment and modifier variables can be incorporated to accommodate potential variations in treatment effects across subgroups. The proposal naturally provides an intuitive causal measure, the survival average treatment effect. Parameter estimation is achieved via an efficient and stable penalized maximum likelihood estimation approach and intervals constructed using related inferential results. We revisit a dataset from the Illinois Reemployment Bonus Experiment to estimate the causal effect of a cash bonus on unemployment duration, unveiling new insights. The modeling framework is incorporated into the R package GJRM, enabling researchers and practitioners to employ the proposed model and ensuring the reproducibility of results.