跳到主要内容

大阪大学 基礎工学研究科 数理科学 (システム創成専攻) 2019年8月実施 数理科学 II [7]

Author​

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

Description​

n,mn,m は2以上の整数、μX,μY\mu_X,\mu_Y は実数とし、X1,…,XnX_1,\ldots,X_n は N(μX,σ2)N(\mu_X,\sigma^2)、Y1,…,YmY_1,\ldots,Y_m は N(μY,σ2)N(\mu_Y,\sigma^2) に従い、全変数は独立とする。σ2>0\sigma^2>0 とし、

X‾=1n∑Xi,Y‾=1m∑Yj,S2=1n∑(Xi−X‾)2,T2=1m∑(Yj−Y‾)2\overline X=\frac1n\sum X_i,\quad\overline Y=\frac1m\sum Y_j,\quad S^2=\frac1n\sum(X_i-\overline X)^2,\quad T^2=\frac1m\sum(Y_j-\overline Y)^2

とおく。帰無仮説 μX=μY\mu_X=\mu_Y の下で次を示せ。

(1) (X‾−Y‾)/σ2/n+σ2/m(\overline X-\overline Y)/\sqrt{\sigma^2/n+\sigma^2/m} は標準正規分布に従う。

(2) Cov⁡(Xi−X‾,X‾)\operatorname{Cov}(X_i-\overline X,\overline X) を求めよ。

(3) X‾−Y‾\overline X-\overline Y と nS2+mT2nS^2+mT^2 は独立である。

(4) (n+m−2)nm(n+m)(nS2+mT2)(X‾−Y‾)\displaystyle\sqrt{\frac{(n+m-2)nm}{(n+m)(nS^2+mT^2)}}(\overline X-\overline Y) は自由度 n+m−2n+m-2 の tt 分布に従う。

Kai​

(1)​

独立な正規変数の線形結合より

X‾−Y‾∼N(0,σ2(1n+1m)).\overline X-\overline Y\sim N\left(0,\sigma^2\left(\frac1n+\frac1m\right)\right).

標準化すればよい。

(2)​

Cov⁡(Xi−X‾,X‾)=σ2n−σ2n=0.\boxed{\operatorname{Cov}(X_i-\overline X,\overline X)=\frac{\sigma^2}{n}-\frac{\sigma^2}{n}=0}.

(3)​

X‾\overline X と残差ベクトル (Xi−X‾)i(X_i-\overline X)_i は同時正規で無相関なので独立。同様に Y‾\overline Y と (Yj−Y‾)j(Y_j-\overline Y)_j も独立で、さらに両標本が独立だから、標本平均の差と残差平方和は独立である。

(4)​

標準正規ベクトルを平均方向とその直交補空間へ直交変換すると

nS2σ2∼χn−12,mT2σ2∼χm−12.\frac{nS^2}{\sigma^2}\sim\chi^2_{n-1},\qquad\frac{mT^2}{\sigma^2}\sim\chi^2_{m-1}.

両者は独立なので U=(nS2+mT2)/σ2∼χn+m−22U=(nS^2+mT^2)/\sigma^2\sim\chi^2_{n+m-2}。 (1) の標準正規変数を ZZ とすれば (3) により Z⊥UZ\perp U。指定の統計量は Z/U/(n+m−2)Z/\sqrt{U/(n+m-2)} であるから結論を得る。