现象学增长模型在流行病预测中的结构性与实用性可识别性研究
Structural and Practical Identifiability of Phenomenological Growth Models for Epidemic Forecasting
摘要 Abstract
现象学模型是利用现实世界数据预测疾病传播动态的强大工具,尤其是在对疾病机制了解有限的情况下。然而,这些模型的可靠性依赖于参数的结构性与实用性可识别性。本研究系统分析了流行病学中常用的六种增长模型的可识别性:广义增长模型、广义Logistic模型、Richards模型、广义Richards模型、Gompertz模型以及具有非均匀混合的改进SEIR模型。为解决这些模型中非整数幂指数带来的挑战,我们通过引入额外的状态变量重新构建模型,并利用JULIA语言中的StructuralIdentifiability.jl软件包进行严格的结构性可识别性分析。我们通过GrowthPredict MATLAB工具箱对参数估计和预测进行验证,该工具箱基于现象学增长模型拟合和预测时间序列轨迹。我们将其应用于三种流行病学数据集:猴痘、COVID-19和埃博拉的每周发病率数据。此外,我们通过蒙特卡洛模拟评估了实际可识别性,以检验不同观测噪声水平下参数估计的稳健性。结果表明,在所提出的重构下,这六种模型均具有结构性可识别性。进一步的实际可识别性分析显示,尽管不同模型和数据集对噪声的敏感性有所差异,但参数估计在各种噪声水平下保持稳健。这些发现为现象学模型刻画流行病传播轨迹提供了重要见解,强调了其适应现实世界挑战的能力及其在公共卫生干预中的作用。
Phenomenological models are highly effective tools for forecasting disease dynamics using real world data, particularly in scenarios where detailed knowledge of disease mechanisms is limited. However, their reliability depends on the model parameters' structural and practical identifiability. In this study, we systematically analyze the identifiability of six commonly used growth models in epidemiology:the generalized growth model, the generalized logistic model, the Richards model, the generalized Richards model, the Gompertz model, and a modified SEIR model with inhomogeneous mixing. To address challenges posed by non-integer power exponents in these models, we reformulate them by introducing additional state variables. This enables rigorous structural identifiability analysis using the StructuralIdentifiability.jl package in JULIA. We validate the structural identifiability results by performing parameter estimation and forecasting using the GrowthPredict MATLAB toolbox. This toolbox is designed to fit and forecast time series trajectories based on phenomenological growth models. We applied it to three epidemiological datasets: weekly incidence data for monkeypox, COVID 19, and Ebola. Additionally, we assess practical identifiability through Monte Carlo simulations to evaluate parameter estimation robustness under varying levels of observational noise. Our results confirm that all six models are structurally identifiable under the proposed reformulation. Furthermore, practical identifiability analyses demonstrate that parameter estimates remain robust across different noise levels, though sensitivity varies by model and dataset. These findings provide critical insights into the strengths and limitations of phenomenological models to characterize epidemic trajectories, emphasizing their adaptability to real world challenges and their role in informing public health interventions.