マルコフ連鎖モンテカルロ法（メトロポリス法）による混合ガウス分布に従う乱数の生成

平均(1,1),(-1,-1)となるような２つのガウシアンを3:7の割合で混ぜ合わせた2次元混合ガウス分布に従う乱数をマルコフ連鎖モンテカルロ法の１つであるメトロポリス法で生成するコードを書いた。

実行すると↓のように混合比率から若干(-1,-1)のガウシアンの方に引き寄せられた形になる。