相关性代码

注意:①xlsl的字体格式-数值②第一行是文本元素,下面的行是数值③更改sheet=④csv格式没有sheet

library('corrplot')
library('eoffice')
library('tidyverse')
dat <- readxl::read_xlsx(file.choose(),sheet=3)#use sheet to select sheet
corr <- cor(dat)
corrplot(corr)
corrplot(corr, tl.col = 'black')
corrplot(corr, tl.col = 'black', order = 'hclust')
corr1 <- cor.mtest(dat)
corr1$p

pic1<-corrplot(corr, tl.col = 'black', order = 'hclust', p.mat = corr1$p,
          insig = 'label_sig', sig.level = c(0.001,0.01,0.05), pch.cex = 0.9,
          pch.col = 'green')
pic2<-corrplot(corr, tl.col = 'black', order = 'hclust', p.mat = corr1$p,
         insig = 'label_sig', sig.level = c(0.001,0.01,0.05), pch.cex = 0.9,
         pch.col = 'green', type = 'lower')
pic3<-corrplot(corr, tl.col = 'black', order = 'hclust', p.mat = corr1$p,
         insig = 'label_sig', sig.level = c(0.001,0.01,0.05), pch.cex = 0.9,
         pch.col = 'green', type = 'lower', method = 'color')

posted @ 2021-12-08 16:32  ^_^Anne  阅读(141)  评论(0)    收藏  举报