ギブスサンプリングに関連する混乱


9

私はこの記事に出くわしましたが、ギブスのサンプリングではすべてのサンプルが受け入れられます。私は少し混乱しています。受け入れたすべてのサンプルが定常分布に収束する場合はどうしてでしょうか。

一般的なMetropolisアルゴリズムでは、min(1、p(x *)/ p(x))として受け入れます。ここで、x *はサンプルポイントです。x *が密度の高い位置を指し示しているため、ターゲットの分布に移動していると想定します。したがって、バーンイン期間の後でターゲットの分布に移動すると思います。

ただし、ギブスサンプリングではすべてを受け入れるため、別の場所に移動する可能性がありますが、定常/ターゲット分布に収束していると言えますか?

分布ます。我々はない計算Zに大都市アルゴリズム我々が使用する用語缶C θ N E W/ C θ O L D分布組み込むC θ を加えた正規化定数Zは相殺し。だから大丈夫p(θ)=c(θ)/Zc(θnew)/c(θold)c(θ)

しかし、ギブスサンプリングでは、分布c θ をどこで使用していc(θ)

たとえば、論文http://books.nips.cc/papers/files/nips25/NIPS2012_0921.pdfで

ですから、サンプリングする正確な条件付き分布はありません。条件付き分布に直接比例するものだけがあります

ここに画像の説明を入力してください


2
p(x)/p(x)

回答:


6

Metropolis-Hastingsアルゴリズムを使用する場合、許容率を計算する必要があります

α=min(1,p(x)p(x))
UUniform(0,1)U<α

α=1Uαp(x)。お役に立てば幸いです。

ここに画像の説明を入力してください


1
p(θ)

変数x1、x2、x3のランダムな値から始めて、その定常分布が必要な分布に収束していると言えますか。その基準は何ですか?
user34790 2013

p(θ)=c(θ)/Zp(θ)

3
なぜ常に1であるかの証明を上に追加しました。ギブスサンプリングを使用するには、完全な条件が何であるかを知る必要があります。

1

受け入れ率がタイプミスとして1に等しいことの証明、つまり中央と3番目の部分の分母では、qの式にz_i素数が必要であるため、最終的にP(z_i素数| z_i素数)が得られます。

アレックス

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