確率母関数

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

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

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| を満たす任意の複素数 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) \isin \mathbb{C}^dであって、 \text{max}\{|z_1|, ..., |z_d|\} \le 1 を満たす任意の d 次元複素数ベクトル z に対して絶対収束する。

## 性質
### 冪級数として
定義式から自明であるが、確率母関数は非負なる係数をもつ冪級数の性質全てに従う。特に、G(1^-) = \lim_{x\to 1, 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</math> により計算される。一方、<math>X</math> の<math>k</math>-次の[[モーメント (数学)|モーメント]]は<math display="block">\operatorname{E}[X^k] = \left(z\frac{\partial}{\partial z}\right)^k G(z) \Big|_{z=1^-}</math>により与えられる。
# <math>M_X(t)</math> を、確率変数 <math>X</math> の[[積率母関数]]としたとき、<math>X</math> の確率母関数と<math>G_X(e^t) = M_X(t)</math>により結びつく。

=== 独立な確率変数に関する関数 ===
確率母関数は、[[独立 (確率論)|独立]]な確率変数に対する関数を考察する際にとりわけ有用である。以下にその例を数点示す:
* <math>X_i (i=1,2,\cdots,N)</math> が独立な(ただし必ずしも同じ分布には属しない)確率変数列とし、それぞれ非負整数を値として取りうるとする。 <math>a_i</math>を整数の定数列としたとき、確率変数 <math>S_N = \sum_{i=1}^N a_i X_i</math>に関する確率母関数は、以下のようにそれぞれの<math>X_i</math>に関する確率母関数を用いて記述できる:
<math>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})</math>
*  特に <math>X,Y</math>が独立かつそれぞれ非負整数値を取る時、以下が成立する:

<math>\begin{align}
G_{X+Y}(z) &= G_X(z) \cdot G_Y(z)\\
G_{X-Y}(z) &= G_X(z) \cdot G_Y(1/z)
\end{align}</math>
* 上述した例では、確率変数列に含まれる確率変数の個数 <math>N</math> は固定されていた。いま、 <math>N</math> もまた非負なる整数を値にとる確率変数であるとし、さらにその分布は <math>X_i</math> とも独立であるとする。<math>X_i</math> が独立同分布であるとすれば、<math>S_N=\sum_{i=1}^NX_i</math> に関する確率母関数は<math display="block">G_{S_N}(z) = G_N(G_X(z))</math>で表せる。ここで、<math>G_N</math> を <math>N</math> に関する確率母関数とした。これは以下のようにして示される:<math display="block">
\begin{align}
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{align}
</math> この性質は、[[ゴルトン・ワトソン過程]]やを考察する際に有用である。
* 再び <math>N</math> が確率変数であるとする。<math>X_i</math> を独立だが必ずしも同一分布とはみなせない確率変数列とする(ただし、<math>N</math> との独立性は保証されているものとする)。この場合、<math>
f_n = \Pr(N=n)
</math>として、<math>
G_{S_N}(z) = \sum_{n \ge 1} f_n \prod_{i=1}^n G_{X_i}(z)
</math> が成立する。独立同分布な <math>X_i</math> に対しては直前の性質を再現するが、一般的な状況で <math>
S_N
</math> を確率母関数の観点で分解する手法があることは時として有用である。

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

## 関連する概念
確率母関数は母関数の一種である。詳細は形式的冪級数を参照せよ。また、確率母関数は確率質量関数に対してZ変換を行ったものと等価であり、そう呼ばれることもある。

確率変数に対する母関数には、他にも積率母関数、特性関数、キュムラント母関数などがある。 確率母関数は、\operatorname{E}\left[z^X\right]で表されるの特別な場合とみなせる。

## 脚注