Bayesian survival analysis with BUGS

Survival analysis is one of the most important fields of statistics in medicine and the biological sciences. In addition, the computational advances in the last decades have favoured the use of Bayesian methods in this context, providing a flexible and powerful alternative to the traditional frequentist approach. The objective of this paper is to summarise some of the most popular Bayesian survival models, such as accelerated failure time, proportional hazards, mixture cure, competing risks, frailty, and joint models of longitudinal and survival data. Moreover, an implementation of each presented model is provided using a BUGS syntax that can be run with JAGS from the R programming language. Reference to other Bayesian R-packages are also discussed.

[1]  D K Dey,et al.  A Weibull Regression Model with Gamma Frailties for Multivariate Survival Data , 1997, Lifetime data analysis.

[2]  Carmen Armero,et al.  Bayesian regularization for flexible baseline hazard functions in Cox survival models , 2020, Biometrical journal. Biometrische Zeitschrift.

[3]  D. Bates,et al.  Output Analysis and Diagnostics for MCMC , 2015 .

[4]  D. Thomas,et al.  Simultaneously modelling censored survival data and repeatedly measured covariates: a Gibbs sampling approach. , 1996, Statistics in medicine.

[5]  Cécile Proust-Lima,et al.  Joint latent class models for longitudinal and time-to-event data: A review , 2014, Statistical methods in medical research.

[6]  James R. Kenyon,et al.  Analysis of Multivariate Survival Data , 2002, Technometrics.

[7]  Geoff Jones,et al.  Book Review: Joint Models for Longitudinal and Time-to-Event Data: With Applications in R, by Dimitris Rizopoulos , 2013 .

[8]  Jiqiang Guo,et al.  Stan: A Probabilistic Programming Language. , 2017, Journal of statistical software.

[9]  W. D. Ray.  4. Modelling Survival Data in Medical Research , 1995 .

[10]  Joseph Berkson,et al.  Survival Curve for Cancer Patients Following Treatment , 1952 .

[11]  H Putter,et al.  Tutorial in biostatistics: competing risks and multi‐state models , 2007, Statistics in medicine.

[12]  Paul J Rathouz,et al.  Bayesian approach for flexible modeling of semicompeting risks data. , 2014, Statistics in medicine.

[13]  C. Armero,et al.  Bayesian joint modeling of bivariate longitudinal and competing risks data: An application to study patient‐ventilator asynchronies in critical care patients , 2017, Biometrical journal. Biometrische Zeitschrift.

[14]  Rizopoulos Dimitris,et al.  Joint Modeling of Longitudinal and Time-to-Event Data , 2014 .

[15]  Carmen Armero,et al.  Approximate Bayesian inference for mixture cure models , 2018, TEST.

[16]  Niels Keiding,et al.  Statistical Models Based on Counting Processes , 1993 .

[17]  Xiaoyan Lin,et al.  A Bayesian proportional hazards model for general interval-censored data , 2015, Lifetime data analysis.

[18]  Yongseok Park,et al.  Real‐Time Individual Predictions of Prostate Cancer Recurrence Using Joint Models , 2013, Biometrics.

[19]  C. Robert,et al.  Bayesian Modeling Using WinBUGS , 2009 .

[20]  D. Kleinbaum,et al.  Survival Analysis: A Self-Learning Text. , 1996 .

[21]  V. De Gruttola,et al.  Modelling progression of CD4-lymphocyte count and its relationship to survival time. , 1994, Biometrics.

[22]  C. Mcgilchrist,et al.  Regression with frailty in survival analysis. , 1991, Biometrics.

[23]  I. Langner.  Survival Analysis: Techniques for Censored and Truncated Data , 2006 .

[24]  David J. Lunn,et al.  The BUGS Book: A Practical Introduction to Bayesian Analysis , 2013 .

[25]  H. Rue,et al.  Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations , 2009 .

[26]  Maja Pohar Perme,et al.  Inference for outcome probabilities in multi-state models , 2008, Lifetime data analysis.

[27]  R Henderson,et al.  Joint modelling of longitudinal measurements and event time data. , 2000, Biostatistics.

[28]  O. Aalen.  Effects of frailty in survival analysis , 1994, Statistical methods in medical research.

[29]  H. Rüden,et al.  Surveillance of nosocomial sepsis and pneumonia in patients with a bone marrow or peripheral blood stem cell transplant: a multicenter project. , 2005, Clinical infectious diseases : an official publication of the Infectious Diseases Society of America.

[30]  E L Hamaker,et al.  A Comparison of Inverse-Wishart Prior Specifications for Covariance Matrices in Multilevel Autoregressive Models , 2016, Multivariate behavioral research.

[31]  M. Wulfsohn,et al.  Modeling the Relationship of Survival to Longitudinal Data Measured with Error. Applications to Survival and CD4 Counts in Patients with AIDS , 1995 .

[32]  D. Rubin,et al.  Inference from Iterative Simulation Using Multiple Sequences , 1992 .

[33]  Kyu Ha Lee,et al.  SemiCompRisks: An R Package for the Analysis of Independent and Cluster-correlated Semi-competing Risks Data , 2019, R J..

[34]  D.,et al.  Regression Models and Life-Tables , 2022 .

[35]  Walter R. Gilks,et al.  A Language and Program for Complex Bayesian Modelling , 1994 .

[36]  McGilchrist Ca,et al.  Regression with frailty in survival analysis. , 1991 .

[37]  J. Ware,et al.  Random-effects models for longitudinal data. , 1982, Biometrics.

[38]  Leo Breiman,et al.  Bagging Predictors , 1996, Machine Learning.

[39]  G. Barrie Wetherill,et al.  Random Effects Models , 1981 .

[40]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[41]  N M Laird,et al.  Increasing efficiency from censored survival data by using random effects to model longitudinal covariates , 1998, Statistical methods in medical research.

[42]  Adrian F. M. Smith,et al.  Bayesian Inference for Generalized Linear and Proportional Hazards Models Via Gibbs Sampling , 1993 .

[43]  Ronald Christensen,et al.  Bayesian Ideas and Data Analysis: An Introduction for Scientists and Statisticians , 2010 .

[44]  Martin Schumacher,et al.  A competing risks analysis of bloodstream infection after stem‐cell transplantation using subdistribution hazards and cause‐specific hazards , 2007, Statistics in medicine.

[45]  Kyu Ha Lee,et al.  Hierarchical Models for Semicompeting Risks Data With Application to Quality of End-of-Life Care for Pancreatic Cancer , 2015, Journal of the American Statistical Association.

[46]  N M Laird,et al.  Mixture models for the joint distribution of repeated measures and event times. , 1997, Statistics in medicine.

[47]  Martyn Plummer,et al.  JAGS: A program for analysis of Bayesian graphical models using Gibbs sampling , 2003 .

[48]  Laurence L. George,et al.  The Statistical Analysis of Failure Time Data , 2003, Technometrics.

[49]  J G Ibrahim,et al.  On Bayesian Inference for Proportional Hazards Models Using Noninformative Priors , 2000, Lifetime data analysis.

[50]  C. Armero.  Bayesian Joint Models for Longitudinal and Survival Data , 2020, 2005.12822.

[51]  Eric R. Ziegel,et al.  The Elements of Statistical Learning , 2003, Technometrics.

[52]  J. Ibrahim,et al.  Handbook of survival analysis , 2014 .

[53]  Dimitris Rizopoulos,et al.  Personalized screening intervals for biomarkers using joint models for longitudinal and survival data. , 2015, Biostatistics.

[54]  Riten Mitra,et al.  Nonparametric Bayesian inference in biostatistics , 2015 .

[55]  N. Breslow.  Covariance analysis of censored survival data. , 1974, Biometrics.

[56]  C Armero,et al.  Bayesian joint ordinal and survival modeling for breast cancer risk assessment , 2016, Statistics in medicine.

[57]  D. Clayton.  A model for association in bivariate life tables and its application in epidemiological studies of familial tendency in chronic disease incidence , 1978 .

[58]  Bradley P Carlin,et al.  Flexible Bayesian survival modeling with semiparametric time-dependent and shape-restricted covariate effects. , 2016, Bayesian analysis.

[59]  Miaomiao Ge,et al.  Bayesian inference of the fully specified subdistribution model for survival data with competing risks , 2012, Lifetime data analysis.

[60]  M. Wulfsohn,et al.  A joint model for survival and longitudinal data measured with error. , 1997, Biometrics.

[61]  Dimitris Rizopoulos,et al.  Joint Models for Longitudinal and Time-to-Event Data: With Applications in R , 2012 .

[62]  N. Keiding,et al.  Multi-state models for event history analysis , 2002, Statistical methods in medical research.

[63]  O. Kardaun,et al.  Statistical survival analysis of male larynx‐cancer patients ‐ a case study , 1983 .

[64]  John K Kruschke,et al.  Bayesian data analysis. , 2010, Wiley interdisciplinary reviews. Cognitive science.

[65]  J. Crowley,et al.  Covariance Analysis of Heart Transplant Survival Data , 1977 .

[66]  Arnošt Komárek,et al.  Survival Analysis with Interval-Censored Data: A Practical Approach with Examples in R, SAS, and BUGS , 2017 .

[67]  Patrick Royston,et al.  Estimating a smooth baseline hazard function for the Cox model , 2012 .

[68]  Anabel Forte,et al.  Bayesian joint modeling for assessing the progression of chronic kidney disease in children , 2018, Statistical methods in medical research.

[69]  A I Goldman,et al.  Comparison of autologous and allogeneic bone marrow transplantation for treatment of high-risk refractory acute lymphoblastic leukemia. , 1987, The New England journal of medicine.

[70]  J. Vaupel,et al.  The impact of heterogeneity in individual frailty on the dynamics of mortality , 1979, Demography.

[71]  Andrew Gelman,et al.  General methods for monitoring convergence of iterative simulations , 1998 .

[72]  Raymond J. Carroll,et al.  Estimation and comparison of changes in the presence of informative right censoring by modeling the censoring process , 1988 .

[73]  H. Borchers.  Practical Numerical Math Functions [R package pracma version 2.2.9] , 2019 .