-------
ベイズ推定の例題とその解 02.10.2013
頻度主義的解法との比較
前回から始めたベイズ推計ですが、そこで提示した例題を、（ａ）ベイズの定理を用いた解法、（ｂ）古典的・頻度主義的な解法の二種類の解を示してみたいと思います。
このような作業をすると、何がベイズ推計の本質なのか、それがやっと見えてくると思うからです。
この例題は、次の書籍から引用したものです。「Ｅｘｃｅｌでスッキリわかる ベイズ統計入門」。涌井良幸、涌井貞美著、日本実業出版社、２０１０年１０月２０日、２２００円。以下、「本書」と記述することにします。
それでは早速取り掛かります。しかし、「誰にも、きちんと読んでもらえない」、という可能性が大きいような気もしますが。。。。。
前回の記事で、問題は５題あった。
２−（１）． 喫煙者の推定の問題
問題：
「男性１０人、女性７人が一室でパーティーを開いた。男子の喫煙者は５人、女性は３人である。部屋に入ったら煙草の吸殻が１本、灰皿の上にあった。このとき、吸った人が女性である確率を求めなさい（煙草の吸い回しはしていないと仮定する）」。
解法：
（ａ）ベイズの定理による解法
吸った人が女性である確率を求めるということは、Ｐ（女性｜喫煙者）を求めることである。
ベイズの定理によって、
Ｐ（喫煙者｜女性）・Ｐ（女性）
Ｐ（女性｜喫煙者）＝---------------------------
Ｐ(喫煙者）
を求めることになる。
Ｐ（喫煙者｜女性）＝３／７
Ｐ（女性）＝７／（７＋１０）＝７／１７ Ｐ（喫煙者）＝（３＋５）／(７＋１０）＝８／１７
上式に代入して、
Ｐ（女性｜喫煙者）＝（３／７ × ７／１７）／（８／１７）＝３／８ （答）
（ｂ）集合論による解法
次の図のような状況であった。吸殻を残したのは喫煙者の誰かであり、喫煙者が女性である確率は、図から明らかに３／８である。
図 状況を示すベン図
解説：この例題から分かるように、この問題をベイズの定理を使って解く理由は特にない。普通にベン図を使った集合論で解けば良い。
それでもベイズの定理を使うことを推奨できるとすれば、定理の各項にかなり機械的に数値を当てはめれば良いことだろうか。「どんな問題でも、ひとつの方法論で解けるならば簡単に思える」、という人であれば、ベイズの定理を使うのも一案かもしれない。
２−（２）． 雨の日のお誘いの問題
問題：
「酒好きのＡさんはＢ氏をよくお酒に誘う。統計をとると、雨の降っていない日に誘うと、Ｂ氏は５回中４回誘いに応じ、雨の降っている日に誘うと、５回中３回誘いに応じることが分かった。Ｂ氏がＡさんの誘いに応じたとき、雨が降っていない確率を求めよ。雨が降った日と降らない日の割合は７：１とする」。
解法：
（ａ）ベイズの定理による解
ベイズの展開公式を使う。
Ｐ（雨ではない｜誘いに乗る）
Ｐ（誘いに乗る｜雨ではない）・Ｐ（雨ではない）
＝------------------------------------------------------------------
Ｐ（誘いに乗る｜雨ではない）・Ｐ（雨ではない）＋ Ｐ（誘いに乗る｜雨）・Ｐ（雨）
Ｐ（誘いに乗る｜雨ではない）＝４／５
Ｐ（誘いに乗る｜雨）＝３／５
Ｐ（雨ではない）＝７／８
Ｐ（雨）＝１／８
を代入する。
４／５ × ７／８
Ｐ（雨ではない｜誘いに乗る）＝---------------------------------------
４／５ × ７／８ ＋ ３／５ × １／８
＝０．９０３
（ｂ）ベン図による解
状況をベン図に書けば、次のようになる。
図 お天気とお誘いのベン図
雨ではないときにお誘いに乗った数＝４×７＝２８
お誘いに乗った総数＝３＋４×７＝３１
お誘いに乗ったときに雨ではない確率は＝２８／３１＝０．９０３
解説：
この例題も（１）の例題と同様で、特に、ベイズの定理を使うことは不必要である。むしろ、使わない方が簡単である。
しかし、この例が何を示唆するか、と言えば、ベイズの定理の左辺をどう解釈するか、というある種の技法が、この問題には潜んでいるからである。
ここで、もう一度、ベイズの定理の図を示したい。
図 ベイズの定理の表現図
左辺Ｐ（A｜X）であるが、ベイズの定理は、集合論から導かれているので、本来の解釈は、単に、ＸであってかつＡでもある集合の大きさをXである集合の大きさで割った数値で表現される確率である。P(A)も、P(X)も、同様で、この図に示されているように、四角（＝全体）を分母とする確率である。
ところが、ベイズ推計では、P(A｜X）ということに、方向性を持たせている。右辺にあるP(X|A)も同様である。
もとに戻って考えなおしてみよう。
「風が吹けば桶屋が儲かる」、の例を再度掲載すれば、
「桶屋が儲かる」ときに「風吹く」確率
「風が吹く」ときに「桶屋が儲かる」確率 × 「風が吹く」確率
＝----------------------------------------------------------
「桶屋が儲かる」確率
である。
ベイズの定理は、もともと集合（ある条件を満たした日数）がその対象であるから、「風が吹いた日」のうち「桶屋が儲かった日」は何日あったか、が正しい表現である。
それが『「風が吹く」ときに「桶屋が儲かる」確率』という表現になっている。
そして、もともとの問題は、「風が吹けば桶屋が儲かる」、という形であった。これは「原因」と「結果」の表現になっているところに特徴がある。
すなわち、ベイズ推計の特徴は何か、と言えば、本来、集合論の問題であったものが、「原因」と「結果」の問題に拡張されていると理解するのが良さそうである。
２−（３）. 取り出した白玉はどちらの壺からの問題
問題：
「外からは区別の付かない壺ａ，ｂがある。壺ａには、白玉が２個、赤玉が３個入っている。壺ｂには、白玉が４個、赤玉が８個入っている。これら２つの壺ａ，ｂの一つを選択し、その壺から一つの玉を取り出したら白玉だった。この玉が壺ａから取り出された確率、壺ｂから取り出された確率を求めよ」。
解法：
（ａ）ベイズの定理による解
まずは、壺ａから白玉が取り出される確率は、５個の玉のうち白玉は２個なので、２／５である。
これは、Ｗを取り出した玉が白玉である確率、Ｈａを壺ａが選択された場合とすれば
「壺ａから取り出された玉が白である確率」＝Ｐ（Ｗ｜Ｈａ）＝２／５
同様に、
「壺ｂから取り出された玉が白である確率」＝Ｐ（Ｗ｜Ｈｂ）＝４／（４＋８）＝４／１２
さて、ここで、壺ａが選択されるか、あるいは、壺ｂが選択されるかについては、なんら情報はないので、「理由不十分の原則」によって、
Ｐ（Ｈａ）＝Ｐ（Ｈｂ）＝１／２
と設定する。
実は、このようなケースについて「理由不十分の原則」の適用を行うとする「本書」の記述は明らかに間違いだと思う。後述する。
Ｐ（Ｈａ｜Ｗ）＝Ｐ（壺ａ｜白玉）＝（２／５ × １／２）／（２／５ × １／２＋４／１２ × １／２）
結局、１／２はすべての項に乗じているので、すべて
＝（２／５）／（２／５ ＋ ４／１２）
＝６／１１
壺ｂから取り出される確率は、
Ｐ（Ｈｂ｜Ｗ）＝Ｐ（壺ｂ｜白玉）＝５／１１
（ｂ）集合論での解法
図 白玉が壺ａから取り出されるベン図
この図より
壺ａが選択されたとき白玉が取り出された確率＝２／５
壺ｂが選択されたとき白玉が取り出される確率＝４／１２
よって、白玉が壺ａから取り出される確率は
＝（２／５）／（２／５ ＋ ４／１２）
＝６／１１
壺ｂから取り出される確率は、
Ｐ（Ｈｂ｜Ｗ）＝Ｐ（壺ｂ｜白玉）＝５／１１
解説：
この例で、「理由不十分の原則」を導入するのは間違だ、と思う。なぜならば、「理由不十分の原則」で導入する１／２という数値は、学問的にきっちりとした頻度主義的な解法では当然だとされる数値だからである。なぜならば、「壺は外からは区別が付かない」ので、頻度主義では１／２になるのが当たり前だからである。すなわち、１／２は「理由充分な数値」なのである。
この例でも、ベイズの原則を用いる妥当性はなくて、むしろ、集合論を用いて解を得ることが簡単かつ明快である。すなわち、ベイズ推定の例としては、不適当である。
２−（４）． 壺から取り出した玉で真贋判定の問題
問題：
「Ａ社の作る本物の壺には、水晶玉とガラス玉が４：１の割合で入っている。Ｂ社の作るにせものの壺には、水晶玉とガラス玉が２：３の割合で入っている。２社の壺は外見から判別できない。また、壺には玉がたくさん入っているものとする。いま、Ａ社製かＢ社製か不明の壺があり、続けて３回玉を取り出したら、順に水晶玉、水晶玉、がらす玉であった。この壺がＡ社製である確率を求めよ」。
解法：
（ａ）ベイズの定理による解
ここでは、ベイズ更新という新しい概念が入ってくる。やっとベイズ推定の例になったという感触である。
Ｐ（Ｈａ｜Ｓ）を求める解とする。これは、水晶球Ｓが取り出されたときに、それが壺ａからである確率である
Ｐ（Ｈａ｜Ｓ）はベイズの定理の展開式によって、
P(S|Ha)・P(Ha)
P(Ha|S)=--------------------------------------
P(S|Ha)・P(Ha) + P(S|Hb)・P(Hb)
ａ→ｂ、Ｓ→Ｇ（ガラス玉）と置き換えれば、同様の式が全部で４式できる。
さて、通常の頻度を考えることで
P(S|Ha)=０．８ なぜならＡ社の壺にはＳ水晶玉が４個、ガラス玉が１個入っているから
P(G|Ha)=１−０．８＝０．２
ｂの壺についても同様に、
P(S|Hb)=０．４
P(G|Hb)=０．６
１）一個目の玉の取り出し
P(Ha|S)を上式で求めようとすると、P(Ha)＝壺ａが選択された場合の値、P(Hb)＝壺ｂが選択された場合の値、が未知である。
そこで、問題の壺が「区別が付かない」と記述されているので、同等にすることが、頻度主義でも当然であるので、いずれも１／２＝０．５とする。
すると、
P(Ha|S)=（０．８×０．５）／（０．８×０．５＋０．４×０．５）＝２／３
P（Hb|S)=１−２／３＝１／３
２）二個目の玉の取り出し
再び、P(Ha)壺ａが選択された場合、P(Hb)壺ｂが選択された場合、が未知である。しかし、一個目の確率が２／３、１／３と求められているので、これを使う方が良さそうである。
「本書」では、「１回目にＡ社の壺が選択されたとするなら、２回目は１回目よりもＡ社の壺が選択される確率は高いはずです。これを学習効果があるからです」、と記述されているが、これは嘘である。なぜならば、題意によれば、２回目も、１回目と同じ製造会社が不明の壺から玉を取り出すのだから。
正確には、１回目の結果により、製造会社がＡ社と推定することが「より正当だ」という結果がでたのだから、２回目には、その結果を使うのが順当な対応だということになる」。これを「ベイズ更新」と呼ぶ。
そこで、２回目の結果は、
０．８×２／３
P(Ha|S)=-------------------------------
０．８×２／３＋０．４×１／３
=４／５
P(Hb|S)=１／５
なお、「本書」では、（注）があって、「尤度（もっともらしさ：A社の壺から水晶が取り出される確率など）は、確率モデルの枠組みを規定するものなので、ベイズ更新でも変化しません、と書かれているが、これも間違い。
問題に記述があるように、「壺には玉がたくさん入っている」と仮定されているので、A社の壺から水晶が取り出される確率は、２回目以降も変わらない、という記述が正しい。
３）三回目の玉の取り出し
ニ回目の結果である
P(Ha)=４／５
P(Hb)=１／５
を使って、３回目を出す。
０．８×４／５
P(Ha|S)＝-------------------------------
０．８×４／５＋０．４×１／５
＝４／７
P(Hb|S)＝１−４／７＝３／７
４）結論
水晶玉、水晶玉、ガラス玉というように玉が取り出された壺がＡ社製である確率は、４／７である。
したがって、Ｂ社製である確率は３／７である。
（ｂ）場合の数による解法
これまで、「集合論による解法」であったが、これは「場合の数による解法」と述べた方が良いと思われる。
ちょっとだけ面倒なことをやらないと答えはでない。
まず、製造社不明の壺から水晶玉、水晶玉、ガラス玉という「順」で玉が取り出されたということであるが、その順番は重要なのだろうか。問題にもあるように、「壺にはたくさんの玉が入っている」ので、一つ取り出したとしても、確率は変わらない。すなわち、ここで重要なことは「順」ではなくて、３個の玉を取り出したところ、水晶玉が２個、ガラス玉が１個取り出されたという個数が重要なのである。
さて、３個の玉を取り出したとき、
場合１： ３個とも水晶玉
場合２： ２個が水晶玉、１個がガラス玉
場合３： １個が水晶玉、２個がガラス玉
場合４： ０個が水晶玉、３個がガラス玉
という４つの場合に別れる。そして、その場合の数は、２の３乗＝８通りである。
水晶＝水、ガラス＝ガと記述すれば、
場合１の数＝１ 水、水、水
場合２の数＝３ 水、水、ガの順番を変えたもの
場合３の数＝３ 水、ガ、ガの順番を変えたもの
場合４の数＝１ ガ、ガ、ガ
もしも壺がＡ社製であれば、
場合１の確率＝０．８×０．８×０．８
場合２の確率＝０．８×０．８×０．２×３
場合３の確率＝０．８×０．２×０．２×３
場合４の確率＝０．２×０．２×０．２
問題で示されているのは場合２なので、０．８×０．８×０．２×３
＝４／５ × ４／５ × １／５ × ３
＝４８／１２５
もしも壺がＢ社製であれば、
場合１の確率＝０．４×０．４×０．４
場合２の確率＝０．４×０．４×０．６×３
場合３の確率＝０．４×０．６×０．６×３
場合４の確率＝０．６×０．６×０．６
問題で示されている場合２では、
０．４×０．４×０．６×３
＝２／５×２／５×３／５×３
＝３６／１２５
したがって、
４８／１２５
壺がＡ社製である確率＝---------------------
４８／１２５＋３６／１２５
＝４８／（４８＋３６）
＝４／７
Ｂ社製である確率＝１−４／７＝３／７
解説：
このように、この問題においても、ベイズ推定の解と場合の数を数える頻度主義的解法による解は一致する。
「理由不十分の原則」や「ベイズ更新」といった仮定を用いたのに、結果が厳密解と一致するのは、奇妙な気がしないでもない。
疑問：「この一致は偶然なのか。それとも必然なのか」。その答えは、「このケースの場合には、必然である」。
まず、１個目の玉の取り出しのとき、「本書」では、「最初の玉を取り出すとき、Ａ社の壺とＢ社の壺の選択確率は不明です。だから、「理由不十分の原則」を用います、と記述されているが、実は、「壺の選択確率」がまず問題で、実は、選択された訳ではない。不明の壺が、最初からそこに１個あったに過ぎない。
１個の玉を取り出したとき、それが水晶玉であったとすると、
Ａ社製の壺であれば、水晶玉である確率は０．８、Ｂ社製の壺であれば０．４である。
だから、不明の壺がＡ社製である確率は、
０．８
＝---------------＝２／３
０．８＋０．４
になるに決まっている。集合論からの結論を得た場合と同じ式を用いているに過ぎない。
二個目の玉の場合は、場合の数を用いた解で説明した三個目のケースを簡単にしたものになって、その答えは、ベイズの定理による解法の二個目の場合と全く同じ式になる。
すなわち、この例題は、「ベイズ更新」といった新しい概念を用いたとは言っているものの、実は、全く頻度主義的な式を変形した形で用いたに過ぎない。
だから、厳密解と同じ解になるのは、当然なのである。
しかし、この方法をベイズ推計と呼ぶのなら、これを使うことに意味があるケースも無い訳ではない、とも言える。なんらかの推定値を求めなければならない場合があったとして、それが、頻度主義的な厳密解を求めることはできないという場合があるからである。
今回は、その例を示す余裕は無いので、次回以降に先送りとしたい。
２−（５）． 病気の診断：陽性と疑陽性の問題
問題：「ある病気を発見する検査Ｔに関して、次のことが知られている。
●病気にかかっている人を検査Ｔを適用すると、９８％の確率で病気であると正しく診断される。
●病気にかかっていない人に検査Ｔを適用すると、５％の確率であやまって病気にかかっていると診断される。
●人全体では、病気にかかっている人と病気にかかっていない人の割合は、それぞれ３％と９７％である。
さて、母集団から無作為に抽出された一人に検査Ｔを適用して、病気にかかっていると診断されたとき（陽性）、この人が本当に病気にかかっている確率を求めよ」。
解法：
（ａ）ベイズの定理による解
Haを病気である人の割合（すべての人が分母）
Hbを病気出ない人＝健康な人の割合（すべての人が分母）
Dを陽性の人（すべての人が分母）
と定義する。
P(病気｜陽性）が求める値で＝P(Ha|D)が求める値だということになる。
ベイズの定理の展開公式（二項の場合）を使う。
P(D|Ha)・P(Ha)
P(Ha|D)=----------------------------------
P(D|Ha)・P(Ha)+P(D|Hb)・P(Hb)
これに次の数値を代入する。
P(D|Ha)は病気である人が陽性と出る確率＝０．９８
P(Ha)はすべての人のうち病気である確率＝０．０３
P(D|Hb)は健康な人が陽性と出る確率＝０．０５
P(Hb)は健康な人の割合＝０．９７
０．９８×０．０３
P(D|Ha)＝----------------------------------
0.98×0.03 ＋ 0.05×0.97
＝０．３７７
陽性なのに本当に病気である確率は３７．７％であるということになる。陽性だからといって、心配しすぎると損をする。
（ｂ）場合の数による解法
このように割合（％）で示された場合、その実態を把握するのは、相当の技量を要する。
そこで、割合で示された場合には、「場合の数にして考える」ことをお薦めする。
この場合、すべてを整数にするために、１００００人を対象に考えるが、別に１０００人でも構わない。
さて、全体として１００００万人が居たとする。色々な条件の人がそれぞれ何人なのだろうか。色々な条件の人といっても、病気である人は３％、病気でない人は９７％とあるから、
病気である人＝１００００×０．０３＝３００人
病気でない人＝９７００人
病気でない人を検査すると５％の人が陽性と判定されるので、
陽性と誤診される人＝９７００×０．０５＝４８５人
病気である人で陽性になる人は９８％＝３００×０．９８＝２９４人
陽性と判定される人は、誤診と本当の病人を合わせて＝４８５＋２９４＝７７９人
求めるのは、陽性の人のうちで、病気である人なので＝２９４／７７９＝０．３７７
解説：
「本書」によれば、「行動経済学で言われることですが、人は事前確率にうといようです。「病気の人は９８％の確率で病気と判断される」という尤度に引きずられて、そのバックグラウンドを見落としてしまう傾向があるのです」と記述されている。
そんなことはない。単に、％などの割合で表現された数字を解析して、実態に変換し理解するのが、意外と厄介だということに過ぎない。
しかし、場合の数にして考えれば、人間の理解力を増強するような直感が利くので、実態を見抜くことが簡単にできる。
今回の最終結論
今回取り上げた確率に関する５つの問題は、実は、ベイズの定理などを使わないで解くことができる問題ばかりであったことがご理解いただけたのではないか、と思う。
それならなぜ「本書」がベイズの定理を使っているのか、そこには、ニつの理由があると思う。
一つ目の理由：式を丸暗記すれば、式に数値を代入するだけで、かなりバラエティーに富んだ問題でも、機械的に解が求まる。
二つ目の理由：この５つの問題については、ベイズの定理を使っても、頻度主義的な厳密解と同じ答えになること、（これは集合論を基本にしている問題である限り、あたり前なのであるが）、これを示すことによって、必ずしも理論的に正しいことが証明できていない部分を含むベイズ推計が正解に近い答えを出すことができることをアッピールしたい。
この二つ目の理由が余りにも目立つと、頻度主義的な統計学者と反りが合わないという結末になるのだろう。
ということで、ここまで紹介してきた「本書」は、入門書としてやさしさを狙っているのは分かるが、初心者を惑わすような決定的な記述ミスが何箇所かあって、お薦めはできない。
ベイズ推計の第二回目の説明としては、こんなところで止めておきたい。
