×

A Gamma-frailty proportional hazards model for bivariate interval-censored data. (English) Zbl 1469.62065

Summary: Correlated survival data naturally arise from many clinical and epidemiological studies. For the analysis of such data, the Gamma-frailty proportional hazards (PH) model is a popular choice because the regression parameters have marginal interpretations and the statistical association between the failure times can be explicitly quantified via Kendall’s tau. Despite their popularity, Gamma-frailty PH models for correlated interval-censored data have not received as much attention as analogous models for right-censored data. A Gamma-frailty PH model for bivariate interval-censored data is presented and an easy to implement expectation-maximization (EM) algorithm for model fitting is developed. The proposed model adopts a monotone spline representation for the purposes of approximating the unknown conditional cumulative baseline hazard functions, significantly reducing the number of unknown parameters while retaining modeling flexibility. The EM algorithm was derived from a data augmentation procedure involving latent Poisson random variables. Extensive numerical studies illustrate that the proposed method can provide reliable estimation and valid inference, and is moreover robust to the misspecification of the frailty distribution. To further illustrate its use, the proposed method is used to analyze data from an epidemiological study of sexually transmitted infections.

MSC:

62-08 Computational methods for problems pertaining to statistics
62P10 Applications of statistics to biology and medical sciences; meta analysis
62N01 Censored data models

Software:

bayesSurv
PDFBibTeX XMLCite
Full Text: DOI Link

References:

[1] Andersen, P.; Klein, J.; Knudsen, K.; Palacios, R., Estimation of variance in cox’s regression model with shared gamma frailties, Biometrics, 53, 1475-1484, (1997) · Zbl 0911.62030
[2] Cai, B.; Lin, X.; Wang, L., Bayesian proportional hazards model for current status data with monotone splines, Comput. Statist. Data Anal., 55, 2644-2651, (2011) · Zbl 1464.62034
[3] Chang, I.; Wen, C.; Wu, Y., A profile likelihood theory for the correlated gamma-frailty model with current status family data, Statist. Sinica, 17, 1023-1046, (2007) · Zbl 1133.62021
[4] Chen, M.; Tong, X.; Sun, J., The proportional odds model for multivariate interval-censored failure time data, Stat. Med., 26, 5147-5161, (2007)
[5] Chen, M.; Tong, X.; Sun, J., A frailty model approach for regression analysis of multivariate current status data, Stat. Med., 28, 3424-3436, (2009)
[6] Chen, M.; Tong, X.; Zhu, L., A linear transformation model for multivariate interval-censored failure time data, Canad. J. Statist., 41, 275-290, (2013) · Zbl 1273.62239
[7] Cui, S.; Sun, Y., Checking for the gamma frailty distribution under the marginal proportional hazards frailty model, Statist. Sinica, 14, 249-267, (2004) · Zbl 1054.62118
[8] Dunson, D.; Dinse, G., Bayesian models for multivariate current status data with informative censoring, Biometrics, 58, 79-88, (2002) · Zbl 1209.62031
[9] Freedman, D., On the so-called huber sandwich estimator and robust standard errors, Amer. Statist., 60, 299-302, (2006)
[10] Ghosh, P.; Tu, W., Assessing sexual attitudes and behaviors of Young women: a joint model with nonlinear time effects, time varying covariates, and dropouts, J. Amer. Statist. Assoc., 104, 474-485, (2009) · Zbl 1388.62254
[11] Goggins, W.; Finkelstein, D., A proportional hazards model for multivariate interval-censored failure time data, Biometrics, 56, 940-943, (2000) · Zbl 1060.62617
[12] Hens, N.; Wienke, A.; Aerts, M.; Molenberghs, G., The correlated and shared gamma frailty model for bivariate current status data: an illustration for cross-sectional serological data, Stat. Med., 28, 2785-2800, (2009)
[13] Henschel, V.; Engel, J.; Hölzel, D.; Mansmann, U., A semiparametric Bayesian proportional hazards model for interval censored data with frailty effects, BMC Med. Res. Methodol., 9, 9, (2009)
[14] Hougaard, P., Analysis of multivariate survival data, (2000), Springer New York · Zbl 0962.62096
[15] Ibrahim, J.; Chen, M.; Sinha, D., Bayesian survival analysis, (2008), Springer New York
[16] Kim, M.; Xue, X., The analysis of multivariate interval-censored survival data, Stat. Med., 21, 3715-3726, (2002)
[17] Klein, J., Semiparametric estimation of random effects using the Cox model based on the EM algorithm, Biometrics, 48, 795-806, (1992)
[18] Komarek, A.; Lessaffre, E., Bayesian accelerated failure time model for correlated interval-censored data with a normal mixture as error distribution, Statist. Sinica, 17, 549-569, (2007) · Zbl 1144.62022
[19] Lam, K.; Xu, Y.; Cheung, T., A multiple imputation approach for clustered interval-censored survival data, Stat. Med., 29, 680-693, (2010)
[20] Lin, X.; Wang, L., A semiparametric probit model for case 2 interval-censored failure time data, Stat. Med., 29, 972-981, (2010)
[21] Lin, X.; Wang, L., Bayesian proportional odds models for analyzing current status data: univariate, clustered, and multivariate, Comm. Statist. Simulation Comput., 40, 1171-1181, (2011) · Zbl 1227.62015
[22] Liu, H.; Shen, Y., A semiparametric regression cure model for interval-censored data, J. Amer. Statist. Assoc., 104, 1168-1178, (2009) · Zbl 1388.62283
[23] Louis, T., Finding the observed information matrix when using the EM algorithm, J. R. Stat. Soc. Ser. B Stat. Methodol., 44, 226-233, (1982) · Zbl 0488.62018
[24] McMahan, C.; Wang, L.; Tebbs, J., Regression analysis for current status data using the EM algorithm, Stat. Med., 32, 4452-4466, (2013)
[25] Ramsay, J., Monotone regression splines in action, Statist. Sci., 3, 425-461, (1988)
[26] Rondeau, V.; Commenges, D.; Joly, P., Maximum penalized likelihood estimation in a gamma-frailty model, Lifetime Data Anal., 9, 139-153, (2003) · Zbl 1116.62408
[27] Shen, P., Additive transformation models for multivariate interval-censored data, Comm. Statist. Theory Methods, 44, 1065-1079, (2015) · Zbl 1325.62182
[28] Tong, X.; Chen, M.; Sun, J., Regression analysis of multivariate interval-censored failure time data with application to tumorigenicity experiments, Biom. J., 50, 364-374, (2008)
[29] Tu, W.; Batteiger, B.; Wiehe, S.; Ofner, S.; Pol, B. V.D.; Katz, B.; Orr, D.; Fortenberry, J., Time from first intercourse to first sexually transmitted infection diagnosis among adolescent women, Arch. Pediatr. Adolesc. Med., 163, 1106-1111, (2009)
[30] Tu, W.; Ghosh, P.; Katz, B., A stochastic model for assessing chlamydia trachomatis transmission risk by using longitudinal observational data, J. R. Stat. Soc. Ser. A, 174, 975-989, (2011)
[31] Turnbull, B., The empirical distribution function with arbitrarily grouped, censored and truncated data, J. R. Stat. Soc. Ser. B Stat. Methodol., 38, 290-295, (1976) · Zbl 0343.62033
[32] Wang, L.; Dunson, D., Semiparametric bayes’ proportional odds models for current status data with underreporting, Biometrics, 67, 1111-1118, (2011) · Zbl 1226.62131
[33] Wang, L.; McMahan, C.; Hudgens, M.; Qureshi, Z., A flexible, computationally efficient method for Fitting the proportional hazards model to interval-censored data, Biometrics, 72, 222-231, (2016) · Zbl 1393.62105
[34] Wang, L.; Sun, L.; Sun, J., A goodness-of-fit test for the marginal Cox model for correlated interval-censored failure time data, Biom. J., 48, 1020-1028, (2006)
[35] Wang, N.; Wang, L.; McMahan, C., Regression analysis of bivariate current status data under the gamma-frailty proportional hazards model using the EM algorithm, Comput. Statist. Data Anal., 83, 140-150, (2015)
[36] Wei, L.; Lin, D.; Weissfeld, L., Regression analysis of multivariate incomplete failure time data by modeling marginal distributions, J. Amer. Statist. Assoc., 84, 1065-1073, (1989)
[37] Wen, C.; Chen, Y., Nonparametric maximum likelihood analysis of clustered current status data with the gamma-frailty Cox model, Comput. Statist. Data Anal., 55, 1053-1060, (2011) · Zbl 1284.62067
[38] Wen, C.; Chen, Y., A frailty model approach for regression analysis of bivariate interval-cenosred survival data, Statist. Sinica, 23, 383-408, (2013) · Zbl 1257.62048
[39] Wienke, A., Frailty models in survival analysis, (2012), Chapman & Hall
[40] Workowski, K.; Bolan, G., Sexually transmitted diseases treatment guidelines, Morb. Mortal. Wkly. Rep. (MMWR), 64, 1-137, (2015)
[41] Yavuz, A.Ç.; Lambert, P., Semi-parametric frailty model for clustered interval-censored data, Stat. Model., 360-391, (2016)
[42] Yin, G.; Ibrahim, J., A class of Bayesian shared gamma frailty models with multivariate failure time data, Biometrics, 61, 208-216, (2005) · Zbl 1077.62088
[43] Yu, Z.; Lin, X.; Tu, W., Semiparametric frailty models for clustered failure time data, Biometrics, 68, 429-436, (2012) · Zbl 1251.62011
[44] Zhang, Z.; Sun, J., Interval censoring, Stat. Methods Med. Res., 19, 53-70, (2010)
[45] Zhou, J.; Zhang, J.; Lu, W., An expectation maximization algorithm for Fitting the generalized odds-rate model to interval censored data, Stat. Med., 36, 7, 1157-1171, (2017)
[46] Zuma, K., A Bayesian analysis of correlated interval-censored data, Comm. Statist. Theory Methods, 36, 725-730, (2007) · Zbl 1109.62113
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. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.