MCMCメソッド-サンプルの書き込み?


12

MCMCの方法、私は約読み続けるburn-in時間やサンプルの数"burn"。これは正確には何ですか?なぜ必要なのですか?

更新:

MCMCが安定したら、安定したままですか?burn-in時間の概念は混合時間のそれとどのように関連していますか?


1
バーンイン」という語句は、サンプルを「バーンイン」することを意味するものではありませんが、システムを「新しい」状態から使用に適した状態に問題なく解決します。(サンプルは落ち着く間に破棄されます、そうですが、このフレーズはその意味で「燃焼」を意味するものではありません。道路に持ち込む前に「エンジンを暖める」のに近いです。)いくつかの議論についてはこの回答を参照してください。それが何のためにあるのか。
Glen_b-モニカを2014

(具体的には、その回答の最初の段落)
Glen_b-モニカを14

1
@Glen_b、私はそれを知っていたと思いますが、私は疑わしいサンプルが煙で上がっていくメンタルなイメージが好きです...
マット・クラウゼ

回答:


17

バーンインは、特に粗い開始点から始まった場合に、マルコフチェーンが平衡分布に到達する時間を与えることを目的としています。チェーンを「バーンイン」するには、ポイントの収集を開始する前に、最初のサンプルを破棄するだけです。n

「悪い」開始点は、それが平衡分布に落ち着く前に、平衡分布の下で実際に非常に低い確率である領域をオーバーサンプリングする可能性があるという考えです。それらのポイントを捨てると、ありそうもないはずのポイントが適切にレアになります。

このページは良い例ですが、バーンインは原理的な手法というよりはハック/アートフォームであることも指摘しています。理論的には、本当に長時間サンプリングするか、まともな開始点を選択する方法を見つけることができます。

編集:混合時間とは、チェーンが定常状態に近づくのにかかる時間を指しますが、直接計算することはしばしば困難です。混合時間を知っている場合は、その数のサンプルを破棄するだけですが、多くの場合、破棄しません。したがって、代わりに十分な大きさのバーンイン時間を選択します。

安定性に関する限り、それは状況によって異なります。チェーンが収束した場合、それは...収束しています。ただし、チェーンが収束したように見えても、実際には状態空間の一部で「ハングアウト」している場合もあります。たとえば、いくつかのモードがあるが、各モードは他のモードとの接続が不十分であると想像してください。サンプラーがそのギャップを通過するまでに非常に長い時間がかかる場合があり、そのジャンプが行われるまで、チェーンが正しく収束したように見えます。

収束の診断はありますが、それらの多くは、真の収束と疑似収束を区別するのに苦労しています。マルコフ連鎖のモンテカルロハンドブックのチャールズガイヤーの章(#1)は、できる限り連鎖を実行することを除いて、すべてについてかなり悲観的です。


ありがとう。コンセプトをさらに明確にすることを期待して、OPに小さな更新を追加しました。
Amelio Vazquez-Reina 14

もう少し追加しました。お役に立てば幸いです。
Matt Krause、2014年

6

Metropolis-Hastingsアルゴリズムは、事後分布からランダムにサンプリングします。マルコフ連鎖は定常分布に対して安定していないため、通常、初期サンプルは完全に有効ではありません。バーンインサンプルを使用すると、まだ静止していない初期サンプルを破棄できます。


1
ありがとう-マルコフ連鎖が安定したら、それは安定したままですか?これは混合時間の概念とどのように関連していますか?同じですか?
Amelio Vazquez-Reina
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.