本文介绍协方差。
协方差
协方差表示的是两个变量的总体的误差,这与只表示一个变量误差的方差不同。 如果两个变量的变化趋势一致,也就是说如果其中一个大于自身的期望值,另外一个也大于自身的期望值,那么两个变量之间的协方差就是正值。 如果两个变量的变化趋势相反,即其中一个大于自身的期望值,另外一个却小于自身的期望值,那么两个变量之间的协方差就是负值。 —— 百度百科
定义
- 在概率论和统计学中,协方差用于衡量两个变量的总体误差。而方差是协方差的一种特殊情况,即当两个变量是相同的情况。
- 期望值分别为 E [ X ] E[X] E[X]与 E [ Y ] E[Y] E[Y]的两个实随机变量 X X X与 Y Y Y之间的**协方差 C o v ( X , Y ) Cov(X,Y) Cov(X,Y)定义为:
{% raw %}
Cov
(
X
,
Y
)
=
E
[
(
X
−
E
[
X
]
)
(
Y
−
E
[
Y
]
)
]
=
E
[
X
Y
]
−
2
E
[
Y
]
E
[
X
]
+
E
[
X
]
E
[
Y
]
=
E
[
X
Y
]
−
E
[
X
]
E
[
Y
]
\begin{aligned} \operatorname{Cov}(X, Y) &=E[(X-E[X])(Y-E[Y])] \\ &=E[X Y]-2 E[Y] E[X]+E[X] E[Y] \\ &=E[X Y]-E[X] E[Y] \end{aligned}
Cov(X,Y)=E[(X−E[X])(Y−E[Y])]=E[XY]−2E[Y]E[X]+E[X]E[Y]=E[XY]−E[X]E[Y]
{% endraw %}
相关系数
- 定义随机变量 X X X与 Y Y Y的相关系数:
ρ X Y = Cov [ X , Y ] Var [ X ] Var [ Y ] \rho_{X Y}=\frac{\operatorname{Cov}[X, Y]}{\sqrt{\operatorname{Var}[X]} \sqrt{\operatorname{Var}[Y]}} ρXY=Var[X]Var[Y]Cov[X,Y]
- 相关系数是协方差的归一化
与方差的关系
- 由定义可知:
{% raw %}
Cov
[
X
,
Y
]
=
Cov
[
Y
,
X
]
Cov
[
X
,
X
]
=
Var
[
X
]
Var
[
X
+
Y
]
=
Var
[
X
]
+
Var
[
Y
]
+
2
Cov
[
X
,
Y
]
\begin{array}{c} \operatorname{Cov}[X, Y]=\operatorname{Cov}[Y, X] \\ \operatorname{Cov}[X, X]=\operatorname{Var}[X] \\ \operatorname{Var}[X+Y]=\operatorname{Var}[X]+\operatorname{Var}[Y]+2 \operatorname{Cov}[X, Y] \end{array}
Cov[X,Y]=Cov[Y,X]Cov[X,X]=Var[X]Var[X+Y]=Var[X]+Var[Y]+2Cov[X,Y]
{% endraw %}
协方差的性质
- a , b a, b a,b 为常数,有:
Cov [ a X , b Y ] = a b Cov [ X , Y ] \operatorname{Cov}[a X, b Y]=a b \operatorname{Cov}[X, Y] Cov[aX,bY]=abCov[X,Y]
- 变量和的协方差:
Cov [ X 1 + X 2 , Y ] = Cov [ X 1 , Y ] + Cov [ X 2 , Y ] \operatorname{Cov}\left[X_{1}+X_{2}, Y\right]=\operatorname{Cov}\left[X_{1}, Y\right]+\operatorname{Cov}\left[X_{2}, Y\right] Cov[X1+X2,Y]=Cov[X1,Y]+Cov[X2,Y]
- 连续变量的协方差:
Cov [ f ( X ) , g ( Y ) ] = E [ ( f ( X ) − E [ f ( X ) ] ) ( g ( Y ) − E [ g ( Y ) ] ) ] \operatorname{Cov}[f(X), g(Y)]=\mathbb{E}[(f(X)-\mathbb{E}[f(X)])(g(Y)-\mathbb{E}[g(Y)])] Cov[f(X),g(Y)]=E[(f(X)−E[f(X)])(g(Y)−E[g(Y)])]
- 连续变量的协方差系数:
ρ [ f ( X ) , g ( Y ) ] = Cov [ f ( X ) , g ( Y ) ] Var [ f ( X ) ] Var [ g ( Y ) ] \rho[f(X), g(Y)]=\frac{\operatorname{Cov}[f(X), g(Y)]}{\sqrt{\operatorname{Var}[f(X)]} \sqrt{\operatorname{Var}[g(Y)]}} ρ[f(X),g(Y)]=Var[f(X)]Var[g(Y)]Cov[f(X),g(Y)]
协方差的物理意义
-
协方差的绝对值越大,说明两个随机变量都远离它们的均值
-
协方差为正,则说明两个随机变量同时趋向于取较大的值或者同时趋向于取较小的值
-
协方差为负,则说明一个随变量趋向于取较大的值,另一个随机变量趋向于取较小的值
-
两个随机变量的独立性可以导出协方差为零,但是两个随机变量的协方差为零无法导出独立性——正交-不相关-独立
因为独立性也包括:没有非线性关系
有可能两个随机变量是非独立的,但是协方差为零
如:假设随机变量 X ∼ U [ − 1 , 1 ] X \sim U[-1,1] X∼U[−1,1]。定义随机变量 的概率分布函数为:
P ( S = 1 ) = 1 2 P ( S = − 1 ) = 1 2 P(S=1)=\frac{1}{2} P(S=-1)=\frac{1}{2} P(S=1)=21P(S=−1)=21
定义随机变量 Y = S X Y=SX Y=SX,则随机变量 X , Y X,Y X,Y是非独立的,但是有:
C o v [ X , Y ] = 0 Cov[X,Y]=0 Cov[X,Y]=0
相关系数的物理意义
考虑以随机变量 X X X的线性函数 a + b X a+bX a+bX 来近似表示 Y Y Y。
- 均方误差, e e e为:
e = E [ ( Y − ( a + b X ) ) 2 ] = E [ Y 2 ] + b 2 E [ X 2 ] + a 2 − 2 b E [ X Y ] + 2 a b E [ X ] − 2 a E [ Y ] e=\mathbb{E}\left[(Y-(a+b X))^{2}\right]=\mathbb{E}\left[Y^{2}\right]+b^{2} \mathbb{E}\left[X^{2}\right]+a^{2}-2 b \mathbb{E}[X Y]+2 a b \mathbb{E}[X]-2 a \mathbb{E}[Y] e=E[(Y−(a+bX))2]=E[Y2]+b2E[X2]+a2−2bE[XY]+2abE[X]−2aE[Y]
- e取极小值,该优化问题为凸优化,直接求导为0即可:
{%raw%}
{
∂
e
∂
a
=
2
a
+
2
b
E
[
X
]
−
2
E
[
Y
]
=
0
∂
e
∂
b
=
2
b
E
[
X
2
]
−
2
E
[
X
Y
]
+
2
a
E
[
X
]
=
0
\left\{ \begin{array}{l} \frac{{\partial e}}{{\partial a}} = 2a + 2b\mathbb{E}[X] - 2\mathbb{E}[Y] = 0\\ \frac{{\partial e}}{{\partial b}} = 2b\mathbb{E}[{X^2}] - 2\mathbb{E}[XY] + 2a\mathbb{E}[X] = 0 \end{array} \right.
{∂a∂e=2a+2bE[X]−2E[Y]=0∂b∂e=2bE[X2]−2E[XY]+2aE[X]=0
{%endraw%}
- 解得:
{%raw%}
KaTeX parse error: Unknown column alignment: * at position 25: …{\begin{array}{*̲{20}{c}} {a = \…
{%endraw%}
- 带入 e e e得:
min ( e ) = E [ ( Y − ( a 0 + b 0 X ) ) 2 ] = ( 1 − ρ X Y 2 ) Var [ Y ] \min (e)=\mathbb{E}\left[\left(Y-\left(a_{0}+b_{0} X\right)\right)^{2}\right]=\left(1-\rho_{X Y}^{2}\right) \operatorname{Var}[Y] min(e)=E[(Y−(a0+b0X))2]=(1−ρXY2)Var[Y]
-
因此有以下定理:
- ∣ ρ X Y ∣ ≤ 1 \left|\rho_{X Y}\right| \leq 1 ∣ρXY∣≤1
- ∣ ρ X Y ∣ = 1 |\rho_{X Y}| = 1 ∣ρXY∣=1的充要条件为:存在常数 a , b a,b a,b使得 Y = a + b X Y=a+bX Y=a+bX
-
当 ∣ ρ X Y ∣ \left|\rho_{X Y}\right| ∣ρXY∣较大时, e e e较小,意味着随机变量 X X X和 Y Y Y联系较紧密。于是 ρ X Y \rho_{X Y} ρXY 是一个表征 X X X、 Y Y Y之间线性关系紧密程度的量
-
当 ρ X Y = 0 \rho_{X Y}=0 ρXY=0时,称 X X X和 Y Y Y 不相关。
- 不相关是就线性关系来讲的,而相互独立是一般关系而言的。
- 相互独立一定不相关;不相关则未必独立。
矩
- 设
X
X
X和
Y
Y
Y是随机变量:
- 若 E [ X k ] , k = 1 , 2 , ⋯ \mathbb{E}\left[X^{k}\right], k=1,2, \cdots E[Xk],k=1,2,⋯存在,则称它为 X X X的 k k k阶原点矩,简称 k k k阶矩。
- 若 E [ ( X − E [ X ] ) k ] , k = 2 , 3 , ⋯ \mathbb{E}\left[(X-\mathbb{E}[X])^{k}\right], k=2,3, \cdots E[(X−E[X])k],k=2,3,⋯存在,则称它为 X X X的 k k k阶中心矩。
- 若 E [ X k Y l ] , k , l = 1 , 2 , ⋯ \mathbb{E}\left[X^{k} Y^{l}\right], k, l=1,2, \cdots E[XkYl],k,l=1,2,⋯存在,则称它为 X X X的 Y Y Y的 k + l k+l k+l阶混合矩。
- 若 E [ ( X − E [ X ] ) k ( Y − E [ Y ] ) l ] , k , l = 1 , 2 , ⋯ \mathbb{E}\left[(X-\mathbb{E}[X])^{k}(Y-\mathbb{E}[Y])^{l}\right], k, l=1,2, \cdots E[(X−E[X])k(Y−E[Y])l],k,l=1,2,⋯存在,则称它为 X X X的 Y Y Y的 k + l k+l k+l阶混合中心矩。
- 期望是一阶原点矩
- 方差是二阶中心矩
- 协方差是二阶混合中心矩
协方差矩阵
- 设 n n n维随机变量 ( X 1 , X 2 , … , X n ) (X_1,X_2, \dots,X_n) (X1,X2,…,Xn)的二阶混合中心矩
c i j = Cov [ X i , X j ] = E [ ( X i − E [ X i ] ) ( X j − E [ X j ] ) ] c_{i j}=\operatorname{Cov}\left[X_{i}, X_{j}\right]=\mathbb{E}\left[\left(X_{i}-\mathbb{E}\left[X_{i}\right]\right)\left(X_{j}-\mathbb{E}\left[X_{j}\right]\right)\right] cij=Cov[Xi,Xj]=E[(Xi−E[Xi])(Xj−E[Xj])]
- 都存在,则称矩阵:
{%raw%}
C
=
[
c
11
c
12
⋯
c
1
n
c
21
c
22
⋯
c
2
n
⋮
⋮
⋱
⋮
c
n
1
c
n
2
⋯
c
n
n
]
\mathbf{C}=\left[\begin{array}{cccc} c_{11} & c_{12} & \cdots & c_{1 n} \\ c_{21} & c_{22} & \cdots & c_{2 n} \\ \vdots & \vdots & \ddots & \vdots \\ c_{n 1} & c_{n 2} & \cdots & c_{n n} \end{array}\right]
C=⎣⎢⎢⎢⎡c11c21⋮cn1c12c22⋮cn2⋯⋯⋱⋯c1nc2n⋮cnn⎦⎥⎥⎥⎤
{%endraw%}
- 为 n n n维随机变量 ( X 1 , X 2 , … , X n ) (X_1,X_2, \dots,X_n) (X1,X2,…,Xn)的协方差矩阵
- 由于 c i j = c j i c_{ij} = c_{ji} cij=cji 因此协方差矩阵是对称阵
- 由于对角线为各个变量的方差,因此对角线非负
- 通常 n n n 维随机变量的分布是不知道的,或者太复杂以致数学上不容易处理,因此实际中协方差矩阵非常重要。
参考资料
-
http://www.huaxiaozhuan.com/%E6%95%B0%E5%AD%A6%E5%9F%BA%E7%A1%80/chapters/2_probability.html
-
https://baike.baidu.com/item/%E5%8D%8F%E6%96%B9%E5%B7%AE/2185936?fr=aladdin