Abstract:Survival models incorporating random effects to account for unmeasured heterogeneity are being increasingly used in biostatistical and applied research. Specifically, unmeasured covariates whose lack of inclusion in the model would lead to biased, inefficient results are commonly modelled by including a subject-specific (or cluster-specific) frailty term that follows a given distribution (e.g. Gamma or log-Normal). Despite that, in the context of parametric frailty models little is known about the impact of misspecifying the baseline hazard, the frailty distribution, or both. Therefore, our aim is to quantify the impact of such misspecification in a wide variety of clinically plausible scenarios via Monte Carlo simulation, using open source software readily available to applied researchers. We generate clustered survival data assuming various baseline hazard functions, including mixture distributions with turning points, and assess the impact of sample size, variance of the frailty, baseline hazard function, and frailty distribution. Models compared include standard parametric distributions and more flexible spline-based approaches; we also included semiparametric Cox models. The resulting bias can be clinically relevant. In conclusion, we highlight the importance of fitting models that are flexible enough and the importance of assessing model fit. We illustrate our conclusions with two applications using data on diabetic retinopathy and bladder cancer. Our results show the importance of assessing model fit with respect to the baseline hazard function and the distribution of the frailty: misspecifying the former leads to biased relative and absolute risk estimates while misspecifying the latter affects absolute risk estimates and measures of heterogeneity.

Fitting conditional survival models to meta-analytic data by using a transformation toward mixed-effects models

Fitting phase-type frailty models

Survival Analysis Using a Scale Change Random Effects Model

Frailty Model with Change Point for Survival Analysis

Bayesian transformation cure frailty models with multivariate failure time data

Multilevel mixed effects parametric survival models using adaptive Gauss-Hermite quadrature with application to recurrent events and individual participant data meta-analysis

Modeling county level breast cancer survival data using a covariate-adjusted frailty proportional hazards model

A Transformation Perspective on Marginal and Conditional Models

Birnbaum–Saunders frailty regression models for clustered survival data

Shared Frailty Methods for Complex Survival Data: A Review of Recent Advances

Impact of model misspecification in shared frailty survival models

Non-proportional hazards model with a PVF frailty term: application with a melanoma dataset

MODELLING COUNTY LEVEL BREAST CANCER SURVIVAL DATA USING A COVARIATE-ADJUSTED FRAILTY PROPORTIONAL HAZARDS MODEL

A framework for meta-analysis of prediction model studies with binary and time-to-event outcomes

Frailty model with change points for survival analysis

On the Addams family of discrete frailty distributions for modeling multivariate case I interval-censored data

On the Addams family of discrete frailty distributions for modelling multivariate case I interval-censored data

Investigating hospital heterogeneity with a competing risks frailty model

Predicting Survival Probabilities with Semiparametric Transformation Models

Analysis of clustered interval‐censored data using a class of semiparametric partly linear frailty transformation models

An alternative estimation method for the accelerated failure time frailty model