確率論における量の集合
確率論 と 統計学 において 、 確率分布 の キュムラント κ nは、分布の モーメント の代わりとなる一連の量です 。モーメントが同一の 2 つの確率分布は、キュムラントも同一であり、その逆も同様です。
最初のキュムラントは平均 、 2 番目のキュムラントは 分散 、3 番目のキュムラントは 3 番目の 中心モーメント と同じです。しかし、4 次以上のキュムラントは中心モーメントと等しくありません。場合によっては、キュムラントによる問題の理論的処理はモーメントを使用するよりも簡単です。特に、2 つ以上のランダム変数が 統計的に独立している 場合、それらの合計の n次キュムラントは、それらの n 次キュムラントの合計に等しくなります 。また、 正規分布 の 3 次以上のキュムラントはゼロであり、この特性を持つ分布は正規分布だけです。
モーメントの場合と同様に、ランダム変数の集合に ジョイントモーメントが使用される場合、 ジョイントキュムラントを 定義することができます 。
意味
ランダム変数X のキュムラントは、 モーメント生成関数 の 自然対数 である キュムラント生成関数 K ( t ) を使用して定義されます 。
け
(
t
)
=
ログ
え
[
e
t
バツ
]
。
{\displaystyle K(t)=\log \operatorname {E} \left[e^{tX}\right].}
キュムラント κn は キュムラント生成関数の
べき級数 展開から得られる。
け
(
t
)
=
∑
ん
=
1
∞
κ
ん
t
ん
ん
!
=
κ
1
t
1
!
+
κ
2
t
2
2
!
+
κ
3
t
3
3
!
+
⋯
=
μ
t
+
σ
2
t
2
2
+
⋯
。
{\displaystyle K(t)=\sum _{n=1}^{\infty }\kappa _{n}{\frac {t^{n}}{n!}}=\kappa _{1}{\frac {t}{1!}}+\kappa _{2}{\frac {t^{2}}{2!}}+\kappa _{3}{\frac {t^{3}}{3!}}+\cdots =\mu t+\sigma ^{2}{\frac {t^{2}}{2}}+\cdots .}
この展開は マクローリン級数 なので、 n番目のキュムラントは上記の展開を n 回微分し、その結果をゼロで評価すること で得られる。 [1]
κ
ん
=
け
(
ん
)
(
0
)
。
{\displaystyle \kappa _{n}=K^{(n)}(0).}
モーメント生成関数が存在しない場合は、キュムラントは、後述するキュムラントとモーメントの関係に基づいて定義できます。
キュムラント生成関数の別の定義
一部の著者 [2] [3] は、キュムラント生成関数を 特性関数 の自然対数として定義することを好む。これは 第2 特性関数 とも呼ばれる。 [4] [5]
H
(
t
)
=
ログ
え
[
e
私
t
バツ
]
=
∑
ん
=
1
∞
κ
ん
(
私
t
)
ん
ん
!
=
μ
私
t
−
σ
2
t
2
2
+
⋯
{\displaystyle H(t)=\log \operatorname {E} \left[e^{itX}\right]=\sum _{n=1}^{\infty }\kappa _{n}{\frac {(it)^{n}}{n!}}=\mu it-\sigma ^{2}{\frac {t^{2}}{2}}+\cdots }
H ( t ) (ある意味で は 関数 K ( t ) を 純虚引数について評価したもの) の利点は、 E[ e itX ] が t の すべて の 実 数値 に対して 明確 に定義されること です。これは、 X が大きな 値 を持つ可能性が「高すぎる」場合に発生することがあります 。関数 H ( t ) は 明確に定義されますが、 マクローリン級数 の長さの点では K ( t )を模倣します。マクローリン級数は引数 t の線形順序を超えて拡張することはできず (まれに線形順序に達することさえありません) 、特に明確に定義されているキュムラントの数は変わりません。ただし、 H ( t ) に長いマクローリン級数がなくても、ランダム変数の解析、特に追加に直接使用できます。コーシー 分布 (ローレンツ分布とも呼ばれる)と、より一般的には 安定分布 (レヴィ分布に関連)は、生成関数の冪級数展開に有限個の明確に定義された項のみが含まれる分布の例です。
いくつかの基本的な特性
ランダム変数(の分布)の 番目 のキュムラントには、 次の特性があります。
ん
{\textstyle n}
κ
ん
(
バツ
)
{\textstyle \kappa _{n}(X)}
バツ
{\textstyle X}
とが 定数(ランダムではない) ならば、 すなわちキュムラントは 並進不変で ある。( ならば、
ん
>
1
{\textstyle n>1}
c
{\textstyle c}
κ
ん
(
バツ
+
c
)
=
κ
ん
(
バツ
)
、
{\textstyle \kappa _{n}(X+c)=\kappa _{n}(X),}
ん
=
1
{\textstyle n=1}
κ
1
(
バツ
+
c
)
=
κ
1
(
バツ
)
+
c
。
)
{\textstyle \kappa _{1}(X+c)=\kappa _{1}(X)+c.)}
が定数(ランダムではない)である 場合、つまり 番目の キュムラントは 次数 の 同次 です 。
c
{\textstyle c}
κ
ん
(
c
バツ
)
=
c
ん
κ
ん
(
バツ
)
、
{\textstyle \kappa _{n}(cX)=c^{n}\kappa _{n}(X),}
ん
{\textstyle n}
ん
{\textstyle n}
ランダム変数 が独立している場合 、つまり、キュムラントは累積的であるため、この名前が付けられます。
バツ
1
、
…
、
バツ
メートル
{\textstyle X_{1},\ldots ,X_{m}}
κ
ん
(
バツ
1
+
⋯
+
バツ
メートル
)
=
κ
ん
(
バツ
1
)
+
⋯
+
κ
ん
(
バツ
メートル
)
。
{\displaystyle \kappa _{n}(X_{1}+\cdots +X_{m})=\kappa _{n}(X_{1})+\cdots +\kappa _{n}(X_{m})\,.}
累積特性は、キュムラント生成関数を考慮することですぐにわかります。
つまり、独立したランダム変数の和の各キュムラントは、 加数 の対応するキュムラントの和です。つまり、加数が統計的に独立している場合、和の平均は平均の合計、和の分散は分散の合計、和の 3 番目のキュムラント (これは 3 番目の中心モーメントになります) は 3 番目のキュムラントの合計であり、キュムラントの各次数について同様に続きます。
け
バツ
1
+
⋯
+
バツ
メートル
(
t
)
=
ログ
え
[
e
t
(
バツ
1
+
⋯
+
バツ
メートル
)
]
=
ログ
(
え
[
e
t
バツ
1
]
⋯
え
[
e
t
バツ
メートル
]
)
=
ログ
え
[
e
t
バツ
1
]
+
⋯
+
ログ
え
[
e
t
バツ
メートル
]
=
け
バツ
1
(
t
)
+
⋯
+
け
バツ
メートル
(
t
)
、
{\displaystyle {\begin{aligned}K_{X_{1}+\cdots +X_{m}}(t)&=\log \operatorname {E} \left[e^{t(X_{1}+\cdots +X_{m})}\right]\\[5pt]&=\log \left(\operatorname {E} \left[e^{tX_{1}}\right]\cdots \operatorname {E} \left[e^{tX_{m}}\right]\right)\\[5pt]&=\log \operatorname {E} \left[e^{tX_{1}}\right]+\cdots +\log \operatorname {E} \left[e^{tX_{m}}\right]\\[5pt]&=K_{X_{1}}(t)+\cdots +K_{X_{m}}(t)、\end{aligned}}}
与えられたキュムラントκ n を持つ分布は、 エッジワース級数 を通じて近似できます 。
モーメントの関数としての最初のいくつかのキュムラント
高次のキュムラントはすべて中心モーメントの多項式関数で、係数は整数ですが、2 次と 3 次だけが実際に中心モーメントとなるキュムラントです。
κ
1
(
バツ
)
=
え
(
バツ
)
=
{\textstyle \kappa _{1}(X)=\オペレータ名 {E} (X)={}}
平均
κ
2
(
バツ
)
=
var
(
バツ
)
=
え
(
(
バツ
−
え
(
バツ
)
)
2
)
=
{\textstyle \kappa _{2}(X)=\オペレーター名 {var} (X)=\オペレーター名 {E} {\big (}(X-\オペレーター名 {E} (X))^{2}{\大)}={}}
分散、または第 2 中心モーメント。
κ
3
(
バツ
)
=
え
(
(
バツ
−
え
(
バツ
)
)
3
)
=
{\textstyle \kappa _{3}(X)=\オペレーター名 {E} {\big (}(X-\オペレーター名 {E} (X))^{3}{\big )}={}}
3番目の中心的な瞬間。
κ
4
(
バツ
)
=
え
(
(
バツ
−
え
(
バツ
)
)
4
)
−
3
(
え
(
(
バツ
−
え
(
バツ
)
)
2
)
)
2
=
{\textstyle \kappa _{4}(X)=\operatorname {E} {\big (}(X-\operatorname {E} (X))^{4}{\big )}-3\left(\operatorname {E} {\big (}(X-\operatorname {E} (X))^{2}{\big )}\right)^{2}={}}
4 番目の中心モーメントから 2 番目の中心モーメントの 2 乗の 3 倍を引いた値です。したがって、これはキュムラントが単なるモーメントや中心モーメントではない最初のケースです。3 次を超える中心モーメントには累積特性がありません。
κ
5
(
X
)
=
E
(
(
X
−
E
(
X
)
)
5
)
−
10
E
(
(
X
−
E
(
X
)
)
3
)
E
(
(
X
−
E
(
X
)
)
2
)
.
{\textstyle \kappa _{5}(X)=\operatorname {E} {\big (}(X-\operatorname {E} (X))^{5}{\big )}-10\operatorname {E} {\big (}(X-\operatorname {E} (X))^{3}{\big )}\operatorname {E} {\big (}(X-\operatorname {E} (X))^{2}{\big )}.}
いくつかの離散確率分布のキュムラント
定数確率変数 X = μ 。キュムラント生成関数は K ( t ) = μt です。最初のキュムラントは κ 1 = K '(0) = μ で、他のキュムラントはゼロ、 κ 2 = κ 3 = κ 4 = ⋅⋅⋅ = 0 です。
ベルヌーイ 分布(成功確率 p での 1 回の試行における成功回数 )。キュムラント生成関数は K ( t ) = log(1 − p + p e t ) です。最初のキュムラントは κ 1 = K '(0) = p および κ 2 = K′′ (0) = p ·(1 − p ) です。キュムラントは再帰式を満たします。
κ
n
+
1
=
p
(
1
−
p
)
d
κ
n
d
p
.
{\displaystyle \kappa _{n+1}=p(1-p){\frac {d\kappa _{n}}{dp}}.}
幾何 分布、(各試行で成功確率 p で1 回の成功までの失敗回数 )。キュムラント生成関数は K ( t ) = log( p / (1 + ( p − 1)et ) ) です。最初のキュムラントは κ 1 = K′ (0) = p −1 − 1 および κ 2 = K′′ (0) = κ 1 p −1 です。p = ( μ + 1) −1を代入すると、 K ( t ) = −log(1 + μ (1−et ) ) および κ 1 = μ と なります 。
ポアソン 分布 。キュムラント生成関数は K ( t ) = μ (e t − 1) です。すべてのキュムラントはパラメータ κ 1 = κ 2 = κ 3 = ... = μ に等しい。
二項分布 、 ( 各試行の成功確率が pである n 回 の独立した 試行 における成功回数)。特殊なケースである n = 1 はベルヌーイ分布です。すべてのキュムラントは、 対応するベルヌーイ分布の対応するキュムラントの n倍です。キュムラント生成関数は K ( t ) = n log(1 − p + p e t ) です。最初のキュムラントは κ 1 = K′ (0) = np および κ 2 = K′′ (0) = κ 1 (1 − p ) です。p = μ· n −1を代入すると、 K '( t ) = ((μ −1 − n −1 )·e − t + n −1 ) −1 および κ 1 = μ が得られます 。極限ケース n −1 = 0 はポアソン分布です 。
負 の二項分布 (各試行で成功確率 p で r 回 の成功を得るまでの失敗回数)。特殊なケースである r = 1 は幾何分布です。すべてのキュムラントは、対応する幾何分布の対応するキュムラントの r 倍です。キュムラント生成関数の導関数は、 K ′( t ) = r ·((1 − p ) −1 ·e − t −1) −1 です。最初のキュムラントは、 κ 1 = K ′(0) = r ·( p −1 −1) および κ 2 = K ''(0) = κ 1 · p −1 です。 p = (μ· r −1 +1) −1を代入すると、 K ′( t ) = (( μ −1 + r −1 ) e − t − r −1 ) −1 および κ 1 = μ となります 。これらの式を二項分布の式と比較すると、「負の二項分布」という名前が説明されます。 極限ケース r −1 = 0 はポアソン分布です。
分散平均比
を導入すると
、上記の確率分布はキュムラント生成関数の導関数の統一された式を得る: [ 要出典 ]
ε
=
μ
−
1
σ
2
=
κ
1
−
1
κ
2
,
{\displaystyle \varepsilon =\mu ^{-1}\sigma ^{2}=\kappa _{1}^{-1}\kappa _{2},}
K
′
(
t
)
=
(
1
+
(
e
−
t
−
1
)
ε
)
−
1
μ
{\displaystyle K'(t)=(1+(e^{-t}-1)\varepsilon )^{-1}\mu }
2番目の導関数は、
最初のキュムラントが κ 1 = K′ (0) = μ であり、2番目のキュムラントが κ 2 = K′ (0) = με で あることを確認します。
K
″
(
t
)
=
(
ε
−
(
ε
−
1
)
e
t
)
−
2
μ
ε
e
t
{\displaystyle K''(t)=(\varepsilon -(\varepsilon -1)e^{t})^{-2}\mu \varepsilon e^{t}}
定数ランダム変数 X = μでは ε = 0 となります 。
二項分布では ε = 1 − p となるので、 0 < ε < 1 となります。
ポアソン分布では ε = 1 です。
負の二項分布では ε = p −1 となるので ε > 1 となる。
離心率 による 円錐断面 の分類との類似点に注意してください : 円 ε = 0 、楕円 0 < ε < 1 、放物線 ε = 1 、双曲線 ε > 1 。
いくつかの連続確率分布のキュムラント
期待値 μ 、 分散 σ 2 の 正規分布 の場合 、キュムラント生成関数は K ( t ) = μt + σ 2 t 2 /2 です。キュムラント生成関数の 1 次導関数と 2 次導関数は、 K ′( t ) = μ + σ 2 · t と K ''( t ) = σ 2 です。キュムラントは、 κ 1 = μ 、 κ 2 = σ 2 、 κ 3 = κ 4 = ⋅⋅⋅ = 0 です。特別な場合 σ 2 = 0 は、定数ランダム変数 X = μ です。
区間 [−1, 0] 上 の一様分布 の キュムラントは κn = Bn / n であり 、ここで Bn は n 番目の ベルヌーイ数 である 。
速度パラメータ λ を持つ指数分布 のキュムラントは κn = λ − n ( n − 1)! です 。
キュムラント生成関数のいくつかの性質
キュムラント生成関数 K ( t ) は 、存在する場合、 無限に微分可能 かつ 凸で 、原点を通ります。その第 1 導関数は 、確率分布のサポートの 最小値から 最大値 までの開区間で単調に変化し、第 2 導関数は 、単一の質点の 退化した分布を除いて、定義されているすべての場所で厳密に正です。キュムラント生成関数が存在するのは、分布の裾が 指数関数的減少 によって主要化される場合、つまり、次の式で表される場合のみです ( Big O 表記を 参照 )
。
ここで、 は 累積分布関数 です 。キュムラント生成関数は、そのような c の 負の上限が存在する場合はその上限に 、 そのような d の上限が存在する場合はその上限に 垂直漸近 線 を持ちます。それ以外の場合は、すべての実数に対して定義されます。
∃
c
>
0
,
F
(
x
)
=
O
(
e
c
x
)
,
x
→
−
∞
;
and
∃
d
>
0
,
1
−
F
(
x
)
=
O
(
e
−
d
x
)
,
x
→
+
∞
;
{\displaystyle {\begin{aligned}&\exists c>0,\,\,F(x)=O(e^{cx}),x\to -\infty ;{\text{ and}}\\[4pt]&\exists d>0,\,\,1-F(x)=O(e^{-dx}),x\to +\infty ;\end{aligned}}}
F
{\textstyle F}
ランダム変数 X のサポート に有限の上限または下限がある 場合、そのキュムラント生成関数 y = K ( t ) が存在する場合は、その関数は、傾きがサポートの上限または下限に等しく、どこでもこれら 2 つの線より上にある漸近線に近づきます。(K ( 0 ) = 0
で ある ため 、 積分 によりこれらの漸近線の y 切片
が得られます 。)
y
=
(
t
+
1
)
inf
supp
X
−
μ
(
X
)
,
and
y
=
(
t
−
1
)
sup
supp
X
+
μ
(
X
)
,
{\displaystyle {\begin{aligned}y&=(t+1)\inf \operatorname {supp} X-\mu (X),{\text{ and}}\\[5pt]y&=(t-1)\sup \operatorname {supp} X+\mu (X),\end{aligned}}}
∫
−
∞
0
[
t
inf
supp
X
−
K
′
(
t
)
]
d
t
,
∫
∞
0
[
t
inf
supp
X
−
K
′
(
t
)
]
d
t
{\displaystyle \int _{-\infty }^{0}\left[t\inf \operatorname {supp} X-K'(t)\right]\,dt,\qquad \int _{\infty }^{0}\left[t\inf \operatorname {supp} X-K'(t)\right]\,dt}
分布が c だけシフトした場合、 c における退化した質点に対して 、キュムラント生成関数は直線であり 、より一般的には、 X と Y が独立であり、それらのキュムラント生成関数が存在する場合のみで ある( 劣独立性 と2次モーメントの存在は独立性を示唆するのに十分である。 [6] )
K
X
+
c
(
t
)
=
K
X
(
t
)
+
c
t
.
{\textstyle K_{X+c}(t)=K_{X}(t)+ct.}
K
c
(
t
)
=
c
t
{\textstyle K_{c}(t)=ct}
K
X
+
Y
=
K
X
+
K
Y
{\textstyle K_{X+Y}=K_{X}+K_{Y}}
分布の自然指数族は、K ( t ) をシフトまたは平行移動し、常に原点を通るように垂直に調整することによって実現できます 。f が キュムラント 生成 関数 を 持つ pdf で 、が その自然指数族である場合 、
K
(
t
)
=
log
M
(
t
)
,
{\textstyle K(t)=\log M(t),}
f
|
θ
{\textstyle f|\theta }
f
(
x
∣
θ
)
=
1
M
(
θ
)
e
θ
x
f
(
x
)
,
{\textstyle f(x\mid \theta )={\frac {1}{M(\theta )}}e^{\theta x}f(x),}
K
(
t
∣
θ
)
=
K
(
t
+
θ
)
−
K
(
θ
)
.
{\textstyle K(t\mid \theta )=K(t+\theta )-K(\theta ).}
K ( t )が t 1 < Re( t ) < t 2 の 範囲で有限である 場合、 t 1 < 0 < t 2 であれば 、 K ( t ) は 解析的であり、 t 1 < Re( t ) < t 2 で無限に微分可能です。さらに、 t が実数で t 1 < t < t 2の場合、 K ( t ) は厳密に凸であり、 K ′( t ) は厳密に増加します。 [ 引用が必要 ]
キュムラントのさらなる性質
否定的な結果
正規分布 のキュムラントの結果を考えると 、
いくつかの m > 3に対して κ m = κ m +1 = ⋯ = 0 となり、低次のキュムラント(3から m − 1 の次数)がゼロでない分布の族が見つかると期待されるかもしれない。そのような分布は存在しない。 [7] ここでの根本的な結果は、キュムラント生成関数は2より大きい次数の有限次多項式にはなり得ないということである。
キュムラントとモーメント
モーメント 生成関数は 次のように表されます。
M
(
t
)
=
1
+
∑
n
=
1
∞
μ
n
′
t
n
n
!
=
exp
(
∑
n
=
1
∞
κ
n
t
n
n
!
)
=
exp
(
K
(
t
)
)
.
{\displaystyle M(t)=1+\sum _{n=1}^{\infty }{\frac {\mu '_{n}t^{n}}{n!}}=\exp \left(\sum _{n=1}^{\infty }{\frac {\kappa _{n}t^{n}}{n!}}\right)=\exp(K(t)).}
つまり、キュムラント母関数はモーメント母関数の対数である。
K
(
t
)
=
log
M
(
t
)
.
{\displaystyle K(t)=\log M(t).}
最初のキュムラントは 期待値 です。2 番目と 3 番目のキュムラントはそれぞれ 2 番目と 3 番目の 中心モーメント です(2 番目中心モーメントは 分散 です)。ただし、高次のキュムラントはモーメントでも中心モーメントでもなく、モーメントのより複雑な多項式関数です。
モーメントは、 におけるの n 次導関数を評価することでキュムラントの観点から復元できる 。
exp
(
K
(
t
)
)
{\textstyle \exp(K(t))}
t
=
0
{\displaystyle t=0}
μ
n
′
=
M
(
n
)
(
0
)
=
d
n
exp
(
K
(
t
)
)
d
t
n
|
t
=
0
.
{\displaystyle \mu '_{n}=M^{(n)}(0)=\left.{\frac {\mathrm {d} ^{n}\exp(K(t))}{\mathrm {d} t^{n}}}\right|_{t=0}.}
同様に、キュムラントは、 におけるの n 次導関数を評価することによってモーメントの観点から回復することができる 。
log
M
(
t
)
{\textstyle \log M(t)}
t
=
0
{\displaystyle t=0}
κ
n
=
K
(
n
)
(
0
)
=
d
n
log
M
(
t
)
d
t
n
|
t
=
0
.
{\displaystyle \kappa _{n}=K^{(n)}(0)=\left.{\frac {\mathrm {d} ^{n}\log M(t)}{\mathrm {d} t^{n}}}\right|_{t=0}.}
最初のn キュムラントによる n 次モーメントの明示的な表現、およびその逆は、合成関数の高次導関数に対する ファディブルーノの公式 を使用することで得られます 。一般に、 となります。
ここで、 は不完全 (または部分) ベル多項式 です。
μ
n
′
=
∑
k
=
1
n
B
n
,
k
(
κ
1
,
…
,
κ
n
−
k
+
1
)
{\displaystyle \mu '_{n}=\sum _{k=1}^{n}B_{n,k}(\kappa _{1},\ldots ,\kappa _{n-k+1})}
κ
n
=
∑
k
=
1
n
(
−
1
)
k
−
1
(
k
−
1
)
!
B
n
,
k
(
μ
1
′
,
…
,
μ
n
−
k
+
1
′
)
,
{\displaystyle \kappa _{n}=\sum _{k=1}^{n}(-1)^{k-1}(k-1)!B_{n,k}(\mu '_{1},\ldots ,\mu '_{n-k+1}),}
B
n
,
k
{\textstyle B_{n,k}}
同様に、平均が で与えられる場合 、中心モーメント生成関数は で与えられ
、 n 番目の中心モーメントはキュムラントで次のように得られる。
μ
{\textstyle \mu }
C
(
t
)
=
E
[
e
t
(
x
−
μ
)
]
=
e
−
μ
t
M
(
t
)
=
exp
(
K
(
t
)
−
μ
t
)
,
{\displaystyle C(t)=\operatorname {E} [e^{t(x-\mu )}]=e^{-\mu t}M(t)=\exp(K(t)-\mu t),}
μ
n
=
C
(
n
)
(
0
)
=
d
n
d
t
n
exp
(
K
(
t
)
−
μ
t
)
|
t
=
0
=
∑
k
=
1
n
B
n
,
k
(
0
,
κ
2
,
…
,
κ
n
−
k
+
1
)
.
{\displaystyle \mu _{n}=C^{(n)}(0)=\left.{\frac {\mathrm {d} ^{n}}{\mathrm {d} t^{n}}}\exp(K(t)-\mu t)\right|_{t=0}=\sum _{k=1}^{n}B_{n,k}(0,\kappa _{2},\ldots ,\kappa _{n-k+1}).}
また、 n > 1 の場合、 中心モーメントによる
n番目のキュムラントは
κ
n
=
K
(
n
)
(
0
)
=
d
n
d
t
n
(
log
C
(
t
)
+
μ
t
)
|
t
=
0
=
∑
k
=
1
n
(
−
1
)
k
−
1
(
k
−
1
)
!
B
n
,
k
(
0
,
μ
2
,
…
,
μ
n
−
k
+
1
)
.
{\displaystyle {\begin{aligned}\kappa _{n}&=K^{(n)}(0)=\left.{\frac {\mathrm {d} ^{n}}{\mathrm {d} t^{n}}}(\log C(t)+\mu t)\right|_{t=0}\\[4pt]&=\sum _{k=1}^{n}(-1)^{k-1}(k-1)!B_{n,k}(0,\mu _{2},\ldots ,\mu _{n-k+1}).\end{aligned}}}
n 次 の モーメント μ ′ n は 、最初のn 個のキュムラントの n 次多項式 です 。最初のいくつかの式は次のとおりです。
μ
1
′
=
κ
1
μ
2
′
=
κ
2
+
κ
1
2
μ
3
′
=
κ
3
+
3
κ
2
κ
1
+
κ
1
3
μ
4
′
=
κ
4
+
4
κ
3
κ
1
+
3
κ
2
2
+
6
κ
2
κ
1
2
+
κ
1
4
μ
5
′
=
κ
5
+
5
κ
4
κ
1
+
10
κ
3
κ
2
+
10
κ
3
κ
1
2
+
15
κ
2
2
κ
1
+
10
κ
2
κ
1
3
+
κ
1
5
μ
6
′
=
κ
6
+
6
κ
5
κ
1
+
15
κ
4
κ
2
+
15
κ
4
κ
1
2
+
10
κ
3
2
+
60
κ
3
κ
2
κ
1
+
20
κ
3
κ
1
3
+
15
κ
2
3
+
45
κ
2
2
κ
1
2
+
15
κ
2
κ
1
4
+
κ
1
6
.
{\displaystyle {\begin{aligned}\mu '_{1}={}&\kappa _{1}\\[5pt]\mu '_{2}={}&\kappa _{2}+\kappa _{1}^{2}\\[5pt]\mu '_{3}={}&\kappa _{3}+3\kappa _{2}\kappa _{1}+\kappa _{1}^{3}\\[5pt]\mu '_{4}={}&\kappa _{4}+4\kappa _{3}\kappa _{1}+3\kappa _{2}^{2}+6\kappa _{2}\kappa _{1}^{2}+\kappa _{1}^{4}\\[5pt]\mu '_{5}={}&\kappa _{5}+5\kappa _{4}\kappa _{1}+10\kappa _{3}\kappa _{2}+10\kappa _{3}\kappa _{1}^{2}+15\kappa _{2}^{2}\kappa _{1}+10\kappa _{2}\kappa _{1}^{3}+\kappa _{1}^{5}\\[5pt]\mu '_{6}={}&\kappa _{6}+6\kappa _{5}\kappa _{1}+15\kappa _{4}\kappa _{2}+15\kappa _{4}\kappa _{1}^{2}+10\kappa _{3}^{2}+60\kappa _{3}\kappa _{2}\kappa _{1}+20\kappa _{3}\kappa _{1}^{3}\\&{}+15\kappa _{2}^{3}+45\kappa _{2}^{2}\kappa _{1}^{2}+15\kappa _{2}\kappa _{1}^{4}+\kappa _{1}^{6}.\end{aligned}}}
「プライム」はモーメント μ ′ n と 中心モーメント μ n を区別します。 中心 モーメントをキュムラントの関数として表すには、これらの多項式から κ 1 が 因子として現れるすべての項を削除するだけです。
μ
1
=
0
μ
2
=
κ
2
μ
3
=
κ
3
μ
4
=
κ
4
+
3
κ
2
2
μ
5
=
κ
5
+
10
κ
3
κ
2
μ
6
=
κ
6
+
15
κ
4
κ
2
+
10
κ
3
2
+
15
κ
2
3
.
{\displaystyle {\begin{aligned}\mu _{1}&=0\\[4pt]\mu _{2}&=\kappa _{2}\\[4pt]\mu _{3}&=\kappa _{3}\\[4pt]\mu _{4}&=\kappa _{4}+3\kappa _{2}^{2}\\[4pt]\mu _{5}&=\kappa _{5}+10\kappa _{3}\kappa _{2}\\[4pt]\mu _{6}&=\kappa _{6}+15\kappa _{4}\kappa _{2}+10\kappa _{3}^{2}+15\kappa _{2}^{3}.\end{aligned}}}
同様に、 n 番目のキュムラント κ n は 、最初のn 個の非中心モーメント
における n 次多項式 です。最初のいくつかの式は次のとおりです。
κ
1
=
μ
1
′
κ
2
=
μ
2
′
−
μ
1
′
2
κ
3
=
μ
3
′
−
3
μ
2
′
μ
1
′
+
2
μ
1
′
3
κ
4
=
μ
4
′
−
4
μ
3
′
μ
1
′
−
3
μ
2
′
2
+
12
μ
2
′
μ
1
′
2
−
6
μ
1
′
4
κ
5
=
μ
5
′
−
5
μ
4
′
μ
1
′
−
10
μ
3
′
μ
2
′
+
20
μ
3
′
μ
1
′
2
+
30
μ
2
′
2
μ
1
′
−
60
μ
2
′
μ
1
′
3
+
24
μ
1
′
5
κ
6
=
μ
6
′
−
6
μ
5
′
μ
1
′
−
15
μ
4
′
μ
2
′
+
30
μ
4
′
μ
1
′
2
−
10
μ
3
′
2
+
120
μ
3
′
μ
2
′
μ
1
′
−
120
μ
3
′
μ
1
′
3
+
30
μ
2
′
3
−
270
μ
2
′
2
μ
1
′
2
+
360
μ
2
′
μ
1
′
4
−
120
μ
1
′
6
.
{\displaystyle {\begin{aligned}\kappa _{1}={}&\mu '_{1}\\[4pt]\kappa _{2}={}&\mu '_{2}-{\mu '_{1}}^{2}\\[4pt]\kappa _{3}={}&\mu '_{3}-3\mu '_{2}\mu '_{1}+2{\mu '_{1}}^{3}\\[4pt]\kappa _{4}={}&\mu '_{4}-4\mu '_{3}\mu '_{1}-3{\mu '_{2}}^{2}+12\mu '_{2}{\mu '_{1}}^{2}-6{\mu '_{1}}^{4}\\[4pt]\kappa _{5}={}&\mu '_{5}-5\mu '_{4}\mu '_{1}-10\mu '_{3}\mu '_{2}+20\mu '_{3}{\mu '_{1}}^{2}+30{\mu '_{2}}^{2}\mu '_{1}-60\mu '_{2}{\mu '_{1}}^{3}+24{\mu '_{1}}^{5}\\[4pt]\kappa _{6}={}&\mu '_{6}-6\mu '_{5}\mu '_{1}-15\mu '_{4}\mu '_{2}+30\mu '_{4}{\mu '_{1}}^{2}-10{\mu '_{3}}^{2}+120\mu '_{3}\mu '_{2}\mu '_{1}\\&{}-120\mu '_{3}{\mu '_{1}}^{3}+30{\mu '_{2}}^{3}-270{\mu '_{2}}^{2}{\mu '_{1}}^{2}+360\mu '_{2}{\mu '_{1}}^{4}-120{\mu '_{1}}^{6}\,.\end{aligned}}}
一般に、 [8] キュムラントは行列の行列式である:
κ
l
=
(
−
1
)
l
+
1
|
μ
1
′
1
0
0
0
0
…
0
μ
2
′
μ
1
′
1
0
0
0
…
0
μ
3
′
μ
2
′
(
2
1
)
μ
1
′
1
0
0
…
0
μ
4
′
μ
3
′
(
3
1
)
μ
2
′
(
3
2
)
μ
1
′
1
0
…
0
μ
5
′
μ
4
′
(
4
1
)
μ
3
′
(
4
2
)
μ
2
′
(
4
3
)
μ
1
′
1
…
0
⋮
⋮
⋮
⋮
⋮
⋱
⋱
⋮
μ
l
−
1
′
μ
l
−
2
′
…
…
…
…
⋱
1
μ
l
′
μ
l
−
1
′
…
…
…
…
…
(
l
−
1
l
−
2
)
μ
1
′
|
{\displaystyle \kappa _{l}=(-1)^{l+1}\left|{\begin{array}{cccccccc}\mu '_{1}&1&0&0&0&0&\ldots &0\\\mu '_{2}&\mu '_{1}&1&0&0&0&\ldots &0\\\mu '_{3}&\mu '_{2}&\left({\begin{array}{l}2\\1\end{array}}\right)\mu '_{1}&1&0&0&\ldots &0\\\mu '_{4}&\mu '_{3}&\left({\begin{array}{l}3\\1\end{array}}\right)\mu '_{2}&\left({\begin{array}{l}3\\2\end{array}}\right)\mu '_{1}&1&0&\ldots &0\\\mu '_{5}&\mu '_{4}&\left({\begin{array}{l}4\\1\end{array}}\right)\mu '_{3}&\left({\begin{array}{l}4\\2\end{array}}\right)\mu '_{2}&\left({\begin{array}{c}4\\3\end{array}}\right)\mu '_{1}&1&\ldots &0\\\vdots &\vdots &\vdots &\vdots &\vdots &\ddots &\ddots &\vdots \\\mu '_{l-1}&\mu '_{l-2}&\ldots &\ldots &\ldots &\ldots &\ddots &1\\\mu '_{l}&\mu '_{l-1}&\ldots &\ldots &\ldots &\ldots &\ldots &\left({\begin{array}{l}l-1\\l-2\end{array}}\right)\mu '_{1}\end{array}}\right|}
n > 1 の キュムラント κ n を 中心モーメントの関数として表すには、これらの多項式から μ' 1 が 因子として現れるすべての項を削除します。
κ
2
=
μ
2
{\displaystyle \kappa _{2}=\mu _{2}\,}
κ
3
=
μ
3
{\displaystyle \kappa _{3}=\mu _{3}\,}
κ
4
=
μ
4
−
3
μ
2
2
{\displaystyle \kappa _{4}=\mu _{4}-3{\mu _{2}}^{2}\,}
κ
5
=
μ
5
−
10
μ
3
μ
2
{\displaystyle \kappa _{5}=\mu _{5}-10\mu _{3}\mu _{2}\,}
κ
6
=
μ
6
−
15
μ
4
μ
2
−
10
μ
3
2
+
30
μ
2
3
.
{\displaystyle \kappa _{6}=\mu _{6}-15\mu _{4}\mu _{2}-10{\mu _{3}}^{2}+30{\mu _{2}}^{3}\,.}
n > 2 の キュムラント κ n を 標準化中心モーメント μ″ n の関数として表すには 、 多項式で
μ' 2 =1も設定します。
κ
3
=
μ
3
″
{\displaystyle \kappa _{3}=\mu ''_{3}\,}
κ
4
=
μ
4
″
−
3
{\displaystyle \kappa _{4}=\mu ''_{4}-3\,}
κ
5
=
μ
5
″
−
10
μ
3
″
{\displaystyle \kappa _{5}=\mu ''_{5}-10\mu ''_{3}\,}
κ
6
=
μ
6
″
−
15
μ
4
″
−
10
μ
3
″
2
+
30
.
{\displaystyle \kappa _{6}=\mu ''_{6}-15\mu ''_{4}-10{\mu ''_{3}}^{2}+30\,.}
キュムラントは、関係 log M ( t ) = K ( t )を t について 微分する ことでモーメントに関連付けることができ 、 M′ ( t ) = K′ ( t ) M ( t ) となります。この式には、便利なことに指数や対数が含まれていません。 t n −1 / ( n −1)! の係数を左辺と右辺で等しくし、 μ′ 0 = 1 を使用すると、 n ≥ 1 について次の式が得られます 。 [9]
これらの式により、低次のキュムラントとモーメントに関する知識を使用して、 またはのいずれかを他方から計算できます。 に対する 中心モーメントに対応する式は、それぞれを に対して に設定し て置き換えること で、これらの式から作成されます 。
μ
1
′
=
κ
1
μ
2
′
=
κ
1
μ
1
′
+
κ
2
μ
3
′
=
κ
1
μ
2
′
+
2
κ
2
μ
1
′
+
κ
3
μ
4
′
=
κ
1
μ
3
′
+
3
κ
2
μ
2
′
+
3
κ
3
μ
1
′
+
κ
4
μ
5
′
=
κ
1
μ
4
′
+
4
κ
2
μ
3
′
+
6
κ
3
μ
2
′
+
4
κ
4
μ
1
′
+
κ
5
μ
6
′
=
κ
1
μ
5
′
+
5
κ
2
μ
4
′
+
10
κ
3
μ
3
′
+
10
κ
4
μ
2
′
+
5
κ
5
μ
1
′
+
κ
6
μ
n
′
=
∑
m
=
1
n
−
1
(
n
−
1
m
−
1
)
κ
m
μ
n
−
m
′
+
κ
n
.
{\displaystyle {\begin{aligned}\mu '_{1}={}&\kappa _{1}\\[1pt]\mu '_{2}={}&\kappa _{1}\mu '_{1}+\kappa _{2}\\[1pt]\mu '_{3}={}&\kappa _{1}\mu '_{2}+2\kappa _{2}\mu '_{1}+\kappa _{3}\\[1pt]\mu '_{4}={}&\kappa _{1}\mu '_{3}+3\kappa _{2}\mu '_{2}+3\kappa _{3}\mu '_{1}+\kappa _{4}\\[1pt]\mu '_{5}={}&\kappa _{1}\mu '_{4}+4\kappa _{2}\mu '_{3}+6\kappa _{3}\mu '_{2}+4\kappa _{4}\mu '_{1}+\kappa _{5}\\[1pt]\mu '_{6}={}&\kappa _{1}\mu '_{5}+5\kappa _{2}\mu '_{4}+10\kappa _{3}\mu '_{3}+10\kappa _{4}\mu '_{2}+5\kappa _{5}\mu '_{1}+\kappa _{6}\\[1pt]\mu '_{n}={}&\sum _{m=1}^{n-1}{n-1 \choose m-1}\kappa _{m}\mu '_{n-m}+\kappa _{n}\,.\end{aligned}}}
κ
n
{\textstyle \kappa _{n}}
μ
n
′
{\textstyle \mu '_{n}}
μ
n
{\textstyle \mu _{n}}
n
≥
2
{\textstyle n\geq 2}
μ
1
′
=
κ
1
=
0
{\textstyle \mu '_{1}=\kappa _{1}=0}
μ
n
′
{\textstyle \mu '_{n}}
μ
n
{\textstyle \mu _{n}}
n
≥
2
{\textstyle n\geq 2}
μ
2
=
κ
2
μ
3
=
κ
3
μ
n
=
∑
m
=
2
n
−
2
(
n
−
1
m
−
1
)
κ
m
μ
n
−
m
+
κ
n
.
{\displaystyle {\begin{aligned}\mu _{2}={}&\kappa _{2}\\[1pt]\mu _{3}={}&\kappa _{3}\\[1pt]\mu _{n}={}&\sum _{m=2}^{n-2}{n-1 \choose m-1}\kappa _{m}\mu _{n-m}+\kappa _{n}\,.\end{aligned}}}
キュムラントと集合分割
これらの多項式は、 係数が 集合の特定の分割を数えるという、注目すべき 組み合わせ論的 解釈を持っています。これらの多項式の一般的な形は
、
μ
n
′
=
∑
π
∈
Π
∏
B
∈
π
κ
|
B
|
{\displaystyle \mu '_{n}=\sum _{\pi \,\in \,\Pi }\prod _{B\,\in \,\pi }\kappa _{|B|}}
π は、サイズ n のセットのすべてのパーティションのリストを実行します 。
「 B ∈ π 」は、 B が集合を分割する「ブロック」の1つであることを意味します。
| B | はセット B のサイズです 。
したがって、各 単項式は、指数の合計が n であるキュムラントの積の定数倍です (たとえば、項 κ 3 κ 2 2 κ 1 では、指数の合計は 3 + 2 + 2 + 1 = 8 です。これは、最初の 8 つのキュムラントの関数として 8 次モーメントを表す多項式に現れます)。 整数 n の分割が 各項に対応します。各項の 係数は、 n メンバーのセットの分割数であり 、セットのメンバーが区別できなくなるときに、その分割は整数 n の分割に縮小されます。
キュムラントと組み合わせ論
キュムラントと組合せ論のさらなる関連は、ジャン・カルロ・ロータ の研究に見出すことができます。ロータ では、 不変理論 、 対称関数 、二項式列との関連が 陰影計算 を介して研究されています。 [10]
結合キュムラント
いくつかのランダム変数 X 1 , ..., X n の結合 キュムラント κは 、多変量キュムラント生成関数のマクローリン級数の係数 κ 1,...,1 ( X 1 , ..., X n )として定義されます 。 [ 11 ] の セクション 3.1 を 参照 し て
ください。
特に、
単一変数の場合と同様に、生成関数とキュムラントは次のように定義できます
。
この場合
、および
G
(
t
1
,
…
,
t
n
)
=
log
E
(
e
∑
j
=
1
n
t
j
X
j
)
=
∑
k
1
,
…
,
k
n
κ
k
1
,
…
,
k
n
t
1
k
1
⋯
t
n
k
n
k
1
!
⋯
k
n
!
.
{\displaystyle G(t_{1},\dots ,t_{n})=\log \mathrm {E} (\mathrm {e} ^{\sum _{j=1}^{n}t_{j}X_{j}})=\sum _{k_{1},\ldots ,k_{n}}\kappa _{k_{1},\ldots ,k_{n}}{\frac {t_{1}^{k_{1}}\cdots t_{n}^{k_{n}}}{k_{1}!\cdots k_{n}!}}\,.}
κ
k
1
,
…
,
k
n
=
(
d
d
t
1
)
k
1
⋯
(
d
d
t
n
)
k
n
G
(
t
1
,
…
,
t
n
)
|
t
1
=
⋯
=
t
n
=
0
,
{\displaystyle \kappa _{k_{1},\dots ,k_{n}}=\left.\left({\frac {\mathrm {d} }{\mathrm {d} t_{1}}}\right)^{k_{1}}\cdots \left({\frac {\mathrm {d} }{\mathrm {d} t_{n}}}\right)^{k_{n}}G(t_{1},\dots ,t_{n})\right|_{t_{1}=\dots =t_{n}=0}\,,}
κ
(
X
1
,
…
,
X
n
)
=
d
n
d
t
1
⋯
d
t
n
G
(
t
1
,
…
,
t
n
)
|
t
1
=
⋯
=
t
n
=
0
.
{\displaystyle \kappa (X_{1},\ldots ,X_{n})=\left.{\frac {\mathrm {d} ^{n}}{\mathrm {d} t_{1}\cdots \mathrm {d} t_{n}}}G(t_{1},\dots ,t_{n})\right|_{t_{1}=\dots =t_{n}=0}\,.}
H
(
t
1
,
…
,
t
n
)
=
log
E
(
e
∑
j
=
1
n
i
t
j
X
j
)
=
∑
k
1
,
…
,
k
n
κ
k
1
,
…
,
k
n
i
k
1
+
⋯
+
k
n
t
1
k
1
⋯
t
n
k
n
k
1
!
⋯
k
n
!
,
{\displaystyle H(t_{1},\dots ,t_{n})=\log \mathrm {E} (\mathrm {e} ^{\sum _{j=1}^{n}it_{j}X_{j}})=\sum _{k_{1},\ldots ,k_{n}}\kappa _{k_{1},\ldots ,k_{n}}i^{k_{1}+\cdots +k_{n}}{\frac {t_{1}^{k_{1}}\cdots t_{n}^{k_{n}}}{k_{1}!\cdots k_{n}!}}\,,}
κ
k
1
,
…
,
k
n
=
(
−
i
)
k
1
+
⋯
+
k
n
(
d
d
t
1
)
k
1
⋯
(
d
d
t
n
)
k
n
H
(
t
1
,
…
,
t
n
)
|
t
1
=
⋯
=
t
n
=
0
,
{\displaystyle \kappa _{k_{1},\dots ,k_{n}}=(-i)^{k_{1}+\cdots +k_{n}}\left.\left({\frac {\mathrm {d} }{\mathrm {d} t_{1}}}\right)^{k_{1}}\cdots \left({\frac {\mathrm {d} }{\mathrm {d} t_{n}}}\right)^{k_{n}}H(t_{1},\dots ,t_{n})\right|_{t_{1}=\dots =t_{n}=0}\,,}
κ
(
X
1
,
…
,
X
n
)
=
(
−
i
)
n
d
n
d
t
1
⋯
d
t
n
H
(
t
1
,
…
,
t
n
)
|
t
1
=
⋯
=
t
n
=
0
.
{\displaystyle \kappa (X_{1},\ldots ,X_{n})=\left.(-i)^{n}{\frac {\mathrm {d} ^{n}}{\mathrm {d} t_{1}\cdots \mathrm {d} t_{n}}}H(t_{1},\dots ,t_{n})\right|_{t_{1}=\dots =t_{n}=0}\,.}
繰り返されるランダム変数と係数の関係 κ k 1 、...、 k n
はとも表記できることに
注目してください。このこと から、
たとえば
および であることが結論付けられます。
特に
、
最後の等式は、単一のランダム変数のキュムラントが、そのランダム変数の複数のコピーの結合キュムラントであることを示しています。
κ
k
1
,
…
,
k
n
(
X
1
,
…
,
X
n
)
{\textstyle \kappa _{k_{1},\dots ,k_{n}}(X_{1},\ldots ,X_{n})}
κ
k
1
,
…
,
k
n
=
d
k
1
d
t
1
,
1
⋯
d
t
1
,
k
1
⋯
d
k
n
d
t
n
,
1
⋯
d
t
n
,
k
n
G
(
∑
j
=
1
k
1
t
1
,
j
,
…
,
∑
j
=
1
k
n
t
n
,
j
)
|
t
i
,
j
=
0
,
{\displaystyle \kappa _{k_{1},\dots ,k_{n}}=\left.{\frac {\mathrm {d} ^{k_{1}}}{\mathrm {d} t_{1,1}\cdots \mathrm {d} t_{1,k_{1}}}}\cdots {\frac {\mathrm {d} ^{k_{n}}}{\mathrm {d} t_{n,1}\cdots \mathrm {d} t_{n,k_{n}}}}G\left(\sum _{j=1}^{k_{1}}t_{1,j},\dots ,\sum _{j=1}^{k_{n}}t_{n,j}\right)\right|_{t_{i,j}=0},}
κ
k
1
,
…
,
k
n
(
X
1
,
…
,
X
n
)
=
κ
1
,
…
,
1
(
X
1
,
…
,
X
1
⏟
k
1
,
…
,
X
n
,
…
,
X
n
⏟
k
n
)
.
{\displaystyle \kappa _{k_{1},\dots ,k_{n}}(X_{1},\ldots ,X_{n})=\kappa _{1,\ldots ,1}(\underbrace {X_{1},\dots ,X_{1}} _{k_{1}},\ldots ,\underbrace {X_{n},\dots ,X_{n}} _{k_{n}}).}
κ
2
,
0
,
1
(
X
,
Y
,
Z
)
=
κ
(
X
,
X
,
Z
)
,
{\displaystyle \kappa _{2,0,1}(X,Y,Z)=\kappa (X,X,Z),\,}
κ
0
,
0
,
n
,
0
(
X
,
Y
,
Z
,
T
)
=
κ
n
(
Z
)
=
κ
(
Z
,
…
,
Z
⏟
n
)
.
{\displaystyle \kappa _{0,0,n,0}(X,Y,Z,T)=\kappa _{n}(Z)=\kappa (\underbrace {Z,\dots ,Z} _{n}).\,}
複雑な瞬間との関係
結合キュムラントまたはランダム変数は、それらの混合モーメント の積の交互和として表現することができます 。 [11]
の式(3.2.7)を参照してください。
ここ で、 πは {1, ..., n } のすべてのパーティションのリストを実行します 。ここで、 Bは パーティションπ のすべてのブロックのリストを実行します 。そして、 | π |は パーティション内の部分の数です。
κ
(
X
1
,
…
,
X
n
)
=
∑
π
(
|
π
|
−
1
)
!
(
−
1
)
|
π
|
−
1
∏
B
∈
π
E
(
∏
i
∈
B
X
i
)
{\displaystyle \kappa (X_{1},\dots ,X_{n})=\sum _{\pi }(|\pi |-1)!(-1)^{|\pi |-1}\prod _{B\in \pi }E\left(\prod _{i\in B}X_{i}\right)}
例えば、
は の期待値 、
は および の 共分散 、そして
κ
(
X
)
=
E
(
X
)
,
{\displaystyle \kappa (X)=\operatorname {E} (X),}
X
{\textstyle X}
κ
(
X
,
Y
)
=
E
(
X
Y
)
−
E
(
X
)
E
(
Y
)
,
{\displaystyle \kappa (X,Y)=\operatorname {E} (XY)-\operatorname {E} (X)\operatorname {E} (Y),}
X
{\textstyle X}
Y
{\textstyle Y}
κ
(
X
,
Y
,
Z
)
=
E
(
X
Y
Z
)
−
E
(
X
Y
)
E
(
Z
)
−
E
(
X
Z
)
E
(
Y
)
−
E
(
Y
Z
)
E
(
X
)
+
2
E
(
X
)
E
(
Y
)
E
(
Z
)
.
{\displaystyle \kappa (X,Y,Z)=\operatorname {E} (XYZ)-\operatorname {E} (XY)\operatorname {E} (Z)-\operatorname {E} (XZ)\operatorname {E} (Y)-\operatorname {E} (YZ)\operatorname {E} (X)+2\operatorname {E} (X)\operatorname {E} (Y)\operatorname {E} (Z).\,}
平均ゼロのランダム変数 の場合、 が の分割で、 その分割が単一の を含む 場合、 形式の混合モーメントは消滅します 。したがって、混合モーメントに関するそれらの結合キュムラントの表現は単純化されます。たとえば、X、Y、Z、W が平均ゼロのランダム変数である場合、次の式が成り立ちます。
X
1
,
…
,
X
n
{\textstyle X_{1},\ldots ,X_{n}}
∏
B
∈
π
E
(
∏
i
∈
B
X
i
)
{\textstyle \prod _{B\in \pi }E\left(\prod _{i\in B}X_{i}\right)}
π
{\textstyle \pi }
{
1
,
…
,
n
}
{\textstyle \{1,\ldots ,n\}}
B
=
{
k
}
{\textstyle B=\{k\}}
κ
(
X
,
Y
,
Z
)
=
E
(
X
Y
Z
)
.
{\displaystyle \kappa (X,Y,Z)=\operatorname {E} (XYZ).\,}
κ
(
X
,
Y
,
Z
,
W
)
=
E
(
X
Y
Z
W
)
−
E
(
X
Y
)
E
(
Z
W
)
−
E
(
X
Z
)
E
(
Y
W
)
−
E
(
X
W
)
E
(
Y
Z
)
.
{\displaystyle \kappa (X,Y,Z,W)=\operatorname {E} (XYZW)-\operatorname {E} (XY)\operatorname {E} (ZW)-\operatorname {E} (XZ)\operatorname {E} (YW)-\operatorname {E} (XW)\operatorname {E} (YZ).\,}
より一般的には、マクローリン級数の係数は混合モーメントで表現することもできるが、簡潔な公式はない。実際、上で述べたように、ランダム変数を適切に繰り返すことで結合キュムラントとして表すことができ、次に上記の公式を適用して混合モーメントで表現することができる。例えば、
κ
201
(
X
,
Y
,
Z
)
=
κ
(
X
,
X
,
Z
)
=
E
(
X
2
Z
)
−
2
E
(
X
Z
)
E
(
X
)
−
E
(
X
2
)
E
(
Z
)
+
2
E
(
X
)
2
E
(
Z
)
.
{\displaystyle \kappa _{201}(X,Y,Z)=\kappa (X,X,Z)=\operatorname {E} (X^{2}Z)-2\operatorname {E} (XZ)\operatorname {E} (X)-\operatorname {E} (X^{2})\operatorname {E} (Z)+2\operatorname {E} (X)^{2}\operatorname {E} (Z).\,}
いくつかのランダム変数が他のすべてのランダム変数から独立している場合、2つ(またはそれ以上)の独立したランダム変数を含むキュムラントはゼロになります。 [ 要出典 ]
キュムラントによる混合モーメントの表現の組み合わせ論的意味は、キュムラントによる混合モーメントの表現の組み合わせ論的意味よりも理解しやすい。 [11]の式(3.2.6)を参照。
E
(
X
1
⋯
X
n
)
=
∑
π
∏
B
∈
π
κ
(
X
i
:
i
∈
B
)
.
{\displaystyle \operatorname {E} (X_{1}\cdots X_{n})=\sum _{\pi }\prod _{B\in \pi }\kappa (X_{i}:i\in B).}
例えば:
E
(
X
Y
Z
)
=
κ
(
X
,
Y
,
Z
)
+
κ
(
X
,
Y
)
κ
(
Z
)
+
κ
(
X
,
Z
)
κ
(
Y
)
+
κ
(
Y
,
Z
)
κ
(
X
)
+
κ
(
X
)
κ
(
Y
)
κ
(
Z
)
.
{\displaystyle \operatorname {E} (XYZ)=\kappa (X,Y,Z)+\kappa (X,Y)\kappa (Z)+\kappa (X,Z)\kappa (Y)+\kappa (Y,Z)\kappa (X)+\kappa (X)\kappa (Y)\kappa (Z).\,}
その他の特性
結合キュムラントのもう一つの重要な特性は多重線性です。
κ
(
X
+
Y
,
Z
1
,
Z
2
,
…
)
=
κ
(
X
,
Z
1
,
Z
2
,
…
)
+
κ
(
Y
,
Z
1
,
Z
2
,
…
)
.
{\displaystyle \kappa (X+Y,Z_{1},Z_{2},\dots )=\kappa (X,Z_{1},Z_{2},\ldots )+\kappa (Y,Z_{1},Z_{2},\ldots ).\,}
2 番目のキュムラントが分散であるのと同様に、2 つのランダム変数の結合キュムラントは共 分散 です。よく知られている恒等式は
キュムラントに一般化されます。
var
(
X
+
Y
)
=
var
(
X
)
+
2
cov
(
X
,
Y
)
+
var
(
Y
)
{\displaystyle \operatorname {var} (X+Y)=\operatorname {var} (X)+2\operatorname {cov} (X,Y)+\operatorname {var} (Y)\,}
κ
n
(
X
+
Y
)
=
∑
j
=
0
n
(
n
j
)
κ
(
X
,
…
,
X
⏟
j
,
Y
,
…
,
Y
⏟
n
−
j
)
.
{\displaystyle \kappa _{n}(X+Y)=\sum _{j=0}^{n}{n \choose j}\kappa (\,\underbrace {X,\dots ,X} _{j},\underbrace {Y,\dots ,Y} _{n-j}\,).\,}
条件付きキュムラントと全キュムラントの法則
全期待値の法則 と 全 分散の法則は、条件付きキュムラントに自然に一般化される 。n = 3 の場合、 キュムラントではなく
(中心) モーメントの言語で表現すると、次のようになる。
μ
3
(
X
)
=
E
(
μ
3
(
X
∣
Y
)
)
+
μ
3
(
E
(
X
∣
Y
)
)
+
3
cov
(
E
(
X
∣
Y
)
,
var
(
X
∣
Y
)
)
.
{\displaystyle \mu _{3}(X)=\operatorname {E} (\mu _{3}(X\mid Y))+\mu _{3}(\operatorname {E} (X\mid Y))+3\operatorname {cov} (\operatorname {E} (X\mid Y),\operatorname {var} (X\mid Y)).}
一般的に [12]
、
κ
(
X
1
,
…
,
X
n
)
=
∑
π
κ
(
κ
(
X
π
1
∣
Y
)
,
…
,
κ
(
X
π
b
∣
Y
)
)
{\displaystyle \kappa (X_{1},\dots ,X_{n})=\sum _{\pi }\kappa (\kappa (X_{\pi _{1}}\mid Y),\dots ,\kappa (X_{\pi _{b}}\mid Y))}
合計はインデックスの 集合 {1, ..., n }のすべての 分割 π にわたっており、
π 1 、...、 π b は、パーティション π のすべての「ブロック」です 。式 κ ( X π m ) は 、インデックスがパーティションのそのブロック内にあるランダム変数の結合キュムラントを示します。
条件付きキュムラントと条件付き期待値
特定の設定では、条件付きキュムラントと条件付き期待値の間に微分恒等式が成立する。例えば、 Y = X + Z ( Zは X に依存しない標準正規分布)とする と、任意の X に対して次の式が成り立つ。 [13]
この結果は指数族にも拡張できる。 [14]
κ
n
+
1
(
X
∣
Y
=
y
)
=
d
n
d
y
n
E
(
X
∣
Y
=
y
)
,
n
∈
N
,
y
∈
R
.
{\displaystyle \kappa _{n+1}(X\mid Y=y)={\frac {\mathrm {d} ^{n}}{\mathrm {d} y^{n}}}\operatorname {E} (X\mid Y=y),\,n\in \mathbb {N} ,\,y\in \mathbb {R} .}
統計物理学との関係
統計物理学 では、 多く の示量 、つまり特定のシステムの体積やサイズに比例する量は、ランダム変数のキュムラントと関連しています。深い関係は、大規模なシステムでは、エネルギーや粒子の数などの示量 (示量) は、(たとえば) ほぼ独立した領域の数に関連付けられたエネルギーの合計として考えることができることです。これらのほぼ独立したランダム変数のキュムラントが (ほぼ) 加算されるという事実から、示量 (示量) がキュムラントと関連していると予想するのは合理的です。
温度T の熱浴と平衡状態にあるシステムには 変動する内部エネルギー E があり、これは分布から抽出されたランダム変数と見なすことができます 。 システムの
分配関数は、 β = 1/( kT )
であり 、 k は ボルツマン定数 であり、エネルギー E との混同を避けるために期待値で はなく 表記法が使用されています。したがって、エネルギー E の 1 番目と 2 番目のキュムラントは、 平均エネルギーと熱容量を与えます。
E
∼
p
(
E
)
{\textstyle E\sim p(E)}
Z
(
β
)
=
⟨
exp
(
−
β
E
)
⟩
,
{\displaystyle Z(\beta )=\langle \exp(-\beta E)\rangle ,}
⟨
A
⟩
{\textstyle \langle A\rangle }
E
[
A
]
{\textstyle \operatorname {E} [A]}
⟨
E
⟩
c
=
∂
log
Z
∂
(
−
β
)
=
⟨
E
⟩
⟨
E
2
⟩
c
=
∂
⟨
E
⟩
c
∂
(
−
β
)
=
k
T
2
∂
⟨
E
⟩
∂
T
=
k
T
2
C
{\displaystyle {\begin{aligned}\langle E\rangle _{c}&={\frac {\partial \log Z}{\partial (-\beta )}}=\langle E\rangle \\[6pt]\langle E^{2}\rangle _{c}&={\frac {\partial \langle E\rangle _{c}}{\partial (-\beta )}}=kT^{2}{\frac {\partial \langle E\rangle }{\partial T}}=kT^{2}C\end{aligned}}}
ヘルムホルツ の自由エネルギー によって表されることにより、
熱力学量がエネルギーのキュムラント生成関数とさらに結び付けられます。 自由エネルギー の導関数である熱力学特性、たとえば 内部エネルギー 、 エントロピー 、 比熱 容量 はすべて、これらのキュムラントによって簡単に表すことができます。 その他の自由エネルギーは、磁場 や化学ポテンシャル などの他の変数の関数になることもあります 。たとえば
、 N は粒子の数、 は総ポテンシャルです。 ここでも、自由エネルギーの定義とキュムラント生成関数の間には密接な関係があるため、この自由エネルギーのさまざまな導関数は、 E と N の結合キュムラントによって表すことができます 。
F
(
β
)
=
−
β
−
1
log
Z
(
β
)
{\displaystyle F(\beta )=-\beta ^{-1}\log Z(\beta )\,}
μ
{\textstyle \mu }
Ω
=
−
β
−
1
log
(
⟨
exp
(
−
β
E
−
β
μ
N
)
⟩
)
,
{\displaystyle \Omega =-\beta ^{-1}\log(\langle \exp(-\beta E-\beta \mu N)\rangle ),\,}
Ω
{\textstyle \Omega }
歴史
キュムラントの歴史については アンダース・ハルト が論じている。 [15] [16]
キュムラントは1889年にソルヴァルド・N・ティーレ によって初めて導入され 、 半不変量 と呼ばれました。 [17] キュムラント と呼ばれたのは、 1932年の ロナルド・フィッシャー と ジョン・ウィシャート による論文においてでした 。 [18] フィッシャーはネイマンによって公にティーレの研究について思い出させられました。ネイマンは、フィッシャーにティーレの以前の引用文献が紹介されていたことにも言及しています。 [19] スティーブン・スティグラーは、 キュムラント という名前は ハロルド・ホテリング からの手紙でフィッシャーに提案されたと 述べています [ 要出典 ] 。1929年に発表された論文では、フィッシャーはそれを 累積モーメント関数 と呼んでいました 。 [20]
統計物理学における分配関数は、 1901 年に ジョサイア ウィラード ギブスによって導入されました。 [ 引用が必要 ] 自由エネルギーは、ギブスの自由エネルギーと呼ばれることがよくあります。 統計力学 では、キュムラントは1927 年の出版物に関連して アーセル関数 としても知られています。 [ 引用が必要 ]
一般化された設定におけるキュムラント
より一般的には、数列{ m n : n = 1, 2, 3, ... } のキュムラント (必ずしも任意の確率分布のモーメントではない) は、定義により、
任意の数列が収束するかどうかの問題を無視して、 n = 1, 2, 3, ... に対する κ n
の値が 形式的に、つまり代数のみで求められるものです。「キュムラントの問題」の困難さはすべて、形式的に作業するとなくなります。最も単純な例は、確率分布の 2 番目のキュムラントは常に非負でなければならず、それより高次のキュムラントがすべてゼロの場合にのみゼロになることです。形式的なキュムラントには、このような制約はありません。
1
+
∑
n
=
1
∞
m
n
t
n
n
!
=
exp
(
∑
n
=
1
∞
κ
n
t
n
n
!
)
,
{\displaystyle 1+\sum _{n=1}^{\infty }{\frac {m_{n}t^{n}}{n!}}=\exp \left(\sum _{n=1}^{\infty }{\frac {\kappa _{n}t^{n}}{n!}}\right),}
ベル番号
組合せ論 では 、 n 番目の ベル数は、サイズ n の集合の分割数です 。ベル数列のキュムラントはすべて 1 に等しくなります 。ベル数は、 期待値が 1 であるポアソン分布のモーメント です。
二項式型の多項式列のキュムラント
特性ゼロの体 における スカラー の任意 のシーケンス { κ n : n = 1, 2, 3, ... } は形式的キュムラントとみなされ、対応する形式的モーメントのシーケンス { μ ′ : n = 1, 2, 3, ... } が存在し、これは上記の多項式によって与えられます。 [ 明確化が必要 ] [ 引用が必要 ] これらの多項式に対して、次の方法で 多項式シーケンスを 構築します。多項式から、
これらと 1 つの追加変数 x の新しい多項式を作成し、
パターンを一般化します。パターンは、前述のパーティション内のブロックの数が x の指数であるというものです。各係数はキュムラントの多項式です。これらは 、 エリック テンプル ベルにちなんで名付けられた ベル多項式 です。 [ 引用が必要 ]
μ
6
′
=
κ
6
+
6
κ
5
κ
1
+
15
κ
4
κ
2
+
15
κ
4
κ
1
2
+
10
κ
3
2
+
60
κ
3
κ
2
κ
1
+
20
κ
3
κ
1
3
+
15
κ
2
3
+
45
κ
2
2
κ
1
2
+
15
κ
2
κ
1
4
+
κ
1
6
{\displaystyle {\begin{aligned}\mu '_{6}={}&\kappa _{6}+6\kappa _{5}\kappa _{1}+15\kappa _{4}\kappa _{2}+15\kappa _{4}\kappa _{1}^{2}+10\kappa _{3}^{2}+60\kappa _{3}\kappa _{2}\kappa _{1}+20\kappa _{3}\kappa _{1}^{3}\\&{}+15\kappa _{2}^{3}+45\kappa _{2}^{2}\kappa _{1}^{2}+15\kappa _{2}\kappa _{1}^{4}+\kappa _{1}^{6}\end{aligned}}}
p
6
(
x
)
=
κ
6
x
+
(
6
κ
5
κ
1
+
15
κ
4
κ
2
+
10
κ
3
2
)
x
2
+
(
15
κ
4
κ
1
2
+
60
κ
3
κ
2
κ
1
+
15
κ
2
3
)
x
3
+
(
45
κ
2
2
κ
1
2
)
x
4
+
(
15
κ
2
κ
1
4
)
x
5
+
(
κ
1
6
)
x
6
,
{\displaystyle {\begin{aligned}p_{6}(x)={}&\kappa _{6}\,x+(6\kappa _{5}\kappa _{1}+15\kappa _{4}\kappa _{2}+10\kappa _{3}^{2})\,x^{2}+(15\kappa _{4}\kappa _{1}^{2}+60\kappa _{3}\kappa _{2}\kappa _{1}+15\kappa _{2}^{3})\,x^{3}\\&{}+(45\kappa _{2}^{2}\kappa _{1}^{2})\,x^{4}+(15\kappa _{2}\kappa _{1}^{4})\,x^{5}+(\kappa _{1}^{6})\,x^{6},\end{aligned}}}
この多項式列は二項型 です 。実際、他の二項型列は存在しません。二項型のすべての多項式列は、その形式キュムラント列によって完全に決定されます。 [ 要出典 ]
自由キュムラント
上記の
結合キュムラントのモーメント-キュムラント公式では、 集合 { 1, ..., n }の すべての 分割について合計します。代わりに、 交差しない分割 についてのみ合計すると、モーメントについて これらの公式を解くことで、 上で扱った従来のキュムラントではなく、 自由キュムラントが得られます。これらの自由キュムラントは、ローランド・シュパイヒャーによって導入され、 自由確率 論で中心的な役割を果たしています。 [21] [22]その理論では、ランダム変数の 代数のテンソル積 で定義される ランダム変数 の 独立性 を考慮するのではなく 、ランダム変数の代数の 自由 積で定義されるランダム変数の 自由独立性 を考慮します。 [22]
E
(
X
1
⋯
X
n
)
=
∑
π
∏
B
∈
π
κ
(
X
i
:
i
∈
B
)
{\displaystyle \operatorname {E} (X_{1}\cdots X_{n})=\sum _{\pi }\prod _{B\,\in \,\pi }\kappa (X_{i}:i\in B)}
κ
{\textstyle \kappa }
正規分布 の2次以上の通常のキュムラント はゼロです。 ウィグナー半円分布 の2次以上の 自由 キュムラントはゼロです。 [22] これは、自由確率論におけるウィグナー分布の役割が従来の確率論における正規分布の役割と類似している点の1つです。
参照
参考文献
^ Weisstein, Eric W. 「Cumulant」。MathWorld より - Wolfram Web リソース。http://mathworld.wolfram.com/Cumulant.html
^ Kendall, MG, Stuart, A. (1969) 『統計学の高度理論 』第 1 巻 (第 3 版)。グリフィン、ロンドン。(セクション 3.12)
^ Lukacs, E. (1970) 特性関数 (第2版)グリフィン、ロンドン(27ページ)
^ Lukacs, E. (1970) 特性関数 (第2版)グリフィン、ロンドン。(セクション2.4)
^ Aapo Hyvarinen、Juha Karhunen、Erkki Oja (2001) 独立成分分析 、 John Wiley & Sons 。 (セクション2.7.2)
^ ハメダニ、GG;フォルクマー、ハンス。 Behboodian、J. (2012-03-01)。 「部分独立確率変数と二変量混合のクラスに関するメモ」。 ハンガリカの数学研究室 。 49 (1): 19-25。 土井 :10.1556/SScMath.2011.1183。
^ Lukacs, E. (1970) 特性関数(第2版)、グリフィン、ロンドン。(定理7.3.5)
^ Bazant, Martin (2005年2月4日). 「MIT 18.366 | 2006年秋 | 大学院 | ランダムウォークと拡散、講義2: モーメント、キュムラント、スケーリング」. MIT OpenCourseWare . 2022年10月7日時点のオリジナルよりアーカイブ 。 2023年9月3日 閲覧。
^スミス、ピーター J. (1995 年 5月)。 「 キュムラントからモーメントを得るという古い問題の再帰的定式化とその逆」。 アメリカ 統計学者 。49 (2): 217–218。doi :10.2307/2684642。JSTOR 2684642。
^ Rota, G.-C.; Shen, J. (2000). 「キュムラントの組合せ論について」. 組合せ理論ジャーナル、シリーズ A. 91 ( 1–2): 283–304. doi : 10.1006/jcta.1999.3017 .
^ abcペッカティ 、 ジョヴァンニ; タク、ムラド S. (2011). 「ウィーナーカオス:モーメント、キュムラント、ダイアグラム」。 ボッコーニ &シュプリンガーシリーズ 。1。doi :10.1007 / 978-88-470-1679-8。ISBN 978-88-470-1678-1 . ISSN 2039-1471.
^ Brillinger, DR (1969). 「条件付けによるキュムラントの計算」. 統計数学研究所紀要 . 21 : 215–218. doi :10.1007/bf02532246. S2CID 122673823.
^ Dytso, Alex; Poor, H. Vincent; Shamai Shitz, Shlomo (2023). 「ガウスノイズにおける条件付き平均推定: 応用を伴うメタ微分恒等式」. IEEE Transactions on Information Theory . 69 (3): 1883–1898. doi :10.1109/TIT.2022.3216012. S2CID 253308274.
^ Dytso, Alex; Cardone, Martina; Zieder, Ian (2023). 「条件付き期待値のメタ微分同一性」. IEEE Transactions on Information Theory . 69 (7): 4284–4302. doi :10.1109/TIT.2023.3249163. S2CID 257247930.
^ Hald, A. (2000)「キュムラントと グラム・シャルリエ級数 の初期の歴史」 国際統計評論 、68(2):137–153。( Lauritzen, Steffen L. 編 (2002) に再録。Thiele : Pioneer in Statistics。Oxford UP ISBN
978-0-19-850972-1 。 )
^ Hald, Anders (1998). 1750年から1930年までの数理統計学の歴史 。ニューヨーク: Wiley。ISBN
978-0-471-17912-2 。
^ H. Cramér (1946) Mathematical Methods of Statistics、プリンストン大学出版局、セクション15.10、p. 186。
^ フィッシャー、RA 、 ジョン・ウィシャート、J. (1932) 2方向分割のパターン式のより単純なパターンからの導出、 ロンドン数学会 紀要、シリーズ2、v. 33、pp. 195–208 doi :10.1112/plms/s2-33.1.195
^ Neyman, J. (1956): 「Sir Ronald Fisher の論文に関する注釈」、 Journal of the Royal Statistical Society 、シリーズ B (方法論)、18、pp. 288–94。
^ Fisher, RA (1929). 「標本分布のモーメントと積モーメント」 (PDF) . ロンドン数学会紀要 . 30 : 199–238. doi :10.1112/plms/s2-30.1.199. hdl : 2440/15200 .
^ Speicher, Roland (1994). 「非交差パーティションの格子上の乗法関数と自由畳み込み」. Mathematische Annalen . 298 (4): 611–628. doi :10.1007/BF01459754. S2CID 123022311.
^ abc Novak, Jonathan; Śniady, Piotr (2011). 「自由キュムラントとは何か?」 アメリカ数学会誌 . 58 (2): 300–301. ISSN 0002-9920.
外部リンク