[论文解读] A Time-dependent SIR model for COVID-19

Yicheng Chen, Ping-En Lu+1 more|arXiv (Cornell University)|Feb 28, 2020

COVID-19 epidemiological studies参考文献 8被引用 52

一句话总结

本文提出了一种时变SIR模型，可动态追踪传播率与恢复率，以高精度预测COVID-19的传播、峰值及结束时间——在中国的每日预测误差约为3%。该模型进一步扩展以评估无症状传播、群体免疫阈值及社交距离措施的有效性，表明拐点后R₀ < 1，且通过控制传播率可降低R₀。

ABSTRACT

In this paper, we conduct mathematical and numerical analyses to address the following crucial questions for COVID-19: (Q1) Is it possible to contain COVID-19? (Q2) When will be the peak and the end of the epidemic? (Q3) How do the asymptomatic infections affect the spread of disease? (Q4) What is the ratio of the population that needs to be infected to achieve herd immunity? (Q5) How effective are the social distancing approaches? (Q6) What is the ratio of the population infected in the long run? For (Q1) and (Q2), we propose a time-dependent susceptible-infected-recovered (SIR) model that tracks 2 time series: (i) the transmission rate at time t and (ii) the recovering rate at time t. Such an approach is more adaptive than traditional static SIR models and more robust than direct estimation methods. Using the data provided by China, we show that the one-day prediction errors for the numbers of confirmed cases are almost in 3%, and the total number of confirmed cases is precisely predicted. Also, the turning point, defined as the day that the transmission rate is less than the recovering rate can be accurately predicted. After that day, the basic reproduction number $R_0$ is less than 1. For (Q3), we extend our SIR model by considering 2 types of infected persons: detectable and undetectable infected persons. Whether there is an outbreak in such a model is characterized by the spectral radius of a 2 by 2 matrix that is closely related to $R_0$. For (Q4), we show that herd immunity can be achieved after at least 1-1/$R_0$ fraction of individuals being infected. For (Q5) and (Q6), we analyze the independent cascade (IC) model for disease propagation in a configuration random graph. By relating the propagation probabilities in the IC model to the transmission rates and recovering rates in the SIR model, we show 2 approaches of social distancing that can lead to a reduction of $R_0$.

研究动机与目标

• 开发一种可动态适应传播率与恢复率变化的SIR模型，以提升流行病预测能力。
• 通过识别传播率低于恢复率的拐点，确定流行病峰值与结束时间。
• 利用双类型感染者模型评估无症状感染在疾病传播中的作用。
• 在时变传播背景下，量化R₀函数的群体免疫阈值。
• 评估通过传播率调控降低R₀的社交距离措施的有效性。

提出的方法

• 构建一种具有独立时变传播率与恢复率的时变SIR模型，以增强对静态模型的适应性。
• 利用中国实时数据标定传播率与恢复率，实现精确的短期预测。
• 将拐点定义为传播率降至恢复率以下的日期，标志着R₀ < 1及流行病开始消退。
• 将SIR模型扩展以包含可检测与不可检测感染者，通过2×2矩阵建模其影响，其谱半径决定疫情爆发潜力。
• 在配置随机图上应用独立级联（IC）模型，模拟疾病传播，并将IC传播概率与SIR参数关联。
• 推导出两种通过降低SIR框架中有效传播率来降低R₀的社交距离策略。

实验结果

研究问题

• RQ1 COVID-19能否被控制？在何种条件下可以实现？
• RQ2 基于时变传播率与恢复率，流行病何时达到峰值并结束？
• RQ3 无症状感染在疾病传播中起何种作用？
• RQ4 需多少比例人口感染才能实现群体免疫？
• RQ5 社交距离措施在降低R₀与控制流行病方面有多有效？

主要发现

• 时变SIR模型对确诊病例的预测误差约为3%，并准确预测了中国的确诊病例总数。
• 当传播率降至恢复率以下时，可精确预测拐点，标志着R₀ < 1及流行病开始消退。
• 若2×2矩阵的谱半径超过1，不可检测感染者的存在可能导致疫情爆发，表明具有类似R₀的临界行为。
• 仅当至少1 - 1/R₀的人口被感染后，群体免疫才可实现，凸显R₀依赖的关键阈值。
• 通过降低传播率来减少R₀的社交距离策略，能有效降低R₀，且IC模型将这些降低与疾病传播动力学的可测量变化相联系。
• 该模型通过在预测精度上优于静态SIR模型与直接估计方法，展现出强鲁棒性与适应性。