ベイズ推定の適用例を公式を使わずに解いてみて、直感的に理解する
ベイズ推定は公式を使わなくても解ける
顧客の指定する商品をミスなく出荷することは倉庫業務における基本ですが、出荷単位が細かくなるほどミスは起こりやすくなります。
そして国によってはミス率が高く、これが結構大きな問題になったりします。
A社は物流業者B社に倉庫業務を委託していますが、誤出荷があまりに多いため、配送トラックに載せる前の出荷前検品をB社に依頼しました。
しかし誤出荷はなかなか減りません。
A社はB社が本当に出荷前検品をしているのか疑念を持っています。
このような状況のもとで、B社の嘘をベイズ推定で暴いてみましょう。
しかし、ただベイズの定理に当てはめるだけでは実感が湧きにくいため、公式はあえて使わずに。
ベイズ推定の考え方は意外と簡単なのだということが分かると思います。
問題設定する
B社が出荷前検品を行う前の誤出荷率は３０％でした。
この誤出荷率は出荷前検品を行うことで５％まで下がることが分かっています。
またA社はB社のことを信頼しており、B社が出荷前検品をやると約束してくれたので、９０％の確率で本当にやってくれているものだと考えています。
このような状況のもと、出荷前検品を開始してから２回連続で誤出荷が発生しました。
A社は、B社が嘘をついている確率を何％だと考えればよいでしょうか？
公式を使わずにベイズ推定してみる
最初から嘘をついていると考えられる初期確率
A社は今までの付き合いから、B社は９０％の確率で出荷前検品をしているだろうと信用していたので、嘘をつく確率は１０％だと考えていました。
１回目の誤出荷を見て嘘をついていると考えられる確率
誤出荷を起こしたということは、
- 出荷前検品をしたのに誤出荷だった
- 出荷前検品をさぼって誤出荷になった
のうちのどちらかです。
ですので、下記のように２つの比率を求めればB社が嘘をついた（出荷前検品をさぼって誤出荷になった）確率が分かります。
２回目の誤出荷を見て嘘をついていると考えられる確率
以上のように１回目の誤出荷を見たA社は、４０％の確率でB社が嘘をついていると思うようになりました。
そうしたら、２回目も誤出荷でした。
先ほどと同じように、誤出荷を起こしたということは、
- 出荷前検品をしたのに誤出荷だった
- 出荷前検品をさぼって誤出荷になった
のうちのどちらかです。
ですので、これも先ほどと同様にB社が嘘をついている確率は次のように求められます。
さすがに２回連続して誤出荷になれば、B社が嘘をついている確率は８０％にまで上がりました。
嘘をついていると考えられる確率の推移
以上をまとめると、B社が嘘をついていると考えられる確率は、誤出荷が起こる度に次のように変化したことが分かります。
最初の初期確率１０％はA社の主観ですが、誤出荷が２回起こったことでベイズ推定によりその確率が４０％、８０％と変化しました。
このように新しい出来事が起こる度（誤出荷が起こる度）、確率が更新されるためベイズ更新とも呼ばれます。
初期確率が１％の場合のベイズ更新
先ほどは初期確率（B社が嘘をついていると考えられる確率）が１０％でしたが、これが１％の場合はどのように確率が更新されるでしょうか？
同じように計算すると次のようになります。
仏の顔も３度まで、B社は嘘をつくなんて９９％あり得ないと思っていたA社でも、さすがに３度も誤出荷が続けば嘘の確率の方が高いとベイズ推定では判断するということです。
ベイズ推定として一般化する
以上がベイズ推定の方法ですが、これを一般化してみましょう。
嘘をついている（検品していない）という事象をAとし、嘘をついていない（検品している）という事象をBとします。
すると、嘘をついていると考えている初期確率が１０％の場合、初期確率P(A)、P(B)は
P(A) = １０％
P(B) = ９０％
と表すことができます。
このP(A)＝１０％がベイズ推定における事前確率です。
１０％というのは、B社が嘘をついているとA社が考えている確率ですね。
あくまでA社の主観でそう思っているだけですので、主観確率です。
次にA社は誤出荷という事実を観察しました。
誤出荷が起こるのは
- 出荷前検品をしたのに誤出荷だった
- 出荷前検品をさぼって誤出荷になった
のうちのどちらかです。
ここで誤出荷という事象をXとすると、誤出荷が起こる確率はP(X)です。
そして問題設定の中にある
「この誤出荷率は出荷前検品を行うことで５％まで下がることが分かっています」
というのは
「検品あり（B）という条件の下で誤出荷（X）が５％の確率で起こる」
ということなので、条件付き確率でP(X|B)で表すことができ、
P(X|B)＝５％
と書けます。
「条件付き確率」の使い方を具体例を使ってわかりやすく解説【初心者向け】
従って、
検品をしたのに誤出荷になった確率は
＝P(X|B)P(B)
＝５％✕９０％
＝４.５％
になります。
これが先の１）の確率です。
次に２）の確率も同じように考えることができます。
問題設定の中にある
「B社が出荷前検品を行う前の誤出荷率は３０％でした」
というのは
「検品なし（A）という条件の下で誤出荷（X）が３０％の確率で起こる」
ということなので、条件付き確率でP(X|A)で表すことができ、
P(X|A)＝３０％
と書けます。
従って、
検品をせずに誤出荷になった確率は
＝P(X|A)P(A)
＝３０％✕１０％
＝３％
になります。
次にP(X|A)P(A)とP(X|B)P(B)を全体とした時のP(X|A)P(A)の割合が、検品なしで誤出荷になった確率なので、
検品なしで誤出荷になった確率
＝B社が嘘をついた確率
＝P(X|A)P(A)／｛P(X|A)P(A)＋P(X|B)P(B)｝
＝３％／（４.５％＋３％）
＝４０％
になります。
図で描くと次のようになります。
これが事後確率です。
A社の主観に過ぎなかった事前確率P(A)がXという事象（誤出荷）が起こったことで、
P(X|A)P(A)／｛P(X|A)P(A)＋P(X|B)P(B)｝
に更新されたということです。
これが正にベイズ推定の公式です。
分母にあるP(X|A)P(A)＋P(X|B)P(B)は周辺確率と呼ばれ、他の事象（AやB）によらずXが起こる確率のことです。
この場合は、B社が嘘をついているか否かにかかわらず、誤出荷が起こる確率です。