186(lab)

2008-07-27

発表の後始末続き 20:11

前回は→g:lab:id:smoking186:20080725:1216995061

google:マルコフ連鎖モンテカルロ法 二項分布で引っかかったので数解研の来嶋秀治さんの博論を軽く読んだ.

n-1次元単体上の分布に関するもの. 近似サンプリングの方だと, 2次元単体 (直線) 上の整数格子点上でのサンプリングをする必要がある. 例の分布を2次元上に拡張してやると上手くいくように見える. (具体的には, xの方を通常の分布に, yの方を一様分布にしてやればよい. 出てきたサンプルのx成分だけを出力すれば所望の分布に従うものが得られる.) 近似の方のマルコフ連鎖を見ると, Step3でLをサンプリングするが, これがまさに例の分布からのサンプリングになってしまうので解決になってない.

次は非負整数の方を探すことになるようだ. こちらを調査しよう.