[论文解读] Mathematical Modeling of Epidemic Diseases; A Case Study of the COVID-19 Coronavirus

Reza Sameni|arXiv (Cornell University)|Mar 25, 2020

COVID-19 epidemiological studies被引用 80

一句话总结

这篇论文将 SEIR 型分区模型扩展用于研究 COVID-19 的传播，展示社交干预如何改变模型参数和结果，MATLAB 代码在线提供。

ABSTRACT

In this research, we study the propagation patterns of epidemic diseases such as the COVID-19 coronavirus, from a mathematical modeling perspective. The study is based on an extensions of the well-known susceptible-infected-recovered (SIR) family of compartmental models. It is shown how social measures such as distancing, regional lockdowns, quarantine and global public health vigilance, influence the model parameters, which can eventually change the mortality rates and active contaminated cases over time, in the real world. As with all mathematical models, the predictive ability of the model is limited by the accuracy of the available data and to the so-called extit{level of abstraction} used for modeling the problem. In order to provide the broader audience of researchers a better understanding of spreading patterns of epidemic diseases, a short introduction on biological systems modeling is also presented and the Matlab source codes for the simulations are provided online.

研究动机与目标

• 动机：利用数学建模来理解疫情传播并为医疗保健规划提供信息。
• 引入一个扩展的分区框架，纳入暴露（Exposed）和致死分区用于 COVID-19 动态。
• 解释非药物干预如何影响模型参数和疫情传播轨迹。
• 强调建模假设、局限性，以及数据质量在参数估计中的作用。

提出的方法

• 将 SEIR 框架扩展为包含暴露和致死分区，形成一个五分区模型（S, E, I, R, P）。
• 推导并给出控制各分区随时间演化的模型方程（Eq. 27）。
• 描述闭系统与测量考虑因素，包括状态方程和测量方程以及可选的降阶（Eq. 28）。
• 讨论分区模型中的繁殖概念与稳定性分析，包括 R0 和一个替代繁殖数 7 0 0 ~（命题 1）。
• 提供关于数理流行病学、分区建模以及对非线性项如 s(t)e(t) 和 s(t)i(t) 的解释的简要入门。
• 请注意，用于复现实验结果的 MATLAB 源代码可在项目仓库的在线版本中获取。

实验结果

研究问题

• RQ1 将 SIR/SEIR 框架扩展如何捕捉包含暴露和致死结果的 COVID-19 进程？
• RQ2 社会措施（社交距离、封锁、警觉性）如何影响模型参数及随时间变化的死亡率和在院/活跃病例？
• RQ3 在爆发稳定性方面，经典的 R0（基于 NGM）与提出的替代繁殖数之间有什么关系？
• RQ4 可观测量有哪些，以及在测量噪声存在的情况下如何进行状态估计？

主要发现

• 提出了一个扩展的 SEIR 型模型，包含暴露（E）和致死（P）分区，以更好地反映 COVID-19 动态。
• 模型显示在不同干预下，传播速率和其他参数如何影响死亡率和活跃病例随时间的变化。
• 提供了一个测量模型，将隐藏分区与可观测量（如感染、治愈和死亡）联系起来。
• 关于繁殖数的讨论强调其与爆发稳定性标准的一致性，并引入一个考虑离散化时间的替代 R0。
• 工作强调由于数据质量和抽象程度，数学模型的局限性，并提供 MATLAB 代码以便复现。