MATLAB 関数で MH sampling

MATLAB の組み込み関数 mhsample
jp.mathworks.com
を使ってみた。目標分布や提案分布の引数を一つしかとれないのが制限きついなあ。データとか、提案分布のパラメータを global で渡すしかないのだろうか。かっこ悪いね。
尤度や提案分布の関数をどうせ自分で定義しなければならないのだから、MH くらいなら自分で書いても手間はほとんど変わらない。ループを書く手間が増えるくらいだ。
まあ、書き方がわかった、ということで。MATLAB の練習にはなりました。
つぎは HMC に取り組んでみたい。
jp.mathworks.com