東京大学 新領域創成科学研究科 メディカル情報生命専攻 2023年8月実施 問題11
Author
zephyr
Description
Let X i X_i X i ( i = 1 , … , n ; n ≥ 2 ) (i = 1, \dots, n; n \geq 2) ( i = 1 , … , n ; n ≥ 2 ) be n n n independent nonnegative real-valued random variables with the same probability density function f ( x ) = e − x f(x) = e^{-x} f ( x ) = e − x . Answer the following questions with mathematical derivation.
(1) Compute the mean and variance of X i X_i X i .
(2) Suppose the value of X i X_i X i is given as x i x_i x i for some i i i . Compute probability P ( X i ≤ X k ∣ X i = x i ) \mathbb{P}(X_i \leq X_k \mid X_i = x_i) P ( X i ≤ X k ∣ X i = x i ) that X i X_i X i is less than or equal to X k X_k X k for a given index k k k with k ≠ i k \neq i k = i .
(3) Compute probability P ( X i ≤ X k ) \mathbb{P}(X_i \leq X_k) P ( X i ≤ X k ) by multiplying the probability of (2) by f ( x i ) f(x_i) f ( x i ) and integrating it with respect to x i x_i x i .
(4) Let X min = min k = 1 , … , n X k X_{\min} = \min_{k=1,\dots,n} X_k X m i n = min k = 1 , … , n X k be the minimum value of set { X i ∣ i = 1 , … , n } \{ X_i \mid i = 1, \dots, n \} { X i ∣ i = 1 , … , n } . Compute the probability P ( x ≤ X min ) \mathbb{P}(x \leq X_{\min}) P ( x ≤ X m i n ) that X min X_{\min} X m i n is greater than or equal to x x x for a given positive real number x x x .
For n n n given positive real numbers λ i \lambda_i λ i ( i = 1 , … , n ) (i = 1, \dots, n) ( i = 1 , … , n ) , define random variables Z i Z_i Z i ( i = 1 , … , n ) (i = 1, \dots, n) ( i = 1 , … , n ) as Z i = X i λ i Z_i = \frac{X_i}{\lambda_i} Z i = λ i X i .
(5) Suppose the value of X i X_i X i is given as x i x_i x i for some i i i . Compute probability P ( Z i ≤ Z k ∣ X i = x i ) \mathbb{P}(Z_i \leq Z_k \mid X_i = x_i) P ( Z i ≤ Z k ∣ X i = x i ) that Z i Z_i Z i is less than or equal to Z k Z_k Z k for a given index k k k with k ≠ i k \neq i k = i .
(6) Suppose the value of X i X_i X i is given as x i x_i x i for some i i i . Compute probability P ( ⋂ k = 1 , … , n ; k ≠ i { Z i ≤ Z k } ∣ X i = x i ) \mathbb{P}\left(\bigcap_{k=1,\dots,n; k \neq i} \{ Z_i \leq Z_k \} \mid X_i = x_i\right) P ( ⋂ k = 1 , … , n ; k = i { Z i ≤ Z k } ∣ X i = x i ) that Z i Z_i Z i is less than or equal to Z k Z_k Z k for all the indices k k k with k ≠ i k \neq i k = i .
(7) Let Z I Z_I Z I be a minimum element in set { Z i ∣ i = 1 , … , n } \{ Z_i \mid i = 1, \dots, n \} { Z i ∣ i = 1 , … , n } . Answer the probability distribution P ( I = i ) \mathbb{P}(I = i) P ( I = i ) ( i = 1 , … , n ) (i = 1, \dots, n) ( i = 1 , … , n ) of index I = arg min k = 1 , … , n Z k I = \mathop{\arg\min}\limits_{k=1,\dots,n} Z_k I = k = 1 , … , n arg min Z k .
设 X i X_i X i ( i = 1 , … , n ; n ≥ 2 ) (i = 1, \dots, n; n \geq 2) ( i = 1 , … , n ; n ≥ 2 ) 是 n n n 个独立的非负实值随机变量,其概率密度函数为 f ( x ) = e − x f(x) = e^{-x} f ( x ) = e − x 。请用数学推导回答以下问题。
(1) 计算 X i X_i X i 的均值和方差。
(2) 假设 X i X_i X i 的值为 x i x_i x i ,计算 P ( X i ≤ X k ∣ X i = x i ) \mathbb{P}(X_i \leq X_k \mid X_i = x_i) P ( X i ≤ X k ∣ X i = x i ) 的概率,即 X i X_i X i 小于或等于 X k X_k X k 的概率,其中 k ≠ i k \neq i k = i 。
(3) 计算 P ( X i ≤ X k ) \mathbb{P}(X_i \leq X_k) P ( X i ≤ X k ) 的概率,通过 (2) 的概率乘以 f ( x i ) f(x_i) f ( x i ) 并对 x i x_i x i 积分。
(4) 令 X min = min k = 1 , … , n X k X_{\min} = \min_{k=1,\dots,n} X_k X m i n = min k = 1 , … , n X k 为集合 { X i ∣ i = 1 , … , n } \{ X_i \mid i = 1, \dots, n \} { X i ∣ i = 1 , … , n } 的最小值。计算 P ( x ≤ X min ) \mathbb{P}(x \leq X_{\min}) P ( x ≤ X m i n ) 的概率,即 X min X_{\min} X m i n 大于或等于给定正实数 x x x 的概率。
对于 n n n 个给定的正实数 λ i \lambda_i λ i ( i = 1 , … , n ) (i = 1, \dots, n) ( i = 1 , … , n ) ,定义随机变量 Z i Z_i Z i ( i = 1 , … , n ) (i = 1, \dots, n) ( i = 1 , … , n ) 为 Z i = X i λ i Z_i = \frac{X_i}{\lambda_i} Z i = λ i X i 。
(5) 假设 X i X_i X i 的值为 x i x_i x i ,计算 P ( Z i ≤ Z k ∣ X i = x i ) \mathbb{P}(Z_i \leq Z_k \mid X_i = x_i) P ( Z i ≤ Z k ∣ X i = x i ) 的概率,即 Z i Z_i Z i 小于或等于 Z k Z_k Z k 的概率,其中 k ≠ i k \neq i k = i 。
(6) 假设 X i X_i X i 的值为 x i x_i x i ,计算 P ( ⋂ k = 1 , … , n ; k ≠ i { Z i ≤ Z k } ∣ X i = x i ) \mathbb{P}\left(\bigcap_{k=1,\dots,n; k \neq i} \{ Z_i \leq Z_k \} \mid X_i = x_i\right) P ( ⋂ k = 1 , … , n ; k = i { Z i ≤ Z k } ∣ X i = x i ) 的概率,即 Z i Z_i Z i 小于或等于 Z k Z_k Z k 的概率,对于所有 k ≠ i k \neq i k = i 。
(7) 令 Z I Z_I Z I 为集合 { Z i ∣ i = 1 , … , n } \{ Z_i \mid i = 1, \dots, n \} { Z i ∣ i = 1 , … , n } 中的最小元素。回答索引 I = arg min k = 1 , … , n Z k I = \mathop{\arg\min}\limits_{k=1,\dots,n} Z_k I = k = 1 , … , n arg min Z k 的概率分布 P ( I = i ) \mathbb{P}(I = i) P ( I = i ) ( i = 1 , … , n ) (i = 1, \dots, n) ( i = 1 , … , n ) 。
题目描述
设 n ≥ 2 n\ge2 n ≥ 2 ,随机变量 X 1 , … , X n X_1,\ldots,X_n X 1 , … , X n 相互独立、均取非负实数,且共同密度为
f ( x ) = e − x ( x ≥ 0 ) . f(x)=e^{-x}\qquad(x\ge0). f ( x ) = e − x ( x ≥ 0 ) .
要求给出数学推导并回答:
求 X i X_i X i 的均值和方差。
对给定 k ≠ i k\ne i k = i ,在 X i = x i X_i=x_i X i = x i 条件下求
P ( X i ≤ X k ∣ X i = x i ) . P(X_i\le X_k\mid X_i=x_i). P ( X i ≤ X k ∣ X i = x i ) .
将第 2 问结果乘以 f ( x i ) f(x_i) f ( x i ) 并对 x i x_i x i 积分,求 P ( X i ≤ X k ) P(X_i\le X_k) P ( X i ≤ X k ) 。
令
X min = min k = 1 , … , n X k . X_{\min}=\min_{k=1,\ldots,n}X_k. X m i n = k = 1 , … , n min X k .
对给定正实数 x x x ,求尾概率 P ( x ≤ X min ) P(x\le X_{\min}) P ( x ≤ X m i n ) 。
给定正实数 λ 1 , … , λ n \lambda_1,\ldots,\lambda_n λ 1 , … , λ n ,定义
Z i = X i λ i . Z_i=\frac{X_i}{\lambda_i}. Z i = λ i X i .
在 X i = x i X_i=x_i X i = x i 条件下,对 k ≠ i k\ne i k = i 求 P ( Z i ≤ Z k ∣ X i = x i ) P(Z_i\le Z_k\mid X_i=x_i) P ( Z i ≤ Z k ∣ X i = x i ) 。
在同一条件下,求
P ( ⋂ k = 1 , … , n k ≠ i { Z i ≤ Z k } | X i = x i ) . P\!\left(\bigcap_{\substack{k=1,\ldots,n\\k\ne i}}
\{Z_i\le Z_k\}\,\middle|\,X_i=x_i\right). P k = 1 , … , n k = i ⋂ { Z i ≤ Z k } X i = x i .
令
I = arg min k = 1 , … , n Z k , I=\underset{k=1,\ldots,n}{\arg\min}\,Z_k, I = k = 1 , … , n arg min Z k ,
求索引 I I I 的分布 P ( I = i ) P(I=i) P ( I = i ) (i = 1 , … , n i=1,\ldots,n i = 1 , … , n )。
指数分布 :由密度积分求均值、方差和生存函数,并利用独立性得到多个变量最小值的分布。
条件概率与全概率公式 :先固定 X i = x i X_i=x_i X i = x i 计算比较事件,再对 x i x_i x i 的密度积分消去条件。
指数竞赛与次序统计量 :识别 Z i = X i / λ i Z_i=X_i/\lambda_i Z i = X i / λ i 具有速率 λ i \lambda_i λ i ,计算每个缩放变量成为总体最小值的概率。
Kai
Written by zephyr
题目背景
令 X i X_i X i ( i = 1 , … , n ; n ≥ 2 ) (i = 1, \dots, n; n \geq 2) ( i = 1 , … , n ; n ≥ 2 ) 为 n n n 个独立的非负实值随机变量,它们具有相同的概率密度函数 f ( x ) = e − x f(x) = e^{-x} f ( x ) = e − x 。
1. Compute the mean and variance of X i X_i X i
For the exponential distribution f ( x ) = e − x f(x) = e^{-x} f ( x ) = e − x , we can calculate the mean and variance as follows:
Mean
E [ X i ] = ∫ 0 ∞ x f ( x ) d x = ∫ 0 ∞ x e − x d x = [ − x e − x ] 0 ∞ + ∫ 0 ∞ e − x d x = 0 + [ − e − x ] 0 ∞ = 1 \begin{align*}
E[X_i] &= \int_0^\infty x f(x) dx \\
&= \int_0^\infty x e^{-x} dx \\
&= [-xe^{-x}]_0^\infty + \int_0^\infty e^{-x} dx \\
&= 0 + [-e^{-x}]_0^\infty \\
&= 1
\end{align*} E [ X i ] = ∫ 0 ∞ x f ( x ) d x = ∫ 0 ∞ x e − x d x = [ − x e − x ] 0 ∞ + ∫ 0 ∞ e − x d x = 0 + [ − e − x ] 0 ∞ = 1
Variance
V a r ( X i ) = E [ X i 2 ] − ( E [ X i ] ) 2 = ∫ 0 ∞ x 2 e − x d x − 1 2 = [ − x 2 e − x ] 0 ∞ + 2 ∫ 0 ∞ x e − x d x − 1 = 0 + 2 − 1 = 1 \begin{align*}
Var(X_i) &= E[X_i^2] - (E[X_i])^2 \\
&= \int_0^\infty x^2 e^{-x} dx - 1^2 \\
&= [-x^2e^{-x}]_0^\infty + 2\int_0^\infty xe^{-x} dx - 1 \\
&= 0 + 2 - 1 \\
&= 1
\end{align*} Va r ( X i ) = E [ X i 2 ] − ( E [ X i ] ) 2 = ∫ 0 ∞ x 2 e − x d x − 1 2 = [ − x 2 e − x ] 0 ∞ + 2 ∫ 0 ∞ x e − x d x − 1 = 0 + 2 − 1 = 1
Therefore, E [ X i ] = 1 E[X_i] = 1 E [ X i ] = 1 and V a r ( X i ) = 1 Var(X_i) = 1 Va r ( X i ) = 1 .
2. Compute probability P ( X i ≤ X k ∣ X i = x i ) \mathbb{P}(X_i \leq X_k \mid X_i = x_i) P ( X i ≤ X k ∣ X i = x i )
Given X i = x i X_i = x_i X i = x i , the distribution of X k X_k X k remains f ( x ) = e − x f(x) = e^{-x} f ( x ) = e − x . Thus:
P ( X i ≤ X k ∣ X i = x i ) = P ( x i ≤ X k ) = 1 − P ( X k < x i ) = 1 − ∫ 0 x i e − x d x = 1 − [ − e − x ] 0 x i = 1 − ( 1 − e − x i ) = e − x i \begin{align*}
\mathbb{P}(X_i \leq X_k \mid X_i = x_i) &= \mathbb{P}(x_i \leq X_k) \\
&= 1 - \mathbb{P}(X_k < x_i) \\
&= 1 - \int_0^{x_i} e^{-x} dx \\
&= 1 - [-e^{-x}]_0^{x_i} \\
&= 1 - (1 - e^{-x_i}) \\
&= e^{-x_i}
\end{align*} P ( X i ≤ X k ∣ X i = x i ) = P ( x i ≤ X k ) = 1 − P ( X k < x i ) = 1 − ∫ 0 x i e − x d x = 1 − [ − e − x ] 0 x i = 1 − ( 1 − e − x i ) = e − x i
3. Compute probability P ( X i ≤ X k ) \mathbb{P}(X_i \leq X_k) P ( X i ≤ X k )
We need to integrate over x i x_i x i :
P ( X i ≤ X k ) = ∫ 0 ∞ P ( X i ≤ X k ∣ X i = x ) f ( x ) d x = ∫ 0 ∞ e − x ⋅ e − x d x = ∫ 0 ∞ e − 2 x d x = [ − 1 2 e − 2 x ] 0 ∞ = 1 2 \begin{align*}
\mathbb{P}(X_i \leq X_k) &= \int_0^\infty \mathbb{P}(X_i \leq X_k \mid X_i = x) f(x) dx \\
&= \int_0^\infty e^{-x} \cdot e^{-x} dx \\
&= \int_0^\infty e^{-2x} dx \\
&= [-\frac{1}{2}e^{-2x}]_0^\infty \\
&= \frac{1}{2}
\end{align*} P ( X i ≤ X k ) = ∫ 0 ∞ P ( X i ≤ X k ∣ X i = x ) f ( x ) d x = ∫ 0 ∞ e − x ⋅ e − x d x = ∫ 0 ∞ e − 2 x d x = [ − 2 1 e − 2 x ] 0 ∞ = 2 1
4. Compute probability P ( x ≤ X min ) \mathbb{P}(x \leq X_{\min}) P ( x ≤ X m i n )
The probability that X min X_{\min} X m i n is greater than or equal to x x x is equal to the probability that all X i X_i X i are greater than or equal to x x x :
P ( x ≤ X min ) = P ( ⋂ i = 1 n { x ≤ X i } ) = ∏ i = 1 n P ( x ≤ X i ) (since X i are independent) = ( P ( x ≤ X 1 ) ) n = ( 1 − P ( X 1 < x ) ) n = ( 1 − ∫ 0 x e − t d t ) n = ( 1 − ( 1 − e − x ) ) n = e − n x \begin{align*}
\mathbb{P}(x \leq X_{\min}) &= \mathbb{P}(\bigcap_{i=1}^n \{x \leq X_i\}) \\
&= \prod_{i=1}^n \mathbb{P}(x \leq X_i) \quad \text{(since $X_i$ are independent)} \\
&= (\mathbb{P}(x \leq X_1))^n \\
&= (1 - \mathbb{P}(X_1 < x))^n \\
&= (1 - \int_0^x e^{-t} dt)^n \\
&= (1 - (1 - e^{-x}))^n \\
&= e^{-nx}
\end{align*} P ( x ≤ X m i n ) = P ( i = 1 ⋂ n { x ≤ X i }) = i = 1 ∏ n P ( x ≤ X i ) (since X i are independent) = ( P ( x ≤ X 1 ) ) n = ( 1 − P ( X 1 < x ) ) n = ( 1 − ∫ 0 x e − t d t ) n = ( 1 − ( 1 − e − x ) ) n = e − n x
5. Compute probability P ( Z i ≤ Z k ∣ X i = x i ) \mathbb{P}(Z_i \leq Z_k \mid X_i = x_i) P ( Z i ≤ Z k ∣ X i = x i )
Given X i = x i X_i = x_i X i = x i , we have:
P ( Z i ≤ Z k ∣ X i = x i ) = P ( x i λ i ≤ X k λ k ) = P ( X k ≥ λ k λ i x i ) = ∫ λ k λ i x i ∞ e − x d x = [ − e − x ] λ k λ i x i ∞ = e − λ k λ i x i \begin{align*}
\mathbb{P}(Z_i \leq Z_k \mid X_i = x_i) &= \mathbb{P}(\frac{x_i}{\lambda_i} \leq \frac{X_k}{\lambda_k}) \\
&= \mathbb{P}(X_k \geq \frac{\lambda_k}{\lambda_i}x_i) \\
&= \int_{\frac{\lambda_k}{\lambda_i}x_i}^\infty e^{-x} dx \\
&= [-e^{-x}]_{\frac{\lambda_k}{\lambda_i}x_i}^\infty \\
&= e^{-\frac{\lambda_k}{\lambda_i}x_i}
\end{align*} P ( Z i ≤ Z k ∣ X i = x i ) = P ( λ i x i ≤ λ k X k ) = P ( X k ≥ λ i λ k x i ) = ∫ λ i λ k x i ∞ e − x d x = [ − e − x ] λ i λ k x i ∞ = e − λ i λ k x i
6. Compute probability P ( ⋂ k = 1 , … , n ; k ≠ i { Z i ≤ Z k } ∣ X i = x i ) \mathbb{P}\left(\bigcap_{k=1,\dots,n; k \neq i} \{ Z_i \leq Z_k \} \mid X_i = x_i\right) P ( ⋂ k = 1 , … , n ; k = i { Z i ≤ Z k } ∣ X i = x i )
This probability is the product of all P ( Z i ≤ Z k ∣ X i = x i ) \mathbb{P}(Z_i \leq Z_k \mid X_i = x_i) P ( Z i ≤ Z k ∣ X i = x i ) for k ≠ i k \neq i k = i :
P ( ⋂ k = 1 , … , n ; k ≠ i { Z i ≤ Z k } ∣ X i = x i ) = ∏ k = 1 , … , n ; k ≠ i P ( Z i ≤ Z k ∣ X i = x i ) = ∏ k = 1 , … , n ; k ≠ i e − λ k λ i x i = exp ( − ∑ k = 1 , … , n ; k ≠ i λ k λ i x i ) = exp ( − x i λ i ∑ k = 1 , … , n ; k ≠ i λ k ) \begin{align*}
\mathbb{P}\left(\bigcap_{k=1,\dots,n; k \neq i} \{ Z_i \leq Z_k \} \mid X_i = x_i\right) &= \prod_{k=1,\dots,n; k \neq i} \mathbb{P}(Z_i \leq Z_k \mid X_i = x_i) \\
&= \prod_{k=1,\dots,n; k \neq i} e^{-\frac{\lambda_k}{\lambda_i}x_i} \\
&= \exp\left(-\sum_{k=1,\dots,n; k \neq i} \frac{\lambda_k}{\lambda_i}x_i\right) \\
&= \exp\left(-\frac{x_i}{\lambda_i}\sum_{k=1,\dots,n; k \neq i} \lambda_k\right)
\end{align*} P k = 1 , … , n ; k = i ⋂ { Z i ≤ Z k } ∣ X i = x i = k = 1 , … , n ; k = i ∏ P ( Z i ≤ Z k ∣ X i = x i ) = k = 1 , … , n ; k = i ∏ e − λ i λ k x i = exp − k = 1 , … , n ; k = i ∑ λ i λ k x i = exp − λ i x i k = 1 , … , n ; k = i ∑ λ k
7. Compute the probability distribution P ( I = i ) \mathbb{P}(I = i) P ( I = i )
To calculate P ( I = i ) \mathbb{P}(I = i) P ( I = i ) , we need to integrate over x i x_i x i :
P ( I = i ) = ∫ 0 ∞ P ( ⋂ k = 1 , … , n ; k ≠ i { Z i ≤ Z k } ∣ X i = x i ) f ( x i ) d x i = ∫ 0 ∞ exp ( − x i λ i ∑ k = 1 , … , n ; k ≠ i λ k ) e − x i d x i = ∫ 0 ∞ exp ( − x i ( 1 λ i ∑ k = 1 , … , n ; k ≠ i λ k + 1 ) ) d x i = [ − 1 1 λ i ∑ k = 1 , … , n ; k ≠ i λ k + 1 exp ( − x i ( 1 λ i ∑ k = 1 , … , n ; k ≠ i λ k + 1 ) ) ] 0 ∞ = 1 1 λ i ∑ k = 1 , … , n ; k ≠ i λ k + 1 = λ i ∑ k = 1 n λ k \begin{align*}
\mathbb{P}(I = i) &= \int_0^\infty \mathbb{P}\left(\bigcap_{k=1,\dots,n; k \neq i} \{ Z_i \leq Z_k \} \mid X_i = x_i\right) f(x_i) dx_i \\
&= \int_0^\infty \exp\left(-\frac{x_i}{\lambda_i}\sum_{k=1,\dots,n; k \neq i} \lambda_k\right) e^{-x_i} dx_i \\
&= \int_0^\infty \exp\left(-x_i\left(\frac{1}{\lambda_i}\sum_{k=1,\dots,n; k \neq i} \lambda_k + 1\right)\right) dx_i \\
&= \left[-\frac{1}{\frac{1}{\lambda_i}\sum_{k=1,\dots,n; k \neq i} \lambda_k + 1} \exp\left(-x_i\left(\frac{1}{\lambda_i}\sum_{k=1,\dots,n; k \neq i} \lambda_k + 1\right)\right)\right]_0^\infty \\
&= \frac{1}{\frac{1}{\lambda_i}\sum_{k=1,\dots,n; k \neq i} \lambda_k + 1} \\
&= \frac{\lambda_i}{\sum_{k=1}^n \lambda_k}
\end{align*} P ( I = i ) = ∫ 0 ∞ P k = 1 , … , n ; k = i ⋂ { Z i ≤ Z k } ∣ X i = x i f ( x i ) d x i = ∫ 0 ∞ exp − λ i x i k = 1 , … , n ; k = i ∑ λ k e − x i d x i = ∫ 0 ∞ exp − x i λ i 1 k = 1 , … , n ; k = i ∑ λ k + 1 d x i = − λ i 1 ∑ k = 1 , … , n ; k = i λ k + 1 1 exp − x i λ i 1 k = 1 , … , n ; k = i ∑ λ k + 1 0 ∞ = λ i 1 ∑ k = 1 , … , n ; k = i λ k + 1 1 = ∑ k = 1 n λ k λ i
Therefore, P ( I = i ) = λ i ∑ k = 1 n λ k \mathbb{P}(I = i) = \frac{\lambda_i}{\sum_{k=1}^n \lambda_k} P ( I = i ) = ∑ k = 1 n λ k λ i .
Knowledge
难点思路
这道题的难点在于处理条件概率和多个随机变量的最小值。特别是在第 6 和第 7 问中,需要仔细处理多个条件的交集概率。
解题技巧和信息
对于指数分布,要熟悉其基本性质,如均值、方差、累积分布函数等。
在处理多个独立随机变量时,要善于利用独立性质简化计算。
在计算条件概率时,要清楚地区分给定条件下的随机变量和非随机变量。
在处理最小值问题时,可以转化为所有变量都大于某个值的概率。
在积分计算中,要善于使用指数函数的性质,如 ∫ e a x d x = 1 a e a x + C \int e^{ax} dx = \frac{1}{a}e^{ax} + C ∫ e a x d x = a 1 e a x + C 。
重点词汇
Probability density function: 概率密度函数
Exponential distribution: 指数分布
Conditional probability: 条件概率
Minimum value: 最小值
Order statistics: 顺序统计量
Independent random variables: 独立随机变量