第 2 次作业:正态抽样分布与 t、F 统计量

题目按源文件 week3.pdf 顺序整理。源文件第 6 题的“最小常数”与不等式方向不相容,解答中如实指出,并给出有实际意义的边界值。

题 1:两个正态样本均值

从 N(100,4)N(100,4) 总体分别独立抽取容量为 1515 和 2020 的样本,样本均值为 Xˉ,Yˉ\bar X,\bar Y。求

P(∣Xˉ−Yˉ∣≥0.2).P(|\bar X-\bar Y|\ge0.2).
展开解答 Xˉ−Yˉ∼N(0,415+420)=N(0,715).\bar X-\bar Y\sim N\left(0,\frac4{15}+\frac4{20}\right) =N\left(0,\frac7{15}\right).

所以

P(∣Xˉ−Yˉ∣≥0.2)=2[1−Φ(0.27/15)].P(|\bar X-\bar Y|\ge0.2) =2\left[ 1-\Phi\left(\frac{0.2}{\sqrt{7/15}}\right) \right].

因标准化值约为 0.2930.293,概率约为

0.770.\boxed{0.770}.

题 2:关于常数 cc 的概率界

设 X1,…,XnX_1,\ldots,X_n 来自 N(μ,1)N(\mu,1)。确定最小的常数 cc,使对任意 μ≥0\mu\ge0,有

P(∣Xˉ∣<c)≤α.P(|\bar X|<c)\le\alpha.
展开解答

对 c>0c>0,

g(μ)=Pμ(∣Xˉ∣<c)=Φ(n(c−μ))−Φ(n(−c−μ)).g(\mu)=P_\mu(|\bar X|<c) =\Phi(\sqrt n(c-\mu)) -\Phi(\sqrt n(-c-\mu)).

在 μ≥0\mu\ge0 上,g(μ)g(\mu) 的最大值出现在 μ=0\mu=0,所以条件等价于

2Φ(n c)−1≤α.2\Phi(\sqrt n\,c)-1\le\alpha.

即

0<c≤Φ−1((1+α)/2)n.0<c\le \frac{\Phi^{-1}((1+\alpha)/2)}{\sqrt n}.

因此原题若坚持问“最小的正数”,并不存在最小值,只有下确界 00;若本意是问最大的允许值,则

cmax⁡=Φ−1((1+α)/2)n.\boxed{ c_{\max}= \frac{\Phi^{-1}((1+\alpha)/2)}{\sqrt n} }.

我当时的手写作答也得到 c≤⋯c\le\cdots,这印证了题面文字存在问题。

题 3:F(n,n)F(n,n) 在 1 两侧的概率

设随机变量 X∼F(n,n)X\sim F(n,n),证明 P(X<1)=0.5P(X<1)=0.5。

展开解答

可写

X=U/nV/n=UV,X=\frac{U/n}{V/n}=\frac UV,

其中 U,VU,V 独立且同服从 χ2(n)\chi^2(n)。交换 U,VU,V 不改变联合分布,因此

1X=VU=dX.\frac1X=\frac VU\overset d=X.

于是

P(X<1)=P(1/X<1)=P(X>1).P(X<1)=P(1/X<1)=P(X>1).

FF 分布连续,P(X=1)=0P(X=1)=0,故

P(X<1)=P(X>1)=12.\boxed{P(X<1)=P(X>1)=\frac12}.

题 4:两个正态线性组合的比值

X1,X2X_1,X_2 来自 N(0,σ2)N(0,\sigma^2),求

Y=(X1+X2X1−X2)2Y=\left(\frac{X_1+X_2}{X_1-X_2}\right)^2

的分布。

展开解答

令

U=X1+X22σ,V=X1−X22σ.U=\frac{X_1+X_2}{\sqrt2\sigma},\qquad V=\frac{X_1-X_2}{\sqrt2\sigma}.

二者均为标准正态,且协方差为零。由于联合正态,所以 U,VU,V 独立。因此

Y=U2V2=χ2(1)/1χ2(1)/1∼F(1,1).Y=\frac{U^2}{V^2} =\frac{\chi^2(1)/1}{\chi^2(1)/1} \sim\boxed{F(1,1)}.

题 5:求临界常数 kk

X1,X2X_1,X_2 来自 N(0,1)N(0,1)。求 kk,使

P{(X1+X2)2(X1−X2)2+(X1+X2)2>k}=0.05.P\left\{ \frac{(X_1+X_2)^2} {(X_1-X_2)^2+(X_1+X_2)^2}>k \right\}=0.05.
展开解答

令 U=(X1+X2)/2U=(X_1+X_2)/\sqrt2、V=(X1−X2)/2V=(X_1-X_2)/\sqrt2,则 U,VU,V 独立标准正态。题中随机变量为

W=U2U2+V2∼Beta⁡(12,12).W=\frac{U^2}{U^2+V^2} \sim\operatorname{Beta}\left(\frac12,\frac12\right).

也可令 F=U2/V2∼F(1,1)F=U^2/V^2\sim F(1,1),则 W=F/(1+F)W=F/(1+F)。要求 P(W>k)=0.05P(W>k)=0.05,所以

k=F0.95(1,1)1+F0.95(1,1)≈0.9938.k=\frac{F_{0.95}(1,1)} {1+F_{0.95}(1,1)} \approx\boxed{0.9938}.

题 6:任意线性组合构造 t 统计量

X1,…,XnX_1,\ldots,X_n 来自 N(μ1,σ2)N(\mu_1,\sigma^2),Y1,…,YmY_1,\ldots,Y_m 来自 N(μ2,σ2)N(\mu_2,\sigma^2),两样本独立。c,dc,d 是任意两个不全为零的常数。证明

t=c(Xˉ−μ1)+d(Yˉ−μ2)Swc2/n+d2/m∼t(n+m−2),t= \frac{c(\bar X-\mu_1)+d(\bar Y-\mu_2)} {S_w\sqrt{c^2/n+d^2/m}} \sim t(n+m-2),

其中

Sw2=(n−1)SX2+(m−1)SY2n+m−2.S_w^2= \frac{(n-1)S_X^2+(m-1)S_Y^2} {n+m-2}.
展开解答

由独立正态样本,

Z=c(Xˉ−μ1)+d(Yˉ−μ2)σc2/n+d2/m∼N(0,1).Z= \frac{c(\bar X-\mu_1)+d(\bar Y-\mu_2)} {\sigma\sqrt{c^2/n+d^2/m}} \sim N(0,1).

又

(n−1)SX2σ2∼χ2(n−1),(m−1)SY2σ2∼χ2(m−1),\frac{(n-1)S_X^2}{\sigma^2}\sim\chi^2(n-1), \qquad \frac{(m-1)S_Y^2}{\sigma^2}\sim\chi^2(m-1),

二者独立,因此

(n+m−2)Sw2σ2∼χ2(n+m−2).\frac{(n+m-2)S_w^2}{\sigma^2} \sim\chi^2(n+m-2).

正态样本中均值与样本方差独立,所以 ZZ 与 Sw2S_w^2 独立。按 tt 分布定义,

Z[(n+m−2)Sw2/σ2]/(n+m−2)∼t(n+m−2),\frac{Z} {\sqrt{[(n+m-2)S_w^2/\sigma^2]/(n+m-2)}} \sim t(n+m-2),

化简即得证。

题 7:成对和的离差平方和

X1,…,X2nX_1,\ldots,X_{2n} 来自 N(μ,σ2)N(\mu,\sigma^2),Xˉ\bar X 是全部 2n2n 个观测的样本均值。令

Γ=∑i=1n(Xi+Xn+i−2Xˉ)2.\Gamma= \sum_{i=1}^n (X_i+X_{n+i}-2\bar X)^2.

求 E(Γ)E(\Gamma)。

展开解答

令 Ai=Xi+Xn+iA_i=X_i+X_{n+i},则

Ai∼i.i.d.N(2μ,2σ2),Aˉ=2Xˉ.A_i\overset{\text{i.i.d.}}\sim N(2\mu,2\sigma^2), \qquad \bar A=2\bar X.

所以

Γ=∑i=1n(Ai−Aˉ)2=(n−1)SA2.\Gamma=\sum_{i=1}^n(A_i-\bar A)^2 =(n-1)S_A^2.

由样本方差无偏,

E(Γ)=2(n−1)σ2.\boxed{E(\Gamma)=2(n-1)\sigma^2}.

题 8:用样本方差控制概率

X1,…,XnX_1,\ldots,X_n 来自正态总体 N(μ,σ2)N(\mu,\sigma^2)。求最小 nn,使

P(S2σ2≤1.5)≥0.95.P\left(\frac{S^2}{\sigma^2}\le1.5\right)\ge0.95.
展开解答

由

(n−1)S2σ2∼χ2(n−1),\frac{(n-1)S^2}{\sigma^2}\sim\chi^2(n-1),

条件等价于

P{χ2(n−1)≤1.5(n−1)}≥0.95.P\{\chi^2(n-1)\le1.5(n-1)\}\ge0.95.

逐个自由度查分布函数:

n=26:P{χ2(25)≤37.5}≈0.9483<0.95,n=27:P{χ2(26)≤39}≈0.9512≥0.95.\begin{aligned} n=26:&\quad P\{\chi^2(25)\le37.5\}\approx0.9483<0.95,\\ n=27:&\quad P\{\chi^2(26)\le39\}\approx0.9512\ge0.95. \end{aligned}

所以

nmin⁡=27.\boxed{n_{\min}=27}.

题 9:内部学生化残差与 t 分布

X1,…,XnX_1,\ldots,X_n 独立同分布于 N(μ,σ2)N(\mu,\sigma^2),令

Xˉ=1n∑Xi,S2=1n−1∑(Xi−Xˉ)2,\bar X=\frac1n\sum X_i,\qquad S^2=\frac1{n-1}\sum(X_i-\bar X)^2, ξ=X1−XˉS.\xi=\frac{X_1-\bar X}{S}.

找出 ξ\xi 与 tt 分布的联系。

展开解答

按题目提示作正交变换,取

Y1=n Xˉ,Y2=nn−1(X1−Xˉ),Y_1=\sqrt n\,\bar X,\qquad Y_2=\sqrt{\frac n{n-1}}(X_1-\bar X),

再补齐正交坐标 Y3,…,YnY_3,\ldots,Y_n。正态向量经正交变换后各坐标独立,且

Y2σ∼N(0,1),1σ2∑i=3nYi2∼χ2(n−2).\frac{Y_2}{\sigma}\sim N(0,1),\qquad \frac1{\sigma^2}\sum_{i=3}^nY_i^2\sim\chi^2(n-2).

因此

T=Y2/σ∑i=3nYi2/[(n−2)σ2]∼t(n−2).T= \frac{Y_2/\sigma} {\sqrt{\sum_{i=3}^nY_i^2/[(n-2)\sigma^2]}} \sim t(n-2).

另一方面,

(n−1)S2=Y22+∑i=3nYi2.(n-1)S^2=Y_2^2+\sum_{i=3}^nY_i^2.

消元得

ξ=(n−1)Tn(n−2+T2),T∼t(n−2).\boxed{ \xi= \frac{(n-1)T} {\sqrt{n(n-2+T^2)}}, \qquad T\sim t(n-2) }.

因此 ∣ξ∣<(n−1)/n|\xi|<(n-1)/\sqrt n。ξ\xi 不是普通的 tt 变量,因为分母 SS 包含了 X1X_1 自己造成的残差。

评论