aidiary https://blog.hatena.ne.jp/aidiary/ 人工知能に関する断創録 https://aidiary.hatenablog.com/ 機械学習 PRML Pythonによるモンテカルロ法入門（2014/6/20）の受理・棄却法（2014/7/12）の2回目。前回は、提案分布に一様分布を使ったためサンプリングの効率が悪かった。今回は、より効率のよい提案分布gとしてベータ分布を使ってみよう。目標分布fは前回と同じベータ分布。下の図を見ればわかるように提案分布では目標分布を覆えていない。すっぽり覆えるように提案分布をM倍する必要がある。Mを求めるには前回と同様にscipy.optimizeを使った。 # f(x)/g(x)を最大化するxoptを求める xopt = scipy.optimize.fmin(lambda x: - f(x) / g(x)… 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Faidiary.hatenablog.com%2Fentry%2F20140713%2F1405255289" title="受理・棄却法 (2) - 人工知能に関する断創録" class="embed-card embed-blogcard" scrolling="no" frameborder="0" style="display: block; width: 100%; height: 190px; max-width: 500px; margin: 10px 0px;"></iframe> http://cdn-ak.f.st-hatena.com/images/fotolife/a/aidiary/20140713/20140713213838.png Hatena Blog https://hatena.blog 2014-07-13 21:41:29 rich https://aidiary.hatenablog.com/entry/20140713/1405255289 1.0 100%