跳到主要内容

電気通信大学 情報理工学研究科 情報・ネットワーク工学専攻 2023年8月実施 選択問題 確率統計

Author​

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

Description​

確率変数 XX の密度関数を

f(x)=1βe−x/β(x≥0, β>0)f(x)=\frac1\beta e^{-x/\beta}\quad(x\geq0,\ \beta>0)

とする。積率母関数、β0>0\beta_0>0 に対する P(X>β0)P(X>\beta_0)、Z=X1{X≥β0}Z=X\mathbf1_{\{X\ge\beta_0\}} の期待値、および Y=XY=\sqrt X の期待値と分散を求めよ。また、独立同分布な標本 x1,…,xnx_1,\ldots,x_n から β\beta を最尤推定し、W=max⁡(X1,…,Xn)W=\max(X_1,\ldots,X_n) の W>β0W>\beta_0 の下での条件付き密度を求めよ。

题目描述​

给定尺度参数为 β\beta 的指数分布,求矩生成函数、尾概率、截断变量和平方根变换的矩;再求 β\beta 的极大似然估计,以及样本最大值在超过阈值条件下的密度。

Kai​

(1)​

φ(θ)=∫0∞eθx1βe−x/β dx=11−βθ(θ<1β).\begin{aligned} \varphi(\theta) &=\int_0^\infty e^{\theta x}\frac1\beta e^{-x/\beta}\,dx\\ &=\boxed{\frac1{1-\beta\theta}} \qquad\left(\theta<\frac1\beta\right). \end{aligned}

(2)​

P(X>β0)=∫β0∞1βe−x/β dx=e−β0/β.\boxed{P(X>\beta_0)=\int_{\beta_0}^{\infty}\frac1\beta e^{-x/\beta}\,dx =e^{-\beta_0/\beta}}.

(3)​

E[Z]=∫β0∞x1βe−x/β dx=(β0+β)e−β0/β.\begin{aligned} E[Z] &=\int_{\beta_0}^{\infty}x\frac1\beta e^{-x/\beta}\,dx\\ &=\boxed{(\beta_0+\beta)e^{-\beta_0/\beta}}. \end{aligned}

(4)​

E[Y]=∫0∞x1βe−x/β dx=πβ2.E[Y] =\int_0^\infty\sqrt{x}\frac1\beta e^{-x/\beta}\,dx =\boxed{\frac{\sqrt{\pi\beta}}2}.

また E[Y2]=E[X]=βE[Y^2]=E[X]=\beta であるから、

V[Y]=β−πβ4=4−π4β.\boxed{V[Y]=\beta-\frac{\pi\beta}{4} =\frac{4-\pi}{4}\beta}.

(5)​

尤度は

L(β)=β−nexp⁡(−1β∑i=1nxi).L(\beta)=\beta^{-n} \exp\left(-\frac1\beta\sum_{i=1}^n x_i\right).

∑ixi>0\sum_i x_i>0 のとき、∂log⁡L/∂β=0\partial\log L/\partial\beta=0 より、

β^=1n∑i=1nxi.\boxed{\widehat\beta=\frac1n\sum_{i=1}^n x_i}.

すべての xix_i が 00 なら、β↓0\beta\downarrow0 で尤度は発散するため、β>0\beta>0 の範囲に最尤推定値は存在しない。

(6)​

w≥0w\geq0 に対して

FW(w)=P(W≤w)=(1−e−w/β)n.F_W(w)=P(W\leq w)=\left(1-e^{-w/\beta}\right)^n.

したがって w≥β0w\geq\beta_0 では、

g(w∣W>β0)=nβe−w/β(1−e−w/β)n−11−(1−e−β0/β)n\boxed{ g(w\mid W>\beta_0)= \frac{\dfrac n\beta e^{-w/\beta} \left(1-e^{-w/\beta}\right)^{n-1}} {1-\left(1-e^{-\beta_0/\beta}\right)^n}}

であり、w<β0w<\beta_0 では g(w∣W>β0)=0g(w\mid W>\beta_0)=0 である。