[論文レビュー] Decay towards the overall-healthy state in SIS epidemics on networks

Piet Van Mieghem|arXiv (Cornell University)|Oct 15, 2013

Complex Network Analysis Techniques参考文献 19被引用数 23

ひとこと要約

本稿は、確率的三重対角行列の第二固有値を代数的手法で計算する方法を開発し、ネットワーク上のSIS疫病における崩壊率の正確な推定を可能にする。大規模Nに対して、疫病閾値を超えた場合の消滅までの期待時間はO(e^{N ln(τ/τ_c)})に比例することを証明し、従来の推定を精緻化し、多様なネットワーク構造において指数的寿命スケーリングが成立することを確認する。

ABSTRACT

The decay rate of SIS epidemics on the complete graph $K_{N}$ is computed analytically, based on a new, algebraic method to compute the second largest eigenvalue of a stochastic three-diagonal matrix up to arbitrary precision. The latter problem has been addressed around 1950, mainly via the theory of orthogonal polynomials and probability theory. The accurate determination of the second largest eigenvalue, also called the \emph{decay parameter}, has been an outstanding problem appearing in general birth-death processes and random walks. Application of our general framework to SIS epidemics shows that the maximum average lifetime of an SIS epidemics in any network with $N$ nodes is not larger (but tight for $K_{N}$) than \[ E\left[ T ight] \sim\frac{1}δ\frac{\fracτ{τ_{c}}\sqrt{2π}% }{\left( \fracτ{τ_{c}}-1 ight) ^{2}}\frac{\exp\left( N\left\{ \log\fracτ{τ_{c}}+\frac{τ_{c}}τ-1 ight\} ight) }{\sqrt {N}}=O\left( e^{N\ln\fracτ{τ_{c}}} ight) \] for large $N$ and for an effective infection rate $τ=\fracβδ$ above the epidemic threshold $τ_{c}$. Our order estimate of $E\left[ T ight] $ sharpens the order estimate $E\left[ T ight] =O\left( e^{bN^{a}} ight) $ of Draief and Massoulié \cite{Draief_Massoulie}. Combining the lower bound results of Mountford \emph{et al.} \cite{Mountford2013} and our upper bound, we conclude that for almost all graphs, the average time to absorption for $τ>τ_{c}$ is $E\left[ T ight] =O\left( e^{c_{G}N} ight) $, where $c_{G}>0$ depends on the topological structure of the graph $G$ and $τ$.

研究の動機と目的

• NノードのネットワークにおけるSIS疫病の吸収（消滅）までの期待時間の鋭い上界を導出すること。
• 大きなネットワークに対して、確率的三重対角行列の第二固有値（崩壊パラメータ）を正確に計算するという長年の課題を解決すること。
• DraiefとMassouliéのO(e^{bN^a})という数量オーダーの推定を、大規模Nに対してより精密なO(e^{N ln(τ/τ_c)})の形に改善すること。
• ほとんどすべてのグラフに対して、平均消滅時間がO(e^{c_G N})に比例し、c_G > 0がネットワークのトポロジーとτに依存することを確立すること。
• 完全グラフK_Nに対してその上界がタイトであることを示し、上界の鋭さを検証すること。

提案手法

• 任意の精度で確率的三重対角行列の第二固有値を計算する代数的手法を適用する。
• 1950年前後ごろに発展した直交多項式および確率論の結果を活用して固有値問題を解く。
• FillとMicloの結果を用い、吸収時間Tを非ゼロ固有値に等しいレートを持つ独立な指数分布の和として表現する。
• ラプラス変換を用いた吸収時間分布の積分に関する境界を、漸近的解析を用いて導出する。
• 階乗のスターリングの近似およびポアソン型の項を含む和の漸近的展開を用いて、最終的な式を簡略化する。
• Mountfordたちはの上界と下界および新たな上界を組み合わせることで、E[T]の指数的スケーリング行動を結論づける。

実験結果

リサーチクエスチョン

• RQ1 有効感染率τが疫病閾値τ_cを超えるとき、Nノードのネットワーク上におけるSIS疫病の期待消滅時間E[T]の正確な数量オーダーは何か？
• RQ2 大規模ネットワークに対して、SIS過程の無限小生成行列の第二固有値を任意の精度で計算する方法は何か？
• RQ3 数量推定E[T] = O(e^{bN^a})を、大規模Nに対してより精密な指数的形に改善できるか？
• RQ4 完全グラフK_Nは、疫病閾値を超えた場合に、期待消滅時間を最大にするネットワークか？
• RQ5 消滅時間のスケーリング定数c_Gは、ネットワークGのトポロジカル構造にどのように依存するか？

主な発見

• Nノードおよびτ > τ_cのネットワーク上におけるSIS疫病の期待消滅時間E[T]は、O(e^{N ln(τ/τ_c)})で上界づけられ、従来のO(e^{bN^a})推定に比べて顕著に精緻化されている。
• 上界は完全グラフK_Nに対してタイトである。これは、K_NがすべてのNノードネットワークの中で最大の期待消滅時間を達成することを意味する。
• ほとんどすべてのグラフに対して、平均消滅時間はE[T] = O(e^{c_G N})に比例し、c_G > 0はネットワークのトポロジーとτに依存する。
• 導出は、ポアソン分布およびガンマ関数を含む積分および和の精密な漸近的解析に依拠している。
• 代数的および特殊関数技法を用いて、生成行列の崩壊パラメータ（第二固有値）を高い精度で計算する手法は、成功裏に実行された。
• E[T]の最終的式はスターリングの近似およびm = N/xを中心にした和の漸近的評価を経て導出され、Nに伴う指数的スケーリングが確認された。