1
|
Liu J, Zhang J, Mitchell A, Fang M, Tian L. Causal inference for longitudinal data based on historical controls. J Biopharm Stat 2022; 33:289-306. [PMID: 36469552 DOI: 10.1080/10543406.2022.2148164] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Use of historical data has become a hot topic recently, considered to provide a way to reduce patient burden, lower drug development cost, and make innovative therapies available to patients earlier. In a single-arm study designed to examine the benefit of an experimental treatment, there is often a desire to compare the outcomes of patients receiving the new intervention with those receiving a control treatment, which can be extracted from sources such as historical trials or electronic medical records. Since the treatment is not randomly assigned, there is a need to adjust for the potential imbalance in key patient characteristics between the current study and historical controls. If the outcome of interest is measured longitudinally and subject to random missing, the required adjustment becomes more complicated. In this paper, we propose a doubly robust adjustment procedure specifically designed for longitudinal data analysis with missing data. The proposed method yields valid analysis results, if either the propensity score model or the mixed effects model for repeated measures (MMRM) regression model is correctly specified. An extensive numerical study is conducted to examine the performance of the proposed method. Data from a real clinical trial comparing with historical data are analyzed as an example applying the proposed procedure.
Collapse
Affiliation(s)
- Jeen Liu
- Data and Statistical Sciences, AbbVie Inc, North Chicago, IL, USA
| | - Jane Zhang
- Data and Statistical Sciences, AbbVie Inc, North Chicago, IL, USA
| | - Alan Mitchell
- Data and Statistical Sciences, AbbVie Inc, North Chicago, IL, USA
| | - Mindy Fang
- Graduate School of Health Innovation, Kanagawa University of Human Services, Yokosuka, Japan
| | - Lu Tian
- Department of Biomedical Data Science, Stanford University, Palo Alto, CA, USA
| |
Collapse
|
2
|
Liu M, Zhao Y. Weighted generalized estimating equations and unified estimation for longitudinal data with nonmonotone missing data patterns. Stat Med 2021; 41:1148-1156. [PMID: 34729797 DOI: 10.1002/sim.9246] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2021] [Revised: 10/07/2021] [Accepted: 10/18/2021] [Indexed: 11/10/2022]
Abstract
Missing data are a major complication in longitudinal data analysis. Weighted generalized estimating equations (WGEEs, Robins et al, J Am Stat Assoc 1995;90:106-121) were developed to deal with missing response data. They have been extended for data with both missing responses and missing covariates (Chen et al, J Am Stat Assoc 2010;105:336-353). However, it may introduce more variability in dealing with the correlation structure of the responses. We propose new WGEEs for missing at random data where both response and (time-dependent) covariates may have values missing in nonmonotone missing data patterns. We also explain how to improve the estimation efficiency of WGEEs using a unified approach (Zhao and Liu, AStA Adv Stat Anal 2021;105(1):87-101). The proposed unified estimator is consistent and more efficient than the regular WGEE estimator. It is computationally simple and can be directly implemented in standard software. Simulation studies for both continuous response and binary response data are provided to examine the performance of the proposed estimators. A clinical trial example investigating the quality of life of women with early-stage breast cancer and the associated factors is analyzed.
Collapse
Affiliation(s)
- Meng Liu
- Finance Business Division, Xiamen ITG Group Corp., Ltd, Fujian, China
| | - Yang Zhao
- Department of Mathematics and Statistics, University of Regina, Regina, Saskatchewan, Canada
| |
Collapse
|
3
|
Josefsson M, Daniels MJ. Bayesian semi-parametric G-computation for causal inference in a cohort study with MNAR dropout and death. J R Stat Soc Ser C Appl Stat 2021; 70:398-414. [PMID: 33692597 PMCID: PMC7939177 DOI: 10.1111/rssc.12464] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
Causal inference with observational longitudinal data and time-varying exposures is often complicated by time-dependent confounding and attrition. The G-computation formula is one approach for estimating a causal effect in this setting. The parametric modeling approach typically used in practice relies on strong modeling assumptions for valid inference, and moreover depends on an assumption of missing at random, which is not appropriate when the missingness is missing not at random (MNAR) or due to death. In this work we develop a flexible Bayesian semi-parametric G-computation approach for assessing the causal effect on the subpopulation that would survive irrespective of exposure, in a setting with MNAR dropout. The approach is to specify models for the observed data using Bayesian additive regression trees, and then use assumptions with embedded sensitivity parameters to identify and estimate the causal effect. The proposed approach is motivated by a longitudinal cohort study on cognition, health, and aging, and we apply our approach to study the effect of becoming a widow on memory. We also compare our approach to several standard methods.
Collapse
Affiliation(s)
- Maria Josefsson
- Centre for Demographic and Ageing Research, Umeå University, Sweden
| | | |
Collapse
|
4
|
Chen C, Shen B, Liu A, Wu R, Wang M. A multiple robust propensity score method for longitudinal analysis with intermittent missing data. Biometrics 2020; 77:519-532. [PMID: 32662124 DOI: 10.1111/biom.13330] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2019] [Revised: 04/15/2020] [Accepted: 06/16/2020] [Indexed: 01/05/2023]
Abstract
Longitudinal data are very popular in practice, but they are often missing in either outcomes or time-dependent risk factors, making them highly unbalanced and complex. Missing data may contain various missing patterns or mechanisms, and how to properly handle it for unbiased and valid inference still presents a significant challenge. Here, we propose a novel semiparametric framework for analyzing longitudinal data with both missing responses and covariates that are missing at random and intermittent, a general and widely encountered situation in observational studies. Within this framework, we consider multiple robust estimation procedures based on innovative calibrated propensity scores, which offers additional relaxation of the misspecification of missing data mechanisms and shows more satisfactory numerical performance. Also, the corresponding robust information criterion on consistent variable selection for our proposed model is developed based on empirical likelihood-based methods. These advocated methods are evaluated in both theory and extensive simulation studies in a variety of situations, showing competing properties and advantages compared to the existing approaches. We illustrate the utility of our approach by analyzing the data from the HIV Epidemiology Research Study.
Collapse
Affiliation(s)
- Chixiang Chen
- Division of Biostatistics and Bioinformatics, Department of Public Health Sciences, Penn State College of Medicine, Hershey, Pennsylvania
| | - Biyi Shen
- Division of Biostatistics and Bioinformatics, Department of Public Health Sciences, Penn State College of Medicine, Hershey, Pennsylvania
| | - Aiyi Liu
- Biostatistics and Bioinformatics Branch, National Institute of Child Health and Human Development, NIH, Bethesda, Maryland
| | - Rongling Wu
- Division of Biostatistics and Bioinformatics, Department of Public Health Sciences, Penn State College of Medicine, Hershey, Pennsylvania
| | - Ming Wang
- Division of Biostatistics and Bioinformatics, Department of Public Health Sciences, Penn State College of Medicine, Hershey, Pennsylvania
| |
Collapse
|
5
|
Cai Q, Wang S. Inferences with generalized partially linear single-index models for longitudinal data. J Stat Plan Inference 2019. [DOI: 10.1016/j.jspi.2018.09.011] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
|
6
|
Xu C, Li Z, Xue Y, Zhang L, Wang M. An R package for model fitting, model selection and the simulation for longitudinal data with dropout missingness. COMMUN STAT-SIMUL C 2018; 48:2812-2829. [PMID: 32346220 PMCID: PMC7188076 DOI: 10.1080/03610918.2018.1468457] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2017] [Revised: 03/08/2018] [Accepted: 04/15/2018] [Indexed: 01/10/2023]
Abstract
Missing data arise frequently in clinical and epidemiological fields, in particular in longitudinal studies. This paper describes the core features of an R package wgeesel, which implements marginal model fitting (i.e., weighted generalized estimating equations, WGEE; doubly robust GEE) for longitudinal data with dropouts under the assumption of missing at random. More importantly, this package comprehensively provide existing information criteria for WGEE model selection on marginal mean or correlation structures. Also, it can serve as a valuable tool for simulating longitudinal data with missing outcomes. Lastly, a real data example and simulations are presented to illustrate and validate our package.
Collapse
Affiliation(s)
- Cong Xu
- Vertex Pharmaceuticals, Boston, Massachusetts, USA
| | - Zheng Li
- Department of Public Health Sciences, Division of Biostatistics and Bioinformatics, College of Medicine, Penn State Hershey Medical Center, Hershey, Pennsylvania, USA
| | - Yuan Xue
- School of Statistics, University of International Business and Economics, Beijing, China
| | - Lijun Zhang
- Department of Biochemistry and Molecular Biology, Institute of Personalized Medicine, Penn State Hershey Medical Center, Hershey, Pennsylvania, USA
| | - Ming Wang
- Department of Public Health Sciences, Division of Biostatistics and Bioinformatics, College of Medicine, Penn State Hershey Medical Center, Hershey, Pennsylvania, USA
| |
Collapse
|
7
|
da Silva JLP, Colosimo EA, Demarqui FN. A general GEE framework for the analysis of longitudinal ordinal missing data and related issues. STAT MODEL 2018. [DOI: 10.1177/1471082x17752753] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Generalized estimating equations (GEEs) are a well-known method for the analysis of categorical longitudinal data. This method presents computational simplicity and provides consistent parameter estimates that have a population-averaged interpretation. However, with missing data, the resulting parameter estimates are consistent only under the strong assumption of missing completely at random (MCAR). Some corrections can be done when the missing data mechanism is missing at random (MAR): inverse probability weighting GEE (WGEE) and multiple imputation GEE (MIGEE). A recent method combining ideas of these two approaches has a doubly robust property in the sense that one only needs to correctly specify the weight or the imputation model in order to obtain consistent estimates for the parameters. In this work, a proportional odds model is assumed and a doubly robust estimator is proposed for the analysis of ordinal longitudinal data with intermittently missing responses and covariates under the MAR mechanism. In addition, the association structure is modelled by means of either the correlation coefficient or local odds ratio. The performance of the proposed method is compared to both WGEE and MIGEE through a simulation study. The method is applied to a dataset related to rheumatic mitral stenosis.
Collapse
Affiliation(s)
- José LP da Silva
- Departamento de Estatística, Universidade Federal de Minas Gerais, Brazil
| | - Enrico A Colosimo
- Departamento de Estatística, Universidade Federal de Minas Gerais, Brazil
| | - Fábio N Demarqui
- Departamento de Estatística, Universidade Federal de Minas Gerais, Brazil
| |
Collapse
|
8
|
Estimation of response from longitudinal binary data with nonignorable missing values in migraine trials. Contemp Clin Trials Commun 2016; 4:90-98. [PMID: 29736472 PMCID: PMC5935886 DOI: 10.1016/j.conctc.2016.06.011] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2015] [Revised: 06/10/2016] [Accepted: 06/17/2016] [Indexed: 11/22/2022] Open
Abstract
In migraine trials pain relief responses from a headache at specific time points and sustained pain relief response over a period of time are important efficacy measures. When there are missing records of individual time point pain scores and/or headache recurrences during a migraine trial, the common approach used in practice to estimate the sustained response is statistically inconsistent even if the data are missing completely at random. Methods dealing with nonignorable longitudinal missing data usually assume certain models for the missing mechanism which can not be checked as they involve unobserved data. Taking advantage of the specific definition of the 'sustained pain relief' response, we propose two estimating methods based on intuitive imputation, which do not require model assumptions on the missing probability or specification of the correlation structure among the longitudinal observations. The consistency of the proposed methods is discussed in theory and their empirical performances are assessed through intensive simulation studies. The simulation results show that the proposed methods perform well in terms of reducing bias and mean square error except in several extreme cases which are unlikely to happen in real trials. The application of the proposed methods is illustrated in a real data analysis.
Collapse
|
9
|
Li H, Yi GY. Missing Data Mechanisms for Analysing Longitudinal Data with Incomplete Observations in Both Responses and Covariates. AUST NZ J STAT 2016. [DOI: 10.1111/anzs.12165] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Haocheng Li
- Departments of Oncology and Community Health Sciences; University of Calgary; Calgary AB Canada
| | - Grace Y. Yi
- Department of Statistics and Actuarial Science; University of Waterloo; Waterloo ON Canada
| |
Collapse
|
10
|
Erler NS, Rizopoulos D, Rosmalen JV, Jaddoe VWV, Franco OH, Lesaffre EMEH. Dealing with missing covariates in epidemiologic studies: a comparison between multiple imputation and a full Bayesian approach. Stat Med 2016; 35:2955-74. [PMID: 27042954 DOI: 10.1002/sim.6944] [Citation(s) in RCA: 48] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2015] [Revised: 03/01/2016] [Accepted: 03/02/2016] [Indexed: 12/19/2022]
Abstract
Incomplete data are generally a challenge to the analysis of most large studies. The current gold standard to account for missing data is multiple imputation, and more specifically multiple imputation with chained equations (MICE). Numerous studies have been conducted to illustrate the performance of MICE for missing covariate data. The results show that the method works well in various situations. However, less is known about its performance in more complex models, specifically when the outcome is multivariate as in longitudinal studies. In current practice, the multivariate nature of the longitudinal outcome is often neglected in the imputation procedure, or only the baseline outcome is used to impute missing covariates. In this work, we evaluate the performance of MICE using different strategies to include a longitudinal outcome into the imputation models and compare it with a fully Bayesian approach that jointly imputes missing values and estimates the parameters of the longitudinal model. Results from simulation and a real data example show that MICE requires the analyst to correctly specify which components of the longitudinal process need to be included in the imputation models in order to obtain unbiased results. The full Bayesian approach, on the other hand, does not require the analyst to explicitly specify how the longitudinal outcome enters the imputation models. It performed well under different scenarios. Copyright © 2016 John Wiley & Sons, Ltd.
Collapse
Affiliation(s)
- Nicole S Erler
- Department of Biostatistics, Erasmus MC, Wytemaweg 80, Rotterdam, 3015CN, The Netherlands.,Department of Epidemiology, Erasmus MC, Wytemaweg 80, Rotterdam, 3015CN, The Netherlands
| | - Dimitris Rizopoulos
- Department of Biostatistics, Erasmus MC, Wytemaweg 80, Rotterdam, 3015CN, The Netherlands
| | - Joost van Rosmalen
- Department of Biostatistics, Erasmus MC, Wytemaweg 80, Rotterdam, 3015CN, The Netherlands
| | - Vincent W V Jaddoe
- Department of Epidemiology, Erasmus MC, Wytemaweg 80, Rotterdam, 3015CN, The Netherlands.,Department of Pediatrics, Erasmus MC, Wytemaweg 80, Rotterdam, 3015CN, The Netherlands
| | - Oscar H Franco
- Department of Epidemiology, Erasmus MC, Wytemaweg 80, Rotterdam, 3015CN, The Netherlands
| | - Emmanuel M E H Lesaffre
- Department of Biostatistics, Erasmus MC, Wytemaweg 80, Rotterdam, 3015CN, The Netherlands.,L-Biostat, KU Leuven, Kapucijnenvoer 35, Leuven, 3000, Belgium
| |
Collapse
|
11
|
Chen B, Qin J. Test the reliability of doubly robust estimation with missing response data. Biometrics 2014; 70:289-98. [PMID: 24571677 DOI: 10.1111/biom.12150] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2012] [Revised: 12/01/2013] [Accepted: 01/01/2014] [Indexed: 11/29/2022]
Abstract
In statistical inference, one has to make sure that the underlying regression model is correctly specified otherwise the resulting estimation may be biased. Model checking is an important method to detect any departure of the regression model from the true one. Missing data are a ubiquitous problem in social and medical studies. If the underlying regression model is correctly specified, recent researches show great popularity of the doubly robust (DR) estimates method for handling missing data because of its robustness to the misspecification of either the missing data model or the conditional mean model, that is, the model for the conditional expectation of true regression model conditioning on the observed quantities. However, little work has been devoted to the goodness of fit test for DR estimates method. In this article, we propose a testing method to assess the reliability of the estimator derived from the DR estimating equation with possibly missing response and always observed auxiliary variables. Numerical studies demonstrate that the proposed test can control type I errors well. Furthermore the proposed method can detect departures from model assumptions in the marginal mean model of interest powerfully. A real dementia data set is used to illustrate the method for the diagnosis of model misspecification in the problem of missing response with an always observed auxiliary variable for cross-sectional data.
Collapse
Affiliation(s)
- Baojiang Chen
- Department of Biostatistics, College of Public Health, University of Nebraska Medical Center, Omaha, Nebraska 68198, U.S.A
| | | |
Collapse
|
12
|
Chen B, Qin J. A new estimation with minimum trace of asymptotic covariance matrix for incomplete longitudinal data with a surrogate process. Stat Med 2013; 32:4763-80. [DOI: 10.1002/sim.5875] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2012] [Accepted: 05/02/2013] [Indexed: 11/08/2022]
Affiliation(s)
- Baojiang Chen
- Department of Biostatistics; University of Nebraska Medical Center; Omaha NE 68198 U.S.A
| | - Jing Qin
- National Institute of Allergy and Infectious Diseases; National Institute of Health; Bethesda MD 20892 U.S.A
| |
Collapse
|
13
|
Chen B, Yi GY, Cook RJ, Zhou XH. Marginal methods for clustered longitudinal binary data with incomplete covariates. J Stat Plan Inference 2013; 142:2819-2831. [PMID: 23805025 DOI: 10.1016/j.jspi.2012.04.006] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
Many analyses for incomplete longitudinal data are directed to examining the impact of covariates on the marginal mean responses. We consider the setting in which longitudinal responses are collected from individuals nested within clusters. We discuss methods for assessing covariate effects on the mean and association parameters when covariates are incompletely observed. Weighted first and second order estimating equations are constructed to obtain consistent estimates of mean and association parameters when covariates are missing at random. Empirical studies demonstrate that estimators from the proposed method have negligible finite sample biases in moderate samples. An application to the National Alzheimer's Coordinating Center (NACC) Uniform Data Set (UDS) demonstrates the utility of the proposed method.
Collapse
Affiliation(s)
- Baojiang Chen
- Department of Biostatistics, University of Nebraska Medical Center, Omaha, NE 68198, USA
| | | | | | | |
Collapse
|
14
|
Chen B, Zhou XH. A latent-variable marginal method for multi-level incomplete binary data. Stat Med 2012; 31:3211-22. [PMID: 22733392 PMCID: PMC3631603 DOI: 10.1002/sim.5394] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2011] [Accepted: 03/13/2012] [Indexed: 12/27/2022]
Abstract
Incomplete multi-level data arise commonly in many clinical trials and observational studies. Because of multi-level variations in this type of data, appropriate data analysis should take these variations into account. A random effects model can allow for the multi-level variations by assuming random effects at each level, but the computation is intensive because high-dimensional integrations are often involved in fitting models. Marginal methods such as the inverse probability weighted generalized estimating equations can involve simple estimation computation, but it is hard to specify the working correlation matrix for multi-level data. In this paper, we introduce a latent variable method to deal with incomplete multi-level data when the missing mechanism is missing at random, which fills the gap between the random effects model and marginal models. Latent variable models are built for both the response and missing data processes to incorporate the variations that arise at each level. Simulation studies demonstrate that this method performs well in various situations. We apply the proposed method to an Alzheimer's disease study.
Collapse
Affiliation(s)
- Baojiang Chen
- Department of Biostatistics, University of Nebraska Medical Center, Omaha, NE 68198, U.S.A.
| | | |
Collapse
|