順序統計量の確率密度関数


シミュレーションを行う際に,事前分布から得た値を小さい順に並び替えて用いるというケースがあります[A]。このような場合には,並び替えた値が従う確率密度関数が必要になります。より具体的に書けば,確率密度関数が ,累積分布関 もっと見る

コルモゴロフ-スミルノフ検定


異なる 2 つのデータセットが異なる確率分布からサンプリングされたものであるかを知りたい場合に,経験分布の検定という統計学的手法があります。検定はいくつか種類がありますが,その 1 つにコルモゴロフ-スミルノフ検定 (K もっと見る

F# でモンテカルロ法


モンテカルロ法は乱数を用いる数値計算手法です。 F# の練習のためにモンテカルロ法を F# で実装してみました。計算するのはおなじみの円周率ですが,よく見かけるモンテカルロ積分[A]ではなく,棄却サンプリング法によりレイ もっと見る

近似ベイズ計算の実装例


『近似ベイズ計算によるベイズ推定』の逐次モンテカルロ法を R で実装してみました。 Sisson et al. (2007) で例示されている 2 つの正規分布 , の混合分布からのサンプリングを行います。

近似ベイズ計算によるベイズ推定


近似ベイズ計算 (approximate Bayesian computation, ABC) と呼ばれるベイズ的手法があります。近似ベイズ計算の「近似」たるゆえんは,複雑な尤度の計算を行わないために,近似計算を行うとい もっと見る