今回はRの話です。 社会心理学会の方法論セミナーでもGLMMをとりあげましたが,階層ベイズの話も久保先生のトークの中にありました。 GLMMでは,変量効果の数が増えると最尤法だと推定が難しくなったりするので,ベイズ推定のほうが向いています。 しかし,GLMMを直接ベイズ推定してくれる商用ソフトもあまりない(あっても機能が部分的)ので,stanなどのフリーのソフトに頼らざるを得ません。しかし,stanは初心者にはなかなかとっつきにくいので,今回はだれでも簡単にGLMMがベイズ推定できる関数を作ってみました。 実は過去に,同様にGLMMを簡単にstanで走らせてくれるglmer2stanというパッケージを紹介したことがあります(こちら)。しかし,glmer2stanはあとで挙げるようにいろいろ使い勝手が悪いところもあり,自分用に使いやすいものを作ろうと思ったのがはじまりです。 追記: Sapp