疫学における区画モデルのソースを表示

'''区画モデル'''（くかくモデル、{{lang-en-short|Compartmental model}}）は、単純かつ一般的な{{仮リンク|感染症の数理モデル化|en|Mathematical modelling of infectious disease}}の方法の総称である。区画モデルでは、集団内の個体は、感染状態に応じて、それぞれ「区画」と呼ばれるラベルを割り当てられる。例えば、'''S'''、'''I'''、'''R'''という区画は、それぞれ'''S'''usceptible〈感受性保持者〉、'''I'''nfectious〈感染者〉、'''R'''ecovered〈免疫保持者〉を表す。個体は感染状態の進行にともなって、区画間を移行する。具体的な区画モデルは、例えばSIRモデル、SIESモデルなどと呼ばれるが、こうしたモデルの名称となるラベルの順番は通常、区画間の移行のパターンを示している。例えば、SEISは、感受性（susceptible）、曝露（exposed）、感染（infectious）、そして再び感受性（susceptible）の順に感染状態が進行することを意味する。

こういったモデルの起源は20世紀初頭であり、1927年の{{仮リンク|ウィリアム・オグルヴィ・カーマック|en|William Ogilvy Kermack|label=カーマック}}と{{仮リンク|アンダーソン・グレイ・マッケンドリック|en|Anderson Gray McKendrick|label=マッケンドリック}}の研究が重要である<ref name="Kermack–McKendrick">{{cite journal |last1=Kermack |first1=W. O. |last2=McKendrick |first2=A. G. |title=A Contribution to the Mathematical Theory of Epidemics |journal=Proceedings of the Royal Society A |volume=115 |issue=772 |pages=700–721 |date=1927 |doi=10.1098/rspa.1927.0118|bibcode=1927RSPSA.115..700K |doi-access=free }}</ref>。

モデルは（決定論的である）[[常微分方程式]]を用いて実行されることが多いが、確率論的（ランダム）な枠組みを用いることも可能である。確率論的モデルはより現実的ではあるが解析がはるかに複雑である。

モデルは、病気がどのように広がるか、感染者の総数、流行の期間などを予測し、[[基本再生産数|再生産数]]などの様々な疫学的パラメータを推定しようとする。このようなモデルは、異なる[[公衆衛生的介入|公衆衛生上の介入]]が流行の結果にどのように影響を与えるかを示すことができる。例えば、与えられた集団に限られた数の[[ワクチン]]を接種するために最も効率的な手法は何かを示すことができる。

== SIRモデル ==
{{Main|SIRモデル}}
'''SIRモデル'''は最も単純な区画モデルの1つであり、多くのモデルはこの基本形から派生している。本モデルは3つの区画で構成されている。
; S
: 感受性（'''s'''usceptible）個体の数。感受性個体と感染個体が「感染性接触」すると、感受性個体が病気に感染し、感染性区画に移行する。
; I
: 感染（'''i'''nfectious）個体数のこと。感染した個体であり、感受性個体を感染させる可能性がある。
; R
: 隔離（'''r'''emoved）（免疫のある）個体、または死亡した個体の数。これらは、感染から回復して隔離区画に入った個体、または死亡した個体である。死亡者数は総人口に対して無視できるほどの数であると仮定している。この区画を「回復（'''r'''ecovered）」または「抵抗性（'''r'''esistant）」と呼ぶこともある。

このモデルは、[[麻疹]]、[[おたふくかぜ]]、[[風疹]]といった、回復が持続的な抵抗性をもたらし、ヒトからヒトへ感染する感染症について合理的に予測可能である{{Citation needed|date=April 2020}}。

[[File:SIR model simulated using python.gif|thumb|空間SIRモデルシミュレーション。個々のセルは隣接した8つの区画を感染させることができる。]]

これらの変数（'''S'''、'''I'''、'''R'''）は、特定の時間に各区画にいる人の数を表す。感受性個体、感染個体、隔離個体の数が（総人口サイズが一定であっても）時間とともに変化する可能性があることを表すために、正確な数をt（時間）の関数'''S'''(''t'')、'''I'''(''t'')、'''R'''(''t'') とする。特定の集団における特定の疾患については、これらの関数は、可能性のある[[アウトブレイク]]を予測し、それらを制御下に置くために働くかもしれない{{Citation needed|date=June 2020}}。

===SIRモデルは3つの意味で動的である===
''t''の変数関数で表わされるように、本モデルは各区画内の数が時間の経過とともに変動するという意味で動的である。この動的な側面の重要性は、1968年にワクチンが導入される前の英国の麻疹のように、感染期間が短い[[風土病]]において最も明白である。このような病気は、感受性保持者数（S(''t'')）が時間の経過とともに変化するため、何度も流行する傾向がある。[[伝染病]]が流行すると、より多くの感受性保持者が感染区画に入るため、感受性保持者の数が急激に減少し、彼らは感染区画と隔離区画に入ることになる。生まれた子どもが感受性区画に入るなどして、感受性個体の数が回復するまでは、病気が再び流行することはできない。

[[File:Graph SIR model without vital dynamics.svg|thumb|黄色=感受性、[[マルーン (色)|えび茶色]]=感染、[[鴨の羽色]]=回復]]

集団の各メンバーは、通常、感受性区画から感染性区画へ、そして隔離区画へと進行する。これは以下のような、箱が異なる区画を表し、矢印が区画間の遷移を表すフロー図として示すことができる。
[[File:SIR Flow Diagram.svg|600px|center|SIR compartment model]]

===遷移速度===
モデルを完全に明確に記述するためには、矢印を区画間の遷移速度でラベル付けするべきである。''S''と''I''の間では、遷移速度を''d(S/N)/dt = -βSI/N<sup>2</sup>''と仮定する。ここで、''N''は総人口、βは時間当たりの1人当たりの平均接触数に感受性保持者と感染者の間の接触における病気の伝播確率を乗じたもの、''SI/N<sup>2</sup>''は感受性保持者と感染者の間の接触のうち、感受性保持者が感染することになる接触の割合である（これは数学的には、分子間のランダムな衝突が化学反応を引き起こし、分率が2つの反応物の濃度に比例するという化学の質量作用の法則に似ている）。

''I''と''R''の間では、遷移速度は感染個体数に比例すると仮定し、これをγ''I''とする。これは、任意の時間間隔''dt''において感染個体が回復する確率が単純にγ''dt''であると仮定していることと等価である。ある個体が平均的な時間間隔''D''で感染している場合、γ = 1/''D''となる。これは、個人が感染状態にある時間の長さが[[指数分布]]を持つランダム変数であるという仮定と等価である。この「古典的」SIRモデルは、I-R遷移速度のためにより複雑で現実的な分布（例えば、[[アーラン分布]]<ref name="Krylova2013">{{cite journal |last1=Krylova |first1=O. |last2=Earn |first2=DJ |date= May 15, 2013|title= Effects of the infectious period distribution on predicted transitions in childhood disease dynamics |url=https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3673147/ |journal=J R Soc Interface|volume=10 |issue=84 |pages= |doi=10.1098/rsif.2013.0098 |access-date=June 13, 2020|doi-access=free }}</ref>）を使用することによって修正することができる。

感染区画からの除去（隔離）がない特別な場合 (γ=0) では、SIRモデルは、すべての個体が最終的に感染する[[ロジスティック分布|ロジスティック]]解を持つ非常に単純なSIモデルに簡約される。

==SIRモデルの生物数理学的・決定論的取り扱い==
===人口動態がないSIRモデル===
伝染病（例えば[[インフルエンザ]]）の動態は、出生と死亡の動態よりもはるかに速いことが多いため、単純な区画モデルでは出生と死亡が省略されることが多い。先に述べたいわゆる人口動態（出生および死亡）を考慮しないSIRシステムは以下の一組の常微分方程式で表わすことができる<ref>{{cite journal |author=Hethcote H |title=The Mathematics of Infectious Diseases |journal=SIAM Review |volume=42 |issue= 4|pages=599–653 |year=2000 |doi=10.1137/s0036144500371907|bibcode=2000SIAMR..42..599H }}</ref>。
:<math>
\begin{align}
& \frac{dS}{dt} = - \frac{\beta I S}{N}, \\[6pt]
& \frac{dI}{dt} = \frac{\beta I S}{N}- \gamma I, \\[6pt]
& \frac{dR}{dt} = \gamma I.
\end{align}
</math>
上式において、<math>S</math>は感受性集団のストック、<math>I</math>は感染集団のストック、<math>R</math>は隔離集団のストック（死亡または回復による）、そして<math>N</math> はこれら3つの和である。

このモデルは、[[ウィリアム・オグルヴィ・カーマック]]と[[アンダーソン・グレイ・マッケンドリック]]によって、現在[[カーマック・マッケンドリック理論]]と呼ばれているものの特別な場合として初めて提案されたものであり、マッケンドリックが[[ロナルド・ロス]]と共に行った研究に続いて行われた。

このシステムは[[非線形]]であるが、陰関数形式でその[[解析解]]を導出することが可能である<ref name="Harko">{{Cite journal|last1=Harko|first1=Tiberiu|last2=Lobo|first2=Francisco S. N.|last3=Mak|first3=M. K. |title=Exact analytical solutions of the Susceptible-Infected-Recovered (SIR) epidemic model and of the SIR model with equal death and birth rates |journal=Applied Mathematics and Computation|language=en|volume=236|pages=184–194|year=2014 |doi=10.1016/j.amc.2014.03.030|bibcode=2014arXiv1403.2160H |arxiv=1403.2160 }}</ref>。他の数値計算ツールとしては、{{仮リンク|ギレスピー・アルゴリズム|en|Gillespie algorithm}}といった[[モンテカルロ法]]がある。

第1に、
:<math> \frac{dS}{dt} + \frac{dI}{dt} + \frac{dR}{dt}  = 0</math>
から、
:<math> S(t) + I(t) + R(t) = \text{一定} = N</math>
となる。この式は数学的に人口<math> N </math>の恒常性を表わしている。

ここで留意すべきは、上記の関係は3つの変数のうち2つに関する式だけを調べる必要があることを示唆している点である。

第2に、伝染病クラスの動態が比率
:<math> R_0 = \frac{\beta}{\gamma},</math>
いわゆる[[基本再生産数]]（基本再生産率とも）に依存することに注目する。この比率は、すべての対象者が感受性である集団における1つの感染からの新規感染（これらの新規感染は二次感染と呼ばれることもある）の予想数として導き出される<ref name=Bailey1975>{{cite book |author=Bailey, Norman T. J. |title=The mathematical theory of infectious diseases and its applications |publisher=Griffin |location=London |year=1975 |isbn=0-85264-231-8 |edition=2nd}}</ref><ref name=nunn2006>{{cite book |author1=Sonia Altizer |author2=Nunn, Charles |title=Infectious diseases in primates: behavior, ecology and evolution |publisher=Oxford University Press |location=Oxford [Oxfordshire] |year=2006 |isbn=0-19-856585-2 |series=Oxford Series in Ecology and Evolution}}</ref>。この考えは、接触の間の典型的な時間は <math>T_{c} = \beta^{-1}</math>であり、隔離までの典型的な時間は <math>T_{r} = \gamma^{-1}</math>であると言えば、おそらくもっとわかりやすいだろう。ここから、平均して、感染者が隔離される「前」に他の人と接触した回数は<math>T_{r}/T_{c}</math>となる。

1つ目の微分方程式を3つ目の微分方程式で割って、[[変数分離|変数を分離]]して[[積分]]すると以下のようになる。
:<math> S(t) = S(0) e^{-R_0(R(t) - R(0))/N}</math>
上式において、<math>S(0)</math>および<math>R(0)</math>はそれぞれ感受性保持者および隔離者の初期数である。感受性保持者の初期割合を<math>s_0 = S(0) / N</math>、極限<math>t \to \infty</math>における感受性保持者および隔離者の割合をそれぞれ<math>s_\infty = S(\infty) / N</math>および<math>r_\infty = R(\infty) / N</math>と書くと。
:<math>s_\infty = 1 - r_\infty = s_0 e^{-R_0(r_\infty - r_0)}</math>
となる（ここで留意すべきは、感染区画はこの極限において空になる点である）。この{{仮リンク|超越方程式|en|Transcendental equation}}は[[ランベルトのW関数|ランベルト{{mvar|W}}関数]]の観点において解を持つ<ref>{{cite web |author1=Wolfram Research, Inc. |title=Mathematica, Version 12.1 |url=https://www.wolfram.com/mathematica |publisher=Champaign IL, 2020|accessdate=2020-06-26}}</ref>、すなわち
:<math>s_\infty = 1-r_\infty = - R_0^{-1}\, W(-s_0 R_0 e^{-R_0(1-r_0)})</math>
となる。これは、[[エピデミック]]の終わりには、<math>s_0=0</math>でない限り、集団のすべての個体が隔離されたわけではないので、一部は感受性のままであることを示している。これは、エピデミックの終わりは、感受性保持者の絶対的な不足ではなく、感染者の数が減少することによって引き起こされることを意味する。

[[基本再生産数]]の役割は極めて重要である。実際、感染者に関する方程式を
:<math> \frac{dI}{dt} = \left(R_0 \frac{S}{N}  - 1\right) \gamma I</math>
と書き換えると、もし
:<math> R_{0} > \frac{N}{S(0)}</math>
ならば、
:<math> \frac{dI}{dt}(0) >0 </math>
となる。すなわち、（人口のかなりの割合に達することができる）感染者数の増加を伴う相応の[[エピデミック]]な[[アウトブレイク]]があるだろう。逆に、もし
:<math> R_{0} < \frac{N}{S(0)}</math>
ならば、
:<math> \frac{dI}{dt}(0) <0 </math>
となる。すなわち、感受性保持者集団の初期の規模とは無関係に、この病気が相応のエピデミックなアウトブレイクを引き起こすことはあり得ない。結果として、基本再生産数が極めて重要であることは明らかである。

===感染力===
ここで留意すべきは、上述したモデルにおいて、関数
:<math> F = \beta I</math>
は、感受性保持者の区画から感染者の区画への遷移速度をモデル化し、これを感染力と呼ぶ。しかしながら、伝染病の多くのグループについては、感染力は感染者の絶対数に依存せず、（一定総人口に対する<math>N</math>に対する）感染者の割合に依存すると考える方がより現実的である。
:<math> F = \beta \frac{I}{N}</math>

Capassoと、その後、他の著者らは、伝染過程をより現実的にモデル化するために、非線形感染力を提案した{{Citation needed|date=June 2020}}。

===SIRモデルの厳密な解析解===
2014年は、Harkoと共著者らは、SIRモデルの厳密な解析解を導出した<ref name="Harko" />。人口動態が設定されていない場合、<math>\mathcal{S}(u)=S(t)</math>などは、
:<math>t= \frac{N}{\beta}\int_u^1 \frac{du^*}{u^*\mathcal{I}(u^*)} , \quad \rho=\frac{\gamma N}{\beta},</math>
に関して、次のような時間パラメータ化に対応する。
:<math>\mathcal{S}(u)= S(0)u </math>
:<math>\mathcal{I}(u)= N -\mathcal{R}(u)-\mathcal{S}(u) </math>
:<math>\mathcal{R}(u)=R(0) -\rho \ln(u)</math>
初期条件は
:<math>t= \frac{N}{\beta}\int_u^1 \frac{du^*}{u^*\mathcal{I}(u^*)} , \quad \rho=\frac{\gamma N}{\beta}</math>
であり、<math>u_T</math>は<math>\mathcal{I}(u_T)=0</math>を満たす。

上述した<math>R_{\infty}</math>に関する超越方程式によって、<math>u_T=e^{-(R_{\infty}-R(0))/\rho}(=S_{\infty}/S(0)</math> if <math>S(0) \neq 0)</math>および<math>I_{\infty}=0</math>となる。

Millerによって見出された等価な解析解<ref>{{cite journal | author = Miller, J.C. | title = A note on the derivation of epidemic final sizes | journal = Bulletin of Mathematical Biology | volume = 74 | issue = 9 | year=2012 | at= section 4.1| doi = 10.1007/s11538-012-9749-6 | pmid = 22829179 | pmc = 3506030 }}</ref><ref>{{cite journal | author = Miller, J.C. | title = Mathematical models of SIR disease spread with combined non-sexual and sexual transmission routes | journal = Infectious Disease Modelling | volume = 2 | issue = 1 | year=2017 | at = section 2.1.3| doi = 10.1016/j.idm.2016.12.003 | pmid = 29928728 | pmc = 5963332 }}</ref>により
:<math>
\begin{align}
S(t) & = S(0) e^{-\xi(t)} \\[8pt]
I(t) & = N-S(t)-R(t) \\[8pt]
R(t) & = R(0) + \rho \xi(t) \\[8pt]
\xi(t) & = \frac{\beta}{N}\int_0^t I(t^*) \, dt^*
\end{align}
</math>
が得られる。

ここで、<math>\xi(t)</math>は、個人が時刻<math>t</math>までに受け取った伝染の期待値と解釈することができる。これらの2つの解は<math>e^{-\xi(t)} = u</math>でつながっている。

事実上、同じ結果はカーマックとマッケンドリックの最初の論文にも見出すことができる<ref name="Kermack–McKendrick" />。

これらの解は、元の微分方程式の右辺のすべての項が<math>I</math>に比例していることに注意すれば、簡単に理解できるかもしれない。したがって、方程式を<math>I</math>で割って、左辺の微分演算子が単純に<math>d/d\tau</math>（<math>d\tau=I dt</math>、すなわち<math>\tau=\int I dt</math>）になるように、時間を再スケーリングすることができる。これで微分方程式はすべて線形となり、<math>dR/d\tau =</math> 定数の形式の第3の方程式は、<math>\tau</math>と<math>R</math>（および上記の<math>\xi</math>）が単純に線形に関係していることを示している。

===人口動態があり人口一定のSIRモデル===
死亡率<math>\mu</math>および出生率<math>\Lambda</math>によって特徴付けられ、伝染病が蔓延している母集団を考えてみる。質量作用伝播を持つモデルは
:<math>
\begin{align}
\frac{dS}{dt} & = \Lambda - \mu S - \frac{\beta I S}{N} \\[8pt]
\frac{dI}{dt} & = \frac{\beta I S}{N} - \gamma I -\mu I \\[8pt]
\frac{dR}{dt} & = \gamma I  - \mu R
\end{align}
</math>
であり、このモデルに対する感染症の無い平衡解（disease-free equilibrium: DFE）は
:<math>\left(S(t),I(t),R(t)\right) =\left(\frac{\Lambda}{\mu},0,0\right)</math>
である。

この場合、[[基本再生産数]]
:<math> R_0 = \frac{ \beta\Lambda }{\mu(\mu+\gamma)} </math>
を導出することができ、これは閾値特性を有さない。実際、生物学的に意味を持つ初期値から独立に、以下を示すことができる。
:<math> R_0 \le 1 \Rightarrow  \lim_{t \to \infty} (S(t),I(t),R(t)) = \textrm{DFE} = \left(\frac{\Lambda}{\mu},0,0\right) </math>
:<math> R_0 > 1 , I(0)> 0 \Rightarrow  \lim_{t \to \infty} (S(t),I(t),R(t)) = \textrm{EE} = \left(\frac{\gamma+\mu}{\beta},\frac{\mu}{\beta}\left(R_0-1\right), \frac{\gamma}{\beta} \left(R_0-1\right)\right) </math>

点EEは[[エンデミック]]な平衡解（Endemic Equilibrium）と呼ばれる（疫病は根絶されず、母集団に留まる）。{{仮リンク|ヒューリスティックな論証|en|Heuristic argument}}によって、<math>R_{0}</math>は完全に感受性な集団の中で1人の感染者によって引き起こされた平均感染数として読めることを示すことができる。上記の関係は、生物学的には、この数が1以下であれば病気は絶滅し、1以上であれば病気は永久に集団の中で流行し続けることを意味する。

==SISモデル==
[[File:SIS System Graph.svg|thumb|黄色=感受性, えび茶色=感染]]
[[風邪]]や[[インフルエンザ]]などの感染症の中には、長期的な[[獲得免疫|免疫]]が得られないものがある。このような感染症は、感染から回復しても免疫力が得られず、再び感受性となる。

[[File:SIS compartment model.svg|400px|center|SIS区画モデル]]

モデルは
:<math>
\begin{align}
\frac{dS}{dt} & = - \frac{\beta S I}{N} + \gamma I \\[6pt]
\frac{dI}{dt} & = \frac{\beta S I}{N} - \gamma I
\end{align}
</math>
である。

ここで留意すべきは、総人口を''N''で表わすと、
:<math>\frac{dS}{dt} + \frac{dI}{dt} = 0 \Rightarrow S(t)+I(t) = N</math>
が成り立つことである。すると、
:<math> \frac{dI}{dt} = (\beta - \gamma)  I - \frac{\beta}{N} I^2 </math>
となる。すなわち、感染の動態は[[ロジスティック関数]]で支配されている。<math>\forall I(0) > 0</math>:
:<math>
\begin{align}
& \frac{\beta}{\gamma} \le 1 \Rightarrow \lim_{t \to +\infty}I(t)=0, \\[6pt]
& \frac{\beta}{\gamma} > 1 \Rightarrow \lim_{t \to +\infty}I(t) = \left(1 - \frac{\gamma}{\beta} \right) N
\end{align}
</math>
基本再生産数が1よりも大きいような、（変数の変換<math>I = y^{-1}</math>を行い、これを平均場方程式に代入することによって）このモデルの解析的解を見つけることができる<ref name=Hethcote1989>{{cite book |last=Hethcote |first=Herbert W. |chapter=Three Basic Epidemiological Models |title=Applied Mathematical Ecology |editor1-last=Levin |editor1-first=Simon A. |editor2-last=Hallam |editor2-first=Thomas G. |editor3-last=Gross |editor3-first=Louis J. |series=Biomathematics |volume=18 |publisher=Springer |location=Berlin |year=1989 |pages=119–144 |isbn=3-540-19465-7 |doi=10.1007/978-3-642-61317-3_5 }}</ref>。解は以下の通りである。
:<math>I(t) = \frac{I_\infty}{1+V e^{-\chi t}}</math>
上式において、<math>I_\infty = (1 -\gamma/\beta)N</math>は伝染感染集団、<math>\chi = \beta-\gamma</math>、<math>V = I_\infty/I_0 - 1</math>である。この系は閉じていると仮定されるため、感受性保持者集団は<math>S(t) = N - I(t)</math>となる。

特殊な場合として、<math>\gamma=0</math>と仮定することによって普通のロジスティック関数が得られる。これは、<math>R=0</math>、すなわち隔離が起こらないSIRモデルにおいても考慮することができる。これは「SIモデル」である<ref>[http://math.colgate.edu/~wweckesser/math312Spring06/handouts/IMM_SI_Model.pdf (p. 19) The SI Model]</ref>。<math>S=N-I</math>を用いる微分方程式はしたがって、
:<math>
\frac{dI}{dt}  \propto  I\cdot (N-I).
</math>
と簡約される。

長期的には、このモデルでは、すべての個体が感染することになる。

==基本的なSIRモデルの精緻化==

===SIRDモデル===
[[File:SIRD.svg|thumb|SIRDモデルの模式図。初期値<math>S(0)=997, I(0)=3, R(0)=0</math>、感染率<math>\beta=0.4</math>、回復率<math>\gamma=0.035</math>、死亡率<math>\mu=0.005</math>。]]
[[File:SIRD model anim.gif|thumb|SIRDモデルのアニメーション。初期値<math display="inline">S(0)=997, I(0)=3, R(0)=0</math>、感染率初期値<math display="inline">\beta=0.5</math>、回復率（一定）<math display="inline">\gamma=0.035</math>、死亡率（一定）<math display="inline">\mu=0.005</math>。利用可能な薬も予防接種もない場合は、適切な対策（例えば「[[社会距離拡大戦略]]（ソーシャル・ディスタンシング）」など）によって感染率を下げる（「カーブを平らにする」としばしば言われる）ことしかできない。このアニメーションは、感染率を76%減少させた場合（<math display="inline">\beta=0.5</math>から<math display="inline">\beta=0.12</math>に低減）の効果を示している。]]

「感受性-感染-回復-死亡モデル」（Susceptible-Infectious-Recovered-Deceased-Model）は、「回復」（具体的には、病気を生き延びて免疫を得た個体を意味する）と「死亡」を区別する<ref>{{cite journal |author=Amenaghawon Osemwinyen, Aboubakary Diakhaby |title=Mathematical Modelling of the Transmission Dynamics of Ebola Virus |url=https://www.researchgate.net/publication/280141961 |journal=Applied and Computational Mathematics |volume=4 |pages=313–320 |year=2015 |issue=4 |doi=10.11648/j.acm.20150404.19 |doi-access=free }}</ref>。このモデルでは、以下の微分方程式の系を使用する。
:<math>
\begin{align}
& \frac{dS}{dt} = - \frac{\beta I S}{N}, \\[6pt]
& \frac{dI}{dt} = \frac{\beta I S}{N} - \gamma I - \mu I, \\[6pt]
& \frac{dR}{dt} = \gamma I, \\[6pt]
& \frac{dD}{dt} = \mu I,
\end{align}
</math>
上式において、<math>\beta, \gamma, \mu</math>はそれぞれ感染率、回復率、死亡率である。

===MSIRモデル===
[[麻疹]]を含む多くの感染症では、新生児は感受性区間へは入らず、生後数か月間は母体由来の抗体（[[胎盤]]を介して、さらに[[初乳]]を介して伝達される）による保護のために、病気に対する免疫を持っている。これは[[受動免疫]]と呼ばれる。この追加の詳細は、モデルの最初にMクラス（母体由来の免疫）を含めることによって示すことができる。

[[File:MSIR.PNG|800px|center|MSIR compartmental model]]

これを数学的に示すために、追加の区間''M''(''t'') が追加される。これにより、次のような微分方程式が得られる。
:<math>
\begin{align}
\frac{dM}{dT} & = \Lambda - \delta M - \mu M\\[8pt]
\frac{dS}{dT} & = \delta M - \frac{\beta SI}{N} - \mu S\\[8pt]
\frac{dI}{dT} & = \frac{\beta SI}{N} - \gamma I - \mu I\\[8pt]
\frac{dR}{dT} & = \gamma I - \mu R
\end{align}
</math> 

===キャリア状態===
[[結核]]といった感染症にかかったことのある人の中には、完全に回復することはなく、病気自体には罹患していないものの、[[無症候性キャリア|感染症を抱え続けている]]人もいる。彼らはその後、再び感染区画に移動し、（結核のように）症状に苦しむかもしれない。または彼らは症状に苦しんでいない間、キャリア状態で他の人に感染し続けるかもしれない。最も有名な例は、22人を[[腸チフス]]に感染させた[[メアリー・マロン]]であろう。キャリア区画はCとラベル付けされる。

[[File:SIR with carrier model.png|A simple modification of previous image by Viki Male to make the word "Carrier" plainly visible.]]

===SEIRモデル===
{{main|SEIRモデル}}
多くの重要な感染症では、個体は感染しているが、まだ発症していない潜伏期間が存在する。この期間中、個体は区画E（exposed）にいる。

[[File:SEIR.PNG|800px|center|SEIR compartmental model]]

潜伏期間がパラメータ<math>a</math>を持つ指数分布の確率変数であると仮定し（すなわち、平均潜伏期間は<math>a^{-1}</math>）、また、出生率<math>\Lambda</math>が死亡率<math>\mu</math>に等しい人口動態の存在を仮定すると、モデルは次のようになる。
:<math>
\begin{align}
\frac{dS}{dt} & = \mu N - \mu S - \beta \frac{I}{N} S \\[8pt]
\frac{dE}{dt} & = \beta \frac{I}{N} S - (\mu +a ) E \\[8pt]
\frac{dI}{dt} & = a E - (\gamma +\mu ) I \\[8pt]
\frac{dR}{dt} & = \gamma I  - \mu R.
\end{align}
</math>
ここで、<math>S+E+I+R=N</math>であるが、 これは、出生率と死亡率が等しいという（縮退した）仮定のために一定であるだけで、一般的に<math>N</math>は変数である。

このモデルについて、基本再生産数は
:<math>R_0 = \frac{a}{\mu+a}\frac{\beta}{\mu+\gamma}</math>
である。

SIRモデルと同様に、この場合においても、感染症の無い平衡解（DFE: ''N'',0,0,0）とエンデミックな平衡解（EE）があり、生物学的に意味を持つ初期条件とは独立に以下を示すことができる。
:<math> \left(S(0),E(0),I(0),R(0)\right) \in \left\{(S,E,I,R)\in [0,N]^4 : S \ge 0, E \ge 0, I\ge 0, R\ge 0, S+E+I+R = N \right\}   </math>

これは
:<math> R_0 \le 1 \Rightarrow \lim_{t \to +\infty} \left(S(t),E(t),I(t),R(t)\right) = DFE = (N,0,0,0),</math>
:<math> R_0 > 1 , I(0)> 0 \Rightarrow \lim_{t \to +\infty} \left(S(t),E(t),I(t),R(t)\right) = EE </math>
を満たす。

接触率<math>\beta(t)</math>が周期的に変化する場合、DFEの大域的誘引性<!--global attractiveness-->の条件は、周期的な係数を持つ以下の線形系
: <math>
\begin{align}
\frac{dE_1}{dt} & = \beta(t) I_1 - (\gamma +a ) E_1 \\[8pt]
\frac{dI_1}{dt} & = a E_1 - (\gamma +\mu ) I_1
\end{align}
</math>
が安定なことである（すなわち、複素平面における単位円内にフロケの固有値を持つ）。

=== SEISモデル ===
SEISモデルは、最後に免疫が獲得できないことを除けば、SEIRモデル（上）と同様である。
:<math>{\color{blue}{\mathcal{S} \to \mathcal{E} \to \mathcal{I} \to \mathcal{S}}}</math>

このモデルでは、感染は免疫を残さないので、回復した個体は再び感受性のある状態に戻り、''S''(''t''') 区画に後退する。以下の微分方程式はこのモデルを記述している。
:<math>
\begin{align}
\frac{dS}{dT} & = \Lambda - \frac{\beta SI}{N} - \mu S + \gamma I \\[6pt]
\frac{dE}{dT} & = \frac{\beta SI}{N} - (\epsilon + \mu)E \\[6pt]
\frac{dI}{dT} & = \varepsilon E - (\gamma + \mu)I
\end{align}
</math>

=== MSEIRモデル ===
受動免疫の要因と潜伏期間を持つ疾患の場合には、MSEIRモデルが存在する。
:<math> \color{blue}{\mathcal{M} \to \mathcal{S} \to \mathcal{E} \to \mathcal{I} \to \mathcal{R}} </math>
:<math>
\begin{align}
\frac{dM}{dT} & = \Lambda - \delta M - \mu M \\[6pt]
\frac{dS}{dT} & = \delta M - \frac{\beta SI}{N} - \mu S \\[6pt]
\frac{dE}{dT} & = \frac{\beta SI}{N} - (\varepsilon + \mu)E \\[6pt]
\frac{dI}{dT} & = \varepsilon E - (\gamma + \mu)I \\[6pt]
\frac{dR}{dT} & = \gamma I - \mu R
\end{align}
</math>

=== MSEIRSモデル ===
MSEIRSモデルはMSEIRに似ているが、Rクラスの免疫は一時的なものであり、一時的な免疫が終了すると個人は再び感受性を取り戻すことになる。
:<math>{\color{blue}{\mathcal{M} \to \mathcal{S} \to \mathcal{E} \to \mathcal{I} \to \mathcal{R} \to \mathcal{S}}}</math>

===変動する接触率、年複数回またはカオス的流行===
病気にかかる確率は時間的に一定ではないことはよく知られている。[[ヒトコロナウイルス|風邪ウイルス]]といった季節性のある病気もあり、これらは冬季に流行しやすくなる。[[麻疹]]、[[おたふくかぜ]]、[[風疹]]などの小児の病気では、[[学校暦]]との相関関係が強く、学校が休みの間はかかる確率が激減する。

結果として、多くの疾患では、周期的に（季節的に）変化する接触率
:<math> F = \beta(t) \frac{I}{N} , \quad \beta(t+T)=\beta(t)</math>
を持つ感染力を考慮する必要がある（周期Tは1年に等しい）。

したがって、モデルは
:<math>
\begin{align}
\frac{dS}{dt} & = \mu N - \mu S - \beta(t) \frac{I}{N} S \\[8pt]
\frac{dI}{dt} & = \beta(t) \frac{I}{N} S - (\gamma +\mu ) I
\end{align}
</math>
となる（回復の動態は簡単に<math>R=N-S-I</math>となる）。すなわち、周期的に変化するパラメータを持つ微分方程式の非線形集合である。このクラスの力学系は非線形パラメトリック共鳴という非常に興味深く複雑な現象を経ることが知られている。

もし
:<math>\frac 1 T \int_0^T \frac{\beta(t)}{\mu+\gamma} \, dt < 1 \Rightarrow \lim_{t \to +\infty} (S(t),I(t)) = DFE = (N,0), </math>
ならば容易に理解できるが、この積分が1よりも大きい場合、病気は消えず、こういった共鳴が存在するかもしれない。例えば、周期的に変化する接触率をシステムの「入力」として考えると、出力は周期が入力の周期の倍数である周期関数であることがわかる。これにより、いくつかの感染症の多年に1回の（通常は2年に1回の）流行の発生を、接触率の振動の周期と、エンデミックな平衡解近傍の減衰振動の疑似周期との間の相互作用として説明することができるようになった。驚くべきことに、いくつかの場合では、その挙動は準周期的なものであったり、カオス的なものであったりすることもある。

==予防接種のモデル化==
SIRモデルは、[[予防接種]]をモデル化するために修正することができる。典型的には、これらはSIRモデルに予防接種を受けた個人のための追加の区画<math>V</math>を導入する。以下、いくつかの例を示す{{cn|date=September 2019}}。

===新生児への予防接種===
伝染病が発生している状況において、主要な課題の1つは、予防策を講じ、可能であれば、集団予防接種体制の確立によって、伝染病を根絶することである。新生児に接種率<math>P \in (0,1)</math>で（生涯免疫を与えるワクチンを）予防接種した場合を考えてみる。

:<math>
\begin{align}
\frac{dS}{dt} & = \nu N (1-P) - \mu S - \beta \frac{I}{N} S \\[8pt]
\frac{dI}{dt} & = \beta \frac{I}{N} S - (\mu+\gamma)  I \\[8pt]
\frac{dV}{dt} & = \nu N P - \mu V
\end{align}
</math>
上式において、<math>V</math>は予防接種を受けた個人のクラスである。
:<math> \lim_{t \to +\infty} V(t)= N P</math>
であることはすぐに分かる、したがって<math>S</math>と<math>I</math>,の長期的な振る舞いを扱うことにする。これらは
:<math> R_0 (1-P) \le 1 \Rightarrow  \lim_{t \to +\infty} \left(S(t),I(t)\right) = DFE = \left(N \left(1-P\right),0\right) </math>
:<math> R_0 (1-P) > 1 , \quad I(0)> 0 \Rightarrow  \lim_{t \to +\infty} \left(S(t),I(t)\right) = EE = \left(\frac{N}{R_0(1-P)},N \left(R_0 (1-P)-1\right)\right) </math>
を満たす。

つまり
:<math> P < P^{*}= 1-\frac{1}{R_0}  </math>
ならば、予防接種計画は病気の根絶には成功せず、それどころか、ワクチン接種をしなかった場合よりは低いレベルではあるが、流行し続けることになる。つまり、この数理モデルは、[[基本再生産数]]が18ほど高い病気を根絶するためには、少なくとも94.4%の新生児にワクチンを接種する必要があることを示唆している。

===予防接種と情報===
現代社会は、「合理的な」免除、すなわち、感染からの知覚されるリスクとワクチンから損害を受けることによるリスクを「合理的」に比較した結果として、子供に予防接種受けさせないという家族の決断に直面している。この行動が本当に合理的であるかどうか、すなわち、それが等しく病気の根絶をもたらすことができるかどうかを評価するため、単に予防接種率が感染者数の増加関数
:<math> P=P(I), \quad P'(I)>0</math>
であると仮定することができる。

このような場合、根絶条件は
:<math> P(0) \ge P^{*}</math>
となる。すなわち、ベースライン予防接種率は、免除の場合は満たされない、「強制接種」の閾値以上でなければならない。このように、「合理的」な免除は、ワクチン接種率が高いために発生率が低い現状のみに基づいており、接種率の低下による将来の感染症再燃を考慮していないため、近視眼的である可能性がある。

===新生児以外への予防接種===
新生児以外にも率ρで予防接種を行う場合、感受性保持者と予防接種済みの人についての方程式は
:<math>
\begin{align}
\frac{dS}{dt} & = \mu N (1-P) - \mu S - \rho S - \beta \frac{I}{N} S \\[8pt]
\frac{dV}{dt} & = \mu N P + \rho S - \mu V
\end{align}
</math>
と修正されなければならず、ここから次のような根絶条件が得られる。
:<math> P \ge 1- \left(1+\frac{\rho}{\mu}\right)\frac{1}{R_0}  </math>

===パルス予防接種戦略===
{{main|パルス予防接種戦略}}

{{仮リンク|パルス予防接種戦略|en|Pulse vaccination strategy}}では、特定の年齢集団（幼児や高齢者など）にワクチンを接種することを繰り返す。この方法を用いると、感受性個体の集団は直ちに除去され、集団全体から（麻疹といった）感染症を排除することが可能となる。T時間単位毎に、感受性保持者の一定の割合pが（病気の動態に対して）比較的短い時間内に接種を受ける。これにより、感受性保持者と予防接種済みの集団について、次のようなインパルス<!--impulsive-->微分方程式が導かれる。
:<math>
\begin{align}
\frac{dS}{dt} & = \mu N  - \mu S - \beta \frac{I}{N} S, \quad S(n T^+) = (1-p) S(n T^-), & & n=0,1,2,\ldots \\[8pt]
\frac{dV}{dt} & = - \mu V, \quad V(n T^+) = V(n T^-) + p S(n T^-), & & n=0,1,2,\ldots
\end{align}
</math>

{{math|1=''I'' = 0}}と設定することで、感受性保持者の動態が
:<math> S^*(t) = 1- \frac{p}{1-(1-p)E^{-\mu T}} E^{-\mu MOD(t,T)}  </math>
で与えられ、根絶条件が
:<math> R_0 \int_0^T S^*(t) \, dt < 1  </math>
であることが容易にわかる。

==年齢の影響: 年齢構造化モデル==
年齢は集団における疾病伝播率、特に接触率に深く影響する。この接触率は、感受性保持者と感染者の接触の有効性を要約する。

:<math>S(t)=\int_0^{a_M} s(t,a)\,da </math>
:<math>I(t)=\int_0^{a_M} i(t,a)\,da</math>
:<math>R(t)=\int_0^{a_M} r(t,a)\,da</math>
（<math>a_M \le +\infty</math>は最大許容年齢）のように伝染病クラス<math>s(t,a),i(t,a),r(t,a)</math>の年齢を考慮に入れる（感受性-感染性-隔離スキームに限定するため）。それらの動態は「単純な」偏微分方程式によって記述されず、以下の[[積分微分方程式]]で記述される。
:<math>\partial_t s(t,a) + \partial_a s(t,a) = -\mu(a) s(a,t) - s(a,t)\int_0^{a_M} k(a,a_1;t)i(a_1,t)\,da_1 </math>
:<math>\partial_t i(t,a) + \partial_a i(t,a) =  s(a,t)\int_{0}^{a_M}{k(a,a_1;t)i(a_1,t)da_1}  -\mu(a) i(a,t) - \gamma(a)i(a,t)  </math>
:<math>\partial_t r(t,a) + \partial_a r(t,a) =   -\mu(a) r(a,t) + \gamma(a)i(a,t)  </math>
上式において、
:<math>F(a,t,i(\cdot,\cdot))=\int_0^{a_M} k(a,a_1;t)i(a_1,t) \, da_1 </math>
は感染力であり、これは、当然、年齢間の相互作用に依存する。

新生児の初期条件（すなわち、a=0の場合）によって複雑さが付加される。この初期条件は感染区分および隔離区分については
:<math>i(t,0)=r(t,0)=0</math>
と明快であるが、感受性を持つ新生児の密度については非局所的である。
:<math>s(t,0)= \int_0^{a_M} \left(\varphi_s(a) s(a,t)+\varphi_i(a) i(a,t)+\varphi_r(a) r(a,t)\right) \, da </math>
上式において、<math>\varphi_j(a), j=s,i,r</math>は大人の出生率である。

さらに、現在の総人口<math>n(t,a)=s(t,a)+i(t,a)+r(t,a)</math>の密度を定義すると、
:<math>\partial_t n(t,a) + \partial_a n(t,a) = -\mu(a) n(a,t) </math>
が得られる。

3つの伝染病クラスにおいて出生率が等しい最も単純な場合において、人口均衡を保つためには、出生率<math>\varphi(.)</math>と死亡率<math>\mu(a)</math>を結び付ける以下の必要十分条件
:<math> 1 = \int_0^{a_M} \varphi(a) \exp\left(- \int_0^a{\mu(q)dq} \right) \, da </math>
と人口均衡
:<math>n^*(a)=C \exp\left(- \int_0^a \mu(q) \, dq \right),</math>
が満たされなけばならず、これによって自動的に感染症の無い解<!--disease-free solution-->
:<math>DFS(a)= (n^*(a),0,0).</math>
の存在が保証される。

基本再生産数は、適切な関数演算子のスペクトル半径として計算することができる。

== 区画感染症モデル内におけるその他の考慮事項 ==

=== 垂直感染 ===
[[AIDS]]や[[B型肝炎]]といった一部の病気の場合、感染した両親の子どもが感染した状態で生まれることがある。このような母親からの伝染は垂直感染と呼ばれる。 感染区画への追加メンバーの流入は、感染区画に新生児の一部を含めることで、モデル内で考慮することができる<ref name=":2">{{cite book |last1=Brauer |first1=F. |last2=Castillo-Chávez |first2=C. |year=2001 |title=Mathematical Models in Population Biology and Epidemiology |location=NY |publisher=Springer |isbn=0-387-98902-1 }}</ref>。

=== 媒介生物感染 ===
ヒトからヒトへ間接的に伝染する病気（例えば、[[マラリア]]は蚊が媒介する）は媒介生物を介して伝播される。これらの場合、感染はヒトから[[昆虫]]へと移行し、伝染モデルは両方の種を含まなければならず、一般的には直接感染のモデルよりも多くの区画を必要とする<ref name=":2" /><ref>For more information on this type of model see {{cite book |editor-last=Anderson |editor-first=R. M. |editor-link=Roy M. Anderson |year=1982 |title=Population Dynamics of Infectious Diseases: Theory and Applications |publisher=Chapman and Hall |location=London-New York |isbn=0-412-21610-8 }}</ref>。

=== その他 ===
伝染病をモデル化する際に考慮する必要があるその他の事象としては、以下のようなものがある<ref name=":2" />。
* 非一様混合
* 変動のある感染力
* 空間的に不均一な分布
* マクロパラサイトによって引き起こされる疾患

==決定論的感染症モデル対確率的感染症モデル==
ここで提示された決定論的モデルは、十分に大きな母集団の場合にのみ有効であることを強調することが重要であり{{Citation needed|date=June 2020}}、そのようなものとして慎重に使用されるべきである<ref>{{cite journal |author=Bartlett MS |s2cid=91114210 |title=Measles periodicity and community size |journal=Journal of the Royal Statistical Society, Series A |volume=120 |pages=48–70 |year=1957 |doi=10.2307/2342553 |issue=1 |jstor=2342553}}</ref>。

より正確には、これらのモデルは[[熱力学的極限]]においてのみ有効であり、ここでは母集団は実質的に無限大である。確率モデルでは、システム内で感染者の数が1人以下になる確率が有限であるため、上述した長期的なエンデミックな平衡解は成り立たない。真のシステムでは、感染する宿主がいないので、病原体は伝播しないかもしれない。しかし、決定論的平均場モデルでは、感染者数は実数、すなわち、感染した宿主の非整数値を取ることができ、モデル内の宿主の数は1以下であっても、0以上であってもよく、それによってモデル内の病原体が伝播することが可能となる。

==出典==
{{reflist}}

==推奨文献==
*{{cite book |last1=May |first1=Robert M. |last2=Anderson |first2=Roy M.|title=Infectious diseases of humans: dynamics and control |publisher=Oxford University Press |location=Oxford |year=1991 |pages= |isbn=0-19-854040-X }}
*{{cite book |first=V. |last=Capasso |title=Mathematical Structure of Epidemic Systems |location=Berlin |publisher=Springer |year=1993 |isbn=3-540-56526-4 }}
*{{cite book |editor1-last=Vynnycky |editor1-first=E. |editor2-last=White |editor2-first=R. G. |title=An Introduction to Infectious Disease Modelling |publisher=Oxford University Press |location=Oxford |year=2010 |isbn=978-0-19-856576-5 }}

==関連項目==
*{{仮リンク|感染症の数理モデル化|en|Mathematical modelling of infectious disease}}
*[[可変単位地区問題]]
*[[次世代行列]]
*[[リスクアセスメント]]
*[[罹患率]]
*[[モデリング (科学的)]]

==外部リンク==
* [http://jsxgraph.uni-bayreuth.de/wiki/index.php/Epidemiology:_The_SIR_model SIR model: Online experiments with JSXGraph]
*{{YouTube|gxAaO2rsdIs|title=Simulating an epidemic}}

{{DEFAULTSORT:えきかくにおけるくかくもてる}}
[[Category:疫学|くかくもてる]]
[[Category:科学モデリング]]