math.luma.dev
検索
/statistics/moment

モーメント

モーメント(moment)の定義

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

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

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

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

  • k次のa周りのモーメント: E[(Xa)k]

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

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

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

MX(θ)=defk=01k!E[Xk]θk

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

証明

M(θ)=k=01k!E[Xk]θk=k=0E[1k!Xkθk]=k=0E[1k!(Xθ)k]=E[k=01k!(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((Xa)θ)]
    • Y=defXa確率変数を新たに定義すれば、 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[X22Xμ+μ2]=E[X2]2E[X]μ+μ2=E[X2]μ2=M(0)(M(0))2

参考