事前確率分布をどのように形式化しますか?経験則や使用すべきヒントはありますか?


9

ベイジアン統計分析と意思決定における事前情報の概念をよく理解していると思いたいのですが、そのアプリケーションに頭を悩ませるのに苦労することがよくあります。私は自分の闘争を例証するいくつかの状況を念頭に置いており、これまでに読んだベイジアン統計教科書ではそれらが適切に扱われていないと感じています。

数年前に私が実施した調査で、68%の人がACME製品の購入に関心があると言ったとします。再度調査を行うことにしました。前回と同じサンプルサイズ(たとえば、n = 400)を使用しますが、それ以降、人々の意見は変わった可能性があります。ただし、以前のバージョンとして、400人中272人が「はい」と回答したベータ版の分布を使用する場合、数年前に行った調査と現在実施している調査に同じ重みを与えます。そのデータが数年前のものであるという理由で、以前に置いておきたいより大きな不確実性を確立するための経験則はありますか?以前のものを272/400から例えば136/200に減らすことができると理解していますが、これは非常に恣意的であり、おそらく文献に何らかの形で正当化があるのだろうかと思います。

別の例として、臨床試験を実行しようとしているとしましょう。試験を開始する前に、専門家の意見、以前の臨床試験(関連性の異なる)の結果、その他の基本的な科学的事実など、予備情報として使用できるいくつかの二次調査を実行します。 (そのうちのいくつかは本質的に非定量的です)事前の確率分布?データに圧倒されることを確実にするためにどの家族を選び、それを普及させるかを決定するだけのケースですか、それともかなり有益な事前配布を確立するために行われる多くの作業がありますか?


回答:


4

400回の試行で272件の成功に関する以前の情報を処理するというあなたの考えは、か​​なり確かなベイズの正当化を持っています。

θ

π(θ)=Γ(α0+β0)Γ(α0)Γ(β0)θα01(1θ)β01
n_=α0+β02n_α01
π(θ)=Γ(α0+β0)Γ(α0)Γ(β0)θα01(1θ)n_(α01)
α0+β02=400α01=272α0=273β0=129α0=137β0=65
μ=αα+βandσ2=αβ(α+β)2(α+β+1)
alpha01 <- 273
beta01 <- 129
(mean01 <- alpha01/(alpha01+beta01))

alpha02 <- 137
beta02 <- 65
(mean02 <- alpha02/(alpha02+beta02))

からの前の分散を増やします

(priorvariance01 <- (alpha01*beta01)/((alpha01+beta01)^2*(alpha01+beta01+1)))
[1] 0.0005407484

(priorvariance02 <- (alpha02*beta02)/((alpha02+beta02)^2*(alpha02+beta02+1)))
[1] 0.001075066

望んだ通りに。

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