跳到主要内容

東京大学 新領域創成科学研究科 メディカル情報生命専攻 2014年8月実施 問題11

Author

zephyr

Description

For an arbitrary random variable XX that takes values in non-negative integers, we define the probability generating function φX(s)\varphi_X(s) of XX as φX(s)=E{sX}=k=0skPr{X=k}\varphi_X(s) = E\{s^X\} = \sum_{k=0}^{\infty} s^k \Pr\{X = k\}. Here, E{A}E\{A\} denotes the expected value of AA, Pr{X=k}\Pr\{X = k\} denotes the probability that XX assumes value kk, and ss represents a real number.

(1) Show the following equalities (a), (b).

  • (a) φX(1)=1\varphi_X(1) = 1
  • (b) dφXds(1)=E{X}\frac{d \varphi_X}{ds}(1) = E\{X\}

In the following, NN and UiU_i (i=1,2,)(i = 1, 2, \ldots) are mutually independent, identically distributed random variables that take values in non-negative integers.

(2) For a positive integer nn, we define a random variable Yn=Yn(U1,U2,)=i=1nUiY_n = Y_n(U_1, U_2, \ldots) = \sum_{i=1}^n U_i. Show φYn(s)=φN(s)n\varphi_{Y_n}(s) = \varphi_N(s)^n.

(3) We define a random variable W=W(N,Y1,Y2,)=n=1YnI(N=n)W = W(N, Y_1, Y_2, \ldots) = \sum_{n=1}^{\infty} Y_n I(N = n). Here,

I(N=n)={1if N=n0if NnI(N = n) = \begin{cases} 1 & \text{if } N = n \\ 0 & \text{if } N \neq n \end{cases}

Show φW(s)=φN(φN(s))\varphi_W(s) = \varphi_N(\varphi_N(s)).

(Hint: Pr{W=k}=n=1Pr{Yn=k}Pr{N=n}\Pr\{W = k\} = \sum_{n=1}^{\infty} \Pr\{Y_n = k\} \Pr\{N = n\} for a positive integer kk)

(4) For Pr{N=k}=qk(1q)\Pr\{N = k\} = q^k (1 - q), (0<q<1)(0 < q < 1), calculate E{N}E\{N\} and E{W}E\{W\}.


对于一个取非负整数值的任意随机变量 XX,我们定义 XX 的概率生成函数 φX(s)\varphi_X(s)φX(s)=E{sX}=k=0skPr{X=k}\varphi_X(s) = E\{s^X\} = \sum_{k=0}^{\infty} s^k \Pr\{X = k\}。其中,E{A}E\{A\} 表示 AA 的期望值,Pr{X=k}\Pr\{X = k\} 表示 XX 取值为 kk 的概率,ss 表示一个实数。

(1) 证明以下等式 (a), (b)。

  • (a) φX(1)=1\varphi_X(1) = 1
  • (b) dφXds(1)=E{X}\frac{d \varphi_X}{ds}(1) = E\{X\}

在下列情形中,NNUiU_i (i=1,2,)(i = 1, 2, \ldots) 是相互独立的同分布随机变量,取非负整数值。

(2) 对于一个正整数 nn,我们定义一个随机变量 Yn=Yn(U1,U2,)=i=1nUiY_n = Y_n(U_1, U_2, \ldots) = \sum_{i=1}^n U_i。证明 φYn(s)=φN(s)n\varphi_{Y_n}(s) = \varphi_N(s)^n

(3) 我们定义一个随机变量 W=W(N,Y1,Y2,)=n=1YnI(N=n)W = W(N, Y_1, Y_2, \ldots) = \sum_{n=1}^{\infty} Y_n I(N = n)。 其中,

I(N=n)={1如果 N=n0如果 Nn I(N = n) = \begin{cases} 1 & \text{如果 } N = n \\ 0 & \text{如果 } N \neq n \end{cases}

证明 φW(s)=φN(φN(s))\varphi_W(s) = \varphi_N(\varphi_N(s))

(提示:Pr{W=k}=n=1Pr{Yn=k}Pr{N=n}\Pr\{W = k\} = \sum_{n=1}^{\infty} \Pr\{Y_n = k\} \Pr\{N = n\} 对于一个正整数 kk

(4) 对于 Pr{N=k}=qk(1q)\Pr\{N = k\} = q^k (1 - q), (0<q<1)(0 < q < 1),计算 E{N}E\{N\}E{W}E\{W\}

Kai

(1)

(a)

φX(1)=E{1X}=k=01kPr{X=k}=k=0Pr{X=k}=1\begin{aligned} \varphi_X(1) &= E\{1^X\} = \sum_{k=0}^{\infty} 1^k \Pr\{X = k\} \\ &= \sum_{k=0}^{\infty} \Pr\{X = k\} \\ &= 1 \end{aligned}

The last step follows from the fact that the sum of probabilities over all possible outcomes is 1.

(b)

dφXds(s)=ddsE{sX}=ddsk=0skPr{X=k}=k=0ksk1Pr{X=k}\begin{aligned} \frac{d \varphi_X}{ds}(s) &= \frac{d}{ds} E\{s^X\} = \frac{d}{ds} \sum_{k=0}^{\infty} s^k \Pr\{X = k\} \\ &= \sum_{k=0}^{\infty} k s^{k-1} \Pr\{X = k\} \end{aligned}

Evaluating at s=1s = 1:

dφXds(1)=k=0k1k1Pr{X=k}=k=0kPr{X=k}=E{X}\begin{aligned} \frac{d \varphi_X}{ds}(1) &= \sum_{k=0}^{\infty} k \cdot 1^{k-1} \Pr\{X = k\} \\ &= \sum_{k=0}^{\infty} k \Pr\{X = k\} \\ &= E\{X\} \end{aligned}

(2)

We need to show φYn(s)=φN(s)n\varphi_{Y_n}(s) = \varphi_N(s)^n where Yn=i=1nUiY_n = \sum_{i=1}^n U_i.

φYn(s)=E{sYn}=E{si=1nUi}=E{sU1sU2...sUn}=E{sU1}E{sU2}...E{sUn}(due to i.i.d)=(E{sN})n(due to i.i.d)=φN(s)n\begin{aligned} \varphi_{Y_n}(s) &= E\{s^{Y_n}\} = E\{s^{\sum_{i=1}^n U_i}\} \\ &= E\{s^{U_1} \cdot s^{U_2} \cdot ... \cdot s^{U_n}\} \\ &= E\{s^{U_1}\} \cdot E\{s^{U_2}\} \cdot ... \cdot E\{s^{U_n}\} \quad \text{(due to i.i.d)} \\ &= (E\{s^{N}\})^n \quad \text{(due to i.i.d)} \\ &= \varphi_N(s)^n \end{aligned}

(3)

We need to show φW(s)=φN(φU(s))\varphi_W(s) = \varphi_N(\varphi_U(s)) where W=n=1YnI(N=n)W = \sum_{n=1}^{\infty} Y_n I(N = n).

Using the hint and the definition of probability generating function:

φW(s)=E{sW}=k=0skPr{W=k}=k=0skn=1Pr{Yn=k}Pr{N=n}=n=1Pr{N=n}k=0skPr{Yn=k}=n=1Pr{N=n}φYn(s)=n=1Pr{N=n}φN(s)n(from result of part 2)=φN(φN(s))\begin{aligned} \varphi_W(s) &= E\{s^W\} = \sum_{k=0}^{\infty} s^k \Pr\{W = k\} \\ &= \sum_{k=0}^{\infty} s^k \sum_{n=1}^{\infty} \Pr\{Y_n = k\} \Pr\{N = n\} \\ &= \sum_{n=1}^{\infty} \Pr\{N = n\} \sum_{k=0}^{\infty} s^k \Pr\{Y_n = k\} \\ &= \sum_{n=1}^{\infty} \Pr\{N = n\} \varphi_{Y_n}(s) \\ &= \sum_{n=1}^{\infty} \Pr\{N = n\} \varphi_N(s)^n \quad \text{(from result of part 2)} \\ &= \varphi_N(\varphi_N(s)) \end{aligned}

(4)

Given Pr{N=k}=qk(1q)\Pr\{N = k\} = q^k (1 - q), (0<q<1)(0 < q < 1), we need to calculate E{N}E\{N\} and E{W}E\{W\}.

First, let's calculate E{N}E\{N\}:

E{N}=k=0kPr{N=k}=k=0kqk(1q)=(1q)k=0kqk=(1q)q(1q)2=q1q\begin{aligned} E\{N\} &= \sum_{k=0}^{\infty} k \Pr\{N = k\} = \sum_{k=0}^{\infty} k q^k (1 - q) \\ &= (1 - q) \sum_{k=0}^{\infty} k q^k = (1 - q) \frac{q}{(1-q)^2} = \frac{q}{1-q} \end{aligned}

Now, for E{W}E\{W\}, we can use the result from part 1(b) and part 3:

E{W}=dφWds(1)=ddsφN(φN(s))s=1=φN(φN(1))φN(1)=E{N}2(using part 1(b))=(q1q)2\begin{aligned} E\{W\} &= \frac{d \varphi_W}{ds}(1) = \frac{d}{ds} \varphi_N(\varphi_N(s)) |_{s=1} \\ &= \varphi_N'(\varphi_N(1)) \cdot \varphi_N'(1) \\ &= E\{N\}^2 \quad \text{(using part 1(b))} \\ &= (\frac{q}{1-q})^2 \end{aligned}

Therefore, E{W}=(q1q)2E\{W\} = (\frac{q}{1-q})^2.

Knowledge

概率论 概率生成函数 条件期望 全期望公式 复合分布

难点思路

  1. 理解概率生成函数的定义和基本性质
  2. 利用独立性推导和的概率生成函数
  3. 使用条件期望和全期望公式推导复合随机变型的概率生成函数
  4. 应用概率生成函数的性质计算具体分布的期望

解题技巧和信息

  1. 概率生成函数的基本性质:
    • φX(1)=1\varphi_X(1) = 1
    • dφXds(1)=E{X}\frac{d \varphi_X}{ds}(1) = E\{X\}
    • φX+Y(s)=φX(s)φY(s)\varphi_{X+Y}(s) = \varphi_X(s) \cdot \varphi_Y(s) (对于独立的 XXYY)
  2. 几何分布的概率生成函数:如果 XGeo(p)X \sim Geo(p),则 φX(s)=p1(1p)s\varphi_X(s) = \frac{p}{1-(1-p)s}
  3. 利用全期望公式:E{Y}=E{E{YX}}E\{Y\} = E\{E\{Y|X\}\}

重点词汇

  • Probability generating function: 概率生成函数
  • Independent and identically distributed (i.i.d.): 独立同分布
  • Compound distribution: 复合分布
  • Conditional expectation: 条件期望
  • Law of total expectation: 全期望公式
  • Geometric distribution: 几何分布