適切な事前および指数化された可能性は、不適切な事後につながる可能性がありますか?


11

(この質問は西安からのこのコメントに触発されています。)

事前分布が適切で、尤度が明確である場合、事後分布はほぼ確実に適切です。π(θ)L(θ|x)π(θ|x)π(θ)L(θ|x)

場合によっては、代わりに、調整された、または指数化された尤度を使用して、疑似事後

π~(θ|x)π(θ)L(θ|x)α
for(たとえば、これには計算上の利点があります)。α>0

この設定では、適切な事前確率を持つことは可能ですが、疑似事後は不適切ですか?


2
実際、数分後、以前のx尤度^α積を考慮すると、以前のx尤度積の発散が減少するため、それはありそうもないと考えます。そして、よりゆっくりとゼロになる条件は、適切な事前計算によって制御されます。したがって、私の賭けは、これは不可能であるということです。(警告:私は間違っていることが知られています!)
Xi'an

1
場合の反例を探すのに役立つ可能性があります。マルコフの不等式から、したがって、に多項式の裾がある場合を見つけることができれば、不適切な疑似事後を構築できる可能性があります。α>1
Eθπ[L(x|θ)α]tαPθπ(L(x|θ)>t)Eθπ[L(x|θ)α]supt>0tαPθπ(L(x|θ)>t)
L(x|θ)
πr8

この引数はでも機能しますか?また、この方法で構築された可能性が適切であることを証明する方法はありますか?α<1
InfProbSciX 2018年

1
実際、場合、であることがわかっているため、RHSの上限は常に有限であり、、あなたはジェンセンの議論を使って同じ控除をする。したがって、その点で議論は失敗します。この引数が成功するためには無限の尤度が必要であること、つまり、すべてのに対してが必要であるという小さな注意。α=1Eπ[L(x|θ)]<α<1LPπ(L(x|θ)>t)>0t
πr8

1
確かに、場合、1つ作成することはできません。私は言わざるを得ません、私は無限の可能性の例を見るために魅了されるでしょう!おそらく、ベータ後方は無制限の可能性の結果でしょう。α=1
InfProbSciX 2018年

回答:


7

以下のために、おそらくこのような事後を構築することは不可能であることを示すための引数のですか?α1

可能かどうかを確認します。π~(θ|x)dθ=

RHSについて:

π(θ)Lα(θ|x)dθ=Eθ(Lα(θ|x))

場合、は凹関数なので、ジェンセンの不等式により、α1xα

Eθ(Lα(θ|x))Eθα(L(θ|x))=m(x)α<

... 西安が指摘したは、正規化定数(証拠)です。m(x)


きちんと、ありがとう。に対して事後が適切であるという事実を使用していることが好きです。α=1
ロビンライダー

1

結果を@InfProbSciXの回答で使用して、結果を一般的に証明することができます。をとして書き換え 場合我々はそれを知っているので、我々は、上記のJensenの不等式ケースを持っている normalisableです。同様に、場合、 を、これも同じケースに該当しますは正規化可能であることがわかっているためです。これで、(強い)帰納法を使用して、ケースを一般的に示すことができます。L(θx)απ(θ)

L(θx)α1L(θx)π(θ).
1α2L(x|θ)π(θ)2α3
L(x|θ)αpL(x|θ)pπ(θ),
1p2L(x|θ)pπ(θ)

古いコメント

これが非常に便利かどうかはわかりませんが、コメントできないので、回答に残しておきます。@InfProbSciXのに関する優れた発言に加えて、であるとさらに仮定すると、適切な事前確率を持つことはできませんが、疑似事後が不適切になります。以下のための。たとえば、の2番目(番目)のモーメントが存在することがわかっている場合、それは()にあることがわかり、したがって、擬事後は。これらのノートのセクション1α1L(θx)Lp1<αppL(θx)L2Lp0α2もう少し詳しく説明しますが、残念ながら、たとえば PDF のクラスがどれほど広いかは明確ではありません。ここで順番を変えて話してしまった場合は申し訳ありません。コメントとして残しておきたかったのです。L10


1
そうです、尤度関数が空間内にある場合、つまり、空間は、事前に誘導された測度に対して、後部はに対して適切です。私はここで完全に推測していますが、空間は私たちが考えることができるほとんどの可能性を包含していると思います-私が何年も前にがリーマン可積分である場合、その正の力も同様であるという証明を読んだかもしれません。は積分可能ですが。参考のために定理1.26L(θ|x)Lp(πθ)Lp1αpffn,nZ+
InfProbSciX

@InfProbSciX、私はここに影に潜んでいる完全な証拠があるかもしれないと思います。私はあなたの答えから、は負になる可能性があると考えています。それが正しい場合、、積分可能関数の逆数が積分可能であるため、疑似尤度が積分可能になることを示すことができます。尤度が可積分である場合、事前分布が有界であり、可積分と有界関数の積が可積分であるため、事後は可積分であると私は主張します(math.stackexchange.com/a/56008/271610)。どう考えているか教えてください。αp>1
Luiz Max Carvalho

1
場合は、質問で明確に想定されているため、無視してかまいません。一般的な場合のの可積分性を示す必要があります。また、事前分布が常に制限されているかどうかもわかりません。たとえば、の密度はそうではありません。α<0LαBeta(0.5,0.5)
InfProbSciX 2018年

@InfProbSciX、私が意味したのは、が問題ではなくても、その条件が成り立つ場合、が可積分であるという事実を利用して可積分性を示すことができるということですはです。あなたが言うように、事前分布が無制限である場合、それらすべてはゼロです。代わりに、可能性を制限することもできます。MLEで使用する可能性は、制限されるか、強い凹状(en.wikipedia.org/wiki/Maximum_likelihood_estimation#Properties)である必要があります。一般的な証拠を構築します。何かご意見は?α<0α>1f1/f
Luiz Max Carvalho

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