尖度

尖度の最新ニュースをまとめて検索!

尖度(せんど、Kurtosis)とは、確率分布や頻度分布のとんがり具合を表す数字。

目次

[編集] 2種類の定義

正規分布の尖度を0とする定義と3とする定義があることに注意。2つの定義の違いは、正規分布との乖離をみるために使われることに起因している。一般には0とすることが多い。Excelの分析ツール等は正規分布の尖度を0としている[1]。英語版のWikipediaでも、正規分布の尖度は0としている。東京大学出版会の「統計学入門」(ISBN 4130420658)やNumerical Recipesなども正規分布の尖度は0である。

[編集] モーメントによる定義

確率変数Xの分布関数をF(X)

\mu=E(X)=\int X dF(x)
\mu_r=E(X-\mu)^r=\int (X-\mu)^r dF(x)(rは正整数)

とする。このとき、分布関数F(X)の尖度β2は以下である。(各積分値が存在すると仮定する)

正規分布の尖度を0とする定義では、

\beta_2=\frac{\mu_4}{{\mu_2}^2} - 3

正規分布の尖度を3とする定義では、

\beta_2=\frac{\mu_4}{{\mu_2}^2}

[編集] キュムラントによる定義

確率変数Xr次のキュムラント[2]κrとすると尖度β2は以下で定義される。

正規分布の尖度を0とする定義では、

\beta_2=\frac{\kappa_4}{{\kappa_2}^2}

正規分布の尖度を3とする定義では、

\beta_2=\frac{\kappa_4}{{\kappa_2}^2}+3

[編集] 計算例

正規分布の尖度。モーメント母関数は

log(MX(t)) = μt + σ2t2 / 2

よりκ1 = μ,κ2 = σ,3次以上のキュムラント\kappa_r=0(r\ge 3)はすべて0であることがわかる。従って正規分布の尖度β2 = 0(または3)となる。

[編集] 尖度の意味

以下に正規分布とそれより尖度が大きくおなじ平均と標準偏差をもつ確率密度関数を示す。

模式的であるが、平均の周りでは尖りが大きく裾を引いた分布であることがわかる。 厳密には、尖度と表現するのは間違いであり、本来は裾の重さというべきであることは明らかであろう。詳細は参考文献を見られたい。

特殊な分布

p(x) =
\frac{1}{\alpha\,\mathrm{\Beta}\!\left(m-\frac12, \frac12\right)}
\left[1 + \left(\frac{x-\lambda}{\alpha}\right)^{\!2\,} \right]^{-m}
\!

では両者(裾の重さと平均の周りでは尖り)の概念は一致するが、一般の分布では一致しないのは明らかである。

[編集] 標本モーメントによる母集団の尖度推定

ここでは、標本の大きさnの標本に基づく母集団の尖度の推定を考える。

厳密な意味での標本における尖度についての不偏推定量の研究はこのウィキペディアの範疇を超えるので省略する。

一般には尖度の定義の分母分子の不偏推定量をもって母集団の尖度の推定量とする方法がもっとも多く使用される。具体的には、母集団のキュムラントの不偏推定量であるk-statistics(k統計量)[3]を使った計算方法である。r次のk統計量をkr、平均周りのr次のモーメントをm_r=\frac{1}{n}\sum_{i=1}^{n}(x_i-\overline{x})^rとすれば、

 k_2=\frac{n}{n-1}m_2
k_4=\frac{n^2}{(n-1)(n-2)(n-3)}\left\{(n+1)m_4-3(n-1){m_2}^2 \right\}

なので、上記2式をキュムラントによる定義に代入して推定量とする方法である。

最終的にはβ2の推定量b2について次を得る。

b_2=\frac{n(n+1)}{(n-1)(n-2)(n-3)} \sum_{i} \left(\frac{x_i-\overline{x}}{s}\right)^4 -\frac{3(3n-5)}{(n-2)(n-3)}
ここにs=\sqrt{\frac{n}{n-1}m_2}(標本標準偏差)

または3を引いた形式では以下となる。

\frac{n(n+1)}{(n-1)(n-2)(n-3)} \sum_{i} \left(\frac{x_i-\overline{x}}{s}\right)^4 -3\frac{(n-1)^2}{(n-2)(n-3)}

[編集] 分布と尖度

分布 尖度
ラプラス分布 -3
双曲線正割分布 2
ロジスティック分布 1.2
正規分布 0
二乗余弦分布 −0.593762…
ウィグナー半円分布 -1
一様分布 -1.2

上記の分布は歪度は全て0である。

[編集] 脚注

  1. ^ 関連項目に挙げたが、歪度(β1 = E(X − μ)3 / [E(X − μ)2](3 / 2))、分散(σ = E(X − μ)2)は不等式による制限関係があるので、注意する必要がある。当然の事ながら標本における歪度、尖度は互いに独立ではない。

    参考までにコーシーシュワルツの不等式を考えれば尖度β2は明らかに1より大きく、3を引いた意味での尖度は-2以上の値をとることは明らかであろう。

  2. ^ 簡単に言えば、モーメント母関数MX(t)の対数をとった関数をt = 0の周りで形式的に展開し、tについて整理し、r次の係数がκrである。
  3. ^ 具体的にはE{kr} = κrとなる性質がある。

[編集] 関連項目

[編集] 参考文献

  • A Dictionary of Statistical Terms, 4th edition(Kendall, Maurice G.;Buckland, William R.)
    • ケンドール 統計学用語辞典 (訳 千葉大学統計グループ) (丸善) ISBN 4621032038

最終更新 2009年11月6日 (金) 15:56 (日時は個人設定で未設定ならばUTC)。
【尖度】変更履歴

ご利用上の注意

もっと調べる!