跳到主要内容

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

Author​

zephyr, 祭音Myyura

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=0∞skPr⁡{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=1∞YnI(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 N≠nI(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=1∞Pr⁡{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(1−q)\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=0∞skPr⁡{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\}

在下列情形中,NN 和 UiU_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=1∞YnI(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如果 N≠n 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=1∞Pr⁡{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(1−q)\Pr\{N = k\} = q^k (1 - q), (0<q<1)(0 < q < 1),计算 E{N}E\{N\} 和 E{W}E\{W\}。

题目描述​

对任意取非负整数值的随机变量 XX,定义其概率生成函数

φX(s)=E{sX}=∑k=0∞skPr⁡{X=k},\varphi_X(s)=E\{s^X\}=\sum_{k=0}^{\infty}s^k\Pr\{X=k\},

其中 ss 为实数。回答下列问题:

  1. 证明 φX(1)=1\varphi_X(1)=1 以及

    dφXds∣s=1=E{X}.\left.\frac{d\varphi_X}{ds}\right|_{s=1}=E\{X\}.
  2. 以下设 N,U1,U2,…N,U_1,U_2,\ldots 为相互独立、同分布且取非负整数值的随机变量。对正整数 nn 定义

    Yn=∑i=1nUi,Y_n=\sum_{i=1}^{n}U_i,

    证明 φYn(s)=φN(s)n\varphi_{Y_n}(s)=\varphi_N(s)^n。

  3. 定义

    W=∑n=1∞YnI(N=n),I(N=n)={1,N=n,0,N≠n,W=\sum_{n=1}^{\infty}Y_n I(N=n),\qquad I(N=n)= \begin{cases} 1,&N=n,\\ 0,&N\ne n, \end{cases}

    证明 φW(s)=φN(φN(s))\varphi_W(s)=\varphi_N(\varphi_N(s))。可使用提示:对正整数 kk,

    Pr⁡{W=k}=∑n=1∞Pr⁡{Yn=k}Pr⁡{N=n}.\Pr\{W=k\}=\sum_{n=1}^{\infty}\Pr\{Y_n=k\}\Pr\{N=n\}.
  4. 当 Pr⁡{N=k}=qk(1−q)\Pr\{N=k\}=q^k(1-q)、0<q<10<q<1 时,计算 E{N}E\{N\} 与 E{W}E\{W\}。

Kai​

We work on 0≤s≤10\le s\le1, where the generating functions are finite. The derivative at s=1s=1 is understood from the left; it may be +∞+\infty if E[X]=+∞E[X]=+\infty. For 0<s<10<s<1, termwise differentiation is valid, and monotone convergence gives the limit as s↑1s\uparrow1.

(1)​

(a)​

φX(1)=E{1X}=∑k=0∞1kPr⁡{X=k}=∑k=0∞Pr⁡{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}=dds∑k=0∞skPr⁡{X=k}=∑k=1∞ksk−1Pr⁡{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=1}^{\infty} k s^{k-1} \Pr\{X = k\} \end{aligned}

Evaluating at s=1s = 1:

dφXds(1)=∑k=1∞kPr⁡{X=k}=E{X}\begin{aligned} \frac{d \varphi_X}{ds}(1) &= \sum_{k=1}^{\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{s∑i=1nUi}=E{sU1⋅sU2⋅...⋅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(φN(s))\varphi_W(s) = \varphi_N(\varphi_N(s)) where W=∑n=1∞YnI(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=0∞skPr⁡{W=k}=Pr⁡{N=0}+∑k=0∞sk∑n=1∞Pr⁡{Yn=k}Pr⁡{N=n}=Pr⁡{N=0}+∑n=1∞Pr⁡{N=n}∑k=0∞skPr⁡{Yn=k}=Pr⁡{N=0}+∑n=1∞Pr⁡{N=n}φYn(s)=∑n=0∞Pr⁡{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\} \\ &=\Pr\{N=0\}+\sum_{k=0}^{\infty}s^k \sum_{n=1}^{\infty}\Pr\{Y_n=k\}\Pr\{N=n\} \\ &=\Pr\{N=0\}+\sum_{n=1}^{\infty}\Pr\{N=n\} \sum_{k=0}^{\infty}s^k\Pr\{Y_n=k\} \\ &=\Pr\{N=0\}+\sum_{n=1}^{\infty}\Pr\{N=n\}\varphi_{Y_n}(s) \\ &= \sum_{n=0}^{\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(1−q)\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=0∞kPr⁡{N=k}=∑k=0∞kqk(1−q)=(1−q)∑k=0∞kqk=(1−q)q(1−q)2=q1−q\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))=(q1−q)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}=(q1−q)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) (对于独立的 XX 和 YY)
  2. 几何分布的概率生成函数:如果 X∼Geo(p)X \sim Geo(p),则 φX(s)=p1−(1−p)s\varphi_X(s) = \frac{p}{1-(1-p)s}
  3. 利用全期望公式:E{Y}=E{E{Y∣X}}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: 几何分布

Reference​