Clustered partly interval-censored survival data naturally arise from many medical and epidemiological studies. We propose a Bayesian semiparametric approach for fitting a mixed effects proportional hazards (PH) model to clustered partly interval-censored data. The proposed method allows for not only a random intercept as most frailty models do for clustered survival data, but also random effects of covariates. We assume a normal prior for each random intercept/random effect, seeing the instability of a gamma prior for a frailty in this situation. Simulation studies with data generated from both mixed effects PH model and mixed effects accelerated failure times model are conducted, to evaluate the performance of the proposed method and compare it with the three methods currently available in the literature. The application of the proposed approach is illustrated through analyzing the progression-free survival data derived from a phase III metastatic colorectal cancer clinical trial.
AslanidouH, DeyDK and SinhaD (1998) Bayesian analysis of multivariate survival data using Monte Carlo methods. The Canadian Journal of Statistics, 26, 33–48.
2.
BellamyS, LiY, RyanLM, LipsitzS, CannerMJ and WrightR (2004) Analysis of clustered and interval censored data from a community-based study in asthma. Statistics in Medicine, 23, 3607–3621
3.
ChenC, ShenP and LeeT (2020) A gamma-frailty model for interval-censored data with dependent examination times: a computationally efficient approach. Communications in Statistics-Simulation and Computation. doi: 10.1080/03610918.2020.1790600.
4.
CoxDR (1972) Regression models and life-tables (with discussion). Journal of the Royal Statistical Society: Series B, 34, 187–220.
5.
DeyDK, ChenM and ChangH (1997) Bayesian approach for nonlinear random effects models. Biometrics, 53, 1239–1252.
6.
GaoF, ZengD and LinDY (2017) Semiparametric estimation of the accelerated failure time model with partly interval-censored data. Biometrics, 73, 1161–1168.
7.
GeisserS and EddyWF (1979) A predictive approach to model selection. Journal of the American Statistical Association, 74, 153–160.
8.
GelfandAE (1992) Model determination using predictive distributions with implementation via sampling-based methods (with discussion). In Bayesian Statistics 4, edited byBernardoJM, BergerJO, DawidAP and SmithAFM, pages 147–167. Oxford:Oxford University Press.
9.
GliddenD and VittinghoffE (2004) Modelling clustered survival data from multicentre clinical trials. Statistics in Medicine, 23, 369–388.
10.
GoethalsK, AmpeB, BerkvensD, LaevensH, JanssenP and DuchateauL (2009) Modeling interval-censored, clustered cow udder quarter infection times through the shared gamma frailty model. Journal of Agricultural, Biological, and Environmental Statistics, 14, 1–14.
11.
GörürD and RasmussenCE (2010) Dirichlet process Gaussian mixture models: Choice of the base distribution. Journal of Computer Science and Technology, 25, 615–626.
12.
HastingsWK (1970) Monte Carlo sampling methods using Markov chains and their applications. Biometrika, 57, 97–109.
13.
HuangJ and WellnerJA (1997) Interval censored survival data: A review of recent progress. In Proceedings of the First Seattle Symposium in Biostatistics, edited byLinDY and FlemingTR, pages 123–169. New York:Springer.
14.
HuangJ (1999) Asymptotic properties of nonparametric estimation based on partly interval-censored data. Statistica Sinica, 9, 501–519.
15.
IshwaranH and ZarepourM (2000) Markov chain Monte Carlo in approximate Dirichlet and beta two-parameter process hierarchical models. Biometrika, 87, 371–390.
16.
IshwaranH and JamesL (2001) Gibbs sampling methods for stick-breaking priors. Journal of the American Statistical Association, 96, 161–173.
17.
KimJS (2003) Maximum likelihood estimation for the proportional hazards model with partly interval-censored data. Journal of the Royal Statistical Society: Series B, 65, 489–502.
18.
KleinJP (1992) Semiparametric estimation of random effects using the Cox model based on the EM algorithm. Biometrics, 48, 795–806.
19.
KomárekA and LesaffreE (2007) Bayesian accelerated failure time model for correlated interval-censored data with a normal mixture as an error distribution. Statistica Sinica, 17, 549–569.
20.
KomárekA, LesaffreE and LegrandC (2007) Baseline and treatment effect heterogeneity for survival times between centers using a random effects accelerated failure time model with flexible error distribution. Statistics in Medicine, 26, 5457–5472.
21.
KomárekA (2018) bayesSurv: Bayesian survival regression with flexible error and random effects distributions. https://cran.r-project.org/package=bayesSurv. R package version 3.2.
22.
PanC, CaiB and WangL (2015) Multiple frailty model for clustered interval-censored data with frailty selection. Statistical Methods in Medical Research, 26, 1308–1322.
23.
PanC, CaiB and WangL (2020) A Bayesian approach for analyzing partly interval-censored data under the proportional hazards model. Statistical Methods in Medical Research. doi: 10.1177/0962280220921552.
24.
PanC and CaiB (2020) A Bayesian model for spatial partly interval-censored data. Communication in Statistics-Simulation and Computation. doi: 10.1080/03610918.2020.1839497.
25.
PanC (2021) PICBayes: Bayesian models for partly interval-censored data. https://cran.r-project.org/package=PICBayes. R package version 1.0.
26.
PeetersM, PriceTJ, CervantesA, SobreroAF, DucreuxM, AndréT, ChanE, LordickF, PuntCJA, StricklandAH, WilsonG, CiuleanuTE, RomanL, CutsemEV, TzekovaV, CollinsS, OlinerKS, RongA and GansertJ (2010) Randomized phase III study of panitumumab with fluorouracil, leucovorin, and irinotecan (FOLFIRI) compared with FOLFIRI alone as second-line treatment in patients with metastatic colorectal cancer. Journal of Clinical Oncology, 28, 4706–4713.
27.
PeetersM, PriceTJ, CervantesA, SobreroAF, DucreuxM, HotkoY, AndréT, ChanE, LordickF, PuntCJA, StricklandAH, WilsonG, CiuleanuTE, RomanL, TianY and SidhuR (2014) Final results from a randomized phase 3 study of FOLFIRI ± panitumumab for second-line treatment of metastatic colorectal cancer. Annals of Oncology, 25, 107–116.
28.
PlummerM, BestN, CowlesK, VinesK, SarkarD, BatesD, AlmondR and MagnussonA (2019) coda: Output analysis and diagnostics for MCMC. https://cran.r-project.org/package=coda. R package version 0.19-3.
RipattiS and PalmgrenJ (2000) Estimation of multivariate frailty models using penalized partial likelihood. Biometrics, 56, 1016–1022.
31.
SargentDJ (1998) A general framework for random effects survival analysis in the Cox proportional hazards setting. Biometrics, 54, 1486–1497.
32.
SpiegelhalterDJ, ThomasA, BestNG and GilksWR (1996) BUGS Examples Volume 1, Version 0.5 (Version ii). Cambridge: Medical Research Council Biostatistics Unit.
33.
SpiegelhalterDJ, BestNG, CarlinBP and van der LindeA (2002) Bayesian measures of model complexity and fit. Journal of the Royal Statistical Society: Series B, 64, 583–639.
34.
SpiegelhalterDJ, BestNG, CarlinBP and van der LindeA (2014) The deviance information criterion: 12 years on. Journal of the Royal Statistical Society: Series B, 76, 458–493.
35.
TherneauTM, LumleyT, AtkinsonE and CrowsonC (2021) survival: Survival analysis. https://cran.r-project.org/package=survival. R package version 3.2-13.
36.
YavuzAC and LambertP (2016) Semi-parametric frailty model for clustered interval-censored data. Statistical Modelling, 16, 360–391.
37.
ZhaoX, ZhaoQ, SunJ and KimJS. (2008) Generalized log-rank tests for partly interval-censored failure time data. Biometrical Journal, 50, 375–385.
38.
ZhouH and HansonT (2018a) A unified framework for fitting Bayesian semiparametric models to arbitrarily censored survival data, including spatially-referenced data. Journal of the American Statistical Association, 113, 571–581.
39.
ZhouH and HansonT (2018b) spBayesSurv: Bayesian modeling and analysis of spatially correlated survival data. https://cran.r-project.org/package=spBayesSurv. R package version 1.1.3.
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
