Aritalab:Lecture/NetworkBiology/Contact Process
m (→ネットワークを考慮しない場合) |
m |
||
Line 95: | Line 95: | ||
<math>S(t) + I(t) + R(t) = 1</math> が成立する。 | <math>S(t) + I(t) + R(t) = 1</math> が成立する。 | ||
− | ==ネットワークを考慮しない場合== | + | ===ネットワークを考慮しない場合=== |
時間<math>\Delta t</math>の間に感染者は<math>\mu</math>の確率で治癒し、健康人が<math>\lambda I(t)</math>の確率で感染するとする。 | 時間<math>\Delta t</math>の間に感染者は<math>\mu</math>の確率で治癒し、健康人が<math>\lambda I(t)</math>の確率で感染するとする。 | ||
Line 106: | Line 106: | ||
</math> | </math> | ||
− | 十分時間が経過した後の定常状態において十分量の<math>R(t)</math> | + | 十分時間が経過した後の定常状態において十分量の<math>R(t)</math>が生じるためには、初期 (<math>t\sim 0</math>) から次第に病人が増えないと駄目。 |
− | + | ||
− | <math> | + | |
<math> | <math> | ||
− | \frac{\lambda}{\mu} = \frac{1}{S(t= 0)} | + | \begin{align} |
+ | \frac{d I(0)}{dt} &= \lambda S(0) I(t) - \mu I(0) > 0 \\ | ||
+ | \frac{\lambda}{\mu} &= \frac{1}{S(t= 0)} | ||
+ | \end{align} | ||
</math> | </math> | ||
Line 118: | Line 119: | ||
SISモデルと同様、今後は一般性を失わずに<math>\mu=1</math>とおく。 | SISモデルと同様、今後は一般性を失わずに<math>\mu=1</math>とおく。 | ||
− | ==一般のネットワークの場合== | + | ===一般のネットワークの場合=== |
− | 次数分布のみを考慮することにし、前出の式を次数 | + | 次数分布のみを考慮することにし、前出の式を次数''k''に限定して考える。 |
<math> | <math> | ||
Line 127: | Line 128: | ||
SISモデルとの違いは、いちど感染した頂点がSに戻らずRになる点である。 | SISモデルとの違いは、いちど感染した頂点がSに戻らずRになる点である。 | ||
− | + | 関数<math>\Theta</math>を求めるには、Sに隣接する頂点がIである割合を考える。 | |
+ | ある頂点がIのとき、必ずその頂点を感染させる原因となる頂点(IかR)が存在する。 | ||
+ | よって、隣接頂点の母関数を''kp(k)''の代わりに、大雑把に''(k-1)p(k)''と見積もろう。 | ||
<math> | <math> | ||
\begin{align} | \begin{align} | ||
− | \Theta(\infty) &= \frac{1}{\langle k \rangle} \ | + | \Theta(\infty) &= \frac{1}{\langle k \rangle} \sum_j (j-1) p(j) I_k (\infty) \\ |
− | &= \frac{1}{\langle k \rangle} \ | + | &= \frac{1}{\langle k \rangle} \sum_j (j-1) p(j) \frac{\lambda k \Theta(\infty)}{1 + \lambda k \Theta(\infty)} |
\end{align} | \end{align} | ||
</math> | </math> | ||
− | |||
この式が<math>\theta = 0</math>以外にも解を持つかどうかの分岐点は | この式が<math>\theta = 0</math>以外にも解を持つかどうかの分岐点は | ||
Revision as of 23:51, 13 June 2010
Wiki Top | Up one level | レポートの書き方 | Arita Laboratory |
|
コンタクトプロセス
感染症モデルでネットワーク構造を考慮したものをコンタクトプロセスと呼ぶ。 それぞれの状態の割合を
- S susceptible (健康状態)
- I infected (感染状態)
- R recovered (治癒状態)
と書く。の場合をSISモデルと呼び、の場合をSIRモデルと呼ぶ。
SIS model
ネットワークを考慮しない場合
時間の間に感染者はの確率で治癒し、健康人がの確率で感染するとする。
十分時間が経過した後の定常状態を考えるとを式変形して
つまり感染率が治癒率を超える場合は健康人の割合が減少し(病人が必ず残る)、感染率のほうが小さい場合はのために病人はゼロになる。
またとは定数倍すれば片方を消せるので、今後は一般性を失わずにとおく。
一般のネットワークの場合
次数分布のみを考慮することにし、前出の式を次数kに限定して考える。
ここでは、次数kの健康人が接続する先(1本あたり)が感染者である期待値を示す。定常状態のとき、であるから
は、次数kの頂点から出る辺が接続する先が感染者である確率を示していた。 辺が接続する先の頂点の次数分布はであるから
これをについて閉じた式にできれば感染率に対する感染者の期待値を解析的に求められるが、それは容易ではない。ここで求めたいのは臨界値を与えるだと考えて、右辺と左辺の関係を考える。
という連立方程式の解を考えよう。下の曲線はで定義され、が解の一つである。またのときにとなり単調増加でもある。
よって連立方程式が以外にも解を持つかどうかの分岐点は
ここから、臨界確率となる。
具体例
- 次数分布がポアソン分布
のとき
- 次数分布が指数分布
のとき
SIR model
が成立する。
ネットワークを考慮しない場合
時間の間に感染者はの確率で治癒し、健康人がの確率で感染するとする。
十分時間が経過した後の定常状態において十分量のが生じるためには、初期 () から次第に病人が増えないと駄目。
つまり感染率が治癒率を超える場合は病人の割合が増加する。逆の場合は病人が減少する一方になる。
SISモデルと同様、今後は一般性を失わずにとおく。
一般のネットワークの場合
次数分布のみを考慮することにし、前出の式を次数kに限定して考える。
SISモデルとの違いは、いちど感染した頂点がSに戻らずRになる点である。 関数を求めるには、Sに隣接する頂点がIである割合を考える。 ある頂点がIのとき、必ずその頂点を感染させる原因となる頂点(IかR)が存在する。 よって、隣接頂点の母関数をkp(k)の代わりに、大雑把に(k-1)p(k)と見積もろう。
この式が以外にも解を持つかどうかの分岐点は
ここから、臨界確率となる。SISモデルとSIRモデルは基本的に結果が変わらない。
SIRモデルとパーコレーション
SIRモデルはいったん感染すると後はRに移動するだけである。したがってにおける最終結果は各辺について一回ずつ感染か非感染かを評価するパーコレーションと変わらない。
このように、頂点ではなく辺に活性、非活性を割り当てる方式をボンド・パーコレーションと呼ぶ。