Goodness-of-fit criteria for survival data

The definition of an appropriate measure for goodness-of-fit in case of survival data comparable to R^2 in linear regression is difficult due to censored observations. In this paper, a variety of answers based on different residuals and variance of survival curves are presented together with a newly introduced criterion. In univariate simulation studies, the presented criteria are examined with respect to their dependence on the value of the coefficient associated with the covariate; underlying covariate distribution and censoring percentage in the data. Investigation of the relations between the values of the different criteria indicates strong dependencies, although the absolute values show high discrepancies and the criteria building processes differ substantially

Using weibull mixture distributions to model heterogeneous survival data

In this article we use Bayesian methods to fit a Weibull mixture model with an unknown number of components to possibly right censored survival data. This is done using the recently developed, birth-death MCMC algorithm. We also show how to estimate the survivor function and the expected hazard rate from the MCMA output

Bayesian Regularisation in Structured Additive Regression Models for Survival Data

During recent years, penalized likelihood approaches have attracted a lot of interest both in the area of semiparametric regression and for the regularization of high-dimensional regression models. In this paper, we introduce a Bayesian formulation that allows to combine both aspects into a joint regression model with a focus on hazard regression for survival times. While Bayesian penalized splines form the basis for estimating nonparametric and flexible time-varying effects, regularization of high-dimensional covariate vectors is based on scale mixture of normals priors. This class of priors allows to keep a (conditional) Gaussian prior for regression coefficients on the predictor stage of the model but introduces suitable mixture distributions for the Gaussian variance to achieve regularization. This scale mixture property allows to device general and adaptive Markov chain Monte Carlo simulation algorithms for fitting a variety of hazard regression models. In particular, unifying algorithms based on iteratively weighted least squares proposals can be employed both for regularization and penalized semiparametric function estimation. Since sampling based estimates do no longer have the variable selection property well-known for the Lasso in frequentist analyses, we additionally consider spike and slab priors that introduce a further mixing stage that allows to separate between influential and redundant parameters. We demonstrate the different shrinkage properties with three simulation settings and apply the methods to the PBC Liver dataset