京都大学 情報学研究科 システム科学専攻 2006年8月実施 確率統計
Author
犬 (finalized by 祭音Myyura with assistance from GPT 6 Astra)
Description
(X1,X2) の結合密度を f(x1,x2)=8x1x2(0<x1<x2<1)、その他では 0 とする。
(i) 各周辺密度を求め、独立かどうか述べよ。
(ii) Y1=X1/X2,Y2=X2 とするとき、結合密度と各周辺密度を求め、独立かどうか述べよ。
自由度 k の χ2 変数 X の密度は
f(x)=2k/2Γ(k/2)xk/2−1e−x/2(x>0),Γ(α)=∫0∞e−ttα−1dt,Γ(α+1)=αΓ(α).
(i) E[X]=k,Var(X)=2k を示せ。
(ii) 正規母集団 N(μx,σ2) からの独立標本 x1,…,xn に対し、xˉ=n−1∑ixi,sx2=(n−1)−1∑i(xi−xˉ)2 と定義する。この定義から E[sx2]=σ2 を示せ。
(iii) sx2=σ2Z/(n−1) とすると Z∼χn−12。別の独立な正規母集団 N(μy,σ2) から独立標本 y1,…,ym をとり、同様に sy2 を定める。sy2=σ2W/(m−1)、W∼χm−12 である。共通の分散の不偏推定量 sp2=λsx2+(1−λ)sy2 の分散を最小にする λ を求めよ。
题目描述
(X1,X2) 在 0<x1<x2<1 上的密度为 8x1x2,其余为 0。
(i) 求边缘密度并判断是否独立。
(ii) 令 Y1=X1/X2,Y2=X2,求联合密度、各边缘密度并判断是否独立。
自由度 k 的卡方变量密度为 xk/2−1e−x/2/[2k/2Γ(k/2)](x>0),可用 Γ(α)=∫0∞e−ttα−1dt 及 Γ(α+1)=αΓ(α)。
(i) 证明 E[X]=k、Var(X)=2k。
(ii) 从 N(μx,σ2) 取 n 个独立样本,定义 xˉ=n−1∑ixi、sx2=(n−1)−1∑i(xi−xˉ)2。从定义证明 E[sx2]=σ2。
(iii) 已知 (n−1)sx2/σ2∼χn−12。从另一个独立总体 N(μy,σ2) 取 m 个独立样本,同样定义 sy2,且 (m−1)sy2/σ2∼χm−12。求使无偏估计 sp2=λsx2+(1−λ)sy2 方差最小的 λ。
Kai
(i) 0<x<1 において
fX1(x)=∫x18xx2dx2=4x(1−x2),fX2(x)=∫0x8x1xdx1=4x3.
その他では両者とも 0。周辺密度の積は三角形の外でも正になり、結合密度と一致しないので独立ではない。
(ii) 逆変換は x1=y1y2,x2=y2、ヤコビアンの絶対値は y2。よって
fY1,Y2(y1,y2)=8y1y23(0<y1,y2<1).
周辺密度は fY1(y)=2y、fY2(y)=4y3(0<y<1、その他は 0)。結合密度はその積に等しく、独立である。
(i) ガンマ積分から E[Xr]=2rΓ(k/2+r)/Γ(k/2)。r=1,2 を代入すると E[X]=k、E[X2]=k(k+2)、したがって Var(X)=2k。
(ii) 恒等式
i∑(xi−xˉ)2=i∑(xi−μx)2−n(xˉ−μx)2
の期待値をとると、右辺は nσ2−n(σ2/n)=(n−1)σ2。よって E[sx2]=σ2。
(iii) (i) から Var(sx2)=2σ4/(n−1)、Var(sy2)=2σ4/(m−1)。標本間の独立性より
Var(sp2)=2σ4(n−1λ2+m−1(1−λ)2).
これは下に凸の二次式であり、微分して λ/(n−1)=(1−λ)/(m−1)。したがって
λ=n+m−2n−1.