syou6162 https://blog.hatena.ne.jp/syou6162/ yasuhisa's blog https://www.yasuhisay.info/ C++ MCMC 飽きもせずにGibbs Samplingとかばっかりやってますが、久しぶりにはまった。離散確率分布とかからサンプリングするのにrandとかを使ってたんだけど、これがとてもとてもとてもいけなかった。CとかC++のrandは線形合同法で実装されているとかで周期性が問題になる場合があるが、これがとてもとてもとても問題になった。対数尤度を見てるとそれなりに収束したのかなと思ってとあるF値を観察していたんだけど、こんな感じの推移を示した。 ちなみに、これは10個の平均をplotしたものである。平均してこれである。死んでしまえ。200くらいで周期があるように見うけられたので、(半日くらい経って)boost… 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Fwww.yasuhisay.info%2Fentry%2F20101218%2F1292666808" title="とあるモデルのMCMC - yasuhisa&#39;s blog" class="embed-card embed-blogcard" scrolling="no" frameborder="0" style="display: block; width: 100%; height: 190px; max-width: 500px; margin: 10px 0px;"></iframe> https://img.skitch.com/20101218-p71p7xyj69yry4ngywtbny17g4.jpg Hatena Blog https://hatena.blog 2010-12-18 19:06:48 rich https://www.yasuhisay.info/entry/20101218/1292666808 1.0 100%