[论文解读] Epidemic analysis of COVID-19 in China by dynamical modeling

Liangrong Peng, Wuyue Yang+3 more|arXiv (Cornell University)|Feb 16, 2020

COVID-19 epidemiological studies参考文献 37被引用 103

一句话总结

该论文将 SEIR 推广为具备 quarantined 状态和自我防护的广义 SEIR 模型，基于早期中国数据估计关键参数，并预测拐点、结束时间以及大陆中国、湖北省和武汉的总感染人数。

ABSTRACT

The outbreak of novel coronavirus-caused pneumonia (COVID-19) in Wuhan has attracted worldwide attention. Here, we propose a generalized SEIR model to analyze this epidemic. Based on the public data of National Health Commission of China from Jan. 20th to Feb. 9th, 2020, we reliably estimate key epidemic parameters and make predictions on the inflection point and possible ending time for 5 different regions. According to optimistic estimation, the epidemics in Beijing and Shanghai will end soon within two weeks, while for most part of China, including the majority of cities in Hubei province, the success of anti-epidemic will be no later than the middle of March. The situation in Wuhan is still very severe, at least based on public data until Feb. 15th. We expect it will end up at the beginning of April. Moreover, by inverse inference, we find the outbreak of COVID-19 in Mainland, Hubei province and Wuhan all can be dated back to the end of December 2019, and the doubling time is around two days at the early stage.

研究动机与目标

• 激发对中国早期疫情爆发期 COVID-19 动态的理解。
• 开发一个考虑自我防护和隔离的广义 SEIR 框架，以分析疫情传播。
• 使用官方数据估计关键参数，并用它们来预测中国各地区的拐点、结束时间和总感染人数。
• 通过逆推推断研究疫情起始日期，以了解早期动态。

提出的方法

• 将 SEIR 模型推广为包含七个状态：S, P, E, I, Q, R, D。
• 引入随时间变化的治愈率和死亡率，以及一个降低有效接触的防护（自我防护）速率。
• 使用带模拟退火的最小二乘回归拟合未知参数和初始条件到日度隔离数据 Q(t)。
• 将潜伏期 gamma^{-1} 固定在合理范围内并探索其对其他参数的影响，以避免过拟合。
• 进行敏感性分析，以评估 β（感染率）和 δ^{-1}（隔离时间）如何影响最终暴发规模和拐点。
• 采用逆推推断（射击法）从早期轨迹估计疫情起始日期。

实验结果

研究问题

• RQ1 最能描述中国早期 COVID-19 动态的关键流行病学参数有哪些（潜伏时间、隔离时间、防护率、感染率）？
• RQ2 自我防护措施和隔离如何影响有效再生数以及疫情拐点的时机？
• RQ3 模型能否预测 Mainland China、Hubei、Wuhan 的拐点、结束时间和总感染人数，并且这些预测与观测数据相比如何？
• RQ4 根据逆推推断，COVID-19 在 Mainland China、Hubei 和 Wuhan 的大致起始日期是什么时候？

主要发现

• 带隔离状态的广义 SEIR 模型能够拟合官方的隔离数据，并捕捉不同地区在防护和隔离动态方面的差异。
• 预测表明北京和上海将很快结束（在 2020-02-16 前后两周内），大陆大部分地区将于三月中旬结束，武汉大约在四月初结束（基于截至 2020-02-16 的公开数据）。
• 逆推推断显示疫情在2020-01-20 前20–25天开始，即大约2019年12月底，早期呈指数增长（翻倍时间约2天）。
• 感染率 β 近似1，表示高度传染性，模型参数显示在更发达地区隔离时间更长、在更严格地区更短。
• 当考虑治愈率和死亡率随时间变化时，预测准确性得到提升，反映医疗与政策应对。
• 北京和上海的总感染病例估计约4百左右，而大陆、中国湖北、武汉等地区的总量预测更大，区域差异可由措施和数据报告变化解释。