Richardson の補外による数値微分


微分は積分に比べて解析的に求まるケースが多いのであまり使われる場面は少ないですが,それでも数値計算的に微分係数が求められると嬉しいです。簡易的に適当に小さい を持ってきて で計算することができるのですが,精度が良くありま もっと見る

F# でデータ拡大サンプリング


データ拡大 (data augmentation) という手法があって,それを MCMC に応用したサンプリング法があります。 分布 からの直接サンプリングが困難である場合, となるような を構築してやります。ここで を もっと見る

F# でハミルトニアンモンテカルロ


ハミルトニアンモンテカルロ (あるいはハイブリッドモンテカルロ) を F# で実装してみたというお話。ハミルトニアンというのはポテンシャルエネルギーと運動エネルギーの和としてあらわされる物理量のことだそうですが,よくわか もっと見る

F# で最高事後密度区間を求める


F# で最高事後密度区間 (highest posterior density interval, HPDI) を求めます。 取りうる区間の中でもっとも幅が狭いものを選択しています。 例として標準正規乱数を 100,00 もっと見る