書籍 「 データ解析のための統計モデリング入門――一般化線形モデル・階層ベイズモデル・MCMC (確率と情報の科学) 」の 3章 「一般化線形モデル(GLM)」 と 9章 「GLMのベイズモデル化と事後分布の推定」 で説明されていたポアソン回帰を下記のような 3通りで試してみました。 (1) GLM によるポアソン回帰 (glm 関数) (2) MCMCpack を使ったベイズ統計によるポアソン回帰1 (MCMCpoisson 関数) (3) MCMCpack を使ったベイズ統計によるポアソン回帰2 (MCMCmetrop1R 関数) 書籍では、R から WinBUGS を呼び出して MCMC サンプリングを行っていましたが、今回は R 上でベイズ統計解析を実施する MCMCpack パッケージを試してみました。 サンプルソースは http://github.com/fits/try_sa
![R でポアソン回帰 - glm, MCMCpack - なんとなくな Developer のメモ](https://cdn-ak-scissors.b.st-hatena.com/image/square/f145bb8710df80312164e860bc23f2e43999faba/height=288;version=1;width=512/http%3A%2F%2Fecx.images-amazon.com%2Fimages%2FI%2F51wEV2xJA1L.jpg)