パッケージを使わないで、一般化線形混合モデルのMCMCアルゴリズムを1から作る.

こんにちは. 金曜日の夜になり、激しめの睡魔に襲われております.先日のこちらの記事で公開したスライドの後半にあるシミュレーションで、地域差を考慮したPoisson - Normalモデルを構築しているのですが、そのコードを載せておきます。ベイズ推論とシミュレーション法の基本について - Data Science by R and Pythontom…