この会の企画・会場設備の提供をして頂きました
㈱ ワークスアプリケーションズ様
にこの場をお借りして御礼申し上げます.
前回は単純な変数消去法による推論アルゴリズムを紹介しましたが, これは効率が悪いです.
本日と来週の冒頭を使ってより効率的な厳密推論アルゴリズムであるジョインツリー法の解説を行います. この辺りもPRMLではあまり解説されていないので以下の参考書を用いました.
以下が前回紹介したアルゴリズムです.
クエリを Q, エビデンスを E とする. S={p(Xi|pa(Xi))}i とする. S 内の各テーブルについて E と一致しない行を 0 にする. Q に含まれない変数 X に対して以下を繰り返す.
S 内のファクターを全て掛け合わせて出力する.
このアルゴリズムは確率密度関数の因数分解 p(X)=∏ip(Xi|pa(Xi))
このアルゴリズムはネットワークのグラフ構造を一切考慮していないという点が問題です. グラフ理論的な考察によって計算量を削減する事が可能となります.
簡単に出来るのが 枝刈り (pruning)です.
前回の例において Q={D}, E={A} である場合を考えましょう.
まず, ノード E は D に関する推論とは全く無関係であるので除去出来ます.
p(A,D)=∑B,C,Ep(E|C)p(D|B,C)p(C|A)p(B|A)p(A)=∑B,Cp(D|B,C)p(C|A)p(B|A)p(A)∑EP(E|C) において ∑EP(E|C)=1 であるからです.
続いてエビデンスノード A も除去してしまう事が出来ます.
例えば, A の観測値が a であるならば p(A=a,D)=∑B,Cp(D|B,C)p(C|A=a)p(B|A=a)p(A=a) であり p(A=a)=1 となる為です.
クエリ集合 Q とエビデンス集合 E が与えられた時,
を除去する事が出来る. これを除去出来るノードがなくなるまで繰り返えす.
また, 枝刈りの結果孤立したQ,Eに含まれないノードも除去する事が出来る.
例えば, 右図のような状況(青がクエリ, 赤がエビデンス)ならば
こうなります.
prog15-1.py に枝刈りの実装を行ってみましたので参考にしてください.
変数消去法の計算量には変数の数 N はもちろんですが, 変数を消去する順番も関係します.
例として, p(A,B,C)=p(C|B)p(B|A)p(A) からA,Bを消去する場合, 消去順によって以下のような違いが出ます.
【A→Bの順に消した場合】
【B→Aの順に消した場合】
B→A の順に消去してしまった場合は φ(A,C) という因子が残り計算量が増加してしまっています.
一般に以下の事が言えます.
従って w が最小となるような変数消去順序を求めれば良いのですが, 残念ながらこの問題はNP困難である事が解っています.
そこで, ヒューリスティックスによる 最小次数法 (minimum degree method) を紹介します.
まず, インタラクショングラフ (interaction graph) というものが必要です. これは以下のように定義されます.
ファクターの集合 S={φ1,…,φn} に対するインタラクショングラフ Gi とは無向グラフであり, 頂点集合 Vと辺集合 E は以下のように定める.
V={φ1,…φnに含まれる各変数}E={(x,y)|x,y が同一のファクターに含まれる}前回やった p(A,B,C,D,E)=p(A)p(B)p(C|A,B)p(D|B)p(E|C,D) を例に考えてみます.
最初の S の状態は以下のようになるのでした S={p(A),p(B),p(C|A,B),p(D|B),p(E|C,D)}
そして, これと対応するインタラクショングラフは下図のようになります.
次にここから B を消去してみると φ(A,C,D)=∑Bp(B)p(C|A,B)p(D|B) より S′={p(A),φ(A,C,D),p(E|C,D)} となります.
これに対応するインタラクショングラフは以下です.
この推移は
という事によって行う事が出来ます.
この時, Xを消去して出来る因子の変数の数=XのGiでの次数 となります.
例えば最初に C を消去すると φ(A,B,D,E)=∑Cp(C|A,B)p(E|C,D) という大きな因子が出来てしまいます.
そこで,
というヒューリスティクス的な手法を考える事が出来ます. これが最小次数法です.
先ほどの問題で最小次数法を適用すると以下のように A→B→C→E などの順番で消去することになり, この場合 w=3 です.
ジョインツリーアルゴリズム (jointree algorithm は O(Nexp(w)) の計算量で厳密な推論を行う事が出来るアルゴリズムです.
ジョインツリーアルゴリズムは ファクター消去法 (factor elimination) というものの特別な場合ですので, まずこれから説明します.
例として再び p(A,B,C,D,E)=p(A)p(B)p(C|A,B)p(D|B)p(E|C,D) において p(D) を計算する問題を考えます.
変数消去法ではD 以外の変数を1つずつ消す事によってこれを計算しました.
ファクター消去法では D を含む適当な1つのファクター(例えば p(D|B)) 以外のファクターを1つずつ消す事によってこれを計算します.
やってみます. S={φ(A),φ(B),φ(A,B,C),φ(B,D),φ(C,D,E)} から φ(B,D) 以外を消去する事にします.
まず φ(A) を消去します. どれか別の適当なファクター(例えば φ(B))に φ(A) を掛けた φ(A,B)=φ(A)φ(B) を計算して φ(B) をこれで置換えます. S={φ(A,B),φ(A,B,C),φ(B,D),φ(C,D,E)}
次に φ(C,D,E) を消去してみましょう. まず E はこのファクターにしか含まれないので 消しちゃいます. ∑Eφ(C,D,E) これを別のどれか(例えば φ(B,D))に掛けた φ(B,C,D)=φ(B,D)∑Eφ(C,D,E) で φ(B,D) を置換えます. S={φ(A,B),φ(A,B,C),φ(B,C,D)}
同様に, 次は φ(A,B) を φ(A,B,C) に掛けて消去しましょう. φ′(A,B,C)=φ(A,B)φ(A,B,C) として S={φ′(A,B,C),φ(B,C,D)} とします.
次は φ′(A,B,C) を消しますが A はこのファクターにしか含まれないので消してしまいます. ∑Aφ′(A,B,C) そして φ(B,C,D) にかけてこのファクターを消します. つまり φ′(B,C,D)=φ(B,C,D)∑Aφ′(A,B,C) として S={φ′(B,C,D)} とします.
最後に残ったファクター φ′(B,C,D) には求める D 以外のファクター B,C が含まれていますがこれを足し合わせれば p(D)=∑B,Cφ′(B,C,D) となります. 以上でファクター消去法による確率計算が終わりました.
ファクター φ の変数 Q 以外を足しあわせて消去するという操作は良く使うので proj(φ,Q)def=∑varsφ−Qφ と書きます. 但し varsφ は φ に含まれる変数の集合です.
まとめると以下の様になります.
クエリ集合を Q とする.
ファクター消去法ではどの順番でファクターを消去するのかが重要となります. これを指定するものが エリミネーションツリー (elimination tree) です.
エリミネーションツリー T とは, 各ファクター φr をノードとする木構造(閉路のない単連結な無向グラフ)である.
エリミネーションツリーを用いたファクター消去法を簡単に書くと以下のようになります.
残すファクター φr をエリミネーションツリー T のルートとみなして以下を繰り返してファクターを消去する.
例えば先ほどの例のファクター消去手順に対応するエリミネーションツリーは下図のようになります.
ファクター消去法はどのようなエリミネーションツリーを用いても必ず正しい答えを得る事が出来ます.
そこで問題となるのが「どのようなエリミネーションツリーを用いれば効率的に計算出来るか?」という事です.
ジョインツリーアルゴリズムとはジョインツリーと呼ばれるものをエリミネーションツリーの構成に利用するものです.
もう少し準備が必要です.
エリミネーションツリー T の辺 (i,j) に対して Sijdef={Tのi 側にある変数}∩{Tのj側にある変数} をこの辺の セパレータ(separator) と呼びます.
ファクター φi を消去し φj にまとめる時には varsφi−Sij が φi のみに含まれる変数となります.
つまり, φj の更新ルールは φj←φjproj(φi,Sij) と書くことが出来ます.
エリミネーションツリー T のノード i に対して Cidef=vars(i)∪⋃jSij をこのノードに対応する クラスター (cluster) と呼びます.
クラスターとは, ファクターを消去していってそのノードが葉になった時に含まれる変数の集合です.
wdef=maxi|Ci|−1 をエリミネーションツリーの幅(width)と言います.
変数消去法の時と同様に w に関して指数的に計算量が増加します. 従って w を出来るだけ小さくする事が目標になります.
ファクター消去法では, エリミネーションツリーの辺に沿ってファクターを流して行くわけですが, これをメッセージ(message) と呼びます.
子ノード i から親ノード j へのメッセージは mij以下の様に定義されます. mijdef=proj(φi∏kmki,Sij) 但し k はノード i の子ノードです.
葉ノード i では子が存在しませんから mijdef=proj(φi,Sij) となります.
このメッセージをルートノード r まで流していった時 p(Cr)=φr∏kmkr となります. つまりクラスターの同時確率が計算されます.
これを周辺化すれば求める確率 p(Q)=proj(p(Cr),Q) を得る事が出来ます.
メッセージパッシングというアルゴリズム, ルートを変えて再計算したい時に, メッセージを再利用出来る というとても良い特徴を持っています.
つまり, 一回ルートまでメッセージを流して, その後ルートから逆にメッセージを再分配すれば一辺に全てのクラスター Ci に関する同時確率を計算してしまうという事が可能です. これによって計算量を削減する事が出来ます.
エリミネーションツリーのルート r を決める.
ノード i から j へのメッセージ mij は以下によって定義される. mijdef=proj(φi∏kmki,Sij) 但し, k はメッセージが流れて来た方向の隣接ノード.
次回前半はジョインツリーの解説を行いジョインツリーアルゴリズムを完成させます. 後半はマルコフ確率場の紹介をします.