跳到主要内容

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

Author

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

Description

問1

(X1,X2)(X_1,X_2) の結合密度を f(x1,x2)=8x1x2f(x_1,x_2)=8x_1x_20<x1<x2<10<x_1<x_2<1)、その他では 00 とする。

(i) 各周辺密度を求め、独立かどうか述べよ。

(ii) Y1=X1/X2,Y2=X2Y_1=X_1/X_2,Y_2=X_2 とするとき、結合密度と各周辺密度を求め、独立かどうか述べよ。

問2

自由度 kkχ2\chi^2 変数 XX の密度は

f(x)=xk/21ex/22k/2Γ(k/2)(x>0),Γ(α)=0ettα1dt,Γ(α+1)=αΓ(α).f(x)=\frac{x^{k/2-1}e^{-x/2}}{2^{k/2}\Gamma(k/2)}\quad(x>0),\qquad \Gamma(\alpha)=\int_0^\infty e^{-t}t^{\alpha-1}\,dt, \quad\Gamma(\alpha+1)=\alpha\Gamma(\alpha).

(i) E[X]=k,Var(X)=2kE[X]=k,\operatorname{Var}(X)=2k を示せ。

(ii) 正規母集団 N(μx,σ2)N(\mu_x,\sigma^2) からの独立標本 x1,,xnx_1,\ldots,x_n に対し、xˉ=n1ixi,sx2=(n1)1i(xixˉ)2\bar x=n^{-1}\sum_ix_i,s_x^2=(n-1)^{-1}\sum_i(x_i-\bar x)^2 と定義する。この定義から E[sx2]=σ2E[s_x^2]=\sigma^2 を示せ。

(iii) sx2=σ2Z/(n1)s_x^2=\sigma^2Z/(n-1) とすると Zχn12Z\sim\chi^2_{n-1}。別の独立な正規母集団 N(μy,σ2)N(\mu_y,\sigma^2) から独立標本 y1,,ymy_1,\ldots,y_m をとり、同様に sy2s_y^2 を定める。sy2=σ2W/(m1)s_y^2=\sigma^2W/(m-1)Wχm12W\sim\chi^2_{m-1} である。共通の分散の不偏推定量 sp2=λsx2+(1λ)sy2s_p^2=\lambda s_x^2+(1-\lambda)s_y^2 の分散を最小にする λ\lambda を求めよ。

题目描述

问1

(X1,X2)(X_1,X_2)0<x1<x2<10<x_1<x_2<1 上的密度为 8x1x28x_1x_2,其余为 00

(i) 求边缘密度并判断是否独立。

(ii) 令 Y1=X1/X2,Y2=X2Y_1=X_1/X_2,Y_2=X_2,求联合密度、各边缘密度并判断是否独立。

问2

自由度 kk 的卡方变量密度为 xk/21ex/2/[2k/2Γ(k/2)]x^{k/2-1}e^{-x/2}/[2^{k/2}\Gamma(k/2)]x>0x>0),可用 Γ(α)=0ettα1dt\Gamma(\alpha)=\int_0^\infty e^{-t}t^{\alpha-1}\,dtΓ(α+1)=αΓ(α)\Gamma(\alpha+1)=\alpha\Gamma(\alpha)

(i) 证明 E[X]=kE[X]=kVar(X)=2k\operatorname{Var}(X)=2k

(ii) 从 N(μx,σ2)N(\mu_x,\sigma^2)nn 个独立样本,定义 xˉ=n1ixi\bar x=n^{-1}\sum_ix_isx2=(n1)1i(xixˉ)2s_x^2=(n-1)^{-1}\sum_i(x_i-\bar x)^2。从定义证明 E[sx2]=σ2E[s_x^2]=\sigma^2

(iii) 已知 (n1)sx2/σ2χn12(n-1)s_x^2/\sigma^2\sim\chi^2_{n-1}。从另一个独立总体 N(μy,σ2)N(\mu_y,\sigma^2)mm 个独立样本,同样定义 sy2s_y^2,且 (m1)sy2/σ2χm12(m-1)s_y^2/\sigma^2\sim\chi^2_{m-1}。求使无偏估计 sp2=λsx2+(1λ)sy2s_p^2=\lambda s_x^2+(1-\lambda)s_y^2 方差最小的 λ\lambda

Kai

問1

(i) 0<x<10<x<1 において

fX1(x)=x18xx2dx2=4x(1x2),fX2(x)=0x8x1xdx1=4x3.f_{X_1}(x)=\int_x^18xx_2\,dx_2=4x(1-x^2),\qquad f_{X_2}(x)=\int_0^x8x_1x\,dx_1=4x^3.

その他では両者とも 00。周辺密度の積は三角形の外でも正になり、結合密度と一致しないので独立ではない。

(ii) 逆変換は x1=y1y2,x2=y2x_1=y_1y_2,x_2=y_2、ヤコビアンの絶対値は y2y_2。よって

fY1,Y2(y1,y2)=8y1y23(0<y1,y2<1).f_{Y_1,Y_2}(y_1,y_2)=8y_1y_2^3\quad(0<y_1,y_2<1).

周辺密度は fY1(y)=2yf_{Y_1}(y)=2yfY2(y)=4y3f_{Y_2}(y)=4y^30<y<10<y<1、その他は 00)。結合密度はその積に等しく、独立である。

問2

(i) ガンマ積分から E[Xr]=2rΓ(k/2+r)/Γ(k/2)E[X^r]=2^r\Gamma(k/2+r)/\Gamma(k/2)r=1,2r=1,2 を代入すると E[X]=kE[X]=kE[X2]=k(k+2)E[X^2]=k(k+2)、したがって Var(X)=2k\operatorname{Var}(X)=2k

(ii) 恒等式

i(xixˉ)2=i(xiμx)2n(xˉμx)2\sum_i(x_i-\bar x)^2=\sum_i(x_i-\mu_x)^2-n(\bar x-\mu_x)^2

の期待値をとると、右辺は nσ2n(σ2/n)=(n1)σ2n\sigma^2-n(\sigma^2/n)=(n-1)\sigma^2。よって E[sx2]=σ2E[s_x^2]=\sigma^2

(iii) (i) から Var(sx2)=2σ4/(n1)\operatorname{Var}(s_x^2)=2\sigma^4/(n-1)Var(sy2)=2σ4/(m1)\operatorname{Var}(s_y^2)=2\sigma^4/(m-1)。標本間の独立性より

Var(sp2)=2σ4(λ2n1+(1λ)2m1).\operatorname{Var}(s_p^2)=2\sigma^4\left(\frac{\lambda^2}{n-1}+\frac{(1-\lambda)^2}{m-1}\right).

これは下に凸の二次式であり、微分して λ/(n1)=(1λ)/(m1)\lambda/(n-1)=(1-\lambda)/(m-1)。したがって

λ=n1n+m2.\boxed{\lambda=\frac{n-1}{n+m-2}}.