SIRモデル

SIRモデル（エスアイアールモデル）は、感染症の短期的な流行過程を決定論的に記述する古典的なモデル方程式である。名称はモデルが対象とする、感受性保持者（Susceptible）、感染者（Infected）、免疫保持者（Recovered、あるいは隔離者 Removed）の頭文字にちなむ。原型となるモデルは、テンプレート:仮リンクとテンプレート:仮リンクの1927年の論文で提案された。単純なSIRモデルであっても、1905–06年のボンベイにおけるペスト流行のデータをうまく再現することが知られている。

概要

SIRモデルの解の挙動例。縦軸は人数、横軸は時間で、青=S, 緑=I, 赤=Rである。

SIRモデルの相平面上の軌道 (S, I)。簡単のため テンプレート:Math とおくと、テンプレート:Math が成り立ち、テンプレート:Mvar は テンプレート:Math のとき最大値 テンプレート:Math をとる。

感染人口密度 テンプレート:Math の時間変化と テンプレート:Math の値の関係。

SIRモデルにおいて、全人口は感受性保持者・感染者・免疫保持者の3つへ分割され、感受性保持者Sは感受性保持者Sと感染者Iの積に比例して定率で感染者Iに移行し、感染者Iは定率で免疫保持者Rに移行する（感染期間は指数分布に従う）と仮定される。この時間発展を非線形常微分方程式で記述される連続力学系として表せば、

$\begin{array}{cc}\frac{dS}{dt}(t)& =-\beta S(t)I(t)\\ \frac{dI}{dt}(t)& =\beta S(t)I(t)-\gamma I(t)\\ \frac{dR}{dt}(t)& =\gamma I(t)\end{array}$

となる。ただし、テンプレート:Math は感染率、テンプレート:Math は回復率（隔離率）を表す（逆数 テンプレート:Math は平均感染期間を表す）。これをフローチャートで

テンプレート:Border$\stackrel{\beta I}{⟶}$テンプレート:Border$\stackrel{\gamma }{⟶}$テンプレート:Border

のように表すこともある。

上記の3式の和を取れば、

$\frac{d}{dt}(S(t)+I(t)+R(t))=0$

であり、これは総人口 テンプレート:Math が一定値をとる保存則（閉鎖人口の仮定）

$S(t)+I(t)+R(t)=\mathrm{c}\mathrm{o}\mathrm{n}\mathrm{s}\mathrm{t}\mathrm{.}$

に対応している。この保存則により、本質的に2変数の方程式である。

簡単のため初期値を テンプレート:Math, テンプレート:Math とおくと

$\frac{dI}{dt}(0)=I_0(\beta S_0-\gamma )>0$

のとき、すなわち

${ℛ}_0=\beta S_0/\gamma >1$

のとき流行が発生する（閾値現象）。この無次元量 テンプレート:Math は基本再生産数と呼ばれる。上のような最も単純なモデルでは $\underset{t\to \mathrm{\infty }}{\lim }I(t)=0$ が成り立ち、エンデミックな定常状態や周期的な流行といった現象は説明できない。

派生モデル

SIRモデルにおいて、出生・死亡などによる人口変動を考慮したモデルや、マスター方程式による確率的モデルが存在する。また、免疫獲得を考慮しないSISモデルや潜伏期間を考慮したSEIRモデルなど色々な区画モデルが知られている。他にも、感染年齢を考慮した偏微分方程式によるモデル（カーマック・マッケンドリック理論）がある。

参考文献

• テンプレート:Cite book
• テンプレート:Cite book
  • テンプレート:Citation（書評）

関連項目

• 疫学における区画モデル
• カーマック・マッケンドリック理論
• SEIRモデル
• 基本再生産数 テンプレート:Math
• テンプレート:仮リンク