キュムラント生成関数のいくつかの性質 キュムラント生成関数K ( t ) は、存在する場合、無限回微分可能 で凸関数 であり、原点を通ります。その 1 階微分は、確率分布のサポートの下限から 上限 までの開区間で単調に変化し、2 階微分は、単一の点質量の退化分布 を除いて、定義されているすべての場所で厳密に正です。キュムラント生成関数が存在するのは、分布の裾が指数関数的減衰によって支配されている 場合のみです 。つまり、(ビッグ O 記法を 参照 ) ∃ c > 0 、 F ( x ) = O ( e c x ) 、 x → − ∞ ; そして ∃ 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{ および }}\\[4pt]&\exists d>0,\,\,1-F(x)=O(e^{-dx}),x\to +\infty ;\end{aligned}}} ここでF {\textstyle F} は累積分布関数 です。キュムラント生成関数は、そのようなc の負の上限 が存在する場合はその上限に、またそのようなd の上限が存在する場合はその上限に 垂直漸近線を 持ち、そうでない場合はすべての実数に対して定義されます。
確率変数X のサポート が有限の上限または下限を持つ場合、そのキュムラント生成関数y = K ( t ) は、存在する場合、サポートの上限または下限に等しい傾きを持つ 漸近線(s) に近づきます。 y = ( t + 1 ) 情報 補足 X − μ ( X ) 、 そして y = ( t − 1 ) すする 補足 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 情報 補足 X − K ′ ( t ) ] d t 、 ∫ ∞ 0 [ t 情報 補足 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} K (0) = 0 なので、これらの漸近線のy 切片 が得られます。
分布をc だけシフトさせると、K X + c ( t ) = K X ( t ) + c t 。 {\textstyle K_{X+c}(t)=K_{X}(t)+ct.} c における縮退点質量の場合、キュムラント生成関数は直線である。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}} X とY が独立であり、かつそれらのキュムラント生成関数が存在する場合に限る。(部分独立性 と2次モーメントの存在は独立性を意味するのに十分である。[ 7 ] )
分布の自然指数族は、K(t)をシフトまたは平行移動し、常に原点を通るように垂直方向に調整することによって実現できます。fが キュムラント 生成関数 を 持つ確率 密度関数である場合K ( t ) = ログ 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 )がt1 < Re( t ) < t2 の 範囲で有限である場合、 t1 < 0 < t2 であれば、K ( t )は t1 < Re( t ) < t2 の範囲で解析的かつ無限回微分可能である。 さらに、tが実数でt1 < t < t2の場合、 K ( t ) は厳密 に凸 で あり 、 K ′ ( t ) は 厳密に 増加 する。
キュムラントのさらなる性質
陰性結果 正規分布 のキュムラントの結果を考慮すると、 低次のキュムラント(3 次からm − 1 次) がゼロでない、あるm > 3に対して κ m = κ m +1 = ⋯ = 0 となる分布の族が見つかることが期待されるかもしれない。そのような分布は存在しない。[ 8 ] ここでの根本的な結果は、キュムラント生成関数は 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 ) = ログ 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 階微分を評価することによってモーメントの観点から復元できます。ログ M ( t ) {\textstyle \log M(t)} で t = 0 {\displaystyle t=0} 、 κ n = K ( n ) ( 0 ) = d n ログ 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}} 不完全(または部分的)ベル多項式 です。
同様に、平均が次のように与えられる場合μ {\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 次中心モーメントはキュムラントを用いて次のように表される。 μ 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 ( ログ 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}}}
一般に、[ 9 ] キュムラントは行列の 行列式である。 κ 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}&{\binom {2}{1}}\mu '_{1}&1&0&0&\ldots &0\\\mu '_{4}&\mu '_{3}&{\binom {3}{1}}\mu '_{2}&{\binom {3}{2}}\mu '_{1}&1&0&\ldots &0\\\mu '_{5}&\mu '_{4}&{\binom {4}{1}}\mu '_{3}&{\binom {4}{2}}\mu '_{2}&{\binom {4}{3}}\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 &{\binom {l-1}{l-2}}\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}\,.}
キュムラントは、関係式K ( t ) = log M ( t )を t に関して微分すること によりモーメントと関連付けることができ、M′ ( t ) = K′ ( t ) M ( t ) が得られます。これは都合よく指数関数や対数を含んでいません。左辺と右辺のt n −1 / ( n −1)!の係数を等しくし、 μ′ 0 = 1を使用すると、 n ≥ 1 に対して次の式が得られます。[ 10 ] μ 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 の分割に縮退する集合の分割の数です。
結合キュムラント 複数の確率変数X 1 、 ...、 X n の同時キュムラント κ は 、多変量キュムラント生成関数のマクローリン級数の係数κ 1,...,1 ( X 1 、 ...、 X n )として定義されます。 [ 12 ] のセクション 3.1 を参照してください。G ( t 1 、 … 、 t n ) = ログ 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 ) = ログ E ( e ∑ j = 1 n 私 t j X j ) = ∑ k 1 、 … 、 k n κ k 1 、 … 、 k n 私 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 = ( − 私 ) 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 ) = ( − 私 ) 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 私 、 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}).\,} 特に、最後の等式は、単一の確率変数のキュムラントが、その確率変数の複数のコピーの同時キュムラントであることを示している。
混合モーメントとの関係 確率変数の同時キュムラントは、それらの混合モーメント の積の交互和として表現できます。 [ 12 ] の式(3.2.7)を参照してください。κ ( X 1 、 … 、 X n ) = ∑ π ( | π | − 1 ) ! ( − 1 ) | π | − 1 ∏ B ∈ π E ( ∏ 私 ∈ B X 私 ) {\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)} ここで、π は {1, ..., n } のすべての分割のリストを走査し、B は 分割π のすべてのブロックのリストを走査し、| π | は分割内の部分の数です。
例えば、 κ ( 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 1 、 … 、 X n {\textstyle X_{1},\ldots ,X_{n}} 、次の形式の混合モーメント∏ B ∈ π E ( ∏ 私 ∈ B X 私 ) {\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、Wが平均ゼロの確率変数である場合、次のようになります。 κ ( 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つ(またはそれ以上)の独立した確率変数を含むキュムラントはゼロになります。
混合モーメントをキュムラントで表した表現の組み合わせ論的な意味は、混合モーメントで表したキュムラントの意味よりも理解しやすい。式(3.2.6)を参照のこと。[ 12 ] E ( X 1 ⋯ X n ) = ∑ π ∏ B ∈ π κ ( X 私 : 私 ∈ 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つの確率変数の同時キュムラントは共分散 である。おなじみの恒等式 変数 ( X + Y ) = 変数 ( X ) + 2 カバー ( X 、 Y ) + 変数 ( 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 カバー ( E ( X ∣ Y ) 、 変数 ( 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)).}
一般的に、[ 13 ] κ ( 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 に対して次の関係が成り立ちます[ 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} .} この結果は指数型分布族にも拡張できる。[ 15 ]