Stan: 2重積分

松浦さんの「階層ベイズモデルとWAIC」にある、シンプソン法による数値積分を、2重積分に拡張してみました。

【追記 やっぱりこんなに単純にはいかないそうです。モンテカルロ積分がよいかも。（松浦さんに多謝。）】
Stanのコードです。

まだじゅうぶんにはテストしていませんが、とりあえずは正しい値がでるようです。

タグ：STAn