n.iter = 9000, n.burnin = 1000, n.thin = 20
でやっていたんだけど (あれ? 400 しかとってなかったか),
n.iter = 26000, n.burnin = 1000, n.thin = 50
としてみる.
3 回ほどくりかえしてみよう.
n.iter = 51000, n.burnin = 1000, n.thin = 50
barchart()
で図示してみる.
統計モデルがうまくいかなくなったら,
とにかく「よけーな加工」をはぶいたシンプルな図示を!
barchart()
のような「データにできるだけ手を加えないで図示」
したものを,
じっと見てると自分の統計モデルがどのあたりでしくじっているか,
よくわかる.