マルコフ連鎖モンテカルロ(MCMC)法#
マルコフ連鎖#
ある状態から、別の状態に推移する確率だけが決まっている
1期前だけで決まる
連鎖はずっと続く(0%/100%がない)
マルコフ過程では、確率変数
となる
定常分布#
各状態にたどり着く確率が変化しなくなってくる。これを定常分布という
任意の状態A,Bについて、Aが生じる確率
を**詳細つり合い条件(detailed balance condition)**という。
詳細釣り合い条件はマルコフ連鎖が定常分布に収束するための十分条件とされる。
MCMCを使う場合は,この条件を満たすように 推移確率
メトロポリス・ヘイスティング(MH)法#
最も単純な方法
だいたいうまくいかない
アルゴリズムの概要:
推移確率行列にまず適当な確率分布を置く
正規分布がよく使われる
その提案分布を補正する処置を施す
提案分布からの採用率を、尤度と事前分布から計算する
正規分布は左右対称→
ハミルトニアン・モンテカルロ(HMC)法#
ハミルトニアン方程式を使って推移確率行列を構成 • 勾配情報をつかい、運動エネルギーと位置を更新していく
NUTS#
not u-turn
HMCの改良版
更新頻度とハイパーパラメータを自動調整
Stanに実装されている