餡子付゛録゛

ソフトウェア開発ツールの便利な使い方を紹介。

2024-10-01から1ヶ月間の記事一覧

最尤法による欠損値の補定と回帰モデルの推定の同時実行

ここ10年間ぐらいで多重代入法を用いた欠損データ処理が一般化しましたが、多重代入法を使わなくても多重代入法と同様の結果が得ることもできます。Fully Bayesian Imputation Methodと呼ばれている手法ですが*1、最尤法による単一代入(imputation model)…

R/FortranでHamiltonian Monte Carlo法によるベイズ推定を書いてみたよ!

マルコフ連鎖モンテカルロ法(MCMC)には様々なアルゴリズムがあり、メトロポリス・ヘイスティングス法は尤度関数があれば推定ができる使い勝手のよい方法だと知られています。しかし、高次元で効率が悪いことが知られており、また多峰の確率分布には上手く…