1
|
Rizopoulos D, Taylor JMG, Papageorgiou G, Morgan TM. Using joint models for longitudinal and time-to-event data to investigate the causal effect of salvage therapy after prostatectomy. Stat Methods Med Res 2024; 33:894-908. [PMID: 38502034 PMCID: PMC11041089 DOI: 10.1177/09622802241239003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/20/2024]
Abstract
Prostate cancer patients who undergo prostatectomy are closely monitored for recurrence and metastasis using routine prostate-specific antigen measurements. When prostate-specific antigen levels rise, salvage therapies are recommended in order to decrease the risk of metastasis. However, due to the side effects of these therapies and to avoid over-treatment, it is important to understand which patients and when to initiate these salvage therapies. In this work, we use the University of Michigan Prostatectomy Registry Data to tackle this question. Due to the observational nature of this data, we face the challenge that prostate-specific antigen is simultaneously a time-varying confounder and an intermediate variable for salvage therapy. We define different causal salvage therapy effects defined conditionally on different specifications of the longitudinal prostate-specific antigen history. We then illustrate how these effects can be estimated using the framework of joint models for longitudinal and time-to-event data. All proposed methodology is implemented in the freely-available R package JMbayes2.
Collapse
Affiliation(s)
- Dimitris Rizopoulos
- Department of Biostatistics, Erasmus University Medical Center, Rotterdam, the Netherlands
- Department of Epidemiology, Erasmus University Medical Center, the Netherlands
| | - Jeremy MG Taylor
- Department of Biostatistics, University of Michigan, Ann Arbor, MI, USA
| | - Grigorios Papageorgiou
- Department of Biostatistics, Erasmus University Medical Center, Rotterdam, the Netherlands
- Department of Epidemiology, Erasmus University Medical Center, the Netherlands
| | - Todd M Morgan
- Department of Urology, University of Michigan, Ann Arbor, MI, USA
| |
Collapse
|
2
|
Shin H, Antonelli J. Improved inference for doubly robust estimators of heterogeneous treatment effects. Biometrics 2023; 79:3140-3152. [PMID: 36745745 DOI: 10.1111/biom.13837] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2021] [Revised: 01/16/2023] [Accepted: 01/30/2023] [Indexed: 02/08/2023]
Abstract
We propose a doubly robust approach to characterizing treatment effect heterogeneity in observational studies. We develop a frequentist inferential procedure that utilizes posterior distributions for both the propensity score and outcome regression models to provide valid inference on the conditional average treatment effect even when high-dimensional or nonparametric models are used. We show that our approach leads to conservative inference in finite samples or under model misspecification and provides a consistent variance estimator when both models are correctly specified. In simulations, we illustrate the utility of these results in difficult settings such as high-dimensional covariate spaces or highly flexible models for the propensity score and outcome regression. Lastly, we analyze environmental exposure data from NHANES to identify how the effects of these exposures vary by subject-level characteristics.
Collapse
Affiliation(s)
- Heejun Shin
- Department of Statistics, University of Florida, Gainesville, Florida, USA
| | - Joseph Antonelli
- Department of Statistics, University of Florida, Gainesville, Florida, USA
| |
Collapse
|
3
|
Li F, Ding P, Mealli F. Bayesian causal inference: a critical review. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2023; 381:20220153. [PMID: 36970828 DOI: 10.1098/rsta.2022.0153] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/28/2022] [Accepted: 10/23/2022] [Indexed: 06/18/2023]
Abstract
This paper provides a critical review of the Bayesian perspective of causal inference based on the potential outcomes framework. We review the causal estimands, assignment mechanism, the general structure of Bayesian inference of causal effects and sensitivity analysis. We highlight issues that are unique to Bayesian causal inference, including the role of the propensity score, the definition of identifiability, the choice of priors in both low- and high-dimensional regimes. We point out the central role of covariate overlap and more generally the design stage in Bayesian causal inference. We extend the discussion to two complex assignment mechanisms: instrumental variable and time-varying treatments. We identify the strengths and weaknesses of the Bayesian approach to causal inference. Throughout, we illustrate the key concepts via examples. This article is part of the theme issue 'Bayesian inference: challenges, perspectives, and prospects'.
Collapse
Affiliation(s)
- Fan Li
- Duke University, Durham, NC, USA
| | - Peng Ding
- University of California, Berkeley, CA, USA
| | | |
Collapse
|
4
|
Josey KP, deSouza P, Wu X, Braun D, Nethery R. Estimating a Causal Exposure Response Function with a Continuous Error-Prone Exposure: A Study of Fine Particulate Matter and All-Cause Mortality. JOURNAL OF AGRICULTURAL, BIOLOGICAL, AND ENVIRONMENTAL STATISTICS 2023; 28:20-41. [PMID: 37063643 PMCID: PMC10103900 DOI: 10.1007/s13253-022-00508-z] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/27/2021] [Revised: 07/08/2022] [Accepted: 07/23/2022] [Indexed: 10/14/2022]
Abstract
Numerous studies have examined the associations between long-term exposure to fine particulate matter (PM2.5) and adverse health outcomes. Recently, many of these studies have begun to employ high-resolution predicted PM2.5 concentrations, which are subject to measurement error. Previous approaches for exposure measurement error correction have either been applied in non-causal settings or have only considered a categorical exposure. Moreover, most procedures have failed to account for uncertainty induced by error correction when fitting an exposure-response function (ERF). To remedy these deficiencies, we develop a multiple imputation framework that combines regression calibration and Bayesian techniques to estimate a causal ERF. We demonstrate how the output of the measurement error correction steps can be seamlessly integrated into a Bayesian additive regression trees (BART) estimator of the causal ERF. We also demonstrate how locally-weighted smoothing of the posterior samples from BART can be used to create a more accurate ERF estimate. Our proposed approach also properly propagates the exposure measurement error uncertainty to yield accurate standard error estimates. We assess the robustness of our proposed approach in an extensive simulation study. We then apply our methodology to estimate the effects of PM2.5 on all-cause mortality among Medicare enrollees in New England from 2000-2012.
Collapse
Affiliation(s)
- Kevin P. Josey
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA
| | - Priyanka deSouza
- Department of Urban and Regional Planning, University of Colorado, Denver, CO
| | - Xiao Wu
- Department of Statistics, Stanford University, Stanford, CA
- Stanford Data Science, Stanford University, Stanford, CA
| | - Danielle Braun
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA
- Department of Data Science, Dana-Farber Cancer Institute, Boston, MA
| | - Rachel Nethery
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA
| |
Collapse
|
5
|
Zhou J, Zhang Y, Tu W. A reference-free R-learner for treatment recommendation. Stat Methods Med Res 2023; 32:404-424. [PMID: 36540907 DOI: 10.1177/09622802221144326] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Assigning optimal treatments to individual patients based on their characteristics is the ultimate goal of precision medicine. Deriving evidence-based recommendations from observational data while considering the causal treatment effects and patient heterogeneity is a challenging task, especially in situations of multiple treatment options. Herein, we propose a reference-free R-learner based on a simplex algorithm for treatment recommendation. We showed through extensive simulation that the proposed method produced accurate recommendations that corresponded to optimal treatment outcomes, regardless of the reference group. We used the method to analyze data from the Systolic Blood Pressure Intervention Trial (SPRINT) and achieved recommendations consistent with the current clinical guidelines.
Collapse
Affiliation(s)
- Junyi Zhou
- Design and Inovation, 7129Amgen Inc., Thousand Oaks, CA, USA
| | - Ying Zhang
- Department of Biostatistics, 12284University of Nebraska Medical Center, Omaha, NE, USA
| | - Wanzhu Tu
- Department of Biostatistics and Health Data Science, Indiana University-School of Medicine and Fairbanks School of Public Health, Indianapolis, IN, USA
| |
Collapse
|
6
|
Gao Q, Zhang Y, Sun H, Wang T. Evaluation of propensity score methods for causal inference with high-dimensional covariates. Brief Bioinform 2022; 23:6603435. [PMID: 35667004 DOI: 10.1093/bib/bbac227] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/26/2021] [Revised: 05/11/2022] [Accepted: 05/17/2022] [Indexed: 11/12/2022] Open
Abstract
In recent work, researchers have paid considerable attention to the estimation of causal effects in observational studies with a large number of covariates, which makes the unconfoundedness assumption plausible. In this paper, we review propensity score (PS) methods developed in high-dimensional settings and broadly group them into model-based methods that extend models for prediction to causal inference and balance-based methods that combine covariate balancing constraints. We conducted systematic simulation experiments to evaluate these two types of methods, and studied whether the use of balancing constraints further improved estimation performance. Our comparison methods were post-double-selection (PDS), double-index PS (DiPS), outcome-adaptive LASSO (OAL), group LASSO and doubly robust estimation (GLiDeR), high-dimensional covariate balancing PS (hdCBPS), regularized calibrated estimators (RCAL) and approximate residual balancing method (balanceHD). For the four model-based methods, simulation studies showed that GLiDeR was the most stable approach, with high estimation accuracy and precision, followed by PDS, OAL and DiPS. For balance-based methods, hdCBPS performed similarly to GLiDeR in terms of accuracy, and outperformed balanceHD and RCAL. These findings imply that PS methods do not benefit appreciably from covariate balancing constraints in high-dimensional settings. In conclusion, we recommend the preferential use of GLiDeR and hdCBPS approaches for estimating causal effects in high-dimensional settings; however, further studies on the construction of valid confidence intervals are required.
Collapse
Affiliation(s)
- Qian Gao
- Department of Health Statistics, School of Public Health, Shanxi Medical University, Taiyuan, China
| | - Yu Zhang
- Department of Health Statistics, School of Public Health, Shanxi Medical University, Taiyuan, China
| | - Hongwei Sun
- Department of Health Statistics, School of Public Health and Management, Binzhou Medical University, Yantai, China
| | - Tong Wang
- Department of Health Statistics, School of Public Health, Shanxi Medical University, Taiyuan, China
| |
Collapse
|
7
|
Linero AR. Simulation-based estimators of analytically intractable causal effects. Biometrics 2021; 78:1001-1017. [PMID: 34051105 DOI: 10.1111/biom.13499] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2020] [Revised: 03/31/2021] [Accepted: 04/28/2021] [Indexed: 11/28/2022]
Abstract
In causal inference problems, one is often tasked with estimating causal effects which are analytically intractable functionals of the data-generating mechanism. Relevant settings include estimating intention-to-treat effects in longitudinal problems with missing data or computing direct and indirect effects in mediation analysis. One approach to computing these effects is to use the g-formula implemented via Monte Carlo integration; when simulation-based methods such as the nonparametric bootstrap or Markov chain Monte Carlo are used for inference, Monte Carlo integration must be nested within an already computationally intensive algorithm. We develop a widely-applicable approach to accelerating this Monte Carlo integration step which greatly reduces the computational burden of existing g-computation algorithms. We refer to our method as accelerated g-computation (AGC). The algorithms we present are similar in spirit to multiple imputation, but require removing within-imputation variance from the standard error rather than adding it. We illustrate the use of AGC on a mediation analysis problem using a beta regression model and in a longitudinal clinical trial subject to nonignorable missingness using a Bayesian additive regression trees model.
Collapse
Affiliation(s)
- Antonio R Linero
- Department of Statistics and Data Science, University of Texas at Austin, Austin, Texas, USA
| |
Collapse
|