A copula-based Markov chain model for serially dependent event times with a dependent terminal event. (English) Zbl 1478.62257

Jpn. J. Stat. Data Sci. 4, No. 2, 917-951 (2021); correction ibid. 4, No. 1, 475 (2021).
Summary: Copula modeling for serial dependence has been extensively discussed in a time series context. However, fitting copula-based Markov models for serially dependent survival data is challenging due to the complex censoring mechanisms. The purpose of this paper is to develop likelihood-based methods for fitting a copula-based Markov chain model to serially dependent event times that are dependently censored by a terminal event, such as death. We propose a novel copula-based Markov chain model for describing serial dependence in recurrent event times. We also apply another copula model for handling dependent censoring. Due to the complex likelihood function with the two copulas, we propose a two-stage estimation method under Weibull distributions for fitting the survival data. The asymptotic normality of the proposed estimator is established through the theory of estimating functions. We propose a jackknife method for interval estimates, which is shown to be asymptotically consistent. To select suitable copulas for a given dataset, we propose a model selection method according to the 2nd stage likelihood. We conduct simulation studies to assess the performance of the proposed methods. For illustration, we analyze survival data from colorectal cancer patients. We implement the proposed methods in our original R package “Copula.Markov.survival” that is made available in CRAN (https://cran.r-project.org/).


62M10 Time series, auto-correlation, regression, etc. in statistics (GARCH)
62H05 Characterization and structure theory for multivariate probability distributions; copulas
60J20 Applications of Markov chains and discrete-time Markov processes on general state spaces (social mobility, learning theory, industrial processes, etc.)
62G05 Nonparametric estimation
62H12 Estimation in multivariate analysis
62E20 Asymptotic distribution theory in statistics
62N01 Censored data models
62N05 Reliability and life testing
62P10 Applications of statistics to biology and medical sciences; meta analysis
Full Text: DOI


[1] Balakrishnan, N.; Chimitova, E., Testing goodness of fit of parametric AFT and PH models with residuals, Communications in Statistics-Simulation and Computation, 42, 6, 1352-1367 (2013) · Zbl 1347.62212
[2] Bedair, K.; Hong, Y., Multivariate frailty models for multi-type recurrent event data and its application to cancer prevention trial, Computational Statistics & Data Analysis, 101, 161-173 (2016) · Zbl 1466.62026
[3] Cai, J.; Schaubel, DE, Marginal means/rates models for multiple type recurrent event data, Lifetime Data Analysis, 10, 121-138 (2004) · Zbl 1058.62098
[4] Chen, X.; Fan, Y., Estimation and model selection of semiparametric copula-based multivariate dynamic models under copula misspecification, Journal of Econometrics, 135, 1-2, 125-154 (2006) · Zbl 1418.62425
[5] Cox, DR; Reid, N., Parameter orthogonality and approximate conditional inference, Journal of the Royal Statistical Society: Series B (Methodological), 49, 1, 1-18 (1987) · Zbl 0616.62006
[6] Darsow, WF; Nguyen, B.; Olsen, ET, Copulas and Markov processes, Illinois Journal of Mathematics, 36, 4, 600-642 (1992) · Zbl 0770.60019
[7] Domma, F.; Giordano, S.; Perri, PF, Statistical modeling of temporal dependence in financial data via a copula function, Communications in Statistics-Simulation and Computation, 38, 4, 703-728 (2009) · Zbl 1290.62096
[8] Duchateau, L.; Janssen, P., The Frailty Model (2007), New York: Springer, New York · Zbl 1210.62153
[9] Durante, F.; Sempi, C., Principles of Copula theory (2016), London: Chapman and Hall/CRC, London · Zbl 1380.62008
[10] Emura, T. (2020). joint.Cox: joint frailty-copula models for tumour progression and death in meta-analysis, CRAN, https://CRAN.R-project.org/package=joint.Cox.
[11] Emura, T.; Chen, Y-H, Analysis of survival data with dependent censoring, copula-based approaches, JSS research series in statistic (2018), Singapore: Springer, Singapore · Zbl 1464.62001
[12] Emura, T.; Hsu, J-H, Estimation of the Mann-Whitney effect in the two-sample problem under dependent censoring, Comptational Statistics & Data Analysis, 150, 106990 (2020) · Zbl 07212330
[13] Emura, T.; Lin, C-W; Wang, W., A goodness-of-fit test for Archimedean copula models in the presence of right censoring, Computational Statistics & Data Analysis, 54, 3033-3043 (2010) · Zbl 1284.62597
[14] Emura, T.; Long, T-H; Sun, L-H, R routines for performing estimation and statistical process control under copula-based time series models, Communications in Statistics-Simulation and Computation, 46, 4, 3067-3087 (2017) · Zbl 1373.62583
[15] Emura, T.; Nakatochi, M.; Murotani, K.; Rondeau, V., A joint frailty-copula model between tumour progression and death for meta-analysis, Statistical Methods in Medical Research, 26, 6, 2649-2666 (2017)
[16] Emura, T.; Matsui, S.; Rondeau, V., Survival analysis with correlated endpoints, joint frailty-copula models, JSS Research Series in Statistics (2019), Singapore: Springer, Singapore · Zbl 1429.62003
[17] Erkal Sonmez, O.; Baray, A.; Calisir, F.; Cevikcan, E.; Camgoz Akdag, H., On copula based serial dependence in statistical process control, Industrial engineering in the big data Era (2019), Cham: Springer, Cham
[18] Gijbels, I.; Veraverbeke, N.; Omelka, M., Conditional copulas, association measures and their applications, Computational Statistics & Data Analysis, 55, 5, 1919-1932 (2011) · Zbl 1328.62366
[19] González, JR; Fernandez, E.; Moreno, V., Sex differences in hospital readmission among colorectal cancer patients, Journal of Epidemiology and Community Health, 59, 6, 506-511 (2005)
[20] Ha, ID; Jeong, JH; Lee, Y., Statistical modelling of survival data with random effects (2017), Singapore: Springer, Singapore · Zbl 1464.62002
[21] Huang, X-W; Emura, T., Model diagnostic procedures for copula-based Markov chain models for statistical process control, Communications in Statistics-Simulation and Computation (2019)
[22] Hung, W-L, Weighted least-squares estimation of the shape parameter of the Weibull distribution, Quality and Reliability Engineering International, 17, 467-469 (2001)
[23] Joe, H., Multivariate models and multivariate dependence concepts (1997), London: Chapman and Hall/CRC, London · Zbl 0990.62517
[24] Joe, H., Asymptotic efficiency of the two-stage estimation method for copula-based models, Journal of Multivariate Analysis, 94, 2, 401-419 (2005) · Zbl 1066.62061
[25] Kim, J-M; Baik, J., Anomaly detection in sensor data, Journal of Applied Reliability, 18, 1, 20-32 (2018)
[26] Kim, J.-M., Baik, J., & Reller, M. (2018). Detecting the change of variance by using conditional distribution with diverse copula functions. In: Paper presented at the Proceedings of the Pacific Rim Statistical Conference for Production Engineering.
[27] Kim, J-M; Baik, J.; Reller, M., Control charts of mean and variance using copula Markov SPC and conditional distribution by copula, Communications in Statistics-Simulation and Computation (2019)
[28] Lawless, JF, Statistical models and methods for lifetime data (2003), Hoboken: Wiley, Hoboken · Zbl 1015.62093
[29] Lawless, JF; Yilmaz, YE, Semiparametric estimation in copula models for bivariate sequential survival times, Biometrical Journal, 53, 5, 779-796 (2011) · Zbl 1226.62118
[30] Li, Z.; Chinchilli, VM; Wang, M., A Bayesian joint model of recurrent events and a terminal event, Biometrical Journal, 61, 1, 187-202 (2019) · Zbl 1412.62166
[31] Li, F.; Tang, Y.; Wang, HJ, Copula-based semiparametric analysis for time series data with detection limits, Canadian Journal of Statistics, 47, 3, 438-454 (2019) · Zbl 07194110
[32] Li, Z.; Chinchilli, VM; Wang, M., A time-varying Bayesian joint hierarchical copula model for analysing recurrent events and a terminal event: an application to the Cardiovascular Health Study, Journal of the Royal Statistical Society: Series C, 69, 1, 151-166 (2020)
[33] Lin, W-C; Emura, T.; Sun, L-H, Estimation under copula-based Markov mixture normal models for serially correlated data, Communications in Statistics-Simulation and Computation (2019)
[34] Liu, X., Planning of accelerated life tests with dependent failure modes based on a gamma frailty model, Technometrics, 54, 4, 398-409 (2012)
[35] Lo, SM; Mammen, E.; Wilke, RA, A nested copula duration model for competing risks with multiple spells, Computational Statistics & Data Analysis, 150, 106986 (2020) · Zbl 07212326
[36] MacDonald, IL, Does Newton-Raphson really fail?, Statistical Methods in Medical Research, 23, 3, 308-311 (2014)
[37] Meyer, R.; Romeo, JS, Bayesian semiparametric analysis of recurrent failure time data using copulas, Biometrical Journal, 57, 6, 982-1001 (2015) · Zbl 1386.62060
[38] Nelsen, RB, Properties of a one-parameter family of bivariate distributions with specified marginals, Communications in Statistics-Theory and Methods, 15, 11, 3277-3285 (1986) · Zbl 0616.62067
[39] Nelsen, RB, An introduction to copulas (2006), New York: Springer, New York · Zbl 1152.62030
[40] Ning, J.; Rahbar, MH, Estimating the ratio of multivariate recurrent event rates with application to a blood transfusion study, Statistical Methods in Medical Research, 26, 4, 1969-1981 (2017)
[41] Rotolo, F.; Paoletti, X.; Michiels, S., surrosurv: An R package for the evaluation of failure time surrogate endpoints in individual patient data meta-analyses of randomized clinical trials, Computer Methods and Programs in Biomedicine, 155, 189-198 (2018)
[42] Schneider, S.; Demarqui, FN; Colosimo, EA; Mayrink, VD, An approach to model clustered survival data with dependent censoring, Biometrical Journal, 62, 1, 157-174 (2020) · Zbl 1441.62482
[43] Shih, JH; Klein, JP; Van Houwelingen, HC; Ibrahim, JG; Scheike, TH, Copula models, Handbook of survival analysis (2014), Boca Raton: CRC Press, Boca Raton
[44] Shinohara, S., Lin, Y. H., Michimae, H., & Emura. T. (2020) Dynamic lifetime prediction using a Weibull-based bivariate failure time model: a meta-analysis of individual-patient data, in reivew.
[45] Su, C.-L., Lin, F.-C. (2020). Analysis of cyclic recurrent event data with multiple event types, Japanese Journal of Statistics and Data Science. doi:10.1007/s42081-020-00088-7.
[46] Sun, L-H; Huang, X-W; Alqawba, M-S; Kim, JM; Emura, T., Copula-based Markov Models for time series-parametric inference and process control, JSS research series in statistics (2020), Berlin: Springer,, Berlin · Zbl 1435.62020
[47] Sun, L-H; Lee, C-S; Emura, T., A Bayesian inference for time series via copula-based Markov chain models, Communications in Statistics-Simulation and Computation (2018)
[48] Valle, LD; Leisen, F.; Rossini, L., Bayesian non-parametric conditional copula estimation of twin data, Journal of the Royal Statistical Society: Series C, 67, 3, 523-548 (2018)
[49] Wang, W.; Emura, T., Comments on inference in multivariate Archimedean copula models by Genest et al, TEST, 20, 276-280 (2011) · Zbl 1367.62171
[50] Wang, K.; Yau, KK; Lee, AH, Multilevel survival modelling of recurrent urinary tract infections, Computer Methods and Programs in Biomedicine, 87, 3, 225-229 (2007)
[51] Wang, YC; Emura, T.; Fan, TH; Lo, SM; Wilke, RA, Likelihood-based inference for a frailty-copula model based on competing risks failure time data, Qual. Reliab. Eng. Int., 36, 5, 1622-1638 (2020)
[52] Weibull, W., Wide applicability, Journal of Applied Mechanics, 103, 730, 293-297 (1951) · Zbl 0042.37903
[53] Wu, BH; Michimae, H.; Emura, T., Meta-analysis of individual patient data with semi-competing risks under the Weibull joint frailty-copula model, Computational Statistics (2020) · Zbl 07311686
[54] Zhang, S.; Zhou, QM; Lin, H., Goodness-of-fit test of copula functions for semi-parametric univariate time series models, Statistical Papers (2020)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.