左手用R右手Python系列11——相关性分析

浏览: 2392

由于最近毕业论文缠身,一直都没有太多时间和精力撰写长篇的干货,但是呢学习的的脚步不能停止,今天跟大家盘点一下R语言与Python中到的相关性分析部分的常用函数。

常用的衡量随机变量相关性的方法主要有三种:


  • pearson相关系数;即皮尔逊相关系数,用于横向两个连续性随机变量间的相关系数。

  • spearman相关系数;即斯皮尔曼相关系数,用于衡量分类定序变量间的相关程度。

  • kendall相关系数;即肯德尔相关系数,也是一种秩相关系数,不过它所计算的对象是分类变量。

R语言:

  • cor

  • cor.test

  • corrplot

cor(x,y=NULL,use="everything",method= c("pearson","kendall","spearman"))

在R语言中,通常使用cor函数进行相关系数分析,可以分别指定向量,也可以指定给cor函数一个数据框。

  • use函数指定处理缺失值的方式

  • method是可选的三种相关系数计算方法。

这里以diamonds数据集为例:

library("ggplot2")

str(diamonds)

image.png

cor(diamonds[,c("carat","depth","price")])

cor(diamonds[,c("carat","depth","price")],method= "pearson")

image.png

默认情况下使用的是pearson相关系数。

corrplot函数可以针对相关系数输出的结果进行可视化:


library("corrplot")

library("dplyr")

cor(diamonds[,c("carat","depth","price")])%>%corrplot()

image.png

使用cor.test函数进行相关性的检验:

cor.test(x, y,    #指定带分析变量

alternative = c("two.sided", "less", "greater"),         

#双侧检验,单侧检验(默认双侧)

method = c("pearson", "kendall", "spearman"), 

 #相关性算法(默认pearson法)

exact = NULL, conf.level = 0.95, continuity = FALSE, ...)

cor.test(diamonds$carat,diamonds$depth)

cor.test(~carat+depth,diamonds)

image.png

以上两种写法都是支持的。

从结果可以看到,两者几乎不相关,pearson相关系数仅有0.02左右。

Pyhton:

import pandas as pd

import numpy as np

diamonds=pd.read_csv('D:/R/File/diamonds.csv',sep = ',',encoding = 'utf-8')

diamonds.info()

image.png

pandas中带有相关系数函数pandas.corr

image.png

mydata=diamonds[["carat","depth","table","price"]]

mydata.info()

image.png

mydata.corr()    

#可以直接给出数据框的相关系数矩阵

mydata.corr()["carat"]    

#给出caret变量与其他变量之间的相关系数

mydata["carat"].corr(mydata["depth"])    

#计算"carat"与"depth"之间的相关系数

与R语言中一样,pandas中内置的相关系数算法也是针对针对数值型变量的pearson法。

mydata.corr(method='pearson')                             

mydata.corr(method='pearson')["carat"]                    

mydata["carat"].corr(method='pearson',mydata["depth"])    

image.png

method也可以指定spearman法和kendall法计算相关系数。

本文小结:

R语言:


  • cor

  • cor.test

  • corplot

Python:

pandas.corr

最近上线的新课程,R语言在商务场景中的应用,内容适合想要通过ggplot2深入学习R语言数据可视化的爱好者,录播已经所有代码和学习资源已经上传了!

image.png

已经有接近200位小伙伴儿们在修习了,报名参加还可以获取小魔方的过往学习笔记、历史案例代码哦!

欢迎关注魔方学院交流群:


推荐 0
本文由 杜雨 创作,采用 知识共享署名-相同方式共享 3.0 中国大陆许可协议 进行许可。
转载、引用前需联系作者,并署名作者且注明文章出处。
本站文章版权归原作者及原出处所有 。内容为作者个人观点, 并不代表本站赞同其观点和对其真实性负责。本站是一个个人学习交流的平台,并不用于任何商业目的,如果有任何问题,请及时联系我们,我们将根据著作权人的要求,立即更正或者删除有关内容。本站拥有对此声明的最终解释权。

1 个评论

python怎么样实现每每两列的相关系数和P值,都显示出来。输出结果为矩阵。

要回复文章请先登录注册