跳到主要内容

東京大学 情報理工学系研究科 コンピュータ科学専攻 2017年8月実施 専門科目II 問題6

Author

祭音Myyura (co-authored with GPT 5.6 SOL)

Description

RNA 分子的分解时间 TT 的概率密度为

fT(t)=λeλt,t0,f_T(t)=\lambda e^{-\lambda t},\qquad t\ge0,

其中 λ>0\lambda>0

(1)求累积分布函数 FT(t)F_T(t)TT 的中位数。

(2)独立同分布地测量 nn 个 RNA 分子的分解时间 T1,,TnT_1,\ldots,T_n。求样本均值 μT=1niTi\mu_T=\frac1n\sum_iT_i 的期望和方差。

(3)令 Tmax=max(T1,,Tn)T_{\max}=\max(T_1,\ldots,T_n)。用 FT(t)F_T(t) 表示 Pr(Tmax>t)\Pr(T_{\max}>t)

(4)求 TmaxT_{\max} 的概率密度和期望。

Kai

(1)

t0t\ge0

FT(t)=0tλeλxdx=1eλt.F_T(t)=\int_0^t\lambda e^{-\lambda x}\,dx=1-e^{-\lambda t}.

FT(m)=1/2F_T(m)=1/2,得中位数

m=ln2λ.m=\frac{\ln 2}{\lambda}.

(2)

指数分布满足 E[Ti]=1/λ\mathbb E[T_i]=1/\lambdaVar(Ti)=1/λ2\operatorname{Var}(T_i)=1/\lambda^2。由独立性,

E[μT]=1λ,Var(μT)=1nλ2.\mathbb E[\mu_T]=\frac1\lambda, \qquad \operatorname{Var}(\mu_T)=\frac1{n\lambda^2}.

(3)

t0t\ge0

Pr(Tmax>t)=1Pr(T1t,,Tnt)=1[FT(t)]n.\Pr(T_{\max}>t) =1-\Pr(T_1\le t,\ldots,T_n\le t) =1-[F_T(t)]^n.

(4)

FTmax(t)=[FT(t)]nF_{T_{\max}}(t)=[F_T(t)]^n

fTmax(t)=n[FT(t)]n1fT(t)=nλeλt(1eλt)n1,t0.f_{T_{\max}}(t) =n[F_T(t)]^{n-1}f_T(t) =n\lambda e^{-\lambda t}(1-e^{-\lambda t})^{n-1}, \qquad t\ge0.

再用尾积分公式,

E[Tmax]=0{1(1eλt)n}dt=1λk=1n1k.\begin{aligned} \mathbb E[T_{\max}] &=\int_0^\infty\{1-(1-e^{-\lambda t})^n\}\,dt\\ &=\frac1\lambda\sum_{k=1}^n\frac1k. \end{aligned}