「キュムラント」の版間の差分
ja>42d8929 細 (数式修正) |
細 (1版 をインポートしました) |
(相違点なし)
|
2018/8/19/ (日) 17:20時点における最新版
確率論や統計学において、キュムラント(きゅむらんと、英: cumulant)は、分布を特徴付ける特性値の一つ。キュムラント母関数を級数展開した際の係数として定義する。その性質を研究したT. N. ティエレに因み[1] 、ティエレの半不変数(英: semi-invariant)とも呼ぶ[2]。
Contents
定義
確率変数Xに対して、
- [math] M(s) := \langle e^{s X} \rangle \, [/math]
で定義されるモーメント母関数の対数log M(s)をキュムラント母関数と呼ぶ。但し、⟨…⟩は期待値を取る操作を表すものとする。キュムラント母関数の級数展開
- [math] \log{M(s)} = \sum_{n=1}^{\infty}\frac{c_n}{n!}s^n [/math]
において、係数cnをn次のキュムラントもしくは、ティエレの半不変数と呼ぶ。この級数展開がn=0の項を含まないことは、M(0) = log 1 = 0よりわかる。キュムラントを表す記号としてcnのほかに、κnやモーメント⟨Xn⟩に対応した⟨Xn⟩cが用いられる。また、確率変数のべき乗Xnに対し
- [math] \langle X^n \rangle_c \, [/math]
を与える操作をキュムラント平均と呼ぶ。
性質
半不変数という呼び名は、確率変数のアフィン変換
- [math] Y=\alpha X + \beta \, [/math]
において、モーメントが
- [math] \langle Y \rangle = \alpha \langle X \rangle + \beta , \,\, \langle Y^2 \rangle = \alpha^2 \langle X^2 \rangle + 2 \alpha \beta \langle X \rangle +\beta^2 , \cdots [/math]
と変換されるのに対し、キュムラントは、
- [math] \langle Y \rangle_c = \alpha \langle X \rangle_c + \beta , \,\, \langle Y^n \rangle_c = \alpha^n \langle X^n \rangle_c \quad n=2,3,\dotsc [/math]
と、ほとんど形を変えないことに因む。
多変数におけるキュムラント
多変数の確率変数X1,..., Xnに対するキュムラントも、そのキュムラント母関数の級数展開の係数として、次式で与えられる。
- [math] \begin{align} \log{M(s_1,\dots,s_n)} & = \sum_{k=1}^n \frac{1}{k!} \sum_{j_1,\dots, j_k=1}^n \langle X_{j_1} \dots X_{j_k} \rangle_c s_{j_1}\dots s_{j_k} \\ & =\sum_{k=1}^n \sum_{k_1+\dots k_n=k}\langle X_1^{\, k_1} \cdots X_n^{\, k_n} \rangle_c \frac{ s_1^{k_1} \dots s_n^{k_n}}{k_1! \dots k_n!} \end{align} [/math]
ここで、第一行の右辺は、j1,..., jkは1,..., nの値を自由に取るものとした、変数の重複も含む表現であり、第二行の右辺は変数の重複は含まない表現である。
多変数のキュムラントは、確率変数X1,..., Xnが二つ以上の互いに独立な組に分かれるとすると、それらのなかで独立な変数にまたがるキュムラントは常に0になるという重要な性質を有する。
- [math] \langle X_{j_1} \cdots X_{j_k} \rangle_c =0 [/math] (独立な組に分けられる場合)
実際、確率変数が例えば、(X1,..., Xm)、(Xm+1,..., Xn)と二つの独立な組に分かれるとすると、モーメント母関数はその性質から二つの積に分解される。
- [math] M(s_1,\dots,s_n) = M(s_1,\dots,s_m) M(s_{m+1},\dots,s_n) [/math]
従って、キュムラント母関数は
- [math] \log{M(s_1,\dots,s_n)} =\log{ M(s_1,\dots,s_m)} + \log{M(s_{m+1},\dots,s_n)} [/math]
となり、左辺の級数展開において、二つの組にまたがる変数の積の項は現れず、対応するキュムラントは0となる。二つ以上の互いに独立な組に分かれる場合も同様である。
例
ポアソン分布
確率分布関数が
- [math] P(x=k) = \frac{\lambda^k e^{-\lambda}}{k!} \, [/math]
で与えられるポアソン分布において、キュムラント母関数は次のように与えられる。
- [math] \log{M(s)} = \lambda(e^s-1)=\sum_{n=1}^{\infty}\frac{\lambda s^n}{n!} [/math]
従って、全てのキュムラント cn は λ となる。
- [math] c_n = \lambda \quad n=1,2,\dots [/math]
ガウス分布
確率密度関数が
- [math] p(x) = \frac{1}{\sqrt{2\pi} \sigma} \exp \left(-\frac{(x-\mu)^2}{2\sigma^2} \right) [/math]
で与えられるガウス分布において、キュムラント母関数は
- [math] \log{M(s)} = \mu s + \frac{\sigma^2}{2!}s^2 [/math]
であり、キュムラントは
- [math] \begin{align} c_1 & = \mu \\ c_2 &= \sigma^2 \\ c_n &= 0 \quad n=3,4,\dots \end{align} [/math]
となる。このように三次以上の高次のキュムラントが全て0になるのが、ガウス分布の特徴である。
物理学への応用
物理学でキュムラント展開が用いられる例として以下のようなものがある。
脚注
- ↑ 伏見康治「確率論及統計論」第III章 記述的統計学 18節 Thieleの半不変数 p.127 ISBN 9784874720127 http://ebsa.ism.ac.jp/ebooks/ebook/204
- ↑ T. N. Thiéle, Almindelig Iagttagelseslaere: Sandsynlighedsregning og mindste Kvadraters Methode, C. A. Reitzel, Copenhagen, 1889.
参考文献