これの続き。多変量の正規分布でギブスサンプリングしたい。 あるデータ点 とそのサイズ 、平均 、分散共分散行列 がある。確率密度関数 は以下となる。 今、データ とそのサイズ がある。ギブスサンプリングでパラメータを推定するには、パラメータ 、 の条…

特に東京に行ったこともないのですが、東京で一番緑が多いのはどこかがふと気になった次第です。将来引っ越しすることがあったら参考になるかもしれないので調べます。 やりかた 行政区域データを集める 土地利用データを集める 各行政区域について土地利用…

階層ベイズモデルで勝敗データからプロ棋士の強さを推定する,StatModeling Memorandumという実験があります。大変面白いので私もやってみました。 データ収集 まず将棋のプロ棋士の対局結果を将棋連盟 棋士別成績一覧から収集しました。こちらのサイトは1964…

ドラゴンボールの勝敗結果から戦闘力を推定する, StatModeling Memorandumという実験記事があります。面白そうなのでやってみました。 準備 RとStanを使いました。tidyrとdplyrはRで使われる変な名前のデータ整形ライブラリです。 set.seed(1) library(rstan…

ギブスサンプリングは確率分布からのサンプリング手法の一つ、各変数について条件付き確率さえわかればサンプリングできる。ベイズ推定での事後分布サンプリングに使われることが有る(Just another Gibbs samplerというそのものズバリなベイズ推定のオープン…