[論文レビュー] A note on Exact solution of SIR and SIS epidemic models

Ghulam Shabbir, Haroon Ahmed Khan+1 more|arXiv (Cornell University)|Dec 22, 2010

Mathematical and Theoretical Epidemiology and Ecology Models参考文献 12被引用数 26

ひとこと要約

本稿では、定常人口を仮定した特定のSIRおよびSIS疫病モデルに対して、直接積分法を用いて正確な解析解を提示している。著者らは、制限のないパrameter条件のもとで、両モデルの閉形式解を導出しており、一般の初期条件のもとでの疾病動態の正確な分析を可能にするという、数学的疫学分野における顕著な進展をもたらしている。

ABSTRACT

In this article we have successfully obtained an exact solution of a particular case of SIR and SIS epidemic models given by Kermack and Mckendrick [1] for constant population, which are described by coupled nonlinear differential equations. Our result has no limiting conditions for any parameter involved in the given models. In epidemiology many researchers believe that it is very hard to get an exact solution for such models. We hope this solution will be an opening window and good addition in the area of epidemiology.

研究の動機と目的

• 定常人口を仮定したSIRおよびSIS疫病モデルの正確な解析解を導出すること。これは、数学的疫学分野における長年の課題を解決することを目的としている。
• 疾病拡散を記述する非線形連立微分方程式における非可積分性という一般的な困難を克服すること。
• モデルパラメータに制限を課さない解法フレームワークを提供することで、モデルの適用可能性と解釈可能性を向上させること。
• 正確な関数形を用いた、疾病動態のより深い定性的分析の基盤を提供すること。

提案手法

• SISモデルに関しては、保存則 s(t) + i(t) = k を用いて、i(t) に関する非線形常微分方程式にシステムを簡略化し、その後 y = i⁻¹ の代入を実行する。
• SISの解は、簡略化された常微分方程式に y = i⁻¹ を代入した後、直接積分を適用することで得られ、積分定数 C を含む指数関数的形をとる。
• SIRモデルに関しては、方程式を加算することで保存則 s(t) + i(t) = 1 + C e⁻ᵘᵗ を導出し、i(t) に関する単一の常微分方程式にシステムを簡略化する。
• SIRの解は、z = i⁻¹ の代入を適用し、e⁻ᵘᵗ の級数展開を用いて高次項を無視することで、積分を簡略化する。
• 両モデルとも、初期条件を用いて積分定数を決定し、s(t) および i(t) の閉形式表現を導出する。
• この手法は、数値的または摂動的アプローチに依存せず、直接積分と代入技術に依存することで、正確性を保証している。

実験結果

リサーチクエスチョン

• RQ1 定常人口のもとで、制限のないパrameter制約のもとでSIRおよびSISモデルの正確な解析解を導出可能か？
• RQ2 これらの古典的疫病モデルの非線形連立微分方程式を、解析的積分法により正確に解く方法は何か？
• RQ3 保存則 s(t) + i(t) = k（または s(t) + i(t) = 1 + C e⁻ᵘᵗ）が正確解の導出を可能にする役割は何か？
• RQ4 正確解は、疫学における疾病動態の理解およびモデル検証の向上に、どの程度寄与するか？

主な発見

• SISモデルの正確解は、i(t) = β / (r + βC e⁻ᵝᵗ) として導出され、ここで β = rk - α であり、C は初期条件 i(0) = i₀ によって決定される。
• SISの解は s(t) + i(t) = k を正確に満たしており、全人口の保存が確認される。
• SIRモデルに関しては、i(t) = λ / (β + λD e⁻ᵘᵗ e^(βC/μ)) として表され、λ = β - μ + βC であり、C = s₀ + i₀ - 1 である。
• SIRの解は、初期条件 s(0) = s₀ および i(0) = i₀ を用いて、積分定数 D および C を完全に決定する。
• SISの場合は近似を用いないが、SIRの場合は簡略化のため、e⁻ᵘᵗ ≈ 1 - μt の級数展開を用い、高次項を無視している。
• 導出された解は、パrameter r, α, β, μ および初期条件 s₀, i₀ > 0 のすべての正の値に対して有効であり、制限条件なしに成立する。