[論文レビュー] Using generalized logistics regression to forecast population infected by Covid-19

Mario Villalobos-Arias|arXiv (Cornell University)|Apr 6, 2020

COVID-19 epidemiological studies参考文献 3被引用数 20

ひとこと要約

本論文では、初期の流行データにS字型成長曲線をフィットさせることで、累積COVID-19感染数を予測するために一般化ロジスティック（GL）およびゴンペルツ回帰モデルの使用を提案する。この手法は非線形最適化を活用して限られたデータ（例：20〜35日分）からパラメータを推定し、複数の国でR² > 0.998を達成した。これにより、最小限のデータ入力で流行のピークサイズと変曲点を信頼性高く予測できる。

ABSTRACT

In this work, a proposal to forecast the populations using generalized logistics regression curve fitting is presented. This type of curve is used to study population growth, in this case population of people infected with the Covid-19 virus; and it can also be used to approximate the survival curve used in actuarial and similar studies.

研究の動機と目的

• 柔軟な成長曲線モデルを用いて、データ効率の良い累積COVID-19症例数の予測手法を開発すること。
• 感染率と回復率を推定する必要があるため、SIRモデルに伴う限界を克服すること。
• 初期流行データ（20〜35日分）が、全流行曲線およびピークサイズを信頼性高く予測できるかどうかを検証すること。
• 一般化ロジスティック関数とゴンペルツ関数が、実世界の感染データに適合する性能を比較すること。
• 公衆衛生当局が流行のピークを予測し、介入の影響を評価するための実用的ツールを提供すること。

提案手法

• 一般化ロジスティック（GL）関数を採用：$ P(t) = \frac{M}{(1 + e^{-at + b})^\alpha} $、ここで$ M $は最大人口、$ a $は成長率を制御、$ b $は曲線のシフトを示し、$ \alpha $は柔軟性を追加する。
• ゴンペルツ関数を適用：$ G(t) = a \cdot e^{-b \cdot e^{-ct}} $、3パラメータモデルで、感染曲線に適した固有の非対称性を有する。
• 線形化による歪みを避けるために、非線形最適化技術を用いてモデルを累積感染データにフィットさせる。
• 決定係数（$ R^2 $）を用いてモデルの適合度を評価し、7%の誤差マージンを含む予測区間を設定する。
• 2020年3月から4月にかけてのデータを用いて、中国、イタリア、スペイン、米国など複数の国でGLおよびゴンペルツモデルの予測性能を比較する。
• 公衆衛生介入の影響を評価するために、時間経過に伴う予測最大症例数（$ M $）の変化を追跡する。

実験結果

リサーチクエスチョン

• RQ1 一般化ロジスティックおよびゴンペルツ回帰モデルは、20〜35日分のデータのみを用いて、累積COVID-19感染数を正確に予測できるか？
• RQ2 中国、イタリア、スペイン、米国などの国々の実世界の流行曲線に適合する際、GLモデルとゴンペルツモデルの予測性能はどのように比較できるか？
• RQ3 公衆衛生介入が実施されるに従い、予測される流行ピークサイズ（$ M $）は時間経過とともに減少するか？また、この変化はモデルで捉えられるか？
• RQ4 日次症例増加が低下し始める時期である変曲点は、初期データからこれらのモデルを用いて信頼性高く推定できるか？
• RQ5 韓国のような異なる流行ダイナミクスを示す国に対しても、モデルの予測はどの程度頑健に保たれるか？

主な発見

• 一般化ロジスティックモデルは、2020年3月31日時点のコスタリカのデータに$ R^2 = 0.99858 $を達成し、予測ピークは886例であった。
• ゴンペルツモデルはコスタリカで$ R^2 = 0.99854 $を達成し、ピークは887例を予測し、GLモデルと強い一致を示した。
• イタリアではゴンペルツモデルが$ R^2 = 0.99968 $を達成し、2020年3月30日までのデータを用いてピーク261,052例を予測した。
• スペインのゴンペルツ適合は$ R^2 = 0.99941 $を達成し、2020年3月31日までのデータでピーク468,495例を予測した。
• 米国のゴンペルツ回帰モデルは$ R^2 = 0.999915 $を達成し、2020年4月5日までのデータでピーク1,486,347例を予測した。
• コスタリカの予測最大症例数（$ M $）は、3月23日の約2,162例から4月5日には743例に減少し、政府の介入の影響が明確に捉えられた。