今日の大台ヶ原MCMC

大台ヶ原階層ベイズモデルMCMC計算の鎖の数を3つに増やし、Gelman-Rubinの収束判定を行なう。うまく1に近い値となった。

なんだか目鼻がつきそうになってきたので、LaTeXで、方法に相当する部分を書いてみる。BUGSのコードだと割とすっきりしているように見えるが、数式にしてみるとなんだか分かりにくい。というか、あっているか不安である。というより、たぶん間違っている気がする。