虫データ階層ベイズ

某所、というか大台ヶ原の虫データの解析。ご存知の方はご存知のように結構複雑な実験系なので、解析もそれなりに複雑になる。ゆーい差決戦主義(Ⓒ久保さん)でやるよりは、階層ベイズにもっていってしまった方がよかろうと、コードを組んでみる。試行錯誤の末、かなり良くなってきたが、それでもgelman.diag()で1.18という微妙な収束具合。トレースを見ると、結構いいような気がするのだが。とりあえずは繰り返し回数を増やしてみるか。

(追記) 回数を増やしたら、最大1.05にまで減少した。

タグ：階層ベイズ