[論文レビュー] The D model for deaths by COVID-19

J. E. Amaro|arXiv (Cornell University)|Mar 30, 2020

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

ひとこと要約

本論文は、SI疫学的モデルから導出された簡素化された解析的フレームワーク、Dモデルを提案する。このモデルは、回復を仮定せず、死亡者数を感染の遅延・スケーリング関数としてモデル化することで、中国、スペイン、イタリアのデータに3つの調整可能なパラメータで適合し、最小限の仮定のもとでパンデミックのピーク時期と plateau 水準を正確に予測できる。

ABSTRACT

We present a simple analytical model to describe the fast increase of deaths produced by the corona virus (COVID-19) infections. The 'D' (deaths) model comes from a simplified version of the SIR (susceptible-infected-recovered) model known as SI model. It assumes that there is no recovery. In that case the dynamical equations can be solved analytically and the result is extended to describe the D-function that depends on three parameters that we can fit to the data. Results for the data from Spain, Italy and China are presented. The model is validated by comparing with the data of deaths in China, which are well described. This allows to make predictions for the development of the disease in Spain and Italy.

研究の動機と目的

• 複雑な数値シミュレーションに依存せずに、COVID-19死亡の推移を迅速に予測できる解析的モデルの開発。
• 回復なしを仮定したSIフレームワークに基づく簡素化モデルを、複数国の実際の死亡データにフィットさせることで、その妥当性を検証すること。
• 進行中の流行がある国での初期段階のパンデミック予測に適した、低パラメータ数かつ解釈可能なツールの提供。
• 中国では既にピークに達していたが、そのデータの適合性を評価し、スペインとイタリアの活発な拡散段階にそのモデルを適用することで、予測性能を検証すること。

提案手法

• モデルはSI方程式に基づき、回復なし（R=0）を仮定することで、感染者数I(t)に関する1つの微分方程式に簡略化される。
• SI方程式の解は、I(t)のロジスティック型関数に変換され、これを利用して死亡者数D(t)をI(t)の遅延およびスケーリング版としてモデル化する。
• D(t)関数は、D(t) = a * exp((t-t₀)/b) / (1 + c * exp((t-t₀)/b)) として表現され、3つのフィットパラメータ（a：振幅、b：時間定数、c：初期条件スケーリング）を有する。
• パラメータa、b、cは、中国、スペイン、イタリアの実際の死亡データに非線形回帰を用いてフィットさせ、感染伝播率や回復率といった生物学的パラメータの知識は不要である。
• 感染から死亡までの固定遅延τと一定の致死率mは、すべてのパラメータに組み込まれており、フィットに必要なパラメータに吸収されている。
• 中国のデータ（すでに plateau に達していた）を用いてモデルを検証した後、スペインとイタリアにおけるパンデミックのピーク時期と最終死亡者数の予測に応用した。

実験結果

リサーチクエスチョン

• RQ1 SIフレームワークに基づく簡素化された解析的モデルは、複数国のCOVID-19死亡推移を正確に記述できるか？
• RQ2 3パラメータのみで構成されるDモデルは、中国、スペイン、イタリアの実際の死亡データにどれほど適合するか？
• RQ3 進行中の流行段階にある国において、モデルはパンデミックのピークと最終死亡者数を信頼性高く予測できるか？
• RQ4 回復なしを仮定することは、モデルの予測精度と適用可能性にどのような影響を与えるか？

主な発見

• Dモデルは中国の死亡データに顕著に適合しており、他の国への予測ツールとしての有効性が裏付けられた。
• スペインについて、モデルはパンデミックが発症後21日目ごろにピークに達すると予測し、2020年4月中旬までに plateau に達すると予想した。
• イタリアでは、データ収集開始から約20日後にパンデミックのピークに達するとモデルが予測した。スペインと同程度の期間であった。
• Dモデルの3パラメータ関数は、3カ国すべての死亡曲線における急激な上昇とその後の plateau 段階を高い忠実度で捉えた。
• パラメータa、b、cのみで、感染伝播率や回復時間といった生物学的パラメータの推定を必要とせずに、死亡曲線全体を記述できた。
• 中国のデータのみを用いても、スペインやイタリアのようなパンデミックの初期段階にある国に対して、Dモデルは強力な予測能力を示した。