Abstract:Background Traditionally, meta‐analysis of time‐to‐event outcomes reports a single pooled hazard ratio assuming proportional hazards (PH). For health technology assessment evaluations, hazard ratios are frequently extrapolated across a lifetime horizon. However, when treatment effects vary over time, an assumption of PH is not always valid. The Royston‐Parmar (RP), piecewise exponential (PE), and fractional polynomial (FP) models can accommodate non‐PH and provide plausible extrapolations of survival curves beyond observed data. Methods Simulation study to assess and compare the performance of RP, PE, and FP models in a Bayesian framework estimating restricted mean survival time difference (RMSTD) at 50 years from a pairwise meta‐analysis with evidence of non‐PH. Individual patient data were generated from a mixture Weibull distribution. Twelve scenarios were considered varying the amount of follow‐up data, number of trials in a meta‐analysis, non‐PH interaction coefficient, and prior distributions. Performance was assessed through bias and mean squared error. Models were applied to a metastatic breast cancer example. Results FP models performed best when the non‐PH interaction coefficient was 0.2. RP models performed best in scenarios with complete follow‐up data. PE models performed well on average across all scenarios. In the metastatic breast cancer example, RMSTD at 50‐years ranged from −14.6 to 8.48 months. Conclusions Synthesis of time‐to‐event outcomes and estimation of RMSTD in the presence of non‐PH can be challenging and computationally intensive. Different approaches make different assumptions regarding extrapolation and sensitivity analyses varying key assumptions are essential to check the robustness of conclusions to different assumptions for the underlying survival function.

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

Nonlinear Multilevel Joint Model for Individual Lesion Kinetics and Survival to Characterize Intra-Individual Heterogeneity in Patients with Advanced Cancer

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

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

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

A Class of Semiparametric Yang and Prentice Frailty Models

PROPORTIONAL HAZARDS MODELS WITH CONTINUOUS MARKS

Impact of model misspecification in shared frailty survival models

A proportional hazards model for interval-censored data subject to instantaneous failures

Semi-parametric Regression Model for Survival Data: Graphical Visualization with R.

Individual frailty excess hazard models in cancer epidemiology

Practical implementation of frailty models in Mendelian risk prediction

Frailty Model with Change Point for Survival Analysis

Functional proportional hazards mixture cure model and its application to modelling the association between cancer mortality and physical activity in NHANES 2003-2006

A Proportional Hazards Model for the Subdistribution of a Competing Risk

Functional Proportional Hazards Mixture Cure Model with Applications in Cancer Mortality in NHANES and Post ICU Recovery.

Generalized accelerated failure time spatial frailty model for arbitrarily censored data

Analysis of survival data with non-proportional hazards: A comparison of propensity score weighted methods

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

Birnbaum–Saunders frailty regression models for clustered survival data

Bayesian pairwise meta‐analysis of time‐to‐event outcomes in the presence of non‐proportional hazards: A simulation study of flexible parametric, piecewise exponential and fractional polynomial models