動き出す次世代サンプラー
区分確定的モンテカルロ
Slide
PDMP
Julia
2025-02-06
PDMP による非絶対連続分布からのサンプリング
司馬博文
5/23/2025
PDMP とはマルコフ連鎖,拡散過程に続いて,近年活発にモンテカルロ法に利用されつつある連続時間マルコフ過程のクラスである.より速く収束するサンプラーが構成しやすいこと,モンテカルロ推定量にバイアスを導入しないようなサブサンプリング(バッチ実行)が可能であることから,近年特に統計・機械学習の分野でも注目が高まっている.
本ポスター発表ではさらに,「 \(\delta_x\) 部分を持った非絶対連続確率分布からも正確なサンプリングが可能」という PDMP に基づくモンテカルロ法の第3の美点に焦点を当て,特にベイズ変数選択を応用として取り上げる.
Date | Location |
---|---|
May 23rd, 13:00-15:00 | ISM 1F |
A Blog Entry on Bayesian Computation by an Applied Mathematician
$$
$$
spike 近傍から脱出するのにかかる時間は \(O_p(1)\) であるのに対して,到着時刻 \(T_\epsilon\) は \(O_p(\epsilon)\) のオーダーである.
従って,\(O(\epsilon^{-1})\) の回数だけ \(T_\epsilon\) のシミュレーションを繰り返すことでやっと \(0\) の近傍から脱出できる.
PDMP とそのシミュレーションに関しては,次のスライドにわかりやすく解説されています:
図は全て発表者開発のパッケージ PDMPFlux.jl
によるものです.
Sticky PDMP (Bierkens et al., 2023) に関するさらに詳しい内容,またはベイズ変数選択一般については,次の記事にまとめています:
@online{2025,
author = {, 司馬博文},
title = {ベイズ変数選択の計算的解決},
date = {2025-05-23},
url = {https://162348.github.io/posts/2025/Posters/OH.html},
langid = {en},
abstract = {PDMP
とはマルコフ連鎖,拡散過程に続いて,近年活発にモンテカルロ法に利用されつつある連続時間マルコフ過程のクラスである.より速く収束するサンプラーが構成しやすいこと,モンテカルロ推定量にバイアスを導入しないようなサブサンプリング(バッチ実行)が可能であることから,近年特に統計・機械学習の分野でも注目が高まっている.
本ポスター発表ではさらに,「 \$\textbackslash delta\_x\$
部分を持った非絶対連続確率分布からも正確なサンプリングが可能」という
PDMP
に基づくモンテカルロ法の第3の美点に焦点を当て,特にベイズ変数選択を応用として取り上げる.}
}