1
|
Tian X, Ciarleglio M, Cai J, Greene EJ, Esserman D, Li F, Zhao Y. Bayesian semi-parametric inference for clustered recurrent events with zero inflation and a terminal event. J R Stat Soc Ser C Appl Stat 2024; 73:598-620. [PMID: 39072299 PMCID: PMC11271983 DOI: 10.1093/jrsssc/qlae003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2022] [Revised: 10/19/2023] [Accepted: 01/05/2024] [Indexed: 07/30/2024]
Abstract
Recurrent events are common in clinical studies and are often subject to terminal events. In pragmatic trials, participants are often nested in clinics and can be susceptible or structurally unsusceptible to the recurrent events. We develop a Bayesian shared random effects model to accommodate this complex data structure. To achieve robustness, we consider the Dirichlet processes to model the residual of the accelerated failure time model for the survival process as well as the cluster-specific shared frailty distribution, along with an efficient sampling algorithm for posterior inference. Our method is applied to a recent cluster randomized trial on fall injury prevention.
Collapse
Affiliation(s)
- Xinyuan Tian
- Department of Biostatistics, Yale University, New Haven, CT, USA
| | - Maria Ciarleglio
- Department of Biostatistics, Yale University, New Haven, CT, USA
| | - Jiachen Cai
- Department of Biostatistics, Yale University, New Haven, CT, USA
| | - Erich J Greene
- Department of Biostatistics, Yale University, New Haven, CT, USA
| | - Denise Esserman
- Department of Biostatistics, Yale University, New Haven, CT, USA
| | - Fan Li
- Department of Biostatistics, Yale University, New Haven, CT, USA
| | - Yize Zhao
- Department of Biostatistics, Yale University, New Haven, CT, USA
| |
Collapse
|
2
|
Huang CH. Determination of correlations in multivariate count data with informative observation times. Stat Methods Med Res 2024; 33:273-294. [PMID: 38297977 DOI: 10.1177/09622802231224632] [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] [Indexed: 02/02/2024]
Abstract
We consider there are various types of recurrent events and the total number of occurrences are collected at the random observation times. It has concerned that the observation process may not be independent to the multivariate event processes, hence the total counts and observation times may be correlated and the dependence may exist among different types of the event processes as well. Many methods have developed nonparametric models to accommodate such unknown structures; however, it is difficult to assess and directly quantify their correlation relationships. A multivariate frailty model is proposed to this study, in which the event and observation processes are linked by frailty variables whose joint distribution can be implicitly specified through the multivariate normal distribution with some unknown covariance matrix. The Bayesian inference method is conducted to obtain the estimates of the regression coefficients and correlation parameters. We use a form of trigonometric functions to represent the covariance matrix, so that it meets the positive-definiteness condition efficiently during the estimation schemes. The simulation studies demonstrate the utility of the proposed models. We apply the model to a skin cancer prevention study, and aim to determine the covariate and association effects. We found treatment is significant in determining the duration of examination times; prior-counts, age and gender are significant variables on the occurrence rates of tumor counts. Using the covariance matrix to access the underlying dependent structure, the mutual correlations among them are all positive, and the basal cell counts are more related to the examination times.
Collapse
Affiliation(s)
- Chia-Hui Huang
- Department of Statistics, National Chengchi University, Taipei, Taiwan
| |
Collapse
|
3
|
Ghosh A, Chan W, Younes N, Davis BR. A Dynamic Risk Model for Multitype Recurrent Events. Am J Epidemiol 2023; 192:621-631. [PMID: 36549905 PMCID: PMC10404068 DOI: 10.1093/aje/kwac213] [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] [Received: 06/28/2021] [Revised: 10/17/2022] [Accepted: 12/14/2022] [Indexed: 12/24/2022] Open
Abstract
Recurrent events can occur more than once in the same individual; such events may be of different types, known as multitype recurrent events. They are very common in longitudinal studies. Often there is a terminating event, after which no further events can occur. The risk of any event, including terminating events such as death or cure, is typically affected by prior events. We propose a flexible joint multitype recurrent-events model that explicitly provides estimates of the change in risk for each event due to subject characteristics, including number and type of prior events and the absolute risk for every event type (terminating and nonterminating), and predicts event-free survival probability over a desired time period. The model is fully parametric, and therefore a standard likelihood function and robust standard errors can be constructed. We illustrate the model with applications to the Antihypertensive and Lipid-Lowering Treatment to Prevent Heart Attack Trial (1994-2002) and provide discussion of the results and model features.
Collapse
Affiliation(s)
- Alokananda Ghosh
- Correspondence to Dr. Alokananda Ghosh, Department of Biostatistics and Bioinformatics, Biostatistics Center, Milken Institute School of Public Health, George Washington University, 6110 Executive Boulevard, Rockville, MD 20852 (e-mail: ); or Dr. Wenyaw Chan, Department of Biostatistics and Data Science, School of Public Health, University of Texas Health Science Center at Houston, 1200 Pressler Street, Houston, TX 77030 (e-mail: )
| | - Wenyaw Chan
- Correspondence to Dr. Alokananda Ghosh, Department of Biostatistics and Bioinformatics, Biostatistics Center, Milken Institute School of Public Health, George Washington University, 6110 Executive Boulevard, Rockville, MD 20852 (e-mail: ); or Dr. Wenyaw Chan, Department of Biostatistics and Data Science, School of Public Health, University of Texas Health Science Center at Houston, 1200 Pressler Street, Houston, TX 77030 (e-mail: )
| | | | | |
Collapse
|
4
|
Ma C, Wang C, Pan J. Multistate modeling and structure selection for multitype recurrent events and terminal event data. Biom J 2023; 65:e2100334. [PMID: 36124712 DOI: 10.1002/bimj.202100334] [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: 10/21/2021] [Revised: 06/19/2022] [Accepted: 07/04/2022] [Indexed: 11/07/2022]
Abstract
In cardiovascular disease studies, a large number of risk factors are measured but it often remains unknown whether all of them are relevant variables and whether the impact of these variables is changing with time or remains constant. In addition, more than one kind of cardiovascular disease events can be observed in the same patient and events of different types are possibly correlated. It is expected that different kinds of events are associated with different covariates and the forms of covariate effects also vary between event types. To tackle these problems, we proposed a multistate modeling framework for the joint analysis of multitype recurrent events and terminal event. Model structure selection is performed to identify covariates with time-varying coefficients, time-independent coefficients, and null effects. This helps in understanding the disease process as it can detect relevant covariates and identify the temporal dynamics of the covariate effects. It also provides a more parsimonious model to achieve better risk prediction. The performance of the proposed model and selection method is evaluated in numerical studies and illustrated on a real dataset from the Atherosclerosis Risk in Communities study.
Collapse
Affiliation(s)
- Chuoxin Ma
- Guangdong Provincial Key Laboratory of Interdisciplinary Research and Application for Data Science, BNU-HKBU United International College, Zhuhai, China
| | - Chunyu Wang
- Department of Mathematics, University of Manchester, Manchester, UK
| | - Jianxin Pan
- Guangdong Provincial Key Laboratory of Interdisciplinary Research and Application for Data Science, BNU-HKBU United International College, Zhuhai, China.,Research Center for Mathematics, Beijing Normal University at Zhuhai, Zhuhai, China
| |
Collapse
|
5
|
Zhang Y, Chen M, Guo F. Bayesian criterion-based assessments of recurrent event models with applications to commercial truck driver behavior studies. Stat Med 2022; 41:4607-4628. [PMID: 35871759 PMCID: PMC9796651 DOI: 10.1002/sim.9528] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2021] [Revised: 03/17/2022] [Accepted: 06/26/2022] [Indexed: 01/01/2023]
Abstract
Multitype recurrent events are commonly observed in transportation studies, since commercial truck drivers may encounter different types of safety critical events (SCEs) and take different lengths of on-duty breaks in a driving shift. Bayesian nonhomogeneous Poisson process models are a flexible approach to jointly model the intensity functions of the multitype recurrent events. For evaluating and comparing these models, the deviance information criterion (DIC) and the logarithm of the pseudo-marginal likelihood (LPML) are studied and Monte Carlo methods are developed for computing these model assessment measures. We also propose a set of new concordance indices (C-indices) to evaluate various discrimination abilities of a Bayesian multitype recurrent event model. Specifically, the within-event C-index quantifies adequacy of a given model in fitting the recurrent event data for each type, the between-event C-index provides an assessment of the model fit between two types of recurrent events, and the overall C-index measures the model's discrimination ability among multiple types of recurrent events simultaneously. Moreover, we jointly model the incidence of SCEs and on-duty breaks with driving behaviors using a Bayesian Poisson process model with time-varying coefficients and time-dependent covariates. An in-depth analysis of a real dataset from the commercial truck driver naturalistic driving study is carried out to demonstrate the usefulness and applicability of the proposed methodology.
Collapse
Affiliation(s)
- Yiming Zhang
- Department of StatisticsUniversity of ConnecticutStorrsConnecticut
| | - Ming‐Hui Chen
- Department of StatisticsUniversity of ConnecticutStorrsConnecticut
| | - Feng Guo
- Department of StatisticsVirginia TechBlacksburgVirginia
| |
Collapse
|
6
|
Rahmati M, Rezanejad Asl P, Mikaeli J, Zeraati H, Rasekhi A. Compound Poisson frailty model with a gamma process prior for the baseline hazard: accounting for a cured fraction. J Appl Stat 2021; 49:3377-3391. [DOI: 10.1080/02664763.2021.1947997] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Maryam Rahmati
- Department of Epidemiology and Biostatistics, School of Public Health, Tehran University of Medical Sciences, Tehran, Iran
- Reproductive Endocrinology Research Center, Research Institute for Endocrine Sciences, Shahid Beheshti University of Medical Sciences, Tehran, Iran
| | - Parisa Rezanejad Asl
- Department of Biostatistics, Faculty of Allied Medical Sciences, Shahid Beheshti University of Medical Sciences, Tehran, Iran
| | - Javad Mikaeli
- Autoimmune and Motility Disorders Research Center, Digestive Diseases Research Institute, Tehran University of Medical Sciences, Tehran, Iran
| | - Hojjat Zeraati
- Department of Epidemiology and Biostatistics, School of Public Health, Tehran University of Medical Sciences, Tehran, Iran
| | - Aliakbar Rasekhi
- Biostatistics Department, Medical Sciences Faculty, Tarbiat Modares University, Tehran, Iran
| |
Collapse
|
7
|
Bedair KF, Hong Y, Al-Khalidi HR. Copula-frailty models for recurrent event data based on Monte Carlo EM algorithm. J STAT COMPUT SIM 2021. [DOI: 10.1080/00949655.2021.1942471] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Affiliation(s)
- Khaled F. Bedair
- Faculty of Commerce, Tanta University, Tanta, Egypt
- School of Medicine, University of Dundee, Dundee, UK
| | - Yili Hong
- Department of Statistics, Virginia Tech, Blacksburg, VA, USA
| | | |
Collapse
|
8
|
Rahmati M, Mahmoudi M, Mohammad K, Mikaeli J, Zeraati H. Bayesian modelling of zero‐inflated recurrent events and dependent termination with compound Poisson frailty model. Stat (Int Stat Inst) 2020. [DOI: 10.1002/sta4.292] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Affiliation(s)
- Maryam Rahmati
- Department of Epidemiology and Biostatistics, School of Public Health Tehran University of Medical Sciences Tehran Iran
| | - Mahmood Mahmoudi
- Department of Epidemiology and Biostatistics, School of Public Health Tehran University of Medical Sciences Tehran Iran
| | - Kazem Mohammad
- Department of Epidemiology and Biostatistics, School of Public Health Tehran University of Medical Sciences Tehran Iran
| | - Javad Mikaeli
- Autoimmune and Motility Disorders Research Center, Digestive Diseases Research Institute Tehran University of Medical Sciences Tehran Iran
| | - Hojjat Zeraati
- Department of Epidemiology and Biostatistics, School of Public Health Tehran University of Medical Sciences Tehran Iran
| |
Collapse
|
9
|
Affiliation(s)
- Yi Liu
- Department of Statistics, Virginia Tech, Blacksburg, VA
- Virginia Tech Transportation Institute, Blacksburg, VA
| | - Feng Guo
- Department of Statistics, Virginia Tech, Blacksburg, VA
- Virginia Tech Transportation Institute, Blacksburg, VA
| |
Collapse
|
10
|
Li Z, Chinchilli VM, Wang M. A Bayesian joint model of recurrent events and a terminal event. Biom J 2018; 61:187-202. [PMID: 30479030 DOI: 10.1002/bimj.201700326] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2017] [Revised: 08/16/2018] [Accepted: 09/07/2018] [Indexed: 11/06/2022]
Abstract
Recurrent events could be stopped by a terminal event, which commonly occurs in biomedical and clinical studies. In this situation, dependent censoring is encountered because of potential dependence between these two event processes, leading to invalid inference if analyzing recurrent events alone. The joint frailty model is one of the widely used approaches to jointly model these two processes by sharing the same frailty term. One important assumption is that recurrent and terminal event processes are conditionally independent given the subject-level frailty; however, this could be violated when the dependency may also depend on time-varying covariates across recurrences. Furthermore, marginal correlation between two event processes based on traditional frailty modeling has no closed form solution for estimation with vague interpretation. In order to fill these gaps, we propose a novel joint frailty-copula approach to model recurrent events and a terminal event with relaxed assumptions. Metropolis-Hastings within the Gibbs Sampler algorithm is used for parameter estimation. Extensive simulation studies are conducted to evaluate the efficiency, robustness, and predictive performance of our proposal. The simulation results show that compared with the joint frailty model, the bias and mean squared error of the proposal is smaller when the conditional independence assumption is violated. Finally, we apply our method into a real example extracted from the MarketScan database to study the association between recurrent strokes and mortality.
Collapse
Affiliation(s)
- Zheng Li
- Department of Public Health Sciences, Penn State College of Medicine, Hershey, Pennslyvania, USA
| | - Vernon M Chinchilli
- Department of Public Health Sciences, Penn State College of Medicine, Hershey, Pennslyvania, USA
| | - Ming Wang
- Department of Public Health Sciences, Penn State College of Medicine, Hershey, Pennslyvania, USA
| |
Collapse
|
11
|
Abstract
Recurrent event outcomes are ubiquitous among clinical trial data which encourages a conventional approach to analysis. Yet a common feature of these data has received less attention, that is, survival times often comprise multiple types of events that may imply a disparity in cost and disease severity. Typically, we neglect this feature of the data by combining event-types or analyzing each type separately, thus ignoring any interdependence among them. This practice may reflect a dearth of readily available methods and software that more appropriately acknowledge the true data structure. We provide a review of the literature on multitype recurrent events and frailty modelling which reflects a renewed interest in the topic over the past decade and the emergence of software for estimation. Thus, a review of available methods seems timely, if not overdue.
Collapse
Affiliation(s)
- Paul M Brown
- Department of Medicine, University of Alberta, Edmonton, Canada
- Canadian VIGOUR Centre, Edmonton, Canada
| | | |
Collapse
|
12
|
Integration of DNA methylation and gene transcription across nineteen cell types reveals cell type-specific and genomic region-dependent regulatory patterns. Sci Rep 2017; 7:3626. [PMID: 28620196 PMCID: PMC5472622 DOI: 10.1038/s41598-017-03837-z] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2017] [Accepted: 05/08/2017] [Indexed: 12/29/2022] Open
Abstract
Despite numerous studies done on understanding the role of DNA methylation, limited work has focused on systems integration of cell type-specific interplay between DNA methylation and gene transcription. Through a genome-wide analysis of DNA methylation across 19 cell types with T-47D as reference, we identified 106,252 cell type-specific differentially-methylated CpGs categorized into 7,537 differentially (46.6% hyper- and 53.4% hypo-) methylated regions. We found 44% promoter regions and 75% CpG islands were T-47D cell type-specific methylated. Pyrosequencing experiments validated the cell type-specific methylation across three benchmark cell lines. Interestingly, these DMRs overlapped with 1,145 known tumor suppressor genes. We then developed a Bayesian Gaussian Regression model to measure the relationship among DNA methylation, genomic segment distribution, differential gene expression and tumor suppressor gene status. The model uncovered that 3′UTR methylation has much less impact on transcriptional activity than other regions. Integration of DNA methylation and 82 transcription factor binding information across the 19 cell types suggested diverse interplay patterns between the two regulators. Our integrative analysis reveals cell type-specific and genomic region-dependent regulatory patterns and provides a perspective for integrating hundreds of various omics-seq data together.
Collapse
|