チラシの裏の落書き日記

統計とか,研究とか,日常についての備忘録的なもの。

2017-12-21から1日間の記事一覧

JAGSで対数尤度を得る方法(Stanも)

R

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