[論文レビュー] Covid-19 spread: Reproduction of data and prediction using a SIR model on Euclidean network

Kathakali Biswas, Abdul Khaleque+1 more|arXiv (Cornell University)|Mar 16, 2020

COVID-19 epidemiological studies参考文献 21被引用数 69

ひとこと要約

本論文は、中国のCOVID-19データをユークリッド網上のSIRモデルで適合させ、時間的・距離的パターンを再現し、局所的な流行結論を予測する。

ABSTRACT

We study the datafor the cumulative as well as daily number of cases in the Covid-19 outbreak in China. The cumulative data can be fit to an empirical form obtained from a Susceptible-Infected-Removed (SIR) model studied on an Euclidean network previously. Plotting the number of cases against the distance from the epicenter for both China and Italy, we find an approximate power law variation with an exponent $\sim 1.85$ showing strongly that the spatial dependence plays a key role, a factor included in the model. We report here that the SIR model on the Eucledean network can reproduce with a high accuracy the data for China for given parameter values, and can also predict when the epidemic, at least locally, can be expected to be over.

研究の動機と目的

• 中国におけるCOVID-19の空間的・時間的拡散を、ユークリッドネットワーク上のSIRモデルを用いて理解する。
• ユークリッドネットワーク上のSIRから導出された経験的形を用いて累積症例データを適合させる。
• 震央からの距離と症例数の関係を検討し、空間的依存性を評価する。

提案手法

• 距離に基づく長距離接続（P(ℓ) ∝ ℓ^{−δ}）と感染確率qを持つ、ユークリッドネットワーク上のSIRモデルを用いる。
• 累積感染データR(t)をR(t) = a exp(t/T) / [1 + c exp(t/T)]に適合させ、a, c, Tを適合パラメータとする。
• R(t)を微分してI(t)を計算し、ピーク時刻t_p = T ln(1/c)を求める。
• 最近傍接続に1本のランダム長距離リンクを加えた2^12ノードネットワークをシミュレートし、100通りの構成で平均結果を比較する。
• モデル出力を現実データの飽和に合わせて再スケーリングし、比較のためにモンテカルロステップを日数に換算する。

実験結果

リサーチクエスチョン

• RQ1 ユークリッドネットワーク上のSIRモデルは、中国におけるCOVID-19の累積症例の時間発展を再現できるか。
• RQ2 震央からの距離に伴う症例の空間的減衰をモデルは捉えられるか、そしてこれは接触範囲δについて何を示唆するか。
• RQ3 データを再現するパラメータ値（δ, q）は何であり、それらが伝播ダイナミクスと流行期間に何を示唆するか。
• RQ4 モデルによる流行終息の推定時期はいつか。

主な発見

• 中国の累積データはR(t) = a exp(t/T) / [1 + c exp(t/T)]で適合でき、a = 1.15698×10^−6, c = 0.0201529, T = 5.22、パラメータ誤差は小さい。
• 距離に対する感染率の減衰はべき乗則R(d) ∝ d^−γに従い、γ = 1.85 ± 0.1で、中国（武漢）およびイタリア（ベルガモ）震央で観測。
• δ ≈ 1.8とq ≈ 0.85のユークリッドネットワーク上のSIRモデルは、適切な再スケーリング（露出割合ρ ≈ 6.0783×10^−5）後に中国データを再現する。
• モデルは、2020年1月21日から約60–78日で中国の流行は終息する可能性があると予測し、感染者0は約82モンテカルロステップ（約78日）で生じる。
• より大きいδは、実質的な接触範囲が短いことを示唆する；比較的高いqは、モデル化された文脈での伝播が強いことを反映している；結果は移動を伴わない抑制解釈と整合する。
• δとqの小さな変動（±0.05）に対して適合は頑健であり、データを再現するために明示的な予防因子を必要としない。