什么是协方差与相关系数?协方差矩阵如何计算?np.cov函数

1 协方差

1.1 定义

协方差(Covariance)在概率论和统计学中用于衡量两个变量的总体误差。而方差是协方差的一种特殊情况,即当两个变量是相同的情况。

协方差的计算公式如下所示:



方差的计算公式如下所示:


可以看到协方差是度量两个变量的总体误差,而方差只考虑单变量的离散程度。

如果连个变量相互独立,则协方差为零。

2 协方差矩阵

假设X是一个n维随机变量

则它的协方差矩阵计算公式为:



其中

我们用一个例子来加以说明,加入我们有三面课程:语文、数学、应用,有5名考生他们的成绩如下:
语文 数学 英语
90 60 90
90 90 30
60 60 60
60 60 90
30 30 30

我们将该矩阵命名为矩阵A,这个矩阵共有三种属性,每种属性有5个变量,我们需要计算学科与学科之间的协方差,综合在一起就构成了协方差矩阵。

我们将语文、数学、英语分别编号为1、2、3,则它们之间的协方差记为E11、E12、E13、E22、E23、E33,最终该矩阵的协方差矩阵为:
\begin{bmatrix} E11 & E12 & E13\\ E21 & E22 & E23 \\ E31 & E32 & E33 \\ \end{bmatrix}
可以根据协方差计算公式进行计算:

首先,我们需要得到这三科的平均成绩:

语文 数学 英语
66 60 60

然后用矩阵A减去平均成绩(三科分别减去各科的均值),得到新的矩阵:
\begin{bmatrix} 24 & 0 & 30\\ 24 & 30 & -30 \\ -6 & 0 & 0 \\ -6 & 0 & 30 \\ -36 & -30 & -30 \\ \end{bmatrix}
E12的计算公式为:
E12=E(语文-语文成绩的均值)(数学-数学成绩的均值)
由于样本减均值刚刚已经计算完成,这里直接进行计算:
E12=24*0+24*30+(-6)*0+(-6)*0+(-36)*(-30)=1800
同理,E13的计算公式为:
E13=24*30+24*(-30)+(-6)*0+(-6)*30+(-36)*(-30)=900

根据Eij=Eji的性质,得到新的矩阵:
\begin{bmatrix} 2520 & 1800 & 900\\ 1800 & 1800 &0 \\ 900 & 0 & 3600 \\ \end{bmatrix}
然后除以样本的个数5,得到最后的协方差矩阵:
\begin{bmatrix} 504 & 360 & 180\\ 360 & 360 &0 \\ 180 & 0 & 720 \end{bmatrix}

知道了协方差矩阵如何计算之后我们来使用numpy内置的函数cov()来计算协方差矩阵。假设有两个变量x0x1,有三组观测(0,2)(1,1)和(2,0)。

import numpy as np

x = np.array([[0, 2], [1, 1], [2, 0]])

值得注意的是,np.cov(x)函数的默认输入矩阵x每一行代表一个特征,每一列代表一个观测,所以在进行协方差矩阵计算的时候需要对输入矩阵进行转置,或者将默认参数设置为False,如np.cov(x,rowvar=False)

import numpy as np
x = np.array([[0, 2], [1, 1], [2, 0]])
np.cov(x.T)

输出:

array([[ 1., -1.],
       [-1.,  1.]])

亦或者:

import numpy as np
x = np.array([[0, 2], [1, 1], [2, 0]])
np.cov(x, rowvar=False)

输出:

array([[ 1., -1.],
       [-1.,  1.]])

3 相关系数

相关系数是研究变量之间线性相关程度的量。

相关系数的计算公式如下所示:


\rho可以表示X和Y线性关系的紧密程度.

  • \rho=0则表示X和Y不线性相关,即\rho=0的充分必要条件是Cov(X,Y)=0,亦即不相关和协方差为零是等价的。
  • |\rho|\leq 1
  • |\rho|较大时,通常说X 和Y相关程度较好;当|\rho|较小时,通常说X和Y相关程度较差;当X和Y不相关,通常认为X和Y之间不存在线性关系,但并不能排除X和Y之间可能存在其他关系

参考:
协方差 - 百度百科
相关系数 - 百度百科
协方差矩阵概念

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 194,491评论 5 459
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 81,856评论 2 371
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 141,745评论 0 319
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 52,196评论 1 263
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 61,073评论 4 355
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 46,112评论 1 272
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 36,531评论 3 381
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 35,215评论 0 253
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 39,485评论 1 290
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 34,578评论 2 309
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 36,356评论 1 326
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 32,215评论 3 312
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 37,583评论 3 299
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 28,898评论 0 17
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 30,174评论 1 250
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 41,497评论 2 341
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 40,697评论 2 335

推荐阅读更多精彩内容