2021-08-01 R语言相关系数

https://www.dxy.cn/bbs/newweb/pc/post/39410770?source=rss

install.packages("Hmisc")

install.packages("ggplot2")

library(Hmisc)

library(ggplot2)

导出：

logout <- rcorr(as.matrix(data))
logout$r

write.csv(logout$r,"logoutr.csv")
logout$P
write.csv(logout$P,"logoutp.csv")

合并R和P

##定义函数
flattenCorrMatrix <- function(cormat, pmat) {
  ut <- upper.tri(cormat)
  data.frame(
    row = rownames(cormat)[row(cormat)[ut]],
    column = rownames(cormat)[col(cormat)[ut]],
    cor  =(cormat)[ut],
    p = pmat[ut]
    )
}

##合并
library(Hmisc)
logout<-rcorr(as.matrix(data[,1:7]))
flattenCorrMatrix(logout$r, logout$P)

##导出

write.csv(flattenCorrMatrix(logout$r, logout$P),"logoutall.csv")