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.

Generalizing the Frailty Assumptions in Survival Analysis

Birnbaum–Saunders frailty regression models for clustered survival data

On the Maximum Penalized Full Likelihood Approach for Cox Model with Extreme Value for Heavily Censored Survival Data

A study on comparisons of additive regression frailty models to counter heterogeneity: Bayesian strategies and case study

Frailty Models for Pneumonia to Death with a Left‐censored Covariate

A simulation study of the estimation quality in the double-Cox model with shared frailty for non-proportional hazards survival analysis

Frailty Model with Change Point for Survival Analysis

CASE-CONTROL SURVIVAL ANALYSIS WITH A GENERAL SEMIPARAMETRIC SHARED FRAILTY MODEL - A PSEUDO FULL LIKELIHOOD APPROACH

A Flexible Adaptive Lasso Cox Frailty Model Based on the Full Likelihood

Generalised Lindley shared additive frailty regression model for bivariate survival data

General Semiparametric Shared Frailty Model: Estimation and Simulation with frailtySurv

Conditional quasi‐likelihood inference for mean residual life regression with clustered failure time data

Estimation of the accelerated failure time frailty model under generalized gamma frailty

Fitting phase-type frailty models

Graphical Tests for the Assumption of Gamma and Inverse Gaussian Frailty Distributions

A Multi-parameter regression model for interval censored survival data

Regression Analysis of Multivariate Incomplete Failure Time Data by Modeling Marginal Distributions

Shared inverse Gaussian frailty models based on additive hazards

Long-Term Dagum-PVF Frailty Regression Model: Application in Health Studies

A General Semiparametric Hazards Regression Model: Efficient Estimation and Structure Selection

Impact of model misspecification in shared frailty survival models