確率母関数

確率母関数（かくりつぼかんすう、英: probability generating function）は、確率論において離散的な確率変数に対し定義される関数の一つ。確率変数の確率質量関数の冪級数表現、すなわち母関数である。確率母関数を用いることで確率分布をより簡潔に表せる状況で利用される。また、非負なる係数をもつ冪級数に関する、よく発展している理論を確率分布に適用する際にも用いられる。

定義

単変量の場合

$X$ を、非負なる整数 $\{0,1,\ldots \}$ を値としてとる離散的確率変数とする。この時、 $X$ の確率母関数は次で定義される: ^[1]

$G(z)=\operatorname {E} (z^{X})=\sum _{x=0}^{\infty }p(x)z^{x}$ ここで、 $p$ は $X$ の確率質量関数であり、 $\operatorname {E} (\cdot )$ は期待値を表す。なお、 $X$ という特定の確率変数に対し定義されるということを強調するために、 $G$ や $p$ が下添字付きの $G_{X}$ や $p_{X}$ として書き表されることも多々ある。この無限級数は少なくとも $|z|<1$ を満たす任意の複素数 $z$ に対して絶対収束するが、実際の収束半径はそれよりも大きい場合も多い。

多変量の場合

$X=(X_{1},\ldots ,X_{d})$ を $d$ 次元の確率変数とし、 $X_{i}$ はそれぞれ非負整数を値として取るとする。この時、 $X$ の確率母関数は次で定義される: $G(z)=G(z_{1},\ldots ,z_{d})=\operatorname {E} {\bigl (}z_{1}^{X_{1}}\cdots z_{d}^{X_{d}}{\bigr )}=\sum _{x_{1},\ldots ,x_{d}=0}^{\infty }p(x_{1},\ldots ,x_{d})z_{1}^{x_{1}}\cdots z_{d}^{x_{d}},$ ここで、再び $p$ は $X$ の確率質量関数であり、 $\operatorname {E} (\cdot )$ は期待値を表す。この無限級数は少なくとも $z=(z_{1},...z_{d})\in \mathbb {C} ^{d}$ であって、 ${\text{max}}\{|z_{1}|,...,|z_{d}|\}\leq 1$ を満たす任意の $d$ 次元複素数ベクトル $z$ に対して絶対収束する。

性質

冪級数として

定義式から自明であるが、確率母関数は非負なる係数をもつ冪級数の性質全てに従う。特に、 $G(1^{-})=\lim _{x\to 1,x<1}G(x)$ を、 $x$ を下側から1に近づける極限として定義した場合、 $G(1^{-})=1$ が成立する。これは、確率分布全体の和が1とならなければならないからである。そのため、非負なる係数をもつべき級数に対するアーベルの連続性定理の帰結として、任意の確率母関数の収束半径は最低でも1となる。

確率と期待値

以下の性質により、 $X$ に関する様々な基本的量を確率母関数から計算できる:

確率母関数の $k$ -次微分を取ることで、 $X=k$ となる確率が得られる: $p(k)=\operatorname {Pr} (X=k)={\frac {G^{(k)}(0)}{k!}}$
確率変数 $X$ 、 $Y$ の確率母関数が一致する( $G_{X}=G_{Y}$ )ならば、確率質量関数も一致する( $p_{X}=p_{Y}$ )。これは1.で述べた性質より導かれる。
確率質量関数が規格化されていることは、確率母関数を用いて次のように表される: $\operatorname {E} [1]=G(1^{-})=\sum _{k=0}^{\infty }p(k)=1$ また、 $X$ の期待値は $\operatorname {E} [X]=G'(1^{-})$ で表されることもわかる。さらに一般的には、 $X$ の $k$ -次の階乗モーメント $\operatorname {E} [X(X-1)\cdots (X-k+1)]$ が $\operatorname {E} \left[{\frac {X!}{(X-k)!}}\right]=G^{(k)}(1^{-}),\quad k\geq 0$ で表される。よって、 $X$ の分散は $\operatorname {Var} (X)=G''(1^{-})+G'(1^{-})-\left[G'(1^{-})\right]^{2}$ により計算される。一方、 $X$ の $k$ -次のモーメントは $\operatorname {E} [X^{k}]=\left(z{\frac {\partial }{\partial z}}\right)^{k}G(z){\Big |}_{z=1^{-}}$ により与えられる。
$M_{X}(t)$ を、確率変数 $X$ の積率母関数としたとき、 $X$ の確率母関数と $G_{X}(e^{t})=M_{X}(t)$ により結びつく。

独立な確率変数に関する関数

確率母関数は、独立な確率変数に対する関数を考察する際にとりわけ有用である。以下にその例を数点示す:

$X_{i}(i=1,2,\cdots ,N)$ が独立な(ただし必ずしも同じ分布には属しない)確率変数列とし、それぞれ非負整数を値として取りうるとする。 $a_{i}$ を整数の定数列としたとき、確率変数 $S_{N}=\sum _{i=1}^{N}a_{i}X_{i}$ に関する確率母関数は、以下のようにそれぞれの $X_{i}$ に関する確率母関数を用いて記述できる:

$G_{S_{N}}(z)=\operatorname {E} \left(z^{\sum _{i=1}^{N}a_{i}X_{i},}\right)=\prod _{k=1}^{N}G_{X_{k}}(z^{a_{k}})$

特に $X,Y$ が独立かつそれぞれ非負整数値を取る時、以下が成立する:

${\begin{aligned}G_{X+Y}(z)&=G_{X}(z)\cdot G_{Y}(z)\\G_{X-Y}(z)&=G_{X}(z)\cdot G_{Y}(1/z)\end{aligned}}$

上述した例では、確率変数列に含まれる確率変数の個数 $N$ は固定されていた。いま、 $N$ もまた非負なる整数を値にとる確率変数であるとし、さらにその分布は $X_{i}$ とも独立であるとする。 $X_{i}$ が独立同分布であるとすれば、 $S_{N}=\sum _{i=1}^{N}X_{i}$ に関する確率母関数は $G_{S_{N}}(z)=G_{N}(G_{X}(z))$ で表せる。ここで、 $G_{N}$ を $N$ に関する確率母関数とした。これは以下のようにして示される: ${\begin{aligned}G_{S_{N}}(z)&=\operatorname {E} (z^{S_{N}})=\operatorname {E} {\big (}\operatorname {E} (z^{\sum _{i=1}^{N}X_{i}}\mid N){\big )}\\[4pt]&=\operatorname {E} {\big (}(G_{X}(z))^{N}{\big )}=G_{N}(G_{X}(z))\end{aligned}}$ この性質は、ゴルトン・ワトソン過程や複合ポアソン過程（英語版）を考察する際に有用である。
再び $N$ が確率変数であるとする。 $X_{i}$ を独立だが必ずしも同一分布とはみなせない確率変数列とする(ただし、 $N$ との独立性は保証されているものとする)。この場合、 $f_{n}=\Pr(N=n)$ として、 $G_{S_{N}}(z)=\sum _{n\geq 1}f_{n}\prod _{i=1}^{n}G_{X_{i}}(z)$ が成立する。独立同分布な $X_{i}$ に対しては直前の性質を再現するが、一般的な状況で $S_{N}$ を確率母関数の観点で分解する手法があることは時として有用である。

いくつかの分布に対する具体例

$X$ が退化分布に従うとき、すなわち $\Pr(X=c)=1$ かつ $\Pr(X\neq c)=0$ であるとき、 $X$ の確率母関数は $G_{X}(z)=z^{c}$ である。
試行回数 $n$ 、成功確率 $p$ を母数に持つ二項分布の確率母関数は $G(z)=\left(1-p+pz\right)^{n}$ である。注: これは、成功確率 $p$ を母数にもつベルヌーイ分布の確率母関数を $n$ 乗したものに一致する。
特に、公正なコインを１回投げる場合の確率母関数は $G(z)={\frac {1}{2}}+{\frac {z}{2}}$ となる。
$\{0,1,2\cdots \}$ を台とする負の二項分布であって、試行終了までの成功回数 $r$ と成功確率 $p$ を母数として持つものの確率母関数は、負の二項定理を用いることで $G(z)=\left({\frac {p}{1-(1-p)z}}\right)^{r}$ により与えられることが示される。なお、 $|z|<{\frac {1}{1-p}}$ であればこの関数は収束する。
注: これは、 $\{0,1,2,\cdots \}$ を台に持ち、成功確率が $1-p$ であるような幾何分布の確率母関数を $n$ 乗したものに一致する。
レートパラメータ $\lambda$ を母数に持つポアソン分布の確率母関数は $G(z)=e^{\lambda (z-1)}$ で与えられる。