階層モデルだとWBICの推定値が上手くいかないようなのですが、
周辺化する以外にいい方法を知っている人がいたら教えてもらえると嬉しいです。
またシミュレーションを追加して資料を修正するかもしれません。
コード
上のスライドで使ったRコードは以下です。
[crayon-60fc2584d42da169876246/]
次に、Stanコードです。
bridgesamplingのため、すべてtarget記法で書いています。なお以下のコードはWBIC用のコードなので、bridgesamplingをする場合はmodelの1.0/log(N)をかけている部分を省いてください。
最初に階層モデルです。
[crayon-60fc2584d42e6837154931/]
次に、周辺化モデルです。
[crayon-60fc2584d42ef920379474/]