Adaptive spline-based weighting functions for blended survival curve extrapolation

Health technology assessment frequently requires survival predictions well beyond the observed trial follow-up, yet single parametric models fitted to short horizons can accumulate long-term bias. Current guidance also encourages the principled use of external evidence, such as registry summaries or expert anchors, while maintaining fidelity to the trial data. We introduce an adaptive spline-weighted blended extrapolation on the cumulative-hazard scale that unifies these aims. The observation side is fitted on the trial window using a piecewise-exponential Cox-type model with a smooth prior-driven continuation beyond the administrative cutoff, implemented via INLA. The external side is an anchored Gompertz tail identified by a prespecified survival level at a clinically relevant time. To blend the two components, we compare their cumulative hazards, learn a monotone P-spline score over time, and pass it through a logistic link to obtain a data-driven, time-varying weight. A simple "temperature" scaling controls the slope of this weight and provides a practical guarantee of non-negative blended hazards on a chosen grid, preserving adaptivity while ensuring feasibility. Across Monte Carlo scenarios spanning multiple tail shapes and censoring levels, the method delivers consistently lower absolute survival error, smaller restricted mean survival time error, and improved stability compared with fixed-schedule blending and single-family parametric models. In a SEER registry study with three-year observation and ten-year extrapolation, the blended curve tracks the Kaplan-Meier estimates within follow-up and transitions smoothly toward the anchored tail, yielding small long-horizon errors across cancer sites and age strata. The framework is modular, interpretable, and easily extended to alternative tails and multiple anchors. An open-source implementation is available in the R package survblendr at https://github.com/haohaostats/survblendr.

Authors

Institutions

Publication Details

Journal
Statistical Methods in Medical Research
Published
2026-07-16
DOI
https://doi.org/10.1177/09622802261458072
Primary Topic
Health Systems, Economic Evaluations, Quality of Life
Type
article
Field-Weighted Citation Impact
0.00
Controls
|||
ALL TIME
JAN
FEB
MAR
APR
MAY
JUN
JUL
AUG
SEP
OCT
article

Adaptive spline-based weighting functions for blended survival curve extrapolation

Clara Grazian, Hao Chen
Statistical Methods in Medical Research
Health Systems, Economic Evaluations, Quality of Life
article

Adaptive spline-based weighting functions for blended survival curve extrapolation

Clara Grazian, Hao Chen
article en

Abstract

Health technology assessment frequently requires survival predictions well beyond the observed trial follow-up, yet single parametric models fitted to short horizons can accumulate long-term bias. Current guidance also encourages the principled use of external evidence, such as registry summaries or expert anchors, while maintaining fidelity to the trial data. We introduce an adaptive spline-weighted blended extrapolation on the cumulative-hazard scale that unifies these aims. The observation side is fitted on the trial window using a piecewise-exponential Cox-type model with a smooth prior-driven continuation beyond the administrative cutoff, implemented via INLA. The external side is an anchored Gompertz tail identified by a prespecified survival level at a clinically relevant time. To blend the two components, we compare their cumulative hazards, learn a monotone P-spline score over time, and pass it through a logistic link to obtain a data-driven, time-varying weight. A simple "temperature" scaling controls the slope of this weight and provides a practical guarantee of non-negative blended hazards on a chosen grid, preserving adaptivity while ensuring feasibility. Across Monte Carlo scenarios spanning multiple tail shapes and censoring levels, the method delivers consistently lower absolute survival error, smaller restricted mean survival time error, and improved stability compared with fixed-schedule blending and single-family parametric models. In a SEER registry study with three-year observation and ten-year extrapolation, the blended curve tracks the Kaplan-Meier estimates within follow-up and transitions smoothly toward the anchored tail, yielding small long-horizon errors across cancer sites and age strata. The framework is modular, interpretable, and easily extended to alternative tails and multiple anchors. An open-source implementation is available in the R package survblendr at https://github.com/haohaostats/survblendr.

Statistical Methods in Medical Research
The University of Sydney (AU), NSW Department of Education (AU)
Openalex Percentile: Top 5%
Health Systems, Economic Evaluations, Quality of Life
AI Navigator

Ask Laika to Summarize, Analyze, and Connect papers live on the map.

Summarize Papers & Methodologies

Extract key findings, datasets, and comparative methods across publications.

Benchmark Rankings & Visual Analytics

Rank top research institutions, authors, funders, topics, and journals by Field-Weighted Citation Impact (FWCI) and paper volume with instant charts.

Connect Distant Disciplines

Bridge topological clusters on the map to find hidden collaborative intersections.