math.luma.dev
検索
/statistics/moment

モーメント

モーメント(moment)の定義

モーメントにはいくつかある。確率変数(random variable) X それぞれに対してモーメントを対応させて考えることができる。

  • k次のゼロモーメント(raw moment): m⁡k=defE⁡[Xk]
  • k次の中心モーメント(central moment): κ⁡k=defE⁡[(X−μ)k]
    • μkと表記されることもあるが、μ≠μ1には注意
  • k次の標準モーメント(standardized moment): μ~⁡k=defE⁡[(X−μσ)k]

それぞれの呼び名「ゼロ」「中心」「標準」にあたる訳語はあまり見当らなかったので、このサイトではそう呼ぶことにする

ゼロモーメントと中心モーメントは、さらに以下のように一般化できる。

モーメント母関数(moment generating function)

モーメント母関数とは、k次のモーメントを数列(number sequence)として生成する指数母関数(exponential generating function)のこと。

ゼロモーメントに関しては特に以下のようにして定義される。

MX(θ)=def∑k=0∞1k!E⁡[Xk]θk

MX の X は文脈から明らかな場合は省略される。
また、 M(θ)=E⁡[exp⁡(Xθ)] というより簡単な表示がある。

証明

M(θ)=∑k=0∞1k!E⁡[Xk]θk=∑k=0∞E⁡[1k!Xkθk]=∑k=0∞E⁡[1k!(Xθ)k]=E⁡[∑k=0∞1k!(Xθ)k](級数と期待値の線形性による入れ替え(注意が必要))=E⁡[exp⁡(Xθ)](指数関数のマクローリン展開)

TODO: 注意の部分

■

期待値(expected value)の定義より、これは以下のように計算される。

  • 離散の場合: M(θ)=∑xexp⁡(xθ)⋅P(X=x)
  • 連続の場合: M(θ)=∫−∞∞exp⁡(xθ)⋅f(x)dx

なお、以降は連続確率分布(probability distribution)の場合のみを記述するが、基本はすべて離散確率分布の場合でも同様。

E⁡[Xk]=((ddθ)kM)(0) として計算できる。 (指数母関数の一般の特徴)

  • 一般のa周りのモーメント母関数: E⁡[exp⁡((X−a)θ)]
    • Y=defX−a で確率変数を新たに定義すれば、 MY(θ) そのものだ。
  • 標準モーメント母関数(standardized moment generating function): E⁡[exp⁡(X−μσθ)]
    • Y=defX−μσ で確率変数を新たに定義すれば、 MY(θ) そのものだ。

モーメント母関数を利用した計算

  • 平均 μ=E⁡[X]=M′(0)
  • 分散(variance) E⁡[X2]−(E⁡[X])2=M′′(0)−(M′(0))2
    • Var⁡(X) とも書き、 κ⁡2 でもある。
    • 通常 σ2 と書かれる。 σ=κ⁡2 のことを標準偏差(standard deviation)と呼ぶ。
    • 非負性: 0≤σ2
  • 歪度 μ~⁡3=E⁡[(X−μ)3]/σ3=κ⁡3(κ⁡2)3/2
    • 負の値をとることがある。
    • 歪度は、確率分布の非対称性(歪み)を表す指標である。
    • 正方向ならば右に歪んでいる、負方向ならば左に歪んでいる。 (TODO: 図解)
  • 尖度 μ~⁡4=E⁡[(X−μ)4]/σ4=κ⁡4(κ⁡2)2
    • 尖度は、確率分布の尖り具合を表す指標である。 (TODO: 図解)
    • 非負性: 0≤μ~⁡4
  • より一般に: TODO

などが得られる。

証明

分散について。

σ2=E⁡[(X−μ)2]=E⁡[X2−2Xμ+μ2]=E⁡[X2]−2E⁡[X]μ+μ2=E⁡[X2]−μ2=M′′(0)−(M′(0))2
■

参考