安定分布

提供: miniwiki
移動先:案内検索
安定分布
確率密度関数
対称安定分布のPDF
特性指数(α)を変えた場合。μ:位置母数, c:尺度母数
非対称安定分布のPDF
歪度指数(β)を変えた場合。μ:位置母数, c:尺度母数
累積分布関数
対称安定分布のCDF
特性指数(α)を変えた場合。μ:位置母数, c:尺度母数
非対称安定分布のCDF
歪度指数(β)を変えた場合。μ:位置母数, c:尺度母数
母数

α ∈ (0, 2] — 特性指数
β ∈ [−1, 1] — 歪度指数
γ ∈ (0, ∞) — 尺度母数

δ ∈ (−∞, ∞) — 位置母数
α < 1 かつ β = 1 の場合
x ∈ [δ, +∞)
α < 1 かつ β = −1 の場合
x ∈ (-∞, δ]
その他の場合、xR(実数)
テンプレート:確率分布/リンク 密度 特別なケースを除き
解析的な数式表現は不可能
累積分布関数 特別なケースを除き
解析的な数式表現は不可能
期待値 α > 1 のとき、δ
その他の場合、平均は存在しない。
中央値 β = 0 のとき、δ
その他の場合、数式表現不可
最頻値 β = 0 のとき、δ
その他の場合、数式表現不可
分散 α = 2 のとき、2γ
その他の場合、∞ (無限大)
歪度 α = 2 のとき、0
尖度 α = 2 のとき、0
エントロピー 特別なケースを除き
解析的な数式表現は不可能
モーメント母関数 undefined
特性関数 本文参照のこと
テンプレートを表示

安定分布(あんていぶんぷ、: stable distribution) は、正規分布コーシー分布を含むより広い概念であり、安定分布に従う確率変数の和は適当な一次変換によって元の分布になる。正規分布やコーシー分布は安定分布の特別な場合である。安定パレート分布 (: stable pareto distribution)、レヴィ分布 (: Lévy distribution) とも呼ばれる。

定義

退化分布を除き、次の性質を満たす分布は安定分布である。 テンプレート:Ref harvard.

X1X2 を確率変数X の独立な複製(copy)とする。確率変数 aX1 + bX2 (ab は定数)が cX + d (cd は定数)と同一な分布であるとき、確率変数 X安定 であると言い、d = 0 のとき、この分布は 厳密に安定 であると言う。
[math] aX_1 + bX_2 \; \overset{\underset{\mathrm{d}}{}}{=} \; cX + d [/math]

確率密度関数

安定分布の確率密度関数を解析的に書くことはできないが、特性関数 ψ(t) を用いて次のように書くことができる。

[math] f(x)=\frac{1}{2\pi}\int_{-\infty}^\infty \varphi(t)e^{-ixt}\,dt [/math]

これを利用して数値計算(数値積分)が可能である。

特性関数

分布の特性関数 ψ(z) は、4つのパラメータ α, β, γ, δ によって以下のように表すことができる。

[math] \varphi(z) = \exp \left[ i \delta z - \gamma |z|^\alpha \left\{ 1 + i \beta \sgn(z) \omega(z, \alpha) \right\} \right] [/math]
[math] \omega(z, \alpha) = \left\{ \begin{matrix} \tan \frac{\pi \alpha}{2} & (\alpha \neq 1) \\ \frac{2}{\pi} \log |z| & (\alpha = 1) \end{matrix} \right. [/math]

ただし、0 < α ≤ 2, -1 ≤ β ≤ 1, γ > 0 、 sgn(x) は x の符号を返す関数。
α は特性指数と呼ばれ、0 < α ≤ 2 の範囲の値をとる安定分布を特徴づける最も重要な量である。安定分布の指数という場合は通常この α のことを指す。α は分布の裾の厚みの尺度であり、小さいほど裾が広い。歪度指数、あるいは非対称パラメータとも呼ばれる β は分布の対称性を支配し -1 ≤ β ≤ 1 の値をとり、β = 0 のときは左右対称な分布となる。位置母数 δ は分布全体を平行移動するパラメータである。規模母数 γ は X の縮尺を変更するパラメータである。

特別なケース

正規分布

α = 2 の場合、(この場合、βは分布に影響を与えない)

[math]\varphi(z) = \exp \left(i \delta z - \gamma z^2 \right)[/math]

となる。これは、平均 δ、分散 2γ の正規分布である。

Holtsmark分布

α = 3/2, β = 0 の場合、Holtsmark distributionEnglish版 となる。この分布は一般超幾何関数English版を使用することにより確率密度関数を記述できる。

コーシー分布

α = 1, β = 0 の場合

[math]\varphi(z) = \exp \left( i \delta z - \gamma |z| \right) [/math]

となる。これは中央値δ、尺度母数γのコーシー分布である。

(狭義)レヴィ分布

α = 0.5, β = 1 の場合

[math] \varphi(z) = \exp \left[ i \delta z - \gamma \sqrt{|z|} \left\{ 1 + i \sgn(z) \right\} \right] [/math]

となる。これは(狭義)レヴィ分布である。

一般化中心極限定理

中心極限定理では、独立で同一の確率分布(ただし分散は有限に限る)に従う確率変数の算術平均の確率分布は、変数の数が多くなるにしたがい正規分布に収束するが、安定分布において 0 < α < 2 の場合は分散が無限大となり、正規分布には収束せず安定分布 [math]\varphi(x;\alpha,0,c,0)[/math] に収束する。 テンプレート:Ref harvard

参考文献

  • Alder, R. J., R. E. Feldman and M. S. Taqqu eds. (1998年) A Practical Guide to Heavy Tails: Statistical Techniques and Applications


外部リンク


関連項目