Bayesian analysis of multivariate survival data using Monte Carlo methods

This paper deals with the analysis of multivariate survival data from a Bayesian perspective using Markov‐chain Monte Carlo methods. The Metropolis along with the Gibbs algorithm is used to calculate some of the marginal posterior distributions. A multivariate survival model is proposed, since survival times within the same group are correlated as a consequence of a frailty random block effect. The conditional proportional‐hazards model of Clayton and Cuzick is used with a martingale structured prior process (Arjas and Gasbarra) for the discretized baseline hazard. Besides the calculation of the marginal posterior distributions of the parameters of interest, this paper presents some Bayesian EDA diagnostic techniques to detect model adequacy. The methodology is exemplified with kidney infection data where the times to infections within the same patients are expected to be correlated.

[1]  D Sinha,et al.  Posterior likelihood methods for multivariate survival data. , 1998, Biometrics.

[2]  Lynn Kuo,et al.  A Bayesian predictive approach to determining the number of components in a mixture distribution , 1995 .

[3]  A. Gelfand,et al.  Bayesian analysis of proportional hazards models built from monotone functions. , 1995, Biometrics.

[4]  Thomas A. Louis,et al.  Assessing gamma frailty models for clustered failure time data , 1995, Lifetime data analysis.

[5]  A. Gelfand,et al.  Bayesian Model Choice: Asymptotics and Exact Calculations , 1994 .

[6]  Bayesian Analysis of Semiparametric Proportional Hazards Models , 1994 .

[7]  Martin Abba Tanner,et al.  Tools for Statistical Inference: Observed Data and Data Augmentation Methods , 1993 .

[8]  McGilchrist Ca REML estimation for survival models with frailty. , 1993 .

[9]  C A McGilchrist REML estimation for survival models with frailty. , 1993, Biometrics.

[10]  Hong Chang,et al.  Model Determination Using Predictive Distributions with Implementation via Sampling-Based Methods , 1992 .

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

[12]  J P Klein,et al.  Semiparametric estimation of random effects using the Cox model based on the EM algorithm. , 1992, Biometrics.

[13]  David Oakes,et al.  Frailty Models For Multiple Event Times , 1992 .

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

[15]  D G Clayton,et al.  A Monte Carlo method for Bayesian inference in frailty models. , 1991, Biometrics.

[16]  D. Gamerman Dynamic Bayesian Models for Survival Data , 1991 .

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

[18]  N. Hjort Nonparametric Bayes Estimators Based on Beta Processes in Models for Life History Data , 1990 .

[19]  Adrian F. M. Smith,et al.  Sampling-Based Approaches to Calculating Marginal Densities , 1990 .

[20]  P. Grambsch,et al.  Martingale-based residuals for survival models , 1990 .

[21]  D. Oakes,et al.  Bivariate survival models induced by frailties , 1989 .

[22]  D. Oakes,et al.  Semiparametric inference in a model for association in bivanate survival data , 1986 .

[23]  D. Clayton,et al.  Multivariate generalizations of the proportional hazards model , 1985 .

[24]  David A. Schoenfeld,et al.  Partial residuals for the proportional hazards regression model , 1982 .

[25]  J. B Urridge,et al.  Empirical Bayes Analysis of Survival Time Data , 1981 .

[26]  J. Kalbfleisch Non‐Parametric Bayesian Analysis of Survival Time Data , 1978 .

[27]  Tom Leonard Density Estimation, Stochastic Processes and Prior Information , 1978 .

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

[29]  Tom Leonard,et al.  A Bayesian method for histograms , 1973 .

[30]  J. Kalbfleisch,et al.  Marginal likelihoods based on Cox's regression and life model , 1973 .

[31]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

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