Aritalab:Lecture/NetworkBiology/Contact Process
m |
m |
||
Line 77: | Line 77: | ||
ここから、臨界確率<math>\lambda_c = \frac{\langle k \rangle}{\langle k^2 \rangle}</math>となる。 | ここから、臨界確率<math>\lambda_c = \frac{\langle k \rangle}{\langle k^2 \rangle}</math>となる。 | ||
+ | |||
+ | ===SIRモデル=== | ||
+ | |||
+ | <math>S(t) + I(t) + R(t) = 1</math> | ||
+ | |||
+ | ===ネットワークを考慮しない場合=== | ||
+ | 時間<math>\delta t</math>の間に感染者は<math>\mu</math>の確率で治癒し、健康人が<math>\lambda I(t)</math>の確率で感染するとする。 | ||
+ | |||
+ | <math> | ||
+ | \begin{align} | ||
+ | \frac{d S(t)}{dt} &= -\lambda I(t)S(t) \\ | ||
+ | \frac{d I(t)}{dt} &= \lambda S(t) I(t) - \mu I(t)\\ | ||
+ | \frac{d R(t)}{dt} &= \mu I(t) | ||
+ | \end{align} | ||
+ | </math> | ||
+ | |||
+ | 十分時間が経過した後の定常状態において十分量の<math>R(t)</math>が生じるためには | ||
+ | |||
+ | <math>\frac{d I(0)}{dt} = \lambda S(0) I(t) - \mu I(0) > 0</math> | ||
+ | |||
+ | <math> | ||
+ | \frac{\lambda}{\mu} = \frac{1}{S(t= 0)} | ||
+ | </math> | ||
+ | |||
+ | つまり感染率<math>\lambda</math>が治癒率<math>\mu</math>を超える場合は病人の割合が増加する。逆の場合は病人が減少する一方になる。 | ||
+ | |||
+ | SISモデルと同様、今後は一般性を失わずに<math>\mu=1</math>とおく。 | ||
+ | |||
+ | ===一般のネットワークの場合=== | ||
+ | |||
+ | 次数分布のみを考慮することにし、前出の式を次数$k$に限定して考える。 | ||
+ | |||
+ | <math> | ||
+ | \frac{d I_k(t)}{dt} = \lambda (1-I_k(t)) k\Theta(t) - I_k(t) | ||
+ | </math> | ||
+ | |||
+ | SISモデルとの違いは、いちど感染した頂点がSに戻らずRになる点である。 | ||
+ | よって | ||
+ | |||
+ | <math> | ||
+ | \begin{align} | ||
+ | \Theta(\infty) &= \frac{1}{\langle k \rangle} \sum_k' (k'-1) p(k') I_k (\infty) \\ | ||
+ | &= \frac{1}{\langle k \rangle} \sum_k' (k'-1) p(k') \frac{\lambda k \Theta(\infty)}{1 + \lambda k \Theta(\infty)} | ||
+ | \end{align} | ||
+ | </math> | ||
+ | |||
+ | ''k' ''が1少なくなっている理由はSに隣接する頂点がIの場合、その頂点を感染させた頂点はIかRであり、二度と感染しないからである。 | ||
+ | この式が<math>\theta = 0</math>以外にも解を持つかどうかの分岐点は | ||
+ | |||
+ | <math> | ||
+ | \frac{d y}{d \Theta} \Bigg|_{\Theta=0} = \frac{1}{\langle k \rangle} \sum_k (k-1)p(k)\lambda_c k = \frac{\langle k^2 \rangle - \langle k \rangle}{\langle k \rangle}\lambda_c = 1 | ||
+ | </math> | ||
+ | |||
+ | ここから、臨界確率<math>\lambda_c = \frac{\langle k \rangle}{\langle k^2 \rangle - \langle k \rangle}</math>となる。SISモデルとSIRモデルは基本的に結果が変わらない。 |
Revision as of 00:59, 10 June 2010
Contents |
コンタクトプロセス
感染症モデルでネットワーク構造を考慮したものをコンタクトプロセスと呼ぶ。 それぞれの状態の割合を
- S susceptible (健康状態)
- I infected (感染状態)
- R recovered (治癒状態)
と書く。の場合をSISモデルと呼び、の場合をSIRモデルと呼ぶ。
SIS model
ネットワークを考慮しない場合
時間の間に感染者はの確率で治癒し、健康人がの確率で感染するとする。
十分時間が経過した後の定常状態を考えるとを式変形して
つまり感染率が治癒率を超える場合は健康人の割合が減少し(病人が必ず残る)、感染率のほうが小さい場合はのために病人はゼロになる。
またとは定数倍すれば片方を消せるので、今後は一般性を失わずにとおく。
一般のネットワークの場合
次数分布のみを考慮することにし、前出の式を次数kに限定して考える。
ここでは、次数kの健康人が接続する先(1本あたり)が感染者である期待値を示す。定常状態のとき、であるから
は、次数kの頂点から出る辺が接続する先が感染者である確率を示していた。 辺が接続する先の頂点の次数分布はであるから
これをについて閉じた式にできれば感染率に対する感染者の期待値を解析的に求められるが、それは容易ではない。ここで求めたいのは臨界値を与えるだと考えて、右辺と左辺の関係を考える。
という連立方程式の解を考えよう。下の曲線はで定義され、が解の一つである。またのときにとなり単調増加でもある。
よって連立方程式が以外にも解を持つかどうかの分岐点は
ここから、臨界確率となる。
SIRモデル
ネットワークを考慮しない場合
時間の間に感染者はの確率で治癒し、健康人がの確率で感染するとする。
十分時間が経過した後の定常状態において十分量のが生じるためには
つまり感染率が治癒率を超える場合は病人の割合が増加する。逆の場合は病人が減少する一方になる。
SISモデルと同様、今後は一般性を失わずにとおく。
一般のネットワークの場合
次数分布のみを考慮することにし、前出の式を次数$k$に限定して考える。
SISモデルとの違いは、いちど感染した頂点がSに戻らずRになる点である。 よって
k' が1少なくなっている理由はSに隣接する頂点がIの場合、その頂点を感染させた頂点はIかRであり、二度と感染しないからである。 この式が以外にも解を持つかどうかの分岐点は
ここから、臨界確率となる。SISモデルとSIRモデルは基本的に結果が変わらない。