Aritalab:Lecture/NetworkBiology/Percolation on Graph

From Metabolomics.JP
< Aritalab:Lecture | NetworkBiology(Difference between revisions)
Jump to: navigation, search
m (占有される頂点の母関数)
m (占有される頂点の母関数)
Line 44: Line 44:
 
\begin{align}
 
\begin{align}
 
F'_v(1) &= q G'_v(1) = q \langle k \rangle \\
 
F'_v(1) &= q G'_v(1) = q \langle k \rangle \\
F''_v(1) &= q^2G''_v(1) = q^2 \langle k^2 \rangle - q \langle k \rangle
+
F''_v(1) &= q^2G''_v(1) - q G'_v(1) = q^2 \langle k^2 \rangle - q \langle k \rangle
 
\end{align}
 
\end{align}
 
</math>
 
</math>

Revision as of 09:33, 17 January 2014

Wiki Top Up one level レポートの書き方 Arita Laboratory

Contents

グラフ上のパーコレーション

グラフをある次数分布に従うツリーとみなすことで、パーコレーションを扱ってみます。 この内容を読む前に母関数の概念を復習しておいてください。

次数分布の母関数

グラフの次数分布 p(k) を与えられたとき、その母関数は

G_v(x) = \sum p(k) x^k

になります。しかし、頂点 v の隣にある頂点 w の次数分布はもはや p(k) ではありません。 選ばれた辺の先には、他の点が等確率でくることはないからです。 辺をランダムに1本選んだとき、もう片方にある端点がどんな次数分布をとるか考えましょう。 次数の大きい点のほうが、次数に比例して存在する確率が高くなります。 つまり、もう片側にある端点の次数分布は、最初の次数分布に辺の数で重みをつけてから正規化した\textstyle kp(k)/\Sigma_j jp(j)となります(分母は正規化定数)。 よってある頂点 v の隣にくる頂点 w の次数分布の母関数は

\sum_{k=0}^{\infty}x^k\frac{kp(k)}{\Sigma_jjp(j)} = x\frac{\Sigma_k kx^{(k-1)}p(k)}{\Sigma_jjp(j)} = x\frac{G'_v(x)}{G'_v(1)}

母関数においては、次数 k に対して x^k が対応します。 頂点 v の先に w がついているとき、w から v 以外に伸びる辺数に注目しましょう。 v から w への辺1本ぶんを除いた w に関する母関数は、上の式をxで割ったものにあたります。

G_w(x) = \frac{G'_v(x)}{G'_v(1)}

占有される頂点の母関数

パーコレーションで相転移が起きるのは、頂点から伸びていく辺の期待値が1をちょうど超えるときです。 ただし、このときの母関数はグラフ全体の次数分布を示す G_v ではなく、占有された頂点のみによる部分ネットワークの母関数です。 したがって、母関数を区別して F_v と記しましょう。関数 FG の関係は、各頂点が確率 q で占有されるため GF の次数平均をそれぞれ\langle k \rangle, \langle l \rangle と書くと 
\begin{align}
\langle l \rangle &= q \langle k \rangle \\
\end{align}
です。

この関係から (式の導出は確率母関数のページを参照)


\begin{align}
F'_v(1) &= q G'_v(1) = q \langle k \rangle \\
F''_v(1) &= q^2G''_v(1) - q G'_v(1) = q^2 \langle k^2 \rangle - q \langle k \rangle
\end{align}

相転移が起きるのはF_w'(1) = 1のときです。


\begin{align}
F'_w(1) &= q G'_w(1) = q \Big(\frac{G'_v(x)}{G'_v(1)}\Big)' \Big|_{x=1}
= q \frac{G''_v(1)}{G'_v(1)}
\end{align}

G''_v(1) = \langle k^2 \rangle, G'_v(1) = \langle k \rangleを代入すると 
q_c = \frac{1}{\frac{\langle k^2 \rangle}{\langle k \rangle} - 1}

すなわち、\langle k^2 \rangleが大きく \langle k \rangle を上回るときに q_c は小さくなります。 ハブの k は大きくなりますから、k^2 は更に大きくなります。次数の大きなハブがあるほど、パーコレーションが起きやすいということです。また全ての次数が同じとき、 \langle k^2 \rangle = \langle k \rangle^2 となりベーテ格子の結果と一致します。

パーコレーションの起こりやすさ

いかなる次数分布でも、\langle k^2 \rangle\langle k \rangleに比して大きいと、低い浸透確率で相転移が起こることをみてきました。

ランダムグラフの場合

次数の分布がポアソン分布になるランダムグラフをかんがえましょう。平均と分散は等しく E[k] = V[k] = \langle k \rangle になります。 このとき\langle k^2 \rangle = \langle k \rangle^2 + \langle k \rangle になるので

q_c = 1 / \langle k \rangle

すなわち、平均次数に反比例して浸透確率が下がっていきます。

スケールフリーネットワークの場合

次数の分布が、べき分布 p(k)=Ck^{-\gamma} のときをかんがえましょう。 ここで \zeta(\gamma)\ はリーマンゼータ関数とします。


\frac{\langle k^2 \rangle}{\langle k \rangle} = \frac{\Sigma_kk^2p(k)}{\Sigma_kkp(k)} = \frac{\Sigma k^{2-\gamma}}{ \Sigma k^{1-\gamma}} = \frac{\zeta(\gamma - 2)}{\zeta(\gamma - 1)}

この値は\gamma \leq 3のときに発散します。また \gamma が大きくなると  k=1 の項が大きな要素を占めるのでほとんど 1 に近づきます。


Personal tools
Namespaces

Variants
Actions
Navigation
metabolites
Toolbox