Statistical model
ベイジアン ネットワーク( ベイズネットワーク 、 ベイズネット 、 ビリーフネットワーク 、 決定ネットワーク とも呼ばれる )は、 一連の変数とその 条件付き依存関係を有向 非巡回グラフ (DAG) で表す 確率的グラフィカルモデルです。 [1] 因果表記法 のいくつかの形式の1つですが 、因果ネットワークはベイジアンネットワークの特殊なケースです。ベイジアンネットワークは、発生したイベントを取得し、考えられる複数の既知の原因のいずれかが要因である可能性を予測するのに最適です。たとえば、ベイジアンネットワークは、病気と症状の確率的な関係を表すことができます。症状が与えられれば、ネットワークを使用してさまざまな病気の存在確率を計算できます。
効率的なアルゴリズムは、ベイジアンネットワークにおける推論 と 学習を 実行できます 。変数のシーケンス(音声信号やタンパク質配列など)をモデル化するベイジアンネットワークは、 動的 ベイジアン ネットワーク と 呼ばれます 。 不確実性の下で意思決定問題を表現および解決できるベイジアンネットワークの一般化は、 影響図 と呼ばれます。
グラフィカルモデル
正式には、ベイジアン ネットワークは、 ノードが ベイジアンセンスの変数を表す 有向非巡回グラフ (DAG) です。変数とは、観測可能な量、 潜在変数 、未知のパラメーター、仮説などです。各エッジは、直接的な条件付き依存関係を表します。接続されていないノードのペア (つまり、1 つのノードを他のノードに接続するパスがない) は、条件付きで互いに独立した変数を表します 。 各 ノードには 、入力としてノードの 親 変数の特定の値セットを受け取り、(出力として) ノードによって表される変数の確率 (または該当する場合は確率分布) を返す確率関数が関連付けられています。たとえば、 親ノードが ブール変数 を表す場合、確率関数は、可能な親の組み合わせごとに 1 つのエントリがあるエントリ のテーブルで表すことができます。同様の考え方は、 マルコフ ネットワーク などの無向で、場合によっては巡回的なグラフにも適用できます 。
m
{\displaystyle m}
m
{\displaystyle m}
2
m
{\displaystyle 2^{m}}
2
m
{\displaystyle 2^{m}}
例
条件付き確率表 を用いた単純なベイジアンネットワーク
3つの変数間の依存関係をモデル化するとします。スプリンクラー(より正確には、スプリンクラーの状態、つまりオンになっているかどうか)、雨の有無、そして芝生が濡れているかどうかです。芝生が濡れる原因となるイベントは2つあります。スプリンクラーの作動と雨です。雨はスプリンクラーの使用に直接影響します(つまり、雨が降っているときはスプリンクラーは通常作動していません)。この状況は、ベイジアンネットワーク(右図)を使用してモデル化できます。各変数には、T(真)とF(偽)の2つの値があります。
確率の連鎖律 によれば、 結合 確率関数 は、
Pr
(
G
,
S
,
R
)
=
Pr
(
G
∣
S
,
R
)
Pr
(
S
∣
R
)
Pr
(
R
)
{\displaystyle \Pr(G,S,R)=\Pr(G\mid S,R)\Pr(S\mid R)\Pr(R)}
ここで、 G = 「芝生が濡れている (true/false)」、 S = 「スプリンクラーがオンになっている (true/false)」、 R = 「雨が降っている (true/false)」です。
このモデルは、条件付き確率の 式 を使用してすべての不要な変数を合計することにより、「草が濡れている場合、雨が降る確率はどれくらいですか?」のような結果の存在を前提とした 原因 の存在に関する質問(いわゆる逆確率)に答えることができます。
Pr
(
R
=
T
∣
G
=
T
)
=
Pr
(
G
=
T
,
R
=
T
)
Pr
(
G
=
T
)
=
∑
x
∈
{
T
,
F
}
Pr
(
G
=
T
,
S
=
x
,
R
=
T
)
∑
x
,
y
∈
{
T
,
F
}
Pr
(
G
=
T
,
S
=
x
,
R
=
y
)
{\displaystyle \Pr(R=T\mid G=T)={\frac {\Pr(G=T,R=T)}{\Pr(G=T)}}={\frac {\sum _{x\in \{T,F\}}\Pr(G=T,S=x,R=T)}{\sum _{x,y\in \{T,F\}}\Pr(G=T,S=x,R=y)}}}
図に示されている条件付き確率表(CPT) の条件付き確率と結合確率関数 の展開を用いて 、分子と分母の和の各項を評価することができます。例えば、
Pr
(
G
,
S
,
R
)
{\displaystyle \Pr(G,S,R)}
Pr
(
G
=
T
,
S
=
T
,
R
=
T
)
=
Pr
(
G
=
T
∣
S
=
T
,
R
=
T
)
Pr
(
S
=
T
∣
R
=
T
)
Pr
(
R
=
T
)
=
0.99
×
0.01
×
0.2
=
0.00198.
{\displaystyle {\begin{aligned}\Pr(G=T,S=T,R=T)&=\Pr(G=T\mid S=T,R=T)\Pr(S=T\mid R=T)\Pr(R=T)\\&=0.99\times 0.01\times 0.2\\&=0.00198.\end{aligned}}}
すると、数値結果(関連する変数値を添え字として)は
Pr
(
R
=
T
∣
G
=
T
)
=
0.00198
T
T
T
+
0.1584
T
F
T
0.00198
T
T
T
+
0.288
T
T
F
+
0.1584
T
F
T
+
0.0
T
F
F
=
891
2491
≈
35.77
%
.
{\displaystyle \Pr(R=T\mid G=T)={\frac {0.00198_{TTT}+0.1584_{TFT}}{0.00198_{TTT}+0.288_{TTF}+0.1584_{TFT}+0.0_{TFF}}}={\frac {891}{2491}}\approx 35.77\%.}
「芝生を濡らした場合、雨が降る確率はどれくらいですか?」といった介入の質問に答えるには、介入後の結合分布関数によって答えが決まります。
Pr
(
S
,
R
∣
do
(
G
=
T
)
)
=
Pr
(
S
∣
R
)
Pr
(
R
)
{\displaystyle \Pr(S,R\mid {\text{do}}(G=T))=\Pr(S\mid R)\Pr(R)}
介入前の分布から因子を除去することによって得られる 。do演算子はGの値を真とするように強制する。降雨確率はアクションの影響を受けない。
Pr
(
G
∣
S
,
R
)
{\displaystyle \Pr(G\mid S,R)}
Pr
(
R
∣
do
(
G
=
T
)
)
=
Pr
(
R
)
.
{\displaystyle \Pr(R\mid {\text{do}}(G=T))=\Pr(R).}
スプリンクラーをオンにした場合の影響を予測するには:
Pr
(
R
,
G
∣
do
(
S
=
T
)
)
=
Pr
(
R
)
Pr
(
G
∣
R
,
S
=
T
)
{\displaystyle \Pr(R,G\mid {\text{do}}(S=T))=\Pr(R)\Pr(G\mid R,S=T)}
用語を 削除すると、その動作は草には影響しますが、雨には影響しないことがわかります。
Pr
(
S
=
T
∣
R
)
{\displaystyle \Pr(S=T\mid R)}
これらの予測は、ほとんどの政策評価問題と同様に、観測されない変数を考慮すると実現可能ではないかもしれない。しかし、バックドア基準が満たされる限り、 行動の効果は予測可能である。 [2] [3] これは、 Xから Y への すべてのバックドアパスをd分離 [4] (またはブロック)するノードの集合 Z が観測できる 場合
、
do
(
x
)
{\displaystyle {\text{do}}(x)}
Pr
(
Y
,
Z
∣
do
(
x
)
)
=
Pr
(
Y
,
Z
,
X
=
x
)
Pr
(
X
=
x
∣
Z
)
.
{\displaystyle \Pr(Y,Z\mid {\text{do}}(x))={\frac {\Pr(Y,Z,X=x)}{\Pr(X=x\mid Z)}}.}
バックドアパスとは、 X への矢印で終わるパスです。バックドア基準を満たす集合は、「十分」または「許容可能」と呼ばれます。例えば、集合 Z = Rは、 S = Tが G に 与える影響を予測するのに許容されます 。なぜなら、 R は(唯一の)バックドアパス S ← R → Gを d 分離するからです。しかし、 S が観測されない場合、このパスを d 分離する他の集合はなく、スプリンクラーをオン( S = T )にした場合の芝生( G )への影響は、受動的な観測から予測できません。その場合、 P ( G | do( S = T ))は「特定」されません。これは、介入データがない場合、 S と G の間に観測される依存関係は因果関係によるものか、あるいは偽の依存関係(共通の原因 R に起因する見かけ上の依存関係 )であるかを反映しています。( シンプソンのパラドックス を参照)
観測されていない変数を持つ任意のベイジアンネットワークから因果関係が特定されるかどうかを判断するには、「 do 計算」の3つのルール [2] [5] を使用し、その関係の表現からすべての do 項を削除できるかどうかをテストすることで、目的の量が頻度データから推定可能であることを確認します。 [6]
ベイジアンネットワークを使用すると、結合分布の依存関係が疎である場合、網羅的な確率テーブルを使用するよりもメモリ使用量を大幅に節約できます。例えば、10個の2値変数の条件付き確率をテーブルとして保存するという単純な方法では、 値を格納するための記憶領域が必要になります。どの変数の局所分布も3つ以上の親変数に依存していない場合、ベイジアンネットワーク表現は最大で以下の 値しか保存しません。
2
10
=
1024
{\displaystyle 2^{10}=1024}
10
⋅
2
3
=
80
{\displaystyle 10\cdot 2^{3}=80}
ベイジアン ネットワークの利点の 1 つは、完全な結合分布よりも、(まばらな) 直接的な依存関係とローカル分布の方が人間にとって直感的に理解しやすいことです。
推論と学習
ベイジアン ネットワークは主に 3 つの推論タスクを実行します。
観測されていない変数の推論
ベイジアンネットワークは変数とその関係性に関する完全なモデルであるため、それらに関する確率的なクエリに答えるために使用できます。例えば、このネットワークは、他の変数( 証拠 変数)が観測されたときに、変数のサブセットの状態に関する知識を更新するために使用できます。証拠を与えられた場合に変数の 事後 分布を計算するこのプロセスは、確率的推論と呼ばれます。事後分布は、検出アプリケーションにおいて、変数サブセットの値を選択する際に、期待損失関数(例えば、意思決定エラーの確率)を最小化するのに 十分な普遍的な統計量を提供します。したがって、ベイジアンネットワークは 、ベイズの定理を 複雑な問題に自動的に適用するためのメカニズムと考えることができます 。
最も一般的な正確な推論方法は、 積分または加算によって、観測されない非クエリ変数を 1 つずつ、その合計を積に分配することによって除去する 変数除去 、一度に多くの変数をクエリして新しい証拠をすばやく伝播できるように計算をキャッシュする クリーク ツリー伝播、および十分なスペースが使用される場合に 空間と時間のトレードオフを 可能にして変数除去の効率に匹敵する再帰条件付けと AND/OR 検索です。これらの方法はすべて、ネットワークの ツリー幅 に比例して複雑になります。最も一般的な 近似推論 アルゴリズムは、 重要度サンプリング 、確率的 MCMC シミュレーション、ミニバケット除去、 ループ状確信伝播 、 一般化確信伝播 、 変分法 です。
パラメータ学習
ベイジアンネットワークを完全に指定し、結合確率分布を 完全に表現するためには、各ノード Xに対して、 X の 親を条件とする X の確率分布を 指定する必要があります。親を条件とする X の分布は、任意の形式を取ることができます。計算を簡素化するため、離散分布または ガウス分布 を使用するのが一般的です。分布に対する制約のみがわかっている場合もあります。その場合は、 最大エントロピーの原理を 使用して、 制約が与えられた場合に エントロピー が最大となる単一の分布を決定できます。(同様に、 動的ベイジアンネットワーク の特定のコンテキストでは、隠れ状態の時間的発展に対する条件付き分布は、通常、暗黙の確率過程の エントロピー率を 最大化するように指定されます。)
多くの場合、これらの条件付き分布には未知のパラメータが含まれており、 最尤法 などを用いてデータから推定する必要があります。尤度(または 事後確率 )を直接最大化することは、観測されていない変数が与えられた場合、しばしば複雑になります。この問題に対する古典的なアプローチは 期待値最大化アルゴリズム です。これは、観測されたデータに基づいて観測されていない変数の期待値を計算することと、以前に計算された期待値が正しいと仮定して完全な尤度(または事後確率)を最大化することを交互に繰り返します。軽度の正則性条件下では、このプロセスはパラメータの最大尤度(または最大事後確率)値に収束します。
パラメータに対するより完全なベイズ的アプローチは、パラメータを追加の観測されていない変数として扱い、観測データに基づいてすべてのノードにわたる完全な事後分布を計算し、その後パラメータを積分して取り除くことです。このアプローチはコストが高く、モデルが大規模になる可能性があるため、従来のパラメータ設定アプローチの方が扱いやすくなります。
構造学習
最も単純なケースでは、ベイジアンネットワークは専門家によって定義され、推論を実行するために使用されます。他のアプリケーションでは、ネットワークを定義する作業は人間にとって複雑すぎるため、ネットワーク構造と局所分布のパラメータをデータから学習する必要があります。
ベイジアンネットワーク(BN)のグラフ構造を自動学習することは、 機械学習における課題の一つです。基本的な考え方は、Rebaneと Pearl [7] によって開発された回復アルゴリズムに遡り 、3ノードDAGで許容される3つのパターンの区別に基づいています。
最初の2つは同じ依存関係( と は与えられた場合独立 )を表しており、したがって区別できません。しかし、 と は 限界独立であり、他のすべてのペアは従属しているため、コライダーは一意に識別できます。したがって、これら3つのトリプレットの スケルトン (矢印を取り除いたグラフ)は同一ですが、矢印の方向性は部分的に識別可能です。 と が 共通の親を持つ場合も同じ区別が適用されますが、まずそれらの親を条件とする必要があります。基になるグラフのスケルトンを体系的に決定し、次に観察された条件付き独立性によって方向性が決定されるすべての矢印の向きを決定するアルゴリズムが開発されています。 [2] [8] [9] [10]
X
{\displaystyle X}
Z
{\displaystyle Z}
Y
{\displaystyle Y}
X
{\displaystyle X}
Z
{\displaystyle Z}
X
{\displaystyle X}
Z
{\displaystyle Z}
構造学習の代替法として、最適化ベースの探索を用いるものがある。この手法では、 スコアリング関数 と探索戦略が必要となる。一般的なスコアリング関数としては、 BIC やBDeuのような、トレーニングデータを与えられた場合の構造の 事後確率 が挙げられる。スコアを最大化する構造を返す 網羅的な探索 に必要な時間は、変数の数に対して 超指数的 である。局所探索戦略では、構造のスコアを向上させることを目的とした段階的な変更が行われる。 マルコフ連鎖モンテカルロ法(MCMC)のようなグローバル探索アルゴリズムでは、 局所的最小値 に陥ることを回避できる。 相互情報 量を最大化する構造を見つけること は、通常、親候補セットを k 個のノードに制限する [11] [12] [13] か、ノードごとに 最適な kを見つける [14] ことによって行われ、ベンチマークデータセットで一貫して高スコアを達成する手法である。
正確なBN学習のための特に高速な方法は、問題を最適化問題として捉え、 整数計画法を用いて解くことです。整数計画法(IP)の解法には、 切断面 の形で非巡回制約が追加されます 。 [15] この方法は、最大100個の変数を持つ問題を扱うことができます。
数千の変数を持つ問題を扱うには、異なるアプローチが必要です。まず、ある順序付けをサンプリングし、その順序付けに関して最適なBN構造を見つけるという方法です。これは、可能な順序付けの探索空間を扱うことを意味しますが、これはネットワーク構造の空間よりも小さいため便利です。次に、複数の順序付けをサンプリングし、評価します。この方法は、変数の数が膨大な場合に、文献で利用可能な最良の方法であることが証明されています。 [16]
もう一つの方法は、 MLEが 閉形式を持つ分解可能なモデルのサブクラスに焦点を当てることです 。これにより、数百の変数に対して一貫した構造を発見することが可能になります。 [17]
木幅が制限されたベイジアンネットワークを学習することは、正確で扱いやすい推論を可能にするために不可欠です。なぜなら、最悪の場合の推論複雑度は木幅kに対して指数関数的になるからです(指数時間仮説の下で)。しかし、グラフのグローバルな特性として、これは学習プロセスの難易度を著しく高めます。この文脈では、 K木を 用いて効果的な学習を行うことが可能です。 [18]
統計入門
データ とパラメータが与えられると 、単純な ベイズ分析は 事前確率 ( 事前 ) と 事後確率 を計算する 尤度 から開始されます 。
x
{\displaystyle x\,\!}
θ
{\displaystyle \theta }
p
(
θ
)
{\displaystyle p(\theta )}
p
(
x
∣
θ
)
{\displaystyle p(x\mid \theta )}
p
(
θ
∣
x
)
∝
p
(
x
∣
θ
)
p
(
θ
)
{\displaystyle p(\theta \mid x)\propto p(x\mid \theta )p(\theta )}
多くの場合、事前確率は 尤度には記載されていない他のパラメータに依存します 。そのため、事前確率は 尤度に置き換えられ 、 新たに導入されたパラメータの事前 確率が必要になります。その結果、事後確率は
θ
{\displaystyle \theta }
φ
{\displaystyle \varphi }
p
(
θ
)
{\displaystyle p(\theta )}
p
(
θ
∣
φ
)
{\displaystyle p(\theta \mid \varphi )}
p
(
φ
)
{\displaystyle p(\varphi )}
φ
{\displaystyle \varphi }
p
(
θ
,
φ
∣
x
)
∝
p
(
x
∣
θ
)
p
(
θ
∣
φ
)
p
(
φ
)
.
{\displaystyle p(\theta ,\varphi \mid x)\propto p(x\mid \theta )p(\theta \mid \varphi )p(\varphi ).}
これは階層的ベイズモデル の最も単純な例です 。
このプロセスは繰り返される場合があります。例えば、パラメータが さらに別のパラメータに依存し 、それらのパラメータにも事前分布が必要となる場合があります。最終的に、このプロセスは、言及されていないパラメータに依存しない事前分布で終了します。
φ
{\displaystyle \varphi }
ψ
{\displaystyle \psi \,\!}
入門例
それぞれの 測定量が既知の 標準偏差の 正規分布 誤差を持つとすると 、
x
1
,
…
,
x
n
{\displaystyle x_{1},\dots ,x_{n}\,\!}
σ
{\displaystyle \sigma \,\!}
x
i
∼
N
(
θ
i
,
σ
2
)
{\displaystyle x_{i}\sim N(\theta _{i},\sigma ^{2})}
を推定したいとします。 最尤推定 法を用いて 推定するというアプローチがあります 。観測値は独立しているので、尤度は因数分解され、最尤推定値は単純に
θ
i
{\displaystyle \theta _{i}}
θ
i
{\displaystyle \theta _{i}}
θ
i
=
x
i
.
{\displaystyle \theta _{i}=x_{i}.}
しかし、もし数量が関連している場合、例えば個体 自体が基礎となる分布から抽出されている場合、この関係は独立性を破壊し、より複雑なモデルを示唆する。例えば、
θ
i
{\displaystyle \theta _{i}}
x
i
∼
N
(
θ
i
,
σ
2
)
,
{\displaystyle x_{i}\sim N(\theta _{i},\sigma ^{2}),}
θ
i
∼
N
(
φ
,
τ
2
)
,
{\displaystyle \theta _{i}\sim N(\varphi ,\tau ^{2}),}
不適切な事前分布 の場合、 となります 。 のとき 、これは 同定されたモデル (つまり、モデルのパラメータに一意の解が存在する)であり、個体の事後分布は 最大尤度推定値から共通平均に向かって 移動する、つまり 縮小する 傾向があります。この縮小は 、階層ベイズモデルにおける典型的な挙動です。
φ
∼
flat
{\displaystyle \varphi \sim {\text{flat}}}
τ
∼
flat
∈
(
0
,
∞
)
{\displaystyle \tau \sim {\text{flat}}\in (0,\infty )}
n
≥
3
{\displaystyle n\geq 3}
θ
i
{\displaystyle \theta _{i}}
事前分布の制限
階層モデルにおいて事前分布を選択する際には、特に例の 変数のように階層の上位レベルにある尺度変数については注意が必要です。Jeffreys 事前分布 のような通常の事前分布は、事後分布が正規化できず、 期待損失を最小化する推定値が許容されないため、多くの場合 機能しませ ん 。
τ
{\displaystyle \tau \,\!}
定義と概念
ベイジアンネットワークには、いくつかの同等の定義が提案されている。以下では、 G = ( V , E ) を 有向非巡回グラフ (DAG) とし、 X = ( X v )、 v ∈ Vを V でインデックス付けされた 確率変数 の集合とする 。
因数分解の定義
Xが G に関してベイジアンネットワークである とは、その結合 確率密度関数( 積測度 に関して )が、親変数を条件として、個々の密度関数の積として表されるときである。
p
(
x
)
=
∏
v
∈
V
p
(
x
v
|
x
pa
(
v
)
)
{\displaystyle p(x)=\prod _{v\in V}p\left(x_{v}\,{\big |}\,x_{\operatorname {pa} (v)}\right)}
ここで、pa( v )は v の親の集合 (つまり、 単一の辺を介して
vに直接指し示す頂点)である。
任意の確率変数の集合について、 結合分布の任意の要素の確率は、 連鎖律 ( X の 位相的順序 が与えられている)を用いた条件付き確率から次のように 計算できる:
P
(
X
1
=
x
1
,
…
,
X
n
=
x
n
)
=
∏
v
=
1
n
P
(
X
v
=
x
v
∣
X
v
+
1
=
x
v
+
1
,
…
,
X
n
=
x
n
)
{\displaystyle \operatorname {P} (X_{1}=x_{1},\ldots ,X_{n}=x_{n})=\prod _{v=1}^{n}\operatorname {P} \left(X_{v}=x_{v}\mid X_{v+1}=x_{v+1},\ldots ,X_{n}=x_{n}\right)}
上記の定義を使用すると、次のように記述できます。
P
(
X
1
=
x
1
,
…
,
X
n
=
x
n
)
=
∏
v
=
1
n
P
(
X
v
=
x
v
∣
X
j
=
x
j
for each
X
j
that is a parent of
X
v
)
{\displaystyle \operatorname {P} (X_{1}=x_{1},\ldots ,X_{n}=x_{n})=\prod _{v=1}^{n}\operatorname {P} (X_{v}=x_{v}\mid X_{j}=x_{j}{\text{ for each }}X_{j}\,{\text{ that is a parent of }}X_{v}\,)}
2 つの式の違いは、 親変数の値が与えられた場合、変数が子孫以外の変数から
条件付きで独立しているかどうかです。
局所マルコフ性
Xが G に関してベイジアンネットワークであるとは、 局所マルコフ性 を満たす場合である 。つまり、各変数は 親変数が与えられた場合に、その非子孫変数から 条件付きで独立している。
X
v
⊥
⊥
X
V
∖
de
(
v
)
∣
X
pa
(
v
)
for all
v
∈
V
{\displaystyle X_{v}\perp \!\!\!\perp X_{V\,\smallsetminus \,\operatorname {de} (v)}\mid X_{\operatorname {pa} (v)}\quad {\text{for all }}v\in V}
ここで、de( v )は子孫の集合であり、 V \de( v )は v の非子孫の集合である 。
これは最初の定義と似た言葉で表現できる。
P
(
X
v
=
x
v
∣
X
i
=
x
i
for each
X
i
that is not a descendant of
X
v
)
=
P
(
X
v
=
x
v
∣
X
j
=
x
j
for each
X
j
that is a parent of
X
v
)
{\displaystyle {\begin{aligned}&\operatorname {P} (X_{v}=x_{v}\mid X_{i}=x_{i}{\text{ for each }}X_{i}{\text{ that is not a descendant of }}X_{v}\,)\\[6pt]={}&P(X_{v}=x_{v}\mid X_{j}=x_{j}{\text{ for each }}X_{j}{\text{ that is a parent of }}X_{v}\,)\end{aligned}}}
グラフは非巡回 なので、親の集合は非子孫の集合のサブセットになります 。
限界独立構造
一般に、データからベイジアンネットワークを学習することは NP困難 であることが知られています。 [21] これは、変数の数が増えるにつれて、 DAGを列挙する際 の 組み合わせ爆発が 一部に起因します。しかし、周辺独立性構造に注目することで、基礎となるベイジアンネットワークについての洞察をデータから多項式時間で得ることができます。 [22] ベイジアンネットワークによってモデル化された分布の条件付き独立性ステートメントは(上記の因数分解とマルコフ特性に従って)DAGによってエンコードされますが、その周辺独立性ステートメント(条件セットが空である条件付き独立性ステートメント)は、等しい 交差 や 独立数 などの特別な特性を持つ 単純な無向グラフ によってエンコードされます。
ベイジアンネットワークの開発
ベイジアンネットワークの開発は、多くの場合、 Xが G に関して局所マルコフ性を満たすよう な DAG G を作成することから始まります。これは 因果 DAGである場合もあります。G内の親変数を与えられた各変数の条件付き確率分布を 評価します。多くの場合、特に変数が離散的な場合、 Xの結合分布 が これらの条件付き分布の積である場合、 Xは G に関してベイジアンネットワークとなります 。 [23]
マルコフブランケット
ノードのマルコフ ブランケット とは、その親、子、そしてその子のさらに別の親からなるノードの集合である。マルコフブランケットは、ノードをネットワークの他の部分から独立させる。ノードのマルコフブランケット内の変数の結合分布は、ノードの分布を計算するのに十分な知識である。X は、その マルコフブランケット が与えられた場合、ネットワーク内の各ノードが他のすべてのノードから条件付きで独立しているとき、 G に関してベイジアンネットワークである 。
d -分離
この定義は、2つのノードの「d」間隔を定義することでより一般化できます。ここで、dは方向性を表します。 [2] まずトレイルの「d」間隔を定義し、次にそれに基づいて2つのノードの「d」間隔を定義します。
P をノード uから v への トレイルとします 。 トレイルとは、2つのノード間のループのない無向パス(つまり、すべてのエッジ方向が無視されるパス)です。以下の条件のいずれかが成立する場合、 P は ノード集合 Zによって d 分離されているとされます 。
P には(完全に有向連鎖である必要はないが)有向連鎖、 つまり が含まれ 、中間ノード mは Z に含まれる 。
u
⋯
←
m
←
⋯
v
{\displaystyle u\cdots \leftarrow m\leftarrow \cdots v}
u
⋯
→
m
→
⋯
v
{\displaystyle u\cdots \rightarrow m\rightarrow \cdots v}
P にはフォークが含まれており 、中間ノード mは Z 内にある 、または
u
⋯
←
m
→
⋯
v
{\displaystyle u\cdots \leftarrow m\rightarrow \cdots v}
P には反転フォーク (またはコライダー) が含まれており、 中間のノード mは Z に存在せず、 m の子孫は Z に存在しません 。
u
⋯
→
m
←
⋯
v
{\displaystyle u\cdots \rightarrow m\leftarrow \cdots v}
ノード u と vは、 それら の間のすべてのトレイルが dで区切られている場合、 Z によって d で区切られています 。u と v がdで区切られていない場合、それらはdで接続されています。
Xが G に関してベイジアンネットワークである とは、任意の 2 つのノード u 、 v に対して次の条件が成立する場合である。
X
u
⊥
⊥
X
v
∣
X
Z
{\displaystyle X_{u}\perp \!\!\!\perp X_{v}\mid X_{Z}}
ここで、 Zは u と vを d 分離する 集合です 。( マルコフブランケット は、ノード v を他のすべてのノードから d 分離するノードの最小集合です 。)
因果ネットワーク
ベイジアンネットワークは因果 関係を表すためによく用いられますが 、必ずしもそうである必要はありません。u から v への有向辺は、 X v が X u に 因果的に依存する ことを要求しません。これは、グラフ上のベイジアンネットワークが以下 の 条件を満たすという事実によって証明されています。
a
→
b
→
c
and
a
←
b
←
c
{\displaystyle a\rightarrow b\rightarrow c\qquad {\text{and}}\qquad a\leftarrow b\leftarrow c}
同等です。つまり、まったく同じ条件付き独立性要件を課します。
因果ネットワークは、関係が因果関係にあるという要件を持つベイジアンネットワークです。因果ネットワークの追加のセマンティクスは、ノード X が能動的に特定の状態 x (do( X = x )と記述されるアクション)に遷移した場合、確率密度関数は、 Xの親から X へのリンクを切断し 、 X をその状態 x に設定することによって得られるネットワークの確率密度関数に変化することを規定します 。 [2] これらのセマンティクスを用いることで、介入前に得られたデータから、外部介入の影響を予測することができます。
推論の複雑さと近似アルゴリズム
1990年、スタンフォード大学で大規模なバイオインフォマティクスアプリケーションに取り組んでいたクーパーは、ベイジアンネットワークにおける正確な推論が NP困難で あることを証明した。 [24] この結果がきっかけとなり、確率的推論の扱いやすい近似を開発することを目的とした近似アルゴリズムの研究が進んだ。1993年、ポール・ダガムと マイケル・ルビーは、 ベイジアンネットワークにおける確率的推論の近似の複雑さに関する2つの驚くべき結果を証明した。 [25] まず、扱いやすい 決定論的アルゴリズムでは確率的推論を 絶対誤差 ɛ < 1/2以内に近似できないこと を証明した。次に、扱いやすい ランダム化アルゴリズムで は信頼確率が 1/2 を超える場合、
確率的推論を絶対誤差 ɛ < 1/2以内に近似できないことを証明した。
ほぼ同じ頃、 ロスは ベイジアンネットワークにおける正確な推論は実際には #P完全(したがって、 連言標準 形式(CNF)の充足割り当ての数を数えるのと同じくらい困難)であり、制限されたアーキテクチャを持つベイジアンネットワークであっても、 すべての ɛ > 0に対して係数2 n 1− ɛ 以内の近似推論はNP困難であることを証明した。 [26] [27]
実用的な観点から言えば、これらの複雑性の結果は、ベイジアンネットワークがAIや機械学習アプリケーションにとって豊富な表現である一方で、実際の大規模なアプリケーションで使用するには、ナイーブベイズネットワークなどの位相的な構造的制約、または条件付き確率に対する制限によって調整する必要があることを示唆しています。 DagumとLubyによって開発された制限付き分散アルゴリズム [28] は、誤差近似を保証したベイジアンネットワークでの確率的推論を効率的に近似する、初めての証明可能な高速近似アルゴリズムでした。 この強力なアルゴリズムでは、ベイジアンネットワークの条件付き確率に対する小さな制限が0と1から制限される必要がありました。 ここで、 はネットワーク内のノードの数の任意の多項式でした 。
1
/
p
(
n
)
{\displaystyle 1/p(n)}
p
(
n
)
{\displaystyle p(n)}
n
{\displaystyle n}
ソフトウェア
ベイジアン ネットワークの注目すべきソフトウェアには次のようなものがあります。
OpenBUGS – WinBUGS のオープンソース開発。
SPSS Modeler – ベイジアン ネットワークの実装を含む商用ソフトウェア。
Stan(ソフトウェア) – Stanは、ハミルトンモンテカルロ法の変種であるNo-U-Turnサンプラー(NUTS) [29] を使用してベイズ推論を行うオープンソースパッケージです。
WinBUGS – MCMCサンプラーの初期の計算実装の一つ。現在はメンテナンスされていません。
歴史
ベイジアンネットワークという用語は、 1985年に ジューディア・パール によって造語され、以下の点を強調した。 [30]
入力情報は主観的であることが多い
情報更新の基礎としてベイズの条件付けに依存すること
因果推論と証拠推論の区別 [31]
1980年代後半には、パールの 「知能システムにおける確率的推論」 [32] と ナポリタン の 「エキスパートシステムにおける確率的推論」 [33]が それらの特性をまとめ、研究分野として確立しました。
参照
注記
^ Ruggeri, Fabrizio; Kenett, Ron S.; Faltin, Frederick W. 編 (2007-12-14). 品質と信頼性に関する統計百科事典 (第1版). Wiley. p. 1. doi :10.1002/9780470061572.eqr089. ISBN 978-0-470-01861-3 。
^ abcde パール、ジューディア (2000年)。『因果律:モデル、推論、推論』 ケンブリッジ大学出版 局 。ISBN 978-0-521-77362-1 . OCLC 42291253。
^ 「バックドア基準」 (PDF) . 2014年9月18日 閲覧 。
^ 「d-分離と涙」 (PDF) . 2014年9月18日 閲覧 。
^ Pearl J (1994). 「行動の確率的計算」. Lopez de Mantaras R, Poole D (編). UAI'94 Proceedings of the Tenth international conference on Uncertainty in Artificial Intelligence . San Mateo CA: Morgan Kaufmann . pp. 454– 462. arXiv : 1302.6835 . Bibcode :2013arXiv1302.6835P. ISBN 1-55860-332-8 。
^ Shpitser I, Pearl J (2023). 「条件付き介入分布の識別」. Dechter R, Richardson TS (編). ベイジアンネットワークの周辺独立性構造に関する組み合わせ論的および代数的視点 . 第14巻. オレゴン州コーバリス: AUAI Press. pp. 437– 444. arXiv : 1206.6876 . doi :10.2140/astat.2023.14.233.
^ Rebane G, Pearl J (1987). 「統計データからの因果多木構造の復元」. 第3回AIにおける不確実性に関するワークショップ議事録 . シアトル, ワシントン州. pp. 222– 228. arXiv : 1304.2736 . {{cite book }}: CS1 maint: location missing publisher (link )
^ Spirtes P, Glymour C (1991). 「疎因果グラフの高速回復アルゴリズム」 (PDF) . Social Science Computer Review . 9 (1): 62– 72. CiteSeerX 10.1.1.650.2922 . doi :10.1177/089443939100900106. S2CID 38398322.
^ Spirtes P、Glymour CN、Scheines R (1993)。因果関係、予測、および検索 (第 1 版)。スプリンガー・フェルラーク。 ISBN 978-0-387-97979-3 。
^ Verma T, Pearl J (1991). 「因果モデルの等価性と合成」. Bonissone P, Henrion M, Kanal LN, Lemmer JF (編). UAI '90 Proceedings of the Sixth Annual Conference on Uncertainty in Artificial Intelligence . Elsevier. pp. 255– 270. ISBN 0-444-89264-8 。
^ Sahami, Mehran (1996-08-02). 「限定依存ベイズ分類器の学習」. 第2回国際知識発見・データマイニング会議議事録 . KDD'96. オレゴン州ポートランド: AAAI Press: 335– 338.
^ Friedman N, Geiger D, Goldszmidt M (1997年11月). 「ベイジアンネットワーク分類器」. 機械学習 . 29 ( 2–3 ): 131–163 . doi : 10.1023/A:1007465528199 .
^ Friedman N, Linial M, Nachman I, Pe'er D (2000年8月). 「ベイジアンネットワークを用いた発現データ解析」. Journal of Computational Biology . 7 ( 3–4 ): 601–20 . CiteSeerX 10.1.1.191.139 . doi :10.1089/106652700750050961. PMID 11108481.
^ Rubio, Arcadio; Gámez, José Antonio (2011-07-12). 「k依存ベイジアンネットワーク分類器の柔軟な学習」. 第13回遺伝的・進化的計算年次会議議事録. GECCO '11. ニューヨーク州ニューヨーク: Association for Computing Machinery. pp. 1219– 1226. doi :10.1145/2001576.2001741. ISBN 978-1-4503-0557-0 。
^ Cussens J (2011). 「切断面を用いたベイジアンネットワーク学習」 (PDF) . Proceedings of the 27th Conference Annual Conference on Uncertainty in Artificial Intelligence : 153– 160. arXiv : 1202.3713 . Bibcode :2012arXiv1202.3713C. 2022年3月27日時点のオリジナルよりアーカイブ。
^ Scanagatta M, de Campos CP, Corani G, Zaffalon M (2015). 「数千の変数を持つベイジアンネットワークの学習」. NIPS-15: ニューラル情報処理システムの進歩 . 第28巻. Curran Associates. pp. 1855– 1863.
^ Petitjean F, Webb GI, Nicholson AE (2013). 高次元データへの対数線形分析のスケーリング (PDF) . 国際データマイニング会議. ダラス, テキサス州, 米国: IEEE.
^ M. Scanagatta, G. Corani, CP de Campos, M. Zaffalon. 数千の変数を持つツリー幅制限ベイジアンネットワークの学習. NIPS-16: Advances in Neural Information Processing Systems 29, 2016.
^ Chickering, David M.; Heckerman, David; Meek, Christopher (2004). 「ベイジアンネットワークの大規模サンプル学習はNP困難である」 (PDF) . Journal of Machine Learning Research . 5 : 1287–1330 .
^ Deligeorgaki, Danai; Markham, Alex; Misra, Pratik; Solus, Liam (2023). 「ベイジアンネットワークの周辺独立性構造に関する組み合わせ的および代数的視点」. 代数統計 . 14 (2): 233– 286. arXiv : 2210.00822 . doi :10.2140/astat.2023.14.233.
^ ナポリ語 RE (2004). ベイジアンネットワークの学習. プレンティス・ホール. ISBN 978-0-13-012534-7 。
^ Cooper GF (1990). 「ベイズ的信念ネットワークを用いた確率的推論の計算複雑性」 (PDF) . 人工知能 . 42 ( 2–3 ): 393–405 . doi :10.1016/0004-3702(90)90060-d. S2CID 43363498.
^ Dagum P, Luby M (1993). 「ベイズ的信念ネットワークにおける確率的推論の近似はNP困難である」. 人工知能 . 60 (1): 141– 153. CiteSeerX 10.1.1.333.1586 . doi :10.1016/0004-3702(93)90036-b.
^ D. Roth, 近似推論の困難さについて, IJCAI (1993)
^ D. Roth, 近似推論の困難性について, 人工知能 (1996)
^ Dagum P, Luby M (1997). 「ベイズ推論のための最適近似アルゴリズム」. 人工知能 . 93 ( 1–2 ): 1–27 . CiteSeerX 10.1.1.36.7946 . doi :10.1016/s0004-3702(97)00013-1. 2017年7月6日時点のオリジナルよりアーカイブ。 2015年12月19日 閲覧 。
^ Hoffman, Matthew D.; Gelman, Andrew (2011). 「No-U-Turn Sampler: Hamiltonian Monte Carlo におけるパス長の適応的設定」 arXiv : 1111.4246 [stat.CO].
^ Pearl J (1985). ベイジアンネットワーク:証拠推論のための自己活性化記憶モデル (UCLA技術レポートCSD-850017) . 認知科学会第7回大会議事録、カリフォルニア大学アーバイン校、カリフォルニア州. pp. 329– 334. 2009年5月1日 閲覧 。
^ ベイズ・T 、プライス( 1763). 「偶然性の教義における問題解決に向けた試論」 王立 協会哲学論文集 53 : 370–418 . doi : 10.1098/rstl.1763.0053 .
^ Pearl J (1988-09-15). 知能システムにおける確率的推論. サンフランシスコ: Morgan Kaufmann . p. 1988. ISBN 978-1-55860-479-7 。
^ ナポリタンRE(1989年)『エキスパートシステムにおける確率的推論:理論とアルゴリズム』Wiley. ISBN 978-0-471-61840-9 。
参考文献
Ben Gal I (2007). 「ベイジアンネットワーク」 (PDF) . Ruggeri F, Kennett RS, Faltin FW (編). サポートページ . 品質と信頼性に関する統計百科事典 . John Wiley & Sons . doi : 10.1002/9780470061572.eqr089 . ISBN 978-0-470-01861-3 . 2016年11月23日時点の オリジナル (PDF)からアーカイブ 。 2007年8月27日 閲覧。
Bertsch McGrayne S (2011). 『死なない理論 』 ニューヘイブン: イェール大学出版局 .
Borgelt C, Kruse R (2002年3月). 『グラフィカルモデル:データ分析とマイニングのための手法』. 英国チチェスター : Wiley . ISBN 978-0-470-84337-6 。
Borsuk ME (2008). 「生態学的情報科学:ベイジアンネットワーク」 Jørgensen, Sven Erik 、 Fath, Brian (編). 『生態学百科事典』 エルゼビア. ISBN 978-0-444-52033-3 。
Castillo E, Gutiérrez JM, Hadi AS (1997). 「ベイジアンネットワークの学習」. エキスパートシステムと確率的ネットワークモデル . コンピュータサイエンスのモノグラフ. ニューヨーク: Springer-Verlag . pp. 481– 528. ISBN 978-0-387-94858-4 。
Comley JW, Dowe DL (2003年6月). 「一般ベイジアンネットワークと非対称言語」. 第2回ハワイ国際統計学および関連分野会議議事録 .
Comley JW, Dowe DL (2005). 「最小メッセージ長と非対称言語における一般化ベイジアンネット」. Grünwald PD, Myung IJ, Pitt MA (編). 最小記述長の進歩:理論と応用 . ニューラル情報処理シリーズ. マサチューセッツ州ケンブリッジ : Bradford Books ( MIT Press ) (2005年4月出版). pp. 265– 294. ISBN 978-0-262-07262-5 。 (本論文では、最小メッセージ長( MML ) を使用して、ベイズネットワークの内部ノードに 決定木を配置します。
Darwiche A (2009). ベイジアンネットワークによるモデリングと推論. ケンブリッジ大学出版局 . ISBN 978-0-521-88438-9 。
Dowe, David L. (2011-05-31). 「ハイブリッド・ベイジアンネットワーク・グラフィカルモデル、統計的一貫性、不変性、そして一意性」 (PDF) . 統計哲学. Elsevier. pp. 901–982. ISBN 978-0-08-093096-1 。
Fenton N, Neil ME (2007年11月). 「現代社会におけるリスク管理:ベイジアンネットワークの応用」 (PDF) . ロンドン数学会と産業数学知識移転ネットワークによる知識移転報告書 . ロンドン(イギリス) : ロンドン数学会 . オリジナル (PDF) から2008年5月14日アーカイブ. 2008年10月29日 閲覧 .
Fenton N, Neil ME (2004年7月23日). 「ベイジアンネットワークを用いたリスク分析における証拠の統合」 (PDF) . Safety Critical Systems Club Newsletter . 第13巻第4号. ニューカッスル・アポン・タイン , イギリス. pp. 8– 13. オリジナル (PDF) から2007年9月27日にアーカイブ。
Gelman A, Carlin JB, Stern HS, Rubin DB (2003). 「第2部 ベイズデータ分析の基礎:第5章 階層モデル」ベイズデータ分析. CRC Press . pp. 120–. ISBN 978-1-58488-388-3 。
ヘッカーマン、デイビッド(1995年3月1日)「ベイジアンネットワークによる学習チュートリアル」。マイケル・アーウィン・ジョーダン(編)『 グラフィカルモデルによる学習 適応型計算と機械学習』、 マサチューセッツ州ケンブリッジ : MIT出版 (1998年出版)。301 ~ 354頁 。ISBN 978-0-262-60032-3 . 2006年7月19日時点のオリジナルよりアーカイブ 。 2006年 9月15日 閲覧。 {{cite book }}: CS1 maint: bot: original URL status unknown (link ):また、 Heckerman, David (1997年3月) としても登場します。「データマイニングのためのベイジアンネットワーク」。 データマイニングと知識発見 。1 ( 1): 79– 119。doi : 10.1023 /A:1009730122752。S2CID 6294315 。
以前のバージョンは、 Microsoft Research 、1995 年 3 月 1 日に掲載されています。この論文は、ベイジアン ネットワークにおけるパラメーターと構造の学習の両方について説明しています。
Jensen FV, Nielsen TD (2007年6月6日). ベイジアンネットワークと決定グラフ. 情報科学・統計シリーズ(第2版). ニューヨーク : Springer-Verlag . ISBN 978-0-387-68281-5 。
Karimi K, Hamilton HJ (2000). 「時間的関係の発見:因果ベイジアンネットワーク vs. C4.5」 (PDF) . 第12回国際知能システム方法論シンポジウム .
Korb KB、Nicholson AE(2010年12月). ベイズ人工知能. CRC Computer Science & Data Analysis(第2版). Chapman & Hall ( CRC Press ). doi :10.1007/s10044-004-0214-5. ISBN 978-1-58488-387-6 . S2CID 22138783。
Lunn D, Spiegelhalter D, Thomas A, Best N (2009年11月). 「BUGSプロジェクト:進化、批判、そして将来の方向性」. Statistics in Medicine . 28 (25): 3049–67 . doi :10.1002/sim.3680. PMID 19630097. S2CID 7717482.
Neil M, Fenton N, Tailor M (2005年8月). Greenberg, Michael R. (編). 「ベイジアンネットワークを用いた期待損失および想定外の運用損失のモデル化」 (PDF) . リスク分析 . 25 (4): 963–72 . Bibcode :2005RiskA..25..963N. doi :10.1111/j.1539-6924.2005.00641.x. PMID 16268944. S2CID 3254505.
Pearl J (1986年9月). 「信念ネットワークにおける融合、伝播、そして構造化」. 人工知能 . 29 (3): 241– 288. doi :10.1016/0004-3702(86)90072-X.
Pearl J (1988). 知能システムにおける確率的推論:尤もらしい推論のネットワーク. 表現と推論シリーズ(第2刷). サンフランシスコ、カリフォルニア州 : Morgan Kaufmann . ISBN 978-0-934613-73-6 。
Pearl J , Russell S (2002年11月). 「ベイジアンネットワーク」. Arbib MA (編). 『脳理論とニューラルネットワークハンドブック』. マサチューセッツ州ケンブリッジ : Bradford Books ( MIT Press ). pp. 157– 160. ISBN 978-0-262-01197-6 。
ラッセル、スチュアート・J. 、 ノーヴィグ、ピーター (2003年)、人工知能:現代的アプローチ(第2版)、アッパーサドルリバー、ニュージャージー:プレンティスホール、 ISBN 0-13-790395-2 。
Zhang NL, Poole D (1994年5月). 「ベイジアンネットワーク計算への シンプルなアプローチ」 (PDF) . 第10回隔年カナダ人工知能会議 (AI-94) 議事録. 171–178 . この論文では、信念ネットワークの変数除去について説明します。
さらに読む
Conrady S, Jouffe L (2015-07-01). ベイジアンネットワークとBayesiaLab ― 研究者のための実践的入門. フランクリン, テネシー州: Bayesian USA. ISBN 978-0-9965333-0-0 。
Charniak E (1991年冬). 「ティアーズのないベイジアンネットワーク」 (PDF) . AI Magazine .
クルーゼ R、ボルゲルト C、クラウォン F、モーヴェス C、スタインブレッチャー M、ヘルド P (2013)。計算知能の方法論的紹介。ロンドン: Springer-Verlag。 ISBN 978-1-4471-5012-1 。
Borgelt C, Steinbrecher M, Kruse R (2009). 『グラフィカルモデル ― 学習、推論、データマイニングのための表現』(第2版). チチェスター: Wiley. ISBN 978-0-470-74956-2 。
外部リンク
ベイジアンネットワークとその現代的応用入門
ベイジアンネットと確率に関するオンラインチュートリアル
ベイジアンネットを作成し、モンテカルロ法で実行するWebアプリ
連続時間ベイジアンネットワーク
ベイジアンネットワーク:説明と類推
ベイジアンネットワークを学ぶためのライブチュートリアル
分類問題におけるサンプルの異質性を処理するための階層的ベイズ モデルは、複製サンプルの測定に関連する不確実性を考慮した分類モデルを提供します。
サンプルの不確実性を処理するための階層的 Naive Bayes モデル ( Wayback Machine に 2007-09-28 でアーカイブ) では、複製された測定による連続変数と離散変数を使用して分類と学習を実行する方法を示します。