積率母関数(moment generating function)とは、実数値の確率変数 $X$ について $M_X(t)=E[e^{tX}]\in(0,\infty]$ で定まる関数で、有限になる $t$ の範囲は $0$ を含む区間であり、$\log M_X$ はその上で凸である。独立な確率変数の和の積率母関数は積になり、$P(X\ge a)\le e^{-ta}M_X(t)$($t\ge0$)という指数型の裾の評価を与える。$0$ の近くで有限なら、すべての積率が有限で $M_X(t)=\sum_kE[X^k]t^k/k!$ と展開され、$0$ の近くで一致する積率母関数をもつ 2 つの確率変数は同じ分布に従う。この仮定がないと、すべての積率が等しくても分布が異なることがある。
前提知識: 期待値, 確率変数, Lebesgue積分, 冪級数
実数値の確率変数 $X$ について、実数 $t$ の関数
$$
M_X(t)=E\bigl[e^{tX}\bigr]
$$
を $X$ の 積率母関数(moment generating function)という。$e^{tX}=\sum_k\frac{t^k}{k!}X^k$ の期待値をとると、$t^k/k!$ の係数に $k$ 次の積率 $E[X^k]$ が現れる。このように積率(モーメント)を係数として「生み出す」関数なので、この名がある。
積率母関数は、独立な確率変数の和を積に変え、指数関数的に小さい裾の確率の評価(Chernoff の評価)を与え、$0$ の近くで有限なら分布を決める。ただし $e^{tX}$ は大きくなりうるので、期待値が有限とは限らない。$0$ 以外のすべての $t$ で $M_X(t)=\infty$ となる分布(Cauchy 分布)もあり、すべての積率が有限でも $t>0$ で $M_X(t)=\infty$ となる分布(対数正規分布)もある。この弱点を持たないのが 特性関数(確率論) で、積率母関数はその「実数版」にあたる。
この記事では、一般の確率空間の上で積率母関数を定め、有限になる範囲が区間であること(prop-mgf-basic)、独立な和が積になること(thm-mgf-sum)、$0$ の近くで有限なら冪級数に展開でき係数が積率になること(thm-mgf-series)、$0$ の近くで一致すれば分布が等しいこと(thm-mgf-unique)を証明する。$0$ の近くで有限という仮定を外すと、積率がすべて等しくても分布が違う例がある(ex-mgf-lognormal)。
以下、$(\Omega,\mathcal F,P)$ を確率空間とし、確率変数は実数値とする。値が $0$ 以上の確率変数 $Y$ の期待値 $E[Y]$ は Lebesgue積分 として $[0,\infty]$ に定まる。
確率変数 $X$ と実数 $t$ について
$$
M_X(t):=E\bigl[e^{tX}\bigr]\in(0,\infty]
$$
を $X$ の 積率母関数 という。$M_X(t)< \infty$ となる $t$ の全体を $D_X$ と書く。自然数 $k$ について $E[X^k]$ を $X$ の $k$ 次の 積率 という($E[\lvert X\rvert^k]< \infty$ のとき定まる)。
$e^{tX}>0$ なので $M_X(t)>0$ であり、$M_X(0)=1$ である。$M_X$ は $X$ の分布だけで決まる。この記事で「積率母関数が $0$ の近くで有限」とは、ある $\delta>0$ について $(-\delta,\delta)\subset D_X$ となることをいう。
$X$ を確率変数とする。
1:$s< v< u$ で $s,u\in D_X$ とする。$x\ge0$ なら $e^{vx}\le e^{ux}$、$x< 0$ なら $e^{vx}\le e^{sx}$ なので、すべての $x$ で $e^{vx}\le e^{sx}+e^{ux}$ である。期待値をとれば $M_X(v)\le M_X(s)+M_X(u)< \infty$ である。$0\in D_X$ と合わせて 1 を得る。
2:$e^{t(aX+b)}=e^{bt}e^{(at)X}$ の期待値をとればよい。
3:$p:=1/\lambda$、$p':=1/(1-\lambda)$ とすると $\frac1p+\frac1{p'}=1$ である。$f:=e^{\lambda sX}$、$g:=e^{(1-\lambda)uX}$ に Hölderの不等式 の定理「Hölderの不等式(積分)」を使うと
$$
E\bigl[e^{(\lambda s+(1-\lambda)u)X}\bigr]=E[fg]\le E[f^p]^{1/p}E[g^{p'}]^{1/p'}=M_X(s)^\lambda M_X(u)^{1-\lambda}
$$
である。両辺の対数をとれば $\log M_X$ の凸性である。$\blacksquare$
$\log M_X$ を キュムラント母関数(対数積率母関数)という。その凸性は大偏差の理論の出発点で、Cramérの定理(大偏差) の命題「対数積率母関数の性質」で詳しく扱われている。
| 分布 | $M_X(t)$ | $D_X$ |
|---|---|---|
| Bernoulli 分布($P(X=1)=p$、$P(X=0)=1-p$) | $1-p+pe^t$ | $\mathbb R$ |
| 二項分布 $B(n,p)$ | $(1-p+pe^t)^n$ | $\mathbb R$ |
| Poisson 分布(平均 $\lambda$) | $\exp(\lambda(e^t-1))$ | $\mathbb R$ |
| 区間 $[0,1]$ 上の一様分布 | $\frac{e^t-1}t$($t\ne0$)、$1$($t=0$) | $\mathbb R$ |
| 指数分布(密度 $\lambda e^{-\lambda x}$、$x\ge0$) | $\frac\lambda{\lambda-t}$ | $(-\infty,\lambda)$ |
| 正規分布 $N(\mu,\sigma^2)$ | $\exp(\mu t+\sigma^2t^2/2)$ | $\mathbb R$ |
| Cauchy 分布(密度 $\frac1{\pi(1+x^2)}$) | $1$($t=0$)、$\infty$($t\ne0$) | $\{0\}$ |
| 対数正規分布($e^Z$、$Z\sim N(0,1)$) | $t\le0$ で有限、$t>0$ で $\infty$ | $(-\infty,0]$ |
Bernoulli 分布と Poisson 分布は定義から計算でき($e^{-\lambda}\sum_k(\lambda e^t)^k/k!$)、二項分布は独立な Bernoulli 分布の和なので thm-mgf-sum から $n$ 乗になる。正規分布は 正規分布 の定理「積率母関数と特性関数」である。指数分布は、$t< \lambda$ なら $\int_0^\infty\lambda e^{-(\lambda-t)x}\,dx=\frac\lambda{\lambda-t}$ で、$t\ge\lambda$ なら被積分関数が $\lambda$ 以上なので積分は $\infty$ である。
Cauchy 分布:$t>0$ なら、$x\ge0$ で $e^{tx}\ge\frac{t^2x^2}2$ なので $M_X(t)\ge\int_0^\infty\frac{t^2x^2}{2\pi(1+x^2)}\,dx=\infty$ である。$t< 0$ は $-X$ が同じ分布に従うことから同じである。
対数正規分布 $X=e^Z$:$t\le0$ なら $0< e^{tX}\le1$ で有限である。$t>0$ なら $M_X(t)=\int_{-\infty}^\infty e^{te^z}\frac{e^{-z^2/2}}{\sqrt{2\pi}}\,dz$ で、$e^z\ge\frac{z^3}6$($z\ge0$)から被積分関数は $z\ge0$ で $\frac1{\sqrt{2\pi}}\exp(\frac{tz^3}6-\frac{z^2}2)$ 以上であり、これは $z\to\infty$ で $\infty$ に発散するので積分は $\infty$ である。一方 $E[X^k]=E[e^{kZ}]=e^{k^2/2}$(正規分布の積率母関数で $t=k$)なので、すべての積率は有限である。
確率変数 $X_1,\dots,X_n$ が 互いに独立 であるとは、すべての Borel集合 $B_1,\dots,B_n\subset\mathbb R$ について $P(X_1\in B_1,\dots,X_n\in B_n)=\prod_iP(X_i\in B_i)$ が成り立つことをいう。
$X_1,\dots,X_n$ が互いに独立で、$f_1,\dots,f_n\colon\mathbb R\to[0,\infty)$ が Borel 可測なら
$$
E\Bigl[\prod_{i=1}^nf_i(X_i)\Bigr]=\prod_{i=1}^nE\bigl[f_i(X_i)\bigr]\qquad(\text{値は }[0,\infty]\text{、}0\cdot\infty=0)
$$
である。
$f_i$ がすべて Borel 集合 $B_i$ の指示関数なら、両辺は $P(X_1\in B_1,\dots,X_n\in B_n)$ と $\prod_iP(X_i\in B_i)$ で、独立性の定義から等しい。
$j=1,\dots,n$ の順に、「$f_1,\dots,f_j$ は $0$ 以上の Borel 可測関数、$f_{j+1},\dots,f_n$ は指示関数」のとき等式が成り立つことを示す。$j-1$ の段階まで示したとし、$f_j$ 以外を固定する。両辺は $f_j$ について、$0$ 以上の係数での 1 次結合を保つので、$f_j$ が $0$ 以上の値の単関数(指示関数の $0$ 以上の係数の 1 次結合)でも成り立つ。一般の $f_j$ は、$0$ 以上の値の単関数の単調増加列 $\varphi_m\uparrow f_j$ の極限である(Lebesgue積分 の補題「単関数による近似」)。$f_j$ を $\varphi_m$ に替えた等式の両辺で $m\to\infty$ とすると、Lebesgue積分 の定理「単調収束定理」により左辺は左辺に、右辺の因子 $E[\varphi_m(X_j)]$ は $E[f_j(X_j)]$ に単調に近づき、等式は極限でも成り立つ($[0,\infty]$ の単調な列の積の極限は極限の積である。いずれかの因子が $0$ のときは、両辺がつねに $0$)。$j=n$ で主張を得る。$\blacksquare$
$X_1,\dots,X_n$ が互いに独立なら、$S:=X_1+\cdots+X_n$ についてすべての実数 $t$ で
$$
M_S(t)=M_{X_1}(t)\,M_{X_2}(t)\cdots M_{X_n}(t)
$$
である($\infty$ を含む等式)。とくに $D_S\supset D_{X_1}\cap\cdots\cap D_{X_n}$ である。
$e^{tS}=\prod_ie^{tX_i}$ であり、$f_i(x):=e^{tx}$ は $0$ 以上の連続関数なので、lem-mgf-indep から従う。各因子は正なので、右辺が有限であることと各因子が有限であることは同値である。$\blacksquare$
$X_1,\dots,X_n$ が独立で同じ分布に従えば $M_S=M_{X_1}^n$ である。表の二項分布はこの形である。正規分布では $\exp(\mu_1t+\sigma_1^2t^2/2)\exp(\mu_2t+\sigma_2^2t^2/2)$ がまた正規分布の積率母関数になるので、thm-mgf-unique と合わせて、独立な正規分布の和は正規分布であることが分かる(正規分布 の定理「独立な正規分布の和」の別証明)。独立性を外すと積の等式は成り立たず、逆に等式から独立性は結論できない(節「条件を外したときの反例」)。
ある $\delta>0$ について $(-\delta,\delta)\subset D_X$ とする。このとき
$\lvert t\rvert< \delta$ とする。$e^{\lvert t\rvert\lvert x\rvert}\le e^{\lvert t\rvert x}+e^{-\lvert t\rvert x}$ なので
$$
E\bigl[e^{\lvert t\rvert\lvert X\rvert}\bigr]\le M_X(\lvert t\rvert)+M_X(-\lvert t\rvert)< \infty\qquad(1)
$$
である。
1:$r:=\delta/2$ とおく。$e^{r\lvert x\rvert}\ge\frac{r^k\lvert x\rvert^k}{k!}$ なので $E[\lvert X\rvert^k]\le\frac{k!}{r^k}E[e^{r\lvert X\rvert}]< \infty$ である。
2:部分和 $S_N:=\sum_{k=0}^N\frac{(tX)^k}{k!}$ は $e^{tX}$ に各点で収束し、$\lvert S_N\rvert\le\sum_{k\le N}\frac{\lvert tX\rvert^k}{k!}\le e^{\lvert t\rvert\lvert X\rvert}$ で、右辺は (1) により可積分である。Lebesgue積分 の定理「優収束定理」と積分の線形性により $M_X(t)=\lim_NE[S_N]=\sum_k\frac{E[X^k]}{k!}t^k$ である。絶対収束は、Lebesgue積分 の定理「単調収束定理」により $\sum_k\frac{E[\lvert X\rvert^k]}{k!}\lvert t\rvert^k=E[e^{\lvert t\rvert\lvert X\rvert}]< \infty$ となることから従う。
3:2 により $M_X$ は $(-\delta,\delta)$ の上で、収束半径が $\delta$ 以上の冪級数の和に等しいので、冪級数 の系「何回でも微分できることと係数の公式」による。$\blacksquare$
3 から、$M_X'(0)=E[X]$、$M_X''(0)=E[X^2]$ であり、キュムラント母関数 $K:=\log M_X$ については $K'(0)=\frac{M_X'(0)}{M_X(0)}=E[X]$、$K''(0)=\frac{M_X''(0)M_X(0)-M_X'(0)^2}{M_X(0)^2}=V[X]$(分散 $E[X^2]-E[X]^2$)である。たとえば Poisson 分布では $K(t)=\lambda(e^t-1)$ なので、期待値と分散はどちらも $K'(0)=K''(0)=\lambda$ である。
指数分布では $\lvert t\rvert< \lambda$ で $M_X(t)=\frac\lambda{\lambda-t}=\sum_{k\ge0}\frac{t^k}{\lambda^k}$ である。thm-mgf-series の 2 と係数を比べて(冪級数 の系「何回でも微分できることと係数の公式」)、$\frac{E[X^k]}{k!}=\frac1{\lambda^k}$、すなわち $E[X^k]=\frac{k!}{\lambda^k}$ を得る。部分積分を $k$ 回繰り返す計算を、等比級数の展開に置き換えたことになる。
$X$ の 特性関数(確率論) を $\varphi_X(t):=E[e^{itX}]$ と書く。特性関数はつねに定まり、特性関数(確率論) の系「特性関数は分布を決める」により分布を決める。
確率変数 $X,Y$ と $\delta>0$ について、$(-\delta,\delta)\subset D_X\cap D_Y$ で、$\lvert t\rvert< \delta$ のとき $M_X(t)=M_Y(t)$ であるとする。このとき $X$ と $Y$ は同じ分布に従う。
段 1(積率の一致):thm-mgf-series の 2 により、$M_X$ と $M_Y$ は $(-\delta,\delta)$ で係数 $E[X^k]/k!$、$E[Y^k]/k!$ の冪級数の和に等しい。冪級数 の系「何回でも微分できることと係数の公式」により、すべての $k$ で $E[X^k]=E[Y^k]$ である。
段 2(特性関数の局所展開):実数 $a$ と $\lvert h\rvert< \delta$ について、$c_k(a):=E\bigl[(iX)^ke^{iaX}\bigr]/k!$ とおくと
$$
\varphi_X(a+h)=\sum_{k=0}^\infty c_k(a)\,h^k\qquad(2)
$$
であり、右辺は絶対収束する。実際、$e^{iaX}\sum_{k\le N}\frac{(ihX)^k}{k!}$ は $e^{i(a+h)X}$ に各点で収束し、絶対値は $e^{\lvert h\rvert\lvert X\rvert}$ 以下で、これは prf-mgf-series の (1) により可積分である。実部と虚部に Lebesgue積分 の定理「優収束定理」を使えば (2) を得る。絶対収束は $\lvert c_k(a)\rvert\le E[\lvert X\rvert^k]/k!$ と (1) による。$Y$ についても同じである。
段 3(局所的な一致は $\delta$ の幅に広がる):$\varphi_X$ と $\varphi_Y$ が $a$ を含むある開区間で一致すれば、(2) の両辺の実部と虚部は $h$ の実係数の冪級数で、$0$ の近くで一致するので、冪級数 の系「何回でも微分できることと係数の公式」により $c_k(a)$ は $X$ と $Y$ で等しい。よって (2) により、$\varphi_X=\varphi_Y$ は $(a-\delta,a+\delta)$ の全体で成り立つ。
段 4(全体へ):$a=0$ では $c_k(0)=i^kE[X^k]/k!$ なので、段 1 と (2) により $\varphi_X=\varphi_Y$ が $(-\delta,\delta)$ で成り立つ。$m\ge2$ について $(-\frac{m\delta}2,\frac{m\delta}2)$ で一致するとすると、$a=\pm\frac{(m-1)\delta}2$ はこの開区間の点なので、段 3 により $(a-\delta,a+\delta)$ でも一致し、$(-\frac{(m+1)\delta}2,\frac{(m+1)\delta}2)$ で一致する。帰納法によりすべての実数で $\varphi_X=\varphi_Y$ であり、特性関数(確率論) の系「特性関数は分布を決める」により $X$ と $Y$ は同じ分布に従う。$\blacksquare$
同じ「特性関数を冪級数でつないでいく」方法は、積率の増え方が $\limsup_k\frac1{2k}E[X^{2k}]^{1/(2k)}< \infty$ を満たすだけで分布が積率で決まることの証明にも使われる(Dur19 Theorem 3.3.25(p. 142))。$0$ の近くで有限という仮定は外せない(ex-mgf-lognormal、反例の表の 2 行目)。
$X_1,\dots,X_n$ が互いに独立で、それぞれ密度 $\lambda e^{-\lambda x}$($x\ge0$)の指数分布に従うとする。thm-mgf-sum により、$t< \lambda$ で $M_S(t)=\bigl(\frac\lambda{\lambda-t}\bigr)^n$ である。密度 $\frac{\lambda^n}{(n-1)!}x^{n-1}e^{-\lambda x}$($x\ge0$)のガンマ分布の積率母関数は、$t< \lambda$ で
$$
\int_0^\infty\frac{\lambda^n}{(n-1)!}x^{n-1}e^{-(\lambda-t)x}\,dx=\frac{\lambda^n}{(\lambda-t)^n}
$$
($\int_0^\infty x^{n-1}e^{-cx}\,dx=(n-1)!/c^n$)なので、thm-mgf-unique により $S$ はこのガンマ分布に従う。畳み込みの積分を繰り返す証明は ガンマ分布 の系「指数分布の和」にある。
確率変数 $X$ と実数 $a$ について、$t\ge0$ なら $P(X\ge a)\le e^{-ta}M_X(t)$ であり、$t\le0$ なら $P(X\le a)\le e^{-ta}M_X(t)$ である。とくに、$X_1,\dots,X_n$ が互いに独立で $X$ と同じ分布に従い、$S_n:=X_1+\cdots+X_n$ なら
$$
P(S_n\ge na)\le\exp\Bigl(-n\sup_{t\ge0}\bigl(ta-\log M_X(t)\bigr)\Bigr)
$$
である($M_X(t)=\infty$ の $t$ では $ta-\log M_X(t)=-\infty$ とする)。
$t\ge0$ なら、$X\ge a$ のとき $e^{t(X-a)}\ge1$ で、ほかでは $e^{t(X-a)}>0$ なので、指示関数について $\mathbf 1_{\{X\ge a\}}\le e^{t(X-a)}$ である。期待値をとれば $P(X\ge a)\le e^{-ta}M_X(t)$ を得る。$t\le0$ も同様である。後半は、$S_n$ に前半を使い、thm-mgf-sum で $M_{S_n}(t)=M_X(t)^n$ とすれば、すべての $t\ge0$ で $P(S_n\ge na)\le\bigl(e^{-ta}M_X(t)\bigr)^n=\exp(-n(ta-\log M_X(t)))$ となり、右辺の下限をとればよい。$\blacksquare$
$X\sim N(0,1)$ なら $M_X(t)=e^{t^2/2}$ で、$a\ge0$ について $t=a$ とすると $P(X\ge a)\le e^{-a^2/2}$ である。指数 $\sup_t(ta-\log M_X(t))$ は率関数と呼ばれ、すべての $t$ で $M_X(t)< \infty$ なら、$a>E[X]$ のとき $P(S_n\ge na)$ が実際にこの速さで $0$ に近づくことが Cramérの定理(大偏差) の定理「Cramér の定理」である。値が有界な独立な和についての指数型の評価も、この命題の $M_X(t)$ を上から評価して得られる。
$f_0(x):=\frac1{\sqrt{2\pi}\,x}\exp\bigl(-\frac{(\log x)^2}2\bigr)$($x>0$、ほかでは $0$)を対数正規分布($Z\sim N(0,1)$ のときの $e^Z$)の密度とし、$-1\le a\le1$ について
$$
f_a(x):=f_0(x)\bigl(1+a\sin(2\pi\log x)\bigr)
$$
とおく。すべての $k=0,1,2,\dots$ について
$$
\int_0^\infty x^kf_0(x)\sin(2\pi\log x)\,dx=0\qquad(3)
$$
が成り立つ。$k=0$ と $f_a\ge0$ から $f_a$ は確率密度であり、$f_a$ の分布の $k$ 次の積率は $f_0$ と同じ $e^{k^2/2}$ である。しかし $a\ne0$ なら、$A:=\{x>0\mid\sin(2\pi\log x)>0\}$ について $\int_Af_a-\int_Af_0=a\int_Af_0(x)\sin(2\pi\log x)\,dx\ne0$ なので、2 つの分布は異なる。どちらの積率母関数も $t>0$ で $\infty$ である(値が正なので $e^{tx}\ge\frac{(tx)^k}{k!}$ から $M(t)\ge\frac{t^ke^{k^2/2}}{k!}$ がすべての $k$ で成り立ち、右辺は $k\to\infty$ で発散する)。
(3) の確かめ:$x=e^z$ と置換すると左辺は $\frac1{\sqrt{2\pi}}\int_{-\infty}^\infty e^{kz-z^2/2}\sin(2\pi z)\,dz$ である。$z=u+k$ とおくと $kz-\frac{z^2}2=\frac{k^2}2-\frac{u^2}2$、$\sin(2\pi(u+k))=\sin(2\pi u)$ なので、これは $\frac{e^{k^2/2}}{\sqrt{2\pi}}\int e^{-u^2/2}\sin(2\pi u)\,du$ で、被積分関数は可積分な奇関数なので $0$ である。この例は Heyde による(Dur19 §3.3.5 の Counterexample 1(p. 140))。
| 外す条件・誤った推論 | 反例 | 成り立たなくなること |
|---|---|---|
| 積率の一致から分布の一致を結論するときの、$0$ の近くで有限 | ex-mgf-lognormal | すべての積率が等しくても、分布が異なる |
| thm-mgf-unique の、$0$ の近くで有限 | Cauchy 分布に従う $X$ と $2X$ | $M_X=M_{2X}$($0$ で $1$、ほかで $\infty$)だが、分布が異なる |
| thm-mgf-series の、$0$ の近くで有限 | 対数正規分布 $X=e^Z$ | 積率はすべて有限だが、$\sum_kE[X^k]t^k/k!$ は $t\ne0$ で発散する |
| thm-mgf-sum の独立性 | $X=Y\sim N(0,1)$ | $M_{X+Y}(t)=e^{2t^2}\ne e^{t^2}=M_X(t)M_Y(t)$ |
| 積の等式から独立を結論する | 確率母関数 の例「反例:確率母関数は積になるが独立でない組」 | $M_{X+Y}=M_XM_Y$ でも独立でない |
2 行目:Cauchy 分布の積率母関数は $t\ne0$ で $\infty$(表の計算)であり、prop-mgf-basic の 2 により $M_{2X}(t)=M_X(2t)$ も $t\ne0$ で $\infty$ である。$P(\lvert X\rvert\le1)=\frac12$、$P(\lvert2X\rvert\le1)=\frac2\pi\arctan\frac12< \frac12$ なので分布は異なる。$D_X=\{0\}$ なので仮定「$(-\delta,\delta)\subset D_X$」だけを破っている。
3 行目:$E[X^k]=e^{k^2/2}$ なので、$t\ne0$ では $\frac{e^{k^2/2}\lvert t\rvert^k}{k!}\ge\frac{e^{k^2/2}\lvert t\rvert^k}{k^k}=\bigl(\frac{e^{k/2}\lvert t\rvert}k\bigr)^k\to\infty$ であり、級数の項が $0$ に近づかない。
4 行目:$X+Y=2X\sim N(0,4)$ なので、正規分布の積率母関数から $M_{2X}(t)=e^{2t^2}$ である。
5 行目:その例の $X,Y$ は $\{0,1,2\}$ に値をとり、$M_X(t)=G_X(e^t)$ なので、確率母関数の等式 $G_{X+Y}=G_XG_Y$ に $s=e^t$ を代入すれば $M_{X+Y}=M_XM_Y$ である($s=e^t>1$ でも、値が有限個なので多項式の等式として成り立つ)。
Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する