[論文レビュー] A graph theoretic interpretation of the mean first passage times

Pavel Chebotarev|arXiv (Cornell University)|Jan 12, 2007

Graph theory and applications参考文献 5被引用数 9

ひとこと要約

本稿は、定常マルコフ連鎖における平均通過時間（MFPT）を、連鎖の有向重み付きグラフ表現における重み付き森の比として表現することにより、グラフ理論的解釈を提供する。$ i \neq j $ に対して $ m_{ij} = f_{ij}/q_j $ が成り立つことが示され、ここで $ f_{ij} $ は $ i $ が1つの木に含まれ、もう一方の木が $ j $ に収束する2木のインフォレストの重み合計であり、$ q_j $ は $ j $ に収束する全域木の重み合計である。また $ m_{jj} = 1/\tilde{q}_j $ である。

ABSTRACT

Let $m_{ij}$ be the mean first passage time from state $i$ to state $j$ in an $n$-state ergodic homogeneous Markov chain with transition matrix $T$. Let $G$ be the weighted digraph without loops whose vertex set coincides with the set of states of the Markov chain and arc weights are equal to the corresponding transition probabilities. We give a graph-theoretic interpretation to $m_{ij}$. Namely, We show that $m_{ij}=f_{ij}/q_j$ if $i e j$ and $m_{ij}=1/ ilde q_j$ if $i=j$, where $f_{ij}$ is the total weight of 2-tree spanning converging forests in $G$ that have one tree containing $i$ and the other tree converging to $j$, $q_j$ is the total weight of spanning trees converging to $j$ in $G$, and $ ilde q_j=q_j/\sum_{k=1}^nq_k$. The result is illustrated by an example. Keywords: Markov chain; Mean first passage time; Spanning rooted forest; Matrix forest theorem; Laplacian matrix

研究の動機と目的

• 定常マルコフ連鎖における平均通過時間（MFPT）の組合せ的・グラフ理論的解釈を確立すること。
• マルコフ連鎖の有向グラフ表現において、MFPTを重み付きインフォレストと全域木を用いて表現すること。
• 特定の2木インフォレストおよび全域木の重み合計を用いて、MFPTの閉形式表現を提供すること。
• 具体的な計算例を通じて理論的枠組みの妥当性を検証すること。

提案手法

• マルコフ連鎖を、状態を頂点とし、弧の重みを遷移確率とする重み付き有向グラフ $ G $ として表現する。
• 2木のインフォレストのうち、1つの木に状態 $ i $ が含まれ、もう一方の木が状態 $ j $ に収束するものの重み合計を $ f_{ij} $ と定義する。
• $ j $ に収束する $ G $ の全域木の重み合計を $ q_j $ と定義する。
• 行列森定理およびラプラシアン行列 $ L = I - T $ の群逆行列を用いて、MFPTの式を導出する。
• 恒等式 $ m_{ij} = f_{ij}/q_j $（$ i \neq j $ の場合）および $ m_{jj} = 1/\tilde{q}_j $（ここで $ \tilde{q}_j = q_j / \sum_k q_k $）を適用する。
• $ n = 4 $ の数値例を用いて、インフォレストおよび全域木からの $ f_{ij} $、$ q_j $、$ m_{ij} $ の計算を通じて結果を検証する。

実験結果

リサーチクエスチョン

• RQ1 定常マルコフ連鎖における平均通過時間は、関連する重み付き有向グラフの組合せ的不変量を用いて表現可能か？
• RQ2 状態 $ i $ から $ j $ へのMFPTに対応する正確なグラフ理論的構造は何か？
• RQ3 2木インフォレストおよび全域木の重みは、MFPT行列とどのように関係するか？
• RQ4 グラフ表現におけるインフォレストおよび木の重みに基づいて、MFPTの閉形式表現が存在するか？
• RQ5 行列逆行列を用いずに、特定のインフォレストおよび全域木の重み合計からMFPT行列を再構成可能か？

主な発見

• $ i \neq j $ の場合、平均通過時間 $ m_{ij} $ は、$ i $ が1つの木に含まれ、もう一方の木が $ j $ に収束する2木インフォレストの重み合計 $ f_{ij} $ と、$ j $ に収束する全域木の重み合計 $ q_j $ の比に等しい。
• 対角成分 $ m_{jj} $ は $ 1/\tilde{q}_j $ で与えられ、ここで $ \tilde{q}_j = q_j / \sum_{k=1}^n q_k $ であり、$ q_j $ は $ j $ に収束する全域木の重み合計である。
• 例として $ n = 4 $ を用いた場合、グラフ理論的公式により計算されたMFPT行列 $ M $ は、群逆行列公式により得られた行列と一致する：$ M = \left[\begin{smallmatrix}12.5&1&6&14.5\\11.5&2.5&5&13.5\\6.5&7.5&5&8.5\\10.5&11.5&4&3.125\end{smallmatrix}\right] $。
• 各状態に収束する全域木の重み合計は $ q = (0.02, 0.1, 0.05, 0.08) $ であり、$ \tilde{q} = (0.08, 0.4, 0.2, 0.32) $ である。これは定常分布 $ \pi $ と一致する。
• $ f_{ij} $ に寄与する2木インフォレストは明示的に列挙されており、それらの重み合計は正しいMFPTをもたらす。例えば $ f_{14} = 1.16 $、$ f_{21} = 0.23 $ などである。
• $ Q_k $ 行列を用いた再帰的手順により、$ f_{ij} $ および $ q_j $ が正しく計算され、グラフ理論的アプローチが行列ベースの手法と整合的であることが確認された。