R语言相关性的度量

R可以计算多种相关系数,包括pearson相关系数、Spearman相关系数、Kendall相关系数、偏相关系数、多分格(polychoric)相关系数和多系列(polyserial)相关系数。

1.pearson、spearson和Kendall相关

pearson积差相关系数衡量了两个定量变量之间的线性相关程度。spearson等级相关系数则衡量分级定序变量之间的相关程度,Kendall's Tau相关系数也是一种非参数的等级相关度量。

cor( )函数可以极端这三种相关系数,而cov()函数可以用来计算协方差。


> states<- state.x77[,1:6]

> cov(states)

              Population      Income  Illiteracy    Life Exp      Murder      HS Grad

Population 19931683.7588 571229.7796  292.8679592 -407.8424612 5663.523714 -3551.509551

Income      571229.7796 377573.3061 -163.7020408  280.6631837 -521.894286  3076.768980

Illiteracy      292.8680  -163.7020    0.3715306  -0.4815122    1.581776    -3.235469

Life Exp      -407.8425    280.6632  -0.4815122    1.8020204  -3.869480    6.312685

Murder        5663.5237  -521.8943    1.5817755  -3.8694804  13.627465  -14.549616

HS Grad      -3551.5096  3076.7690  -3.2354694    6.3126849  -14.549616    65.237894

> #可以看到,收入和高中毕业率之间存在显著的正先关性,文盲和预期寿命之间存在很强的负相关

> cor(states)

            Population    Income Illiteracy    Life Exp    Murder    HS Grad

Population  1.00000000  0.2082276  0.1076224 -0.06805195  0.3436428 -0.09848975

Income      0.20822756  1.0000000 -0.4370752  0.34025534 -0.2300776  0.61993232

Illiteracy  0.10762237 -0.4370752  1.0000000 -0.58847793  0.7029752 -0.65718861

Life Exp  -0.06805195  0.3402553 -0.5884779  1.00000000 -0.7808458  0.58221620

Murder      0.34364275 -0.2300776  0.7029752 -0.78084575  1.0000000 -0.48797102

HS Grad    -0.09848975  0.6199323 -0.6571886  0.58221620 -0.4879710  1.00000000

> #可以看到,pearson积差相关系数结果中,收入依然和高中毕业率有着很强的正相关;人口和文盲率之间有着很强的负相关性

> cor(states, method="pearson")

            Population    Income Illiteracy    Life Exp    Murder    HS Grad

Population  1.00000000  0.2082276  0.1076224 -0.06805195  0.3436428 -0.09848975

Income      0.20822756  1.0000000 -0.4370752  0.34025534 -0.2300776  0.61993232

Illiteracy  0.10762237 -0.4370752  1.0000000 -0.58847793  0.7029752 -0.65718861

Life Exp  -0.06805195  0.3402553 -0.5884779  1.00000000 -0.7808458  0.58221620

Murder      0.34364275 -0.2300776  0.7029752 -0.78084575  1.0000000 -0.48797102

HS Grad    -0.09848975  0.6199323 -0.6571886  0.58221620 -0.4879710  1.00000000

> cor(states, method="spearman")

          Population    Income Illiteracy  Life Exp    Murder    HS Grad

Population  1.0000000  0.1246098  0.3130496 -0.1040171  0.3457401 -0.3833649

Income      0.1246098  1.0000000 -0.3145948  0.3241050 -0.2174623  0.5104809

Illiteracy  0.3130496 -0.3145948  1.0000000 -0.5553735  0.6723592 -0.6545396

Life Exp  -0.1040171  0.3241050 -0.5553735  1.0000000 -0.7802406  0.5239410

Murder      0.3457401 -0.2174623  0.6723592 -0.7802406  1.0000000 -0.4367330

HS Grad    -0.3833649  0.5104809 -0.6545396  0.5239410 -0.4367330  1.0000000

> #举例分析,可以看到,文盲率和谋杀犯是紧密正相关的,收入和谋杀犯是强负相关的。


> x<- states[,c("Population","Income","Illiteracy","HS Grad")]

> y<- states[,c("Life Exp","Murder")]

> cor(x,y)

              Life Exp    Murder

Population -0.06805195  0.3436428

Income      0.34025534 -0.2300776

Illiteracy -0.58847793  0.7029752

HS Grad    0.58221620 -0.4879710

> #也可以有选择的对变量之间的相关性进行分析,结果表明:人口和预期寿命之间存在显著的负相关;文盲和犯罪率之间存在显著的正相关。

> library(ggm)

> colnames(states)

[1] "Population" "Income"    "Illiteracy" "Life Exp"  "Murder"    "HS Grad" 

> pcor(c(1,5,2,3,6),cov(states))

[1] 0.3462724

> #偏相关是指在控制一个或多个变量定量变量时,另外两个定量变量之间的相互关系。可以使用ggm包中的pcor()函数计算偏相关系数。格式为pcor()

> #pcor(μ,S)

> #其中μ是一个数值向量,前两个数值表示要计算相关系数的变量下标,其余的数值为条件变量的下标。s为变量的协方差阵。

> #这里的ggm包需要安装

> #分析可知,在控制了收入、文盲率和高中毕业率的影响时,人口和谋杀率之间的相关系数为0.346.偏相关系数常用于社会科学的研究中。


> cor.test(states[,3],states[,5])

        Pearson's product-moment correlation

data:  states[, 3] and states[, 5]

t = 6.8479, df = 48, p-value = 1.258e-08

alternative hypothesis: true correlation is not equal to 0

95 percent confidence interval:

0.5279280 0.8207295

sample estimates:

      cor

0.7029752

> #上述检验的是文盲率和谋杀率之间的关系,明显是显著相关的,相关系数为0.703

> library(psych)

> corr.test(states, use="complete")

Call:corr.test(x = states, use = "complete")

Correlation matrix

          Population Income Illiteracy Life Exp Murder HS Grad

Population      1.00  0.21      0.11    -0.07  0.34  -0.10

Income          0.21  1.00      -0.44    0.34  -0.23    0.62

Illiteracy      0.11  -0.44      1.00    -0.59  0.70  -0.66

Life Exp        -0.07  0.34      -0.59    1.00  -0.78    0.58

Murder          0.34  -0.23      0.70    -0.78  1.00  -0.49

HS Grad        -0.10  0.62      -0.66    0.58  -0.49    1.00

Sample Size

[1] 50

Probability values (Entries above the diagonal are adjusted for multiple tests.)

          Population Income Illiteracy Life Exp Murder HS Grad

Population      0.00  0.59      1.00      1.0  0.10      1

Income          0.15  0.00      0.01      0.1  0.54      0

Illiteracy      0.46  0.00      0.00      0.0  0.00      0

Life Exp        0.64  0.02      0.00      0.0  0.00      0

Murder          0.01  0.11      0.00      0.0  0.00      0

HS Grad          0.50  0.00      0.00      0.0  0.00      0

To see confidence intervals of the correlations, print with the short=FALSE option

> #虽然毕业率和人口数的相关系数为0.10,但是并不显著,其中p=0.50;而收入和毕业率的相关系数为0.62,显示显著。


相关性度量的基本知识到这就结束了,咱们下期再见!O(∩_∩)O哈哈~

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

推荐阅读更多精彩内容

  • rljs by sennchi Timeline of History Part One The Cognitiv...
    sennchi阅读 7,279评论 0 10
  • R中的线性回归函数比较简单,就是lm(),比较复杂的是对线性模型的诊断和调整。这里结合Statistical Le...
    真依然很拉风阅读 65,126评论 1 64
  • 时间久了,什么都会过期的,谁见过一直保鲜的爱情,我们村有人离婚了,也许是感情破碎了,也许是没钱买米了,也许是彼此越...
    自寅日记阅读 666评论 0 15
  • 这些话是在月考完的晚自习上突然想的,因为做的卷子感觉还是有些难度,也怕自己没有考好,考完学校安排还是要补课一天,晚...
    虾塘阿熙阅读 149评论 0 0
  • 前言:朋友说她今天遇见了一些不开心的人与事,我说那都不是事,睡一觉起来都忘了,明天依旧可以很快乐(๑>؂<๑)。 ...
    时遇小姐阅读 613评论 0 0