Abstract:Chronic disease microsimulation models often simulate disease incidence as a function of risk factors that evolve over time (e.g., blood pressure increasing with age) in order to facilitate decision analyses of different disease screening and prevention strategies. Existing models typically rely on incidence rates estimated with standard survival analysis techniques (e.g., proportional hazards from baseline data) that are not designed to be continually updated each model cycle. We introduce the use of joint longitudinal and time-to-event to parameterize microsimulations to avoid potential issues from using these existing methods. These joint models include random effects regressions to estimate the risk factor trajectories and a survival model to predict disease risk based on those estimated trajectories. In a case study on cardiovascular disease (CVD), we compare the validity of microsimulation models parameterized with this joint model approach to those parameterized with the standard approaches. A CVD microsimulation model was constructed that modeled the trajectory of seven CVD risk factors/predictors as a function of age (smoking, diabetes, systolic blood pressure, antihypertensive medication use, total cholesterol, HDL, and statin use) and predicted yearly CVD incidence as a function of these predictors, plus age, sex, and race. We parameterized the model using data from the Atherosclerosis Risk in Communities study (ARIC). The risk of CVD in the microsimulation was parameterized with three approaches: (1) joint longitudinal and time-to-event model, (2) proportional hazards model estimated using baseline data, and (3) proportional hazards model estimated using time-varying data. We accounted for non-CVD mortality across all the parameterization approaches. We simulated risk factor trajectories and CVD incidence from age 70y to 85y for an external test set comprised of individuals from the Multi-Ethnic Study of Atherosclerosis (MESA). We compared the simulated to observed incidence using both average survival curves and the E50 and E90 calibration metrics (the median and 90th percentile absolute difference between observed and predicted incidence) to measure the validity of each parameterization approach. The average CVD survival curve estimated by the microsimulation model parameterized with the joint model approach matched the observed curve from the test set relatively closely. The other parameterization methods generally performed worse, especially the proportional hazards model estimated using baseline data. Similar results were observed for the calibration metrics, with the joint model performing particularly well on the E90 metric compared to the other models. Using a joint longitudinal and time-to-event model to parameterize a CVD simulation model produced incidence predictions that more accurately reflected observed data than a model parameterized with standard approaches. This parameterization approach could lead to more reliable microsimulation models, especially for models that evaluate policies which depend on tracking dynamic risk factors over time. Beyond this single case study, more work is needed to identify the specific circumstances where the joint model approach will outperform existing methods.

Using Age-Specific Rates for Parametric Survival Function Estimation in Simulation Models

Parmsurv: a SAS Macro for Flexible Parametric Survival Analysis with Long-Term Predictions

A Fast Nonparametric Sampling (NPS) Method for Time-to-Event in Individual-Level Simulation Models.

Modelling multiple time-scales with flexible parametric survival models

Parametric Regression Model for Survival Data: Weibull Regression Model As an Example

The Impact of the Approach to Accounting for Age and Sex in Economic Models on Predicted Quality-Adjusted Life-Years

Nonlinear multiparametric modeling of life-table data with adaptive distributions: time evolution of hazard ratios

A flexible parametric modelling framework for survival analysis

Impact of limited sample size and follow-up on single event survival extrapolation for health technology assessment: a simulation study

Bayesian parametric models for survival prediction in medical applications

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

Beyond standardized mortality ratios; some uses of smoothed age-specific mortality rates on small areas studies

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

The use of restricted mean survival time to estimate treatment effect under model misspecification, a simulation study

Probabilistic mortality forecasting with varying age-specific survival improvements

Comparing Survival Extrapolation within All-Cause and Relative Survival Frameworks by Standard Parametric Models and Flexible Parametric Spline Models Using the Swedish Cancer Registry

Grouped Functional Time Series Forecasting: An Application to Age-Specific Mortality Rates

Using Joint Longitudinal and Time-to-Event Models to Improve the Parameterization of Chronic Disease Microsimulation Models: an Application to Cardiovascular Disease

paramix : An R package for parameter discretisation in compartmental models, with application to calculating years of life lost

Impact of model misspecification in shared frailty survival models

SurvInt: a simple tool to obtain precise parametric survival extrapolations