R计算两列数据的相关系数_R语言入门之相关性

今天这一期的内容主要是如何在R中进行数据之间的相关性分析,其实这一部分的内容和独立性检验的有点类似,大家可以对比着学习!

1. 相关性度量的assocstats()函数

在这里,我想和大家简单介绍一下如何度量列联表里分类变量之间的相关性。我们可以使用“vcd”包里的assocstats(x) 函数,这里x是一个列联表,示例如下:

A <- c(rep("male",15),rep("female",20),rep("male",15))# 创建变量A
B <- c(rep("healthy",4),rep("sick",35),rep("healthy",11)) # 创建变量B
C <- c(rep("smoker",26), rep("nonsmoker",24)) # 创建变量C
mydata <- data.frame(A,B,C) # 利用以创建的变量构建数据框
library(vcd) #加载R包
mytable <- xtabs(~C+B+A, data=mydata) #生成列联表
assocstats(mytable) #计算关联度

R计算两列数据的相关系数_R语言入门之相关性_第1张图片

从结果中可以看到,男性中的吸烟和患病有一定相关性(Phi-Coefficient=0.467 > 0.3,P值小于0.05)。由于数据的问题,女性没有计算出结果来,因为表格里有数据是0。

2. 协方差与相关系数

在R中你可以使用基础函数cor()来计算相关系数,用cov()函数来计算协方差。

其简单使用如下:

cor(x, use=, method= ),这里x是矩阵或者数据框,参数use=是用来指定缺失值的处理方法,而method=则是用来指定计算方法,默认的是计算Pearson相关系数。

# 计算mtcars数据集里变量之间的相关系数,删除缺失值 
cor(mtcars, use="complete.obs", method="kendall") #计算kendall相关系数
cov(mtcars, use="complete.obs") #计算协方差

不过,cor()和cov()不能进行显著性检验,虽然可以用corr.test()去进行单次的显著性检验,但这略显麻烦。

这里我想大家推荐使用”psych”包的corr.test()函数,它不仅给出相关系数,也给出各个相关系数的p值,使用很方便。

library(psych) #加载R包
corr.test(mtcars) #计算相关系数及其显著性
#另外我们也可以绘制相关系数图
data <- corr.test(mtcars)
r <- data$r
p <- data$p
corPlot(r,pval=p,numbers=TRUE,diag=FALSE,stars=TRUE)

R计算两列数据的相关系数_R语言入门之相关性_第2张图片

上面这幅图很好地展示了各个变量之间的相关系数及显著性,数字部分代表相关系数,星号代表显著性,星号越多表明越显著。

好了,今天的内容就分享到这里,最后提醒一下大家,“psych”包真的很强大,大家学有余力的时候一定要好好掌握它!

你可能感兴趣的:(R计算两列数据的相关系数)