统计量的分布

常用的抽样分布

以下分布是后续常见的分布,分别是 χ2分布,t分布,F 分布。

χ2 分布

χ2 分布
n 个随机变量 X1,X2,,Xn 相互独立,都服从标准正态分布 N(0,1),则称 χ2=i=1nXi2 服从自由度为 nχ2 分布,记为 χ2χ2(n)。(读作“卡方”分布)
自由度指 包含的独立变量的个数

χ2(n) 分布的概率密度为(不用记):

fn(x)={12Γ(n/2)(x2)2n2ex2x>00x0Γ(α)=0+xα1exdx

性质:

  1. E(χ2)=nD(χ2)=2n
  2. 可加性:Y1χ2(n1)Y2χ2(n2),且 Y1,Y2 相互独立,则 Y1+Y2χ2(n1+n2)
α 分位数
给定 0<α<1,称满足条件 P(χ2>χ2_α(n))=α 的点 χα2(n)χ2 分布的 α 分位数
可以通过查表找到不同自由度 n 的上 α 分位数。

chi2distrubution

Tip

例题:XN(0,4),抽取样本 x1,x2,x3,x4,要使 a(x1+x2)2+b(3x34x4)2χ2(2),求 a,b

Note

解:

xiN(0,4)1i4x1+x2N(0,8)x1+x222N(0,1)3x34x4N(0,100)3x34x410N(0,1)18(x1+x2)2+1100(3x34x4)2χ2(2)a=18b=1100

t 分布

t 分布
XN(0,1)Yχ2(n) ,且 X,Y 相互独立,则称 T=XY/n 服从自由度为 nt分布,记为 Tt(n) (也称为学生分布)

t 分布的概率密度为(不用记):

f(x;n)=Γ(n+12)nπΓ(n2)(1+x2n)n+12

特别地,n=1t 分布就是柯西分布:f(x;1)=1π(1+x2)nt 分布就是标准正态分布:f(x;n)12πex2/2

α 分位数
给定 0<α<1,称满足 tα(n)f(x;n)dx=α 的点 tα(n)t(n) 分布的 α 分位数
同样可查表得出。

image_1bc2de1u91i7gd9q1shjm7spnam

[!TIP] 例题:XN(0,σ2)YN(0,σ2)X,Y 独立。取样 (X1,X2,X3,X4)(Y1,Y2,Y3,Y4),构造 U=X1+X2+X3+X4Y12+Y22+Y32+Y42,问 U 服从什么分布?

[!NOTE] 解:分子部分:

$$

X_1+X_2+X_3+X_4 \sim N(0,4\sigma^2)\ \frac{X_1+X_2+X_3+X_4}{2\sigma} \sim N(0,1) $$

分母部分:

YiσN(0,1)1σ2(Y12+Y22+Y32+Y42)χ2(4)


X1+X2+X3+X42σ/1σ2(Y12+Y22+Y32+Y42)/4t(4)U=X1+X2+X3+X4Y12+Y22+Y32+Y42t(4)

[!TIP] 例题:XN(0,σ2),抽取样本 X1,X2,问 X1|X2| 服从什么分布

[!NOTE] 解:分子部分:X1N(0,σ2)X1σN(0,1)
  分母部分:X2N(0,σ2)X2σN(0,1)
  X22σ2χ2(1)X1σX22σ2/1t(1)
  即 X1|X2|t(1)

F 分布

F 分布
Xχ2(n1)Yχ2(n2),且 X,Y 独立,则称随机变量 F=X/n1Y/n2 服从自由度为 (n1,n2) 的 F 分布,记作 FF(n1,n2)
n1 是第一自由度,n2 是第二自由度
F 分布的倒数依然服从 F 分布,即 1FF(n2,n1)

F 分布的密度函数为(不用记):

f(x;n1,n2)={1B(n1/2,n2/2)n1n1/2n2n2/2xn1/21(n2+n1x)(n1+n2)/2x>00x0B(a,b)=01xa1(1x)b1dx
α 分位数
给定 0<α<1,称满足 Fα(n1,n2)f(x;n1,x2)dx=α 的点 tα(n)t(n) 分布的 α 分位数
同样可查表得出。
查表可能用到的性质:F1α(n1,n2)=1/Fα(n2,n1)

image_1b2b9718u1s8hg3h1m6v10fqs8l83

[!TIP] 例题:已知 Xt(n),求 X2 ?

[!NOTE] 解:Xt(n),则存在 UN(0,1)Vχ2(n)U,V 独立,使得 X=UV/n
   X2=U2V/n
   UN(0,1)U2χ2(1)
   X2=U2/1V/nF(1,n)

单个正态总体的抽样分布

请拉到下面看修改后的内容。

**📚定理一:**

设总体 XN(μ,σ2)X1,X2,,Xn 是样本,样本均值 X=1ni=1nXi,样本方差 S2=1n1i=1n(XiX)2,则:

  1. XN(μ,σ2n) ;建议记 Xμσ/nN(0,1)
  2. (n1)S2σ2χ2(n1)
  3. XS2 相互独立。

🔍证明:

E(X)=E(1ni=1nXi)=1ni=1nE(Xi)=μD(X)=D(1ni=1nXi)=1n2i=1nD(Xi)=σ2nX线XN(μ,σ2n)

第二个结论的证明很复杂,在此不给出。

[!TIP] 例题:设总体 XN(μ,σ2)X1,X2,,Xn 是样本,求证

1i=1n(XiX)2σ2χ2(n1)2i=1n(Xiμ)2σ2χ2(n)

Note

解:对于(1),因为有一个约束条件:i=1n(XiX)=0,所以自由度为 n1;而对于(2),这 n 个随机变量之间相互独立,所以自由度为 nn 个标准正态分布的平方和服从 χ2 分布。

请记住上面例题中的两个结论。(上面的公式 (1)其实就是定理一中的公式(2))


我们可以对定理一中的 X 进行标准化:

Xμσ/nN(0,1)

如果不知道 σ,可以用 S 来代替 σ,此时有:

📚定理二:

设总体 XN(μ,σ2)X1,X2,,Xn 是样本,样本均值 X=1ni=1nXi,样本方差 S2=1n1i=1n(XiX)2,则:

XμS/nt(n1)

🔍证明:

Xμσ/nN(0,1)(n1)S2σ2χ2(n1)XμS/n=Xμσ/n(n1)S2σ2t(n1)

分割线(上面写的乱得一批,下面重新梳理过)


已知 XN(μ,σ2),从中抽取简单随机样本 (X1,X2,,Xn),样本的统计量满足如下分布。

结论1{:.success}

E(X)=E(1ni=1nXi)=1ni=1nE(Xi)=μD(X)=D(1ni=1nXi)=1n2i=1nD(Xi)=σ2nX线XN(μ,σ2n)Xμσ/nN(0,1)

结论2{:.success}

样本方差 S2=1n1i=1n(XiX)2 应该与总体方差 σ2 相近,故作替换:

XμS/nt(n1)

n 时,S2σ2t(n1)N(0,1)

结论3{:.success}

简单随机样本与总体同分布,故 XiN(μ,σ2),标准化得到:XiμσN(0,1),根据 χ2 分布的定义,有:

i=1n(Xiμσ)2χ2(n)i=1n(Xiμ)2σ2χ2(n)

结论4{:.success}

总体均值为 EX=μ,样本均值为 X=1ni=1nXi,这两者的值应该近似,故用样本均值替换结论3 中的总体均值,得到:

1σ2i=1n(XiX)2χ2(n1)n1σ21n1i=1n(XiX)2χ2(n1)(n1)S2σ2χ2(n1)

因为存在约束条件:i=1n(XiX)=0,故自由度少一个。

结论5{:.success}

(n1)S2σ2χ2(n1)E[(n1)S2σ2]=n1n1σ2ES2=n1ES2=σ2

注:这个结论不仅仅适用于正态总体,而是对于任意分布的总体,样本方差的均值都等于总体方差。

结论6{:.success}

XS2 相互独立。

两个正态总体的抽样分布

Xi,(i=1n1)Yj,(j=1n2) 分别来自总体 N(μ,σ12)N(μ2,σ22),并且它们相互独立。记样本均值为 X,Y,样本方差为 S12,S22,则有:

定理1{:.success}

F=S12σ12/S22σ22=S12S22/σ12σ22F(n11,n21)

🔍证明:

F=X/n1Y/n2,(n1)S2σ2χ2(n1)F=S12σ12/S22σ22=χ12/(n11)χ22/(n21)F(n11,n21)

定理2{:.success}

(XY)(μ1μ2)σ12n1+σ22n2N(0,1)

🔍证明:

XN(μ1,σ12n1),YN(μ2,σ22n2)XYN(μ1μ2,σ12n1+σ22n2)

定理3{:.success}

σ1=σ2=σ(XY)(μ1μ2)Sw1n1+1n2t(n1+n22)Sw2=(n11)S12+(n21)S22n1+n22

🔍证明:

(XY)(μ1μ2)σ12n1+σ22n2=(XY)(μ1μ2)σ1n1+1n2N(0,1)(n11)S12σ12χ2(n11),(n21)S22σ22χ2(n21)(n11)S12+(n21)S22σ2χ2(n1+n22)XS12,YS22t(XY)(μ1μ2)σ12n1+σ22n2/(n11)S12+(n21)S22σ2=(XY)(μ1μ2)Sw1n1+1n2t(n1+n22)

为什么用 Sw 不用 S1,S2

批注 2020-05-13 232004


6.2小结:

  • 对于单个正态总体,得到了 X,S2 的分布,用于对 μ,σ 进行推断。
  • 对于两个独立正态总体,得到了 XY,S12/S22 的分布,用于对 μ1μ2,σ12/σ22 进行推断。

参考资料

分位数的图片取自:吴良超的学习笔记:概率论与数理统计知识整理(6)-参数估计,符合 CC BY 4.0 许可协议。