ベイジアン更新-コイン投げの例


7

ベイジアン更新について質問があります。一般に、ベイジアン更新は、以前の信念分布から事後を取得するプロセスを指します。

別の方法として、最初のステップの後を、さらなる計算のための事前の入力として使用することで、用語を理解できます。

以下は簡単な計算例です。方法aは標準計算です。方法bでは、次の事後を計算する前に、事後出力を入力として使用します。

メソッドaを使用して、P(F | HH)= 0.2を取得します。メソッドbを使用すると、P(F | HH)= 0.05が得られます。私の質問は、方法bが有効なアプローチであるかどうかです。


問題:あなたはコインを2回投げ、2つのヘッドを獲得します。コインが公正である確率はどのくらいですか、すなわちPr(Fair coin|HH)

では、最初のトスについて: Pr(Fair coin|H)=Pr(Head|Fair)P(Fair)Pr(Head|Fair)P(Fair)+Pr(Head|Biased)P(Biased)=Pr(H|F)P(F)P(H)(1)

事前信念P(Fair)= 0.5を開始すると仮定して、最初のトスのP(F | H)を見つけたい

以下は、中間ステップの計算です。

P(H|F)=(nx)θx(1θ)nx=(11)0.51(0.5)0=0.5

P(H)=P(H|F)P(F)+P(H|Biased)P(Biased)=(0.50.5)+(10.5)=0.75

(注:P(H |バイアス)= 1なので、コインの両側にヘッドがある極端な例を想定すると、バイアスされたコインでヘッドが得られる確率= 1(計算が簡単になります))

したがって、(1)に接続すると、次のようになります。

Pr(F|H)=Pr(H|F)P(F)P(H)=0.50.50.75=0.33


ここで、もう一度コインを投げて、別のHを得ます。計算するには Pr(F|HH) 、 我々

a)P(Fair)= 0.5を使い続ける

Pr(F|HH)=Pr(HH|F)P(F)Pr(HH|F)P(F)+Pr(HH|Biased)P(Biased)=Pr(HH|F)P(F)P(HH)(2)

P(HH|F)=(nx)θx(1θ)nx=(22)0.52(0.5)0=0.25

P(HH)=P(HH|F)P(F)+P(HH|Biased)P(Biased)=(0.250.5)+(10.5)=0.625

したがって、(2)に接続すると、 Pr(F|HH)=Pr(HH|F)P(F)P(HH)=0.250.50.625=0.2


または、計算するとどうなるか Pr(F|HH) を使用して

b)最初のステップでPr(F | H)から得た更新された信念P(Fair)= 0.33

この場合、

P(HH|F)=(nx)θx(1θ)nx=(22)0.332(10.33)0=0.1089

P(HH)=P(HH|F)P(F)+P(HH|Biased)P(Biased)=(0.10890.33)+(10.67)=0.705937

したがって、(2)に接続すると、 Pr(F|HH)=Pr(HH|F)P(F)P(HH)=0.10890.330.705937=0.05091


メソッドaを使用して、P(F | HH)= 0.2を取得します。メソッドbを使用すると、P(F | HH)= 0.05が得られます。私の質問は、方法bが有効なアプローチであるかどうかです。


バイアスされたコインが与えられた場合に表が出る確率は1であると、どのように考えましたか?とにかく、コインを何回投げても、公平である確率はゼロです。
Neil G、

コインがバイアスされていて、ヘッドが表示されている場合は、コインの両面にヘッドがあることを意味します。私たちは常にヘッドを見るので、バイアスされたコイン= 1でヘッドを得る確率。
TinaW 2016年

1
通常、バイアスされたコインは公正でないことを意味しますが、バイアスがある可能性があります。質問では、コインが完全に公正であるか、常に頭に浮かぶ可能性のみを検討していることを明確にする必要があります。これは非常に現実的なコインではないので、骨壷の問題としてこれを作り直すことをお勧めします。
Neil G

1
「確かに、コインは完全に公正であるか、完全に偏っている。」- あんまり。完全に公正な、または「完全に偏った」硬貨はほとんどありません。
Neil G

1
両側首
TinaW 2016年

回答:


11

あなたのアプローチb)は間違っています:すべてのデータが一緒に使用されて前のものを更新して事後に到達するシングルステップ更新データが1つずつ使用されるベイズ逐次再帰とも呼ばれる)更新の両方後続の反復の前になる事後を取得するには、まったく同じ結果を与える必要があります。これは、ベイジアン統計の柱の1つである一貫性です。

エラーは単純です。最初のサンプル(最初の「Head」)で以前のサンプルを更新すると、新しい以前のサンプルを更新するために、見込みに含めるサンプルは残り1つだけになります。数式では:

P(F|HH)=P(H|H,F)P(F|H)P(H|H)

この式はベイズの定理にすぎず、最初のイベント「Head」がすでに発生した後で適用されます。条件付き確率はそれ自体確率であるため、ベイズの定理はイベント「Head」に条件付けされた確率にも有効であり、実際に証明するものはありません。 。しかし、この結果が自明ではない場合があることに気づいたため、少し長めの証拠を提示しました。

P(F|HH)=P(HH|F)P(F)P(HH)=P(H|H,F)P(H|F)P(F)P(HH)

条件付き確率の連鎖規則による。次に、分子と分母にP(H)、あなたは得る

P(H|H,F)P(H|F)P(F)P(HH)=P(H|H,F)P(H|F)P(F)P(H)P(HH)P(H)=P(H|H,F)P(H)P(HH)P(H|F)P(F)P(H)=P(H|H,F)P(H|H)P(H|F)P(F)P(H)=P(H|H,F)P(F|H)P(H|H)

最後のステップでは、ベイズの定理を適用しました。今:

P(H|H,F)=P(H|F)=0.5

これは明白です:条件付きでコインが公正(またはバイアス)であることを条件として、コイントスをiidとしてモデル化しています。この同じアイデアを分母に適用すると、次のようになります。

P(H|H)=P(H|F,H)P(F|H)+P(H|B,H)P(B|H)=P(H|F)P(F|H)+P(H|B)P(B|H)=0.50.3¯+10.6¯

最終的に:

P(F|HH)=P(H|H,F)P(F|H)P(H|H)=0.50.3¯0.50.3¯+10.6¯=0.2

QED


それだけです。ベイジアン逐次更新を使用して楽しんでください。多くの状況で非常に役立ちます。詳細を知りたい場合は、インターネット上に多くのリソースがあります。これは非常に便利です。

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.