非対称の事前知識を持つ多項分布のベイズ推定?


8

二項分布からいくつかのサンプルを取得するとします。私の以前の知識をモデル化する1つの方法は、パラメーターおよびベータ分布を使用することです。私が理解しているように、これは試験で「頭」を回見たのと同じです。そのため、本格的なベイジアン推論を行うための良い近道は、回の試行で頭を見た後の「頭」の確率の新しい平均としてを使用することです。β α α + β H + ααβαα+β HのNh+α+α+βh

ここで、3つ以上の状態があると仮定します。そのため、多項分布からいくつかのサンプルを取得します。事前分布としてパラメーターを使用したディリクレ分布を使用するとします。ここでもショートカットとして、これをイベントの確率の事前知識としてと同等に扱うことができます回の試行でイベント回を目撃した場合、私の事後なります。。I α Iα IHNIH+αIαΣαj hh+α+Σαj

今二項の場合には、それは「頭」の事前知識が発生していることをうまくいくで時間を裁判起こる「尾」に相当しますで時間をトライアル。論理的には、「尾」よりも「頭」の可能性についてより強い知識を持つことができるとは思いません。ただし、これは2つ以上の結果を伴ってより興味深いものになります。私が6面ダイスと言った場合、50トライアルではサイド1の事前知識は10に相当し、100トライアルではサイド2の事前知識は15 2に相当すると想像できます。α + β β α + βαα+ββα+β

それで、すべての紹介の後で、私の質問は、多項式の場合にそのような非対称の事前知識を適切にモデル化する方法ですか?注意しないと、合計確率/尤度が1にならないため、簡単に非論理的な結果が得られるようです。ディリクレショートカットを引き続き使用できる方法はありますか、これを完全に犠牲にして使用する必要がありますか。他の事前配布は完全に?

上記の表記法や用語の乱用によって引き起こされる混乱を許してください。


If I have say a 6-sided die, I can imagine my prior knowledge of side "1" being equivalent to 10 ones in 50 trials and my prior knowledge of side "2" as being equivalent to 15 twos in 100 trials.最小公倍数を使用して「」と「2」(つまり50100)の「同等のサンプルサイズ」を統一したくない特別な理由はありますか?この場合にはα 1 = 20 / 100およびα 2 = 15 / 1001250100α1=20/100α2=15/100
Zhubarb 2013年

それは確保しながら、あなたが非対称事前分布を表すことができますΣ=1=6α=1
Zhubarb

最後に、ケンプらによるこの論文。とても関連性があり、読むのが楽しいです。
Zhubarb 2013年

@Berkan私はまだ紙を読んでいませんが(ありがとう!)サンプルサイズを統一しなかった理由は、20/100が10/50と同じではないためです。その事実は、この質問が尋ねられた理由。
Michael McGowan 2013年

同等性が高いため、前者の方が自信があるからです。サンプルサイズ?
Zhubarb 2013年

回答:


2

あなたはあなたの質問を非常にうまく組み立てました。

ここで探しているのは階層モデリングの場合だと思います。また、階層の複数の層をモデル化することもできます(現時点では、事前についてのみ話します)。ハイパーパラメーターのハイパー優先度の別のレイヤーがあると、ハイパーパラメーターの追加の変動性をモデル化できます(ハイパーパラメーターの変動性の問題が懸念されるため)。また、モデリングが柔軟で堅牢になります(遅くなる場合があります)。

特にあなたのケースでは、ディリクレ分布パラメーターの事前情報を取得することでメリットが得られます(ベータは特殊なケースです)。Gelmanによるこの投稿では、ディリクレ分布のパラメーターに事前分布を課す方法について説明しています。彼はまた、毒物学のジャーナルで自分の論文を引用しています。


もっとモデリングに特化することもできますが、ここでは推測したくありません。質問の変更を選択した場合は、詳細を追加します。
suncoolsu
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.