协方差与相关系数
层级:A|必学
1. 两个变量是否一起变化
方差看单变量波动,协方差看两个变量偏离各自均值时是否同向。相关系数把协方差标准化到无单位范围。它们只捕捉线性共同变化,不等同独立或因果。
2. 协方差定义
Cov(X,Y)=E[(X−E[X])(Y−E[Y])].
展开:
Cov(X,Y)=E[XY]−E[X]E[Y].
- 正:常同向偏离均值;
- 负:常反向偏离;
- 0:无线性共同变化,但可能强非线性依赖。
Cov(X,X)=Var(X).
3. 尺度性质
Cov(aX+b,cY+d)=acCov(X,Y).
平移不影响协方差,缩放会改变数值和单位。若温度从摄氏换华氏,协方差随尺度变,因此跨变量难直接比较。
4. 相关系数
若 σX,σY>0:
ρXY=σXσYCov(X,Y).
Cauchy–Schwarz 保证
−1≤ρ≤1.
ρ=1 或 −1(有限方差且非退化)表示几乎确定的正/负线性关系。ρ=0 只表示不相关。
5. 独立与不相关
若独立且二阶矩存在:
E[XY]=E[X]E[Y]⇒Cov(X,Y)=0.
反例:X 关于 0 对称,Y=X2:
Cov(X,Y)=E[X3]−E[X]E[X2]=0,
但 Y 完全由 X 决定。
联合高斯是重要例外:零协方差等价于独立。
6. 协方差矩阵
随机向量 X∈Rd,均值 μ:
Σ=E[(X−μ)(X−μ)T].
Σij=Cov(Xi,Xj),对角是方差。它对称半正定:
aTΣa=Var(aTX)≥0.
线性变换 Y=AX+b:
E[Y]=Aμ+b,
Cov(Y)=AΣAT.
7. 样本协方差
S=n−11i=1∑n(xi−xˉ)(xi−xˉ)T.
若中心化数据矩阵 Xc 每行一个样本:
S=n−11XcTXc.
因此 S PSD,秩最多 min(d,n−1)。当 d≥n,样本协方差必奇异。
8. Pearson、Spearman 与 Kendall
- Pearson:标准化线性协方差,关注线性关系;
- Spearman:对秩做 Pearson,关注单调关系、对离群更稳健;
- Kendall tau:基于样本对的一致/不一致排序。
非线性非单调关系三者都可能接近 0。散点图与领域知识仍重要。
9. 偏相关与混杂
两个变量的相关可能由第三变量驱动。偏相关尝试在线性意义下控制其他变量后测量剩余关系,例如先分别回归掉 Z,再相关残差。
偏相关仍不是因果保证;测量误差、未观测混杂、非线性和 collider 条件化都可能误导。
10. 多重共线性
特征间高度线性相关使 XTX 条件差,回归系数不稳定、标准误增大。预测可能仍好,但单个系数解释困难。处理:
- 收集更多独立信息;
- 删除/合并特征;
- PCA;
- 岭/Elastic Net;
- 依据业务而非只看相关阈值。
11. 相关不等于因果
可能原因:
- X→Y;
- Y→X;
- 共同原因 Z;
- 选择偏差;
- 时间趋势;
- 巧合/多重检验。
因果推断需要干预或明确结构假设,不靠一个相关系数解决。
12. 相关矩阵的数值检查
应对称、对角为 1(非零方差变量)、PSD。成对删除缺失值可能让每个相关系数用不同样本,拼出的矩阵不再 PSD,影响分解。需要一致缺失处理或最近相关矩阵修正。
易错点
- 零相关不等于独立。
- 高相关不等于因果,也不说明斜率为 1。
- 协方差有单位,相关系数无单位。
- 相关只衡量特定关系,Pearson 主要线性。
- 标准化不等于让特征独立,只让边缘均值/方差改变。
常见问答
Q1:相关系数 0.8 是否一定“强”?
没有脱离领域的统一阈值。样本量、测量误差、分布、非线性和决策成本都影响解释。
Q2:两个 one-hot 列为何负相关?
互斥类别中一个为 1 会迫使另一个为 0;完整 one-hot 列还与截距线性相关。
Q3:协方差矩阵为何常加 shrinkage?
高维小样本下样本协方差噪声大或奇异,向对角/单位矩阵收缩可改善估计与求逆稳定性。
练习
- 证明 Cov(X,X)=Var(X)。
- 若 Cov(X,Y)=2,求 Cov(3X+1,−2Y+5)。
- 独立是否推出零相关?反向呢?
- n=20,d=50 时样本协方差最大秩是多少?
- 为什么 Pearson 可能漏掉 U 形关系?
答案与提示
- 直接代定义。
- 3(−2)2=−12。
- 有限二阶矩下推出;反向一般不成立。
- n−1=19。
- 正负线性趋势在两侧抵消,尽管存在确定非线性依赖。