Posterior likelihood methods for multivariate survival data.

This article deals with the semiparametric analysis of multivariate survival data with random block (group) effects. Survival times within the same group are correlated as a consequence of a frailty random block effect. The standard approaches assume either a parametric or a completely unknown baseline hazard function. This paper considers an intermediate solution, that is, a nonparametric function that is reasonably smooth. This is accomplished by a Bayesian model in which the conditional proportional hazards model is used with a correlated prior process for the baseline hazard. The posterior likelihood based on data, as well as the prior process, is similar to the discretized penalized likelihood for the frailty model. The methodology is exemplified with the recurrent kidney infections data of McGilchrist and Aisbett (1991, Biometrics 47, 461-466), in which the times to infections within the same patients are expected to be correlated. The reanalysis of the data has shown that the estimates of the parameters of interest and the associated standard errors depend on the prior knowledge about the smoothness of the baseline hazard.

[1]  D. Cox,et al.  Analysis of Survival Data. , 1985 .

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

[3]  Nozer D. Singpurwalla,et al.  Understanding the Kalman Filter , 1983 .

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

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

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

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

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

[9]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[10]  I. Good.  Nonparametric roughness penalties for probability densities , 1971 .

[11]  Debajyoti Sinha,et al.  Semiparametric Bayesian Analysis of Multiple Event Time Data , 1993 .

[12]  P. Hougaard.  A class of multivanate failure time distributions , 1986 .

[13]  William Feller,et al.  An Introduction to Probability Theory and Its Applications , 1951 .

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

[15]  Ludwig Fahrmeir,et al.  Dynamic modelling and penalized likelihood estimation for discrete time survival data , 1994 .

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

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

[18]  K S Chan,et al.  A dynamic frailty model for multivariate survival data. , 1997, Biometrics.

[19]  M. A. Tanner,et al.  Tools for Statistical Inference: Methods for the Exploration of Posterior Distributions and Likelihood Functions, 3rd Edition , 1998 .

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

[21]  J. Anderson,et al.  Penalized maximum likelihood estimation in logistic regression and discrimination , 1982 .

[22]  C. Mcgilchrist,et al.  ML and REML estimation in survival analysis with time dependent correlated frailty. , 1998, Statistics in medicine.

[23]  G. Wahba.  Bayesian "Confidence Intervals" for the Cross-validated Smoothing Spline , 1983 .

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

[25]  D. Cox.  Regression Models and Life-Tables , 1972 .

[26]  D. W. Scott,et al.  Nonparametric Probability Density Estimation by Discrete Maximum Penalized- Likelihood Criteria , 1980 .