跳到正文
格致开物MATHWIKI

卡方分布

卡方分布(chi-squared distribution)描述独立标准正态变量的平方和,是正态样本方差和许多检验统计量的参照分布。自由度记录独立平方项的数量;其密度也可推广到任意正实数自由度。记号 χν2 表示分布,ν>0 表示自由度。

若两个独立测量误差都已经除以各自已知标准差,得到标准正态变量 Z1,Z2,则 Q=Z12+Z22 衡量总的平方偏离。一次实现为 (1,2) 时,Q=5;要判断它在该误差模型下是否常见,需要计算整个平方和分布的尾概率。

从一个平方到独立平方和

设 Z∼N(0,1)。对 q>0,变换 Q=Z2 有两个原像 q 与 −q,每支的雅可比绝对值都是 1/(2q)。因此 fZ2(q)=2e−q/22π2q=q−1/2e−q/221/2Γ(1/2). 这就是形状 1/2、尺度二的伽马分布。独立性使卷积适用;同尺度伽马变量相加时形状相加,所以对正整数 ν, Q=∑j=1νZj2∼χν2,fν(q)=qν/2−1e−q/22ν/2Γ(ν/2),q>0. 这里 Zj 必须相互独立、均值为零、方差为一。密度在负半轴为零,闭支撑为 [0,∞),在零点没有原子。由于它是伽马密度,归一化也适用于非整数 ν>0,此时不再解释为非整数个平方项。[1]

自由度二四六的卡方密度比较
每条曲线均未额外缩放;自由度改变位置和相对分散程度。

由伽马矩公式立即得到 𝔼Q=ν、Var⁡(Q)=2ν。也可从独立平方项核对:𝔼Z2=1、𝔼Z4=3,故每项方差为二。自由度增大时绝对方差增大,而相对标准差 2/ν 减小。

一个尾概率的完整计算

当 ν=2 时,密度简化为 f(q)=12e−q/2。于是开头观测的右尾概率为 Pr⁡(Q≥5)=∫5∞12e−q/2dq=e−5/2≈0.082085. 从五直到无穷远的面积正是此值,图只截取其中一段。它表示在声明的独立正态误差模型下,平方偏离至少这么大的概率,不是“模型为真的概率”。

自由度二时五右侧的尾部及阈值线
阴影显示右尾在图窗内的部分;连同十四以右的部分,总面积为 P(Q≥5)=e^(−5/2)。竖线处的密度高度不等于尾概率。

若允许右尾概率为 0.05,解 e−c/2=0.05 得阈值 c=−2log⁡0.05≈5.991465。作为变式,把误差实现改为 (1,1),则 Q=2、右尾为 e−1≈0.367879,尚未超过该阈值。

样本均值带走的一个自由度

设 X1,…,Xn 独立同分布于 N(μ,σ2),n≥2、σ>0,样本方差定义为 S2=(n−1)−1∑i(Xi−X¯)2。令标准化向量 Zi=(Xi−μ)/σ。在 ℝn 中选第一轴为 (1,…,1)/n,其余轴组成与它正交的标准正交基。正交变换后的联合密度仍为 (2π)−n/2exp⁡(−12∑j=1nwj2), 因为长度和雅可比绝对值一都保持不变。因此新坐标仍是独立标准正态。第一坐标为 W1=n(X¯−μ)/σ;减去样本均值恰好删去这个方向,于是 (n−1)S2σ2=∑j=2nWj2∼χn−12. 这解释了自由度的损失,也说明为何“标准化后有 n 个残差”不等于“有 n 个独立残差”:它们的和为零。该结论的精确分布依赖正态抽样,不能仅凭有限方差推广到任意总体。

独立条件也不能删掉。若令 Z2=Z1,两者各自都是标准正态,但 Q=2Z12 的方差为八,而 χ22 的方差为四。非零均值平方和涉及非中心卡方分布,单纯平移一个中心卡方变量不能替代它。

参考文献

  1. ↑ SciPy,scipy.stats.chi2,Notes:密度、自由度及 Gamma(df/2, scale=2) 关系,访问日期:2026年10月25日。