概率与统计 · 抽样分布

抽样分布

统计量的抽样分布,是从总体里反复独立地抽大小为n的样本时,这个统计量取值的分布:每个样本贡献一个值、一个点。它的中心和真值之差是偏差,它的宽度是标准误。

抽样分布与本页画面

统计量是由样本算出来的一个数,比如均值、中位数。从同一个总体里一次又一次独立地抽大小为nn的样本,每次算一遍统计量,这些值的分布就是它的抽样分布。它的期望减去真值是偏差,它的标准差是标准误。画面的三层把三种分布分开:总体,一个样本的nn个值,以及每个样本贡献一个点的抽样分布。

为什么一个样本始终散得很开,均值却越收越拢?样本是总体的一份粗糙复制,nn变大时,它的标准差ss稳定在总体的σ\sigma附近,并不缩小。nn个独立值的均值就不同了:偏高和偏低的值互相抵消一部分,只要总体方差有限,就有Var⁡(xˉ)=σ2/n\operatorname{Var}(\bar x) = \sigma^2/n,所以均值的标准误是σ/n\sigma/\sqrt{n}。样本量变成四倍,标准误减半。可见标准误不是一个样本的散布,它说的是换一个样本时均值会移动多少。均值为什么还会排成钟形,见中心极限定理。

样本方差为什么除以n−1n - 1?每个样本都满足∑(x−μ)2=∑(x−xˉ)2+n(xˉ−μ)2\sum (x - \mu)^2 = \sum (x - \bar x)^2 + n(\bar x - \mu)^2。从样本自己的均值量起,离差平方和最小;它比从总体均值量起的平方和少了n(xˉ−μ)2n(\bar x - \mu)^2,平均起来,这一块是xˉ\bar x的方差的nn倍,恰好是σ2\sigma^2。所以∑(x−xˉ)2\sum (x - \bar x)^2平均是(n−1)σ2(n - 1)\sigma^2:除以n−1n - 1,s2s^2就是无偏的;除以nn,平均少σ2/n\sigma^2/n。这个修正只管方差:开平方是往下弯的,ss本身平均仍比σ\sigma小(正态数据、n=5n = 5时约小6%)。最大值偏小的原因更直接:均匀总体里每个数都小于100,样本的最大值平均起来也就比100小。

估计中心,均值总是更好吗?要看总体。对称的总体里,均值和中位数都以真正的中心为中心,标准误小的那个更好。正态数据里中位数的标准误更大,nn大时约是均值的π/2≈1.25\sqrt{\pi/2} \approx 1.25倍;尾巴更厚的拉普拉斯数据里它反而更小,约是1/2≈0.711/\sqrt{2} \approx 0.71倍。一个极端值能把均值拉走,却几乎挪不动中位数,而厚尾的总体常出极端值。偏斜的总体里,两者估计的可能不是同一个数:本页的指数总体均值是50,中位数是40+10ln⁡2≈46.9340 + 10\ln 2 \approx 46.93。

值得注意

  • 黄色的点和柱是一次固定种子的模拟,墨色曲线是精确分布(细虚线表示数值算出)。“模拟显示偏移”只描述这一次模拟;统计量有没有偏,看理论偏差。

操作

  • 换统计量或除数,只是对同一批样本重算,回放不中断;改nn或换总体会清空。

相关概念