数据格式要求模板
读取数据并计算相关性
library(tidyverse)raw_data<-read_csv("data/20231110/Source data/Fig2c&d/environment.csv") %>% select(-c(1,2,3))
corrmatrix <- cor(raw_data, method = "spearman")corrmatrix
res1 <-corrplot::cor.mtest(corrmatrix, conf.level= .95)res1$pres1$lowCIres1$uppCI
绘制相关性分析热图
col3 <- grDevices::colorRampPalette(c("#2082DD", "white", "#FF3F3F")) col3corrplot::corrplot.mixed(corr=corrmatrix, lower="number", upper="circle", diag="u", upper.col =col3(20), lower.col = col3(20), number.cex=0.9, p.mat= res1$p, sig.level= 0.05, bg = "white", is.corr = TRUE, outline = FALSE, mar = c(0,0,3,0), addCoef.col = NULL, addCoefasPercent = FALSE, order = c("original"), rect.col = "black", rect.lwd = 1, tl.cex = 1.2, tl.col = "black", tl.offset = 0.4, tl.srt = 90, cl.cex = 1.1, cl.ratio = 0.2, tl.pos="lt", cl.offset = 0.5 )title(main="Correlation coefficient of 19 environmental variables",cex.main=2.1)
