logo

ビンガム-マルディア 分布 📂確率分布論

ビンガム-マルディア 分布

定義 1

ユニークモードunique mode μSp1\mu \in S^{p-1}集中concentration κ>0\kappa > 0 そして半径 ν>0\nu > 0 について、次のような確率密度関数を持つ多変量分布 BMp(μ,κ,ν)\text{BM}_{p} \left( \mu , \kappa, \nu \right)ビンガム・マーディア分布bingham-Mardia distributionと言う。 f(x)=1α(κ,ν)exp(κ(μTxν)2),xSp1 f \left( \mathbf{x} \right) = {{ 1 } \over { \alpha \left( \kappa , \nu \right) }} \exp \left( - \kappa \left( \mu^{T} \mathbf{x} - \nu \right)^{2} \right) \qquad , \mathbf{x} \in S^{p-1} ここで、α(κ,ν)>0\alpha \left( \kappa , \nu \right) > 0Sp1f(x)dx=1\int_{S^{p-1}} f(\mathbf{x}) d \mathbf{x} = 1 となるようにする正規化定数normalizing constantだ。


説明

ビンガム・マーディア分布は、球面上で小さい円small Circle形のクラスターを成す確率分布だ。

フォンミーゼス・フィッシャー分布の確率密度関数: f(x)=(κ2)p/211γ(p/2)Ip/21(κ)exp(κμTx),xSp1 f \left( \mathbf{x} \right) = \left( {{ \kappa } \over { 2 }} \right)^{p/2-1} {{ 1 } \over { \gamma \left( p/2 \right) I_{p/2-1} \left( \kappa \right) }} \exp \left( \kappa \mu^{T} \mathbf{x} \right) \qquad , \mathbf{x} \in S^{p-1}

球面上の正規分布のような感じだったフォンミーゼス・フィッシャー分布と比較してみれば、ビンガム・マーディア分布の確率密度関数では、κ(μTxν)2\kappa \left( \mu^{T} \mathbf{x} - \nu \right)^{2} が円形を成す役割をするという点が容易に理解できるだろう。


  1. Kim. (2019). Small sphere distributions for directional data with application to medical imaging. https://doi.org/10.1111/sjos.12381 ↩︎