[論文レビュー] A simple mathematical model for the evolution of the corona virus

Stefan Tappe|arXiv (Cornell University)|Mar 20, 2020

COVID-19 epidemiological studies参考文献 6被引用数 4

ひとこと要約

本稿では、対数変換と区分的関数を用いたシンプルな2パラメータの数学的モデルを提案し、COVID-19の死亡者数を予測する。政府の対策前には凸的成長、その後には凹的安定化を示す。モデルは中国のデータに良好に適合し、イタリアの死亡者トレンドを予測する。早期の対策が死亡者数を半減させられることを示しており、学生や公衆衛生計画者にとっても容易に利用可能である。

ABSTRACT

The goal of this note is to present a simple mathematical model with two parameters for the number of deaths due to the corona (COVID-19) virus. The model only requires basic knowledge in differential calculus, and can also be understood by pupils attending secondary school. The model can easily be implemented on a computer, and we will illustrate it on the basis of case studies for different countries.

研究の動機と目的

• 基本的な微積分のみを用いて、COVID-19の死亡者数の推移を数学的にシンプルかつ効果的に予測するモデルを開発すること。
• 最小限の数学的訓練しか持たない実務家や学生が、このモデルを理解し、応用できるようにすること。
• 死亡者数の指数的成長から飽和成長への遷移をモデル化することで、政府の対策の影響を推定すること。
• 既知の致死率と報告率を用いて、感染人数と感染確率を推定するためのモデルの拡張。
• 中国およびイタリア、ドイツ、アメリカの実世界データを用いて、モデルの妥当性を検証し、将来の結果を予測すること。

提案手法

• 死亡者数 D(t) を対数関数 L(t) の指数関数としてモデル化し、D(t) = exp(L(t)) とすることで、成長ダイナミクスを線形化する。
• L(t) を区分的に定義する：[t₁, T₁] では、凸性パラメータ β ∈ (0,1] を持つべきべき関数を用い、加速的成長を捉える。
• (T₁, ∞) では、L(t) = L(T₁) + λ(1 − exp(−ν(t − T₁))) とし、ロジスティック型関数を用いて減速的成長と安定化をモデル化する。
• 初期の死亡者データを用いて β と λ を推定し、同様の対策の厳しさを持つ国々では λ が共通であると仮定する（例：中国・イタリアでは λ = ln(3.25)）。
• T₁ を政府対策開始後17日目とする。これは平均潜伏期間から死亡までの期間に基づく。
• I(t) = D(t+17)/(μ·κ) を用いて、感染済み人数 I(t) を推定する。μ = 0.01 および κ = 5 を標準値として用いる。

実験結果

リサーチクエスチョン

• RQ1 2つのパラメータのみを有するシンプルな数学的モデルが、異なる国におけるCOVID-19死亡者数の推移をどのように捉えることができるか？
• RQ2 政府の対策のタイミングが、モデル予測の最終死亡者数に与える影響は何か？
• RQ3 本モデルは中国の歴史的死亡者データに良好に適合し、イタリアのような国における将来の傾向を信頼性高く予測できるか？
• RQ4 凸性パラメータ β は、対策前後における死亡者数の増加パターンとどのように関係しているか？
• RQ5 死亡者データのみを用いて、真の感染人数や感染確率をどの程度正確に推定できるか？

主な発見

• 中国の死亡者データに良好に適合し、β = 0.6 および λ = ln(3.25) を用いて、2020年2月10日頃に凸的成長から凹的成長への遷移が明確に現れる。
• イタリアでは β = 0.65 と推定され、λ = ln(3.25) を同じく使用し、2020年4月下旬までに死亡者数のピークが約13,000人になると予測される。
• シナリオ分析により、イタリアの政府対策を3日間早める（3月28日から3月25日へ）ことで、合計死亡者数が約半減すると示された。
• β > 0.6 を示す国（例：ドイツ、イギリス）では、より少ない凹性成長が予測され、死亡者数の増加がより持続的になると示唆される。
• μ = 0.01 および κ = 5 を用いた推定により、イタリアの真の感染人数は報告済み症例の最大50倍に達する可能性があると示された。