vlambda博客
学习文章列表

R语言相关系数、显著性检验及可视化的尝试


最近在练习R实现相关分析,不试一下不知道,用R去做一个相关系数矩阵以及显著性检验,还这么麻烦,输出的结果也并不能灵活的展示,如果和SPSS或JASP相比的话,的确不怎么方便。


案例数据:犯罪数据state.x77

后台回复【犯罪数据】下载案例数据自由练习。


R语言相关系数、显著性检验及可视化的尝试


犯罪率及相关指标共8个变量,现在小兵希望考察8个连续数据间的相关系数,并做显著性检验。


相关系数及显著性检验


使用R的Hmisc包中rcorr函数实现,建议不要使用R内置的cor函数。


library(Hmisc)mycor <- rcorr(state.x77)


相关系数矩阵和p值结果


R语言相关系数、显著性检验及可视化的尝试


此时相关系数和显著性p值的结果都有了,不过可视化程度较低。所以接下来我还需要找一个可视化的函数来帮忙。


这就是R,所有的事情都要用户自己来决定,包括使用哪个函数以及何种展示效果。


好吧,我决定使用ggcorrplot包的函数实现相关系数可视化。


library(ggcorrplot)ggcorrplot(mycor$r, hc.order = TRUE, type = "lower", lab = TRUE)


可视化结果



这个还不错的吧。红色表示正相关,蓝色表示负相关,颜色的深浅代表相关程度的大小。


本文完

文/图=数据小兵


历史文章 ↓