研究を行う上で非常に重要なのは，やはりどのような問題に決着をつけるべきなのかその問題を見極めることだと思う。問題を設定するために，はじめは指導教員や先輩の研究者などに教えを請うのが望ましい。それによって，何が研究としてあるべき形なのかを身…

IRTモデルで予測分布や対数尤度の計算などをおこなってみる。 対数密度の記述は以下の頁の質問に少し書かれている。 https://sourceforge.net/p/mcmc-jags/discussion/610037/thread/c728e87b/ The general rule is that if y ~ dfoo(alpha, beta, gamma, ..…

やり方は簡単で，予測分布を未知の変数としてモデルに書けばいい。 Stanではgenerated quantitiyブロックでrng系の関数を使えばいいが，JAGSではそのまま"x_pred ~ dnorm(mu, sigma)"などとするとRでいうところのx_pred sigma)のように乱数を生成してくれる…

WAICなどを計算するには，対数尤度のMCMCサンプルが必要。 個々のMCMCサンプルからがんばって計算する方法もあるが，より簡単にはMCMCサンプラーのなかで対数尤度を得ると楽。 Stanではgenerated quantityが容易されているが，JAGSでは少し工夫が必要。参考 …