跳到主要内容

京都大学 情報学研究科 システム科学専攻 2009年8月実施 確率統計

Author

犬 (finalized by 祭音Myyura with assistance from GPT 6 Astra)

Description

問1

(X1,X2)(X_1,X_2) の結合確率密度が x1>0,x2>0x_1>0,x_2>0f(x1,x2)=14e(x1+x2)/2f(x_1,x_2)=\frac14e^{-(x_1+x_2)/2}、その他で 00 である。Y=(X1X2)/2Y=(X_1-X_2)/2 の確率密度を求めよ。

問2

X1,,XnX_1,\ldots,X_n は独立に、母数 θ>0\theta>0 の指数分布 f(x)=θ1ex/θf(x)=\theta^{-1}e^{-x/\theta}x>0x>0、その他は 00)に従う。

(i) θ^=n1i=1nXi\widehat\theta=n^{-1}\sum_{i=1}^nX_i が不偏、すなわち E[θ^]=θE[\widehat\theta]=\theta であることを示せ。

(ii) E[etXi]E[e^{tX_i}] を求め、Y=(2/θ)i=1nXiY=(2/\theta)\sum_{i=1}^nX_i のモーメント母関数を求めよ。tt は期待値の存在する範囲とする。

(iii) 自由度 kkχ2\chi^2 分布に従う ZZZ=i=1kWi2Z=\sum_{i=1}^kW_i^2 と表される。WiW_i は互いに独立な平均 00、分散 11 の正規変数である。E[etWi2]E[e^{tW_i^2}] および E[etZ]E[e^{tZ}] を求めよ。

(iv) YYχ2\chi^2 分布に従うことを示し、自由度を求めよ。

(v) Pr(Z>χα,k2)=α\Pr(Z>\chi^2_{\alpha,k})=\alpha と定義する。このとき

Pr(χ1α/2,k2<Z<χα/2,k2)=1α.\Pr\bigl(\chi^2_{1-\alpha/2,k}<Z<\chi^2_{\alpha/2,k}\bigr)=1-\alpha.

これを用い、θ^\widehat\theta に基づく θ\theta の信頼係数 1α1-\alpha の信頼区間を示せ。

题目描述

问1

(X1,X2)(X_1,X_2) 的联合密度在 x1,x2>0x_1,x_2>0 时为 14e(x1+x2)/2\frac14e^{-(x_1+x_2)/2},其余为 00。求 Y=(X1X2)/2Y=(X_1-X_2)/2 的密度。

问2

X1,,XnX_1,\ldots,X_n 独立服从尺度参数 θ>0\theta>0 的指数分布,密度为 θ1ex/θ\theta^{-1}e^{-x/\theta}x>0x>0,其他为 00)。

(i) 证明 θ^=n1iXi\widehat\theta=n^{-1}\sum_iX_i 无偏。

(ii) 求 XiX_iY=(2/θ)iXiY=(2/\theta)\sum_iX_i 的矩母函数,tt 限于相应期望存在的范围。

(iii) 若 Z=i=1kWi2Z=\sum_{i=1}^kW_i^2,其中 WiW_i 独立服从标准正态分布,求 Wi2W_i^2ZZ 的矩母函数。

(iv) 证明 YY 服从卡方分布并求自由度。

(v) 定义上尾分位数 Pr(Z>χα,k2)=α\Pr(Z>\chi^2_{\alpha,k})=\alpha。利用 Pr(χ1α/2,k2<Z<χα/2,k2)=1α\Pr(\chi^2_{1-\alpha/2,k}<Z<\chi^2_{\alpha/2,k})=1-\alpha,给出基于 θ^\widehat\thetaθ\theta1α1-\alpha 置信区间。

Kai

問1

Z=X1Z=X_1 とおけば X1=Z,X2=Z2YX_1=Z,X_2=Z-2Y、ヤコビアンの絶対値は 22。したがって

fY(y)=max(0,2y)12eyzdz=12ey(yR).f_Y(y)=\int_{\max(0,2y)}^\infty\frac12e^{y-z}\,dz =\boxed{\frac12e^{-|y|}}\qquad(y\in\mathbb R).

問2

(i) E[Xi]=0(x/θ)ex/θdx=θE[X_i]=\int_0^\infty(x/\theta)e^{-x/\theta}\,dx=\theta より、E[θ^]=θE[\widehat\theta]=\theta

(ii) 直接積分と独立性から

MXi(t)=11θt(t<1/θ),MY(t)=i=1nMXi(2t/θ)=(12t)n(t<1/2).M_{X_i}(t)=\frac1{1-\theta t}\quad(t<1/\theta),\qquad M_Y(t)=\prod_{i=1}^nM_{X_i}(2t/\theta)=(1-2t)^{-n}\quad(t<1/2).

(iii) ガウス積分より

E[etWi2]=12πe(12t)w2/2dw=(12t)1/2.E[e^{tW_i^2}]=\frac1{\sqrt{2\pi}}\int_{-\infty}^\infty e^{-(1-2t)w^2/2}\,dw=(1-2t)^{-1/2}.

よって MZ(t)=(12t)k/2M_Z(t)=(1-2t)^{-k/2}t<1/2t<1/2)。

(iv) モーメント母関数の一致より Yχ2n2Y\sim\chi^2_{2n}

(v) Y=2nθ^/θY=2n\widehat\theta/\theta を分位点の間にはさみ、正数の逆数をとると

(2nθ^χα/2,2n2, 2nθ^χ1α/2,2n2).\boxed{\left(\frac{2n\widehat\theta}{\chi^2_{\alpha/2,2n}},\ \frac{2n\widehat\theta}{\chi^2_{1-\alpha/2,2n}}\right)}.