Skip to content

方法Method

Hotelling T² 检验

Hotelling T-squared test · Hotelling T2 statistic

用样本协方差逆度量均值向量偏离,并在独立正态样本下通过精确 F 分布统一校准全部方向。

形式陈述 ​

多个响应彼此相关时,均值向量离原假设“多远”才算异常?设 X1,…,Xn 独立同分布于 Np(μ,Σ),Σ 正定,且 n>p。记样本均值为 X¯、使用分母 n−1 的样本协方差为 S。检验 H0:μ=μ0 的 Hotelling 统计量为

T2=n(X¯−μ0)TS−1(X¯−μ0).

Wishart 样本协方差理论保证 S 几乎必然可逆,并给出它与 X¯ 的独立性。原假设下有精确分布

n−pp(n−1)T2∼Fp,n−p.

因此一个 level-α 的联合均值检验在

T2>cα,cα=p(n−1)n−pFp,n−p;1−α

时拒绝。Fd1,d2;q 表示相应 F 分布的 q 分位数。

正态性、固定样本量、正定总体协方差和 n>p 都是这里有限样本精确性的条件;它不同于一般估计量的大样本 Wald 卡方近似。

直觉

如果两个响应经常一起上升,那么“两个都略高”不一定罕见;但“一个高、另一个低”可能违背常见共同波动,值得更大警觉。S−1 正是在按各方向的噪声大小重新衡量偏差。

逐坐标做两个 t 检验只看坐标轴。Hotelling 检验则考虑所有线性方向中最异常的那个,同时为这种选择付出正确的校准代价。这个最异常方向本身会随数据变化,不能把选中之后的一维 t 值当作事先固定方向的统计量。

例子与边界

同样欧氏距离,联合异常程度不同 ​

取 n=10,p=2,观察到

S=(10.80.81).

若均值差为 d=(0.5,0.5)T,它沿着方差较大的同向方向,得到

T2=102(0.5)21.8=259≈2.7778.

若均值差改为 d=(0.5,−0.5)T,欧氏长度完全相同,却沿着方差只有 0.2 的反向方向,得到 T2=25。转换成 F2,8 统计量,分别为 100/81≈1.2346 与 100/9≈11.1111。

5% 临界值为 F2,8;0.95≈4.45897,等价的 T2 阈值为 c0.05≈10.03268。因此前一个偏差不拒绝,后一个拒绝。

维数过大与正态假设失效 ​

若 n≤p,中心化数据至多有 n−1 个行方向,S 必然奇异。用广义逆或加上岭项可以构造别的统计量,但不会自动保留本页的精确 F 分布。

非正态样本下,即使 S 可逆,均值与协方差一般不独立,精确 F 校准也就没有依据。固定维数的大样本近似可能成立,但要另述矩条件和渐近论证;维数随样本量增长时还需不同分析。

推论与应用

精确 F 分布的核心比值 ​

先证明一个标准化结论:令 Z∼Np(0,Ip),W∼Wp(ν,Ip),两者独立,ν≥p,则

ν−p+1pZTW−1Z∼Fp,ν−p+1.

写 Z=RD,其中 R2∼χp2,D 为单位方向。W 的旋转不变性意味着,给定任意 D,DTW−1D 的分布都等于 e1TW−1e1,且不依赖 R,D。

为找这个分布,将生成 W 的标准正态数据矩阵按列分成 (a,B),其中 a∈Rν、B∈Rν×(p−1)。分块求逆给出

1e1TW−1e1=aT[I−B(BTB)−1BT]a.

给定 B,括号内是秩 ν−p+1 的正交投影,且 a 是独立标准正态向量,所以右侧服从 χν−p+12。p=1 时把 B 视为空矩阵,投影就是 I。

因此 ZTW−1Z 等于两个独立卡方变量的比值,除以各自自由度后正是 F 分布。这里的独立性来自条件于 Z 后,分母分布仍相同,而不只是两者边缘各自为卡方。

把样本均值和协方差代入 ​

原假设下令

Z=nΣ−1/2(X¯−μ0),W=Σ−1/2(n−1)SΣ−1/2.

它们分别是标准正态和自由度 ν=n−1 的标准 Wishart,且独立。又有 T2=(n−1)ZTW−1Z,代入上段便得到形式陈述中的系数和自由度。

当 p=1 时,T2 就是一元单样本 t 统计量的平方,F1,n−1=tn−12,作为边界核对。

椭球同时控制所有线性方向 ​

由检验反演得到精确覆盖 1−α 的均值置信椭球

C={μ:n(X¯−μ)TS−1(X¯−μ)≤cα}.

对任意 d,Cauchy–Schwarz 不等式及取等方向给出

supa≠0n(aTd)2aTSa=ndTS−1d.

所以同一个覆盖事件等价于对所有 a 同时满足

|aT(X¯−μ)|≤cαnaTSa.

这比“每个固定坐标分别有 95% 区间”包含更强的同时量词。它允许在观察后探查方向,而校准仍来自一开始就控制了所有方向的椭球。

参考资料
  • Stanford STATS305C,《Normal Theory》,Hotelling’s T2 与 Theorem 3.4.7 的正态—Wishart 比值。
  • Stanford STATS305C,《Hotelling T2》,Largest t statistic interpretation 与联合方向校准。
关系图谱14 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系