其他
R绘图模板——热图+网络图展示mantel test相关性！！！
点击上方
“科研后花园”
关注我们
图片灵感来源：
通过R语言和AI复现得到的效果图：
更多精彩欢迎订阅并关注公众号：
1、设置工作环境并加载所需R包：
rm(list=ls())#clear Global Environmentsetwd('D:/桌面/mantel test')#设置工作路径
#加载包library(vegan) # Community Ecology Packagelibrary(dplyr) # A Grammar of Data Manipulationlibrary(ggcor) # Extended tools for correlation analysis and visualizationlibrary(ggplot2) # Create Elegant Data Visualisations Using the Grammar of Graphics2、加载绘图数据并进行格式转换（根据原文图片随机编写数据，无实际意义）：
#加载数据#OTU表格，并根据绘图需求对其进行编辑和转换df <- read.table("otu.txt",header = 1,row.names = 1,check.names = F,sep = "\t")df <-data.frame(t(df))df <- df[c(1,3,4,2,8,10,6,9,7,11,5,12),]#环境因子数据env <- read.table("env.txt",sep="\t",header = T,row.names = 1,check.names = F)3、绘制环境因子间相关性热图：
#绘制环境因子间相关性热图，pearson方法p <- quickcor(env,#数据 method = "pearson",#方法 show.diag = T) +#是否显示对角线 geom_square()+#类型 scale_fill_gradient2( high = '#fe6263', mid = 'white',low = '#7caaff') #颜色设置p4、mantel test分析及自定义标签：
#计算OTU与环境因子之间的mantel test的r值和p值并定义OTU数据分组df_mantel <- mantel_test(df, env, #数据 mantel.fun = 'mantel',#方法 spec.dist.method = 'bray', env.dist.method = 'euclidean', spec.select = list("Yield components" = 1:4, "Wheat growth indicators" = 5:8, "Yield" = 9:12))#将群落数据按组进行分开#定义标签,即对mantel分析后的数据按照个人需求进行分割df_mantel <- df_mantel %>% mutate(df_r = cut(r, breaks = c(-Inf, 0.25, 0.5, Inf), labels = c("< 0.25", "0.25 - 0.5", ">= 0.5")),#定义Mantel的R值范围标签 df_p = cut(p.value, breaks = c(-Inf, 0.01, 0.05, Inf), labels = c("< 0.01", "0.01 - 0.05", ">= 0.05")))#定义Mantel的P值范围标签#自定义连线类型df_mantel$linetype <- ifelse(df_mantel$p.value>=0.05,2,1)df_mantel$linetype <- factor(df_mantel$linetype,levels = c("1","2"))5、在热图基础上绘制图形展示mantel test分析结果：
quickcor(env, type = "upper",method = "pearson",show.diag = T,cor.test = T) + geom_square()+ scale_fill_gradient2( high = '#fe6263', mid = 'white',low = '#7caaff')+ geom_square() + geom_mark(r = NA,sig.thres = 0.05, size = 6, color = "black")+#显著性标签 anno_link(df_mantel, aes(color = df_p, size = df_r, linetype = linetype), label.size = 4, label.fontface = 1, curvature = 0.2,#连接线变为曲线 nudge_x =0.2)+#标签位置 scale_size_manual(values = c(0.8, 1.4, 2))+#连线粗细设置 scale_color_manual(values = c("#ffa83a","#dc4fff","#cacdd2"))+#线条颜色设置 scale_linetype_manual(values = c(1,2))+ guides(fill = guide_colorbar(title = "Pearson's r", order = 1),#图例相关设置 size = guide_legend(title = "Mantel's r",order = 2), color = guide_legend(title = "p-value", order = 3), linetype = "none") # 设置线条粗细###最后用AI进行微调即可
PS: 以上内容是小编个人学习代码笔记分享，仅供参考学习，欢迎大家一起交流学习。
完整代码（附带注释）及原始数据获取：
夸克网盘（手机端下载夸克APP进行搜索提取）：
链接：https://pan.quark.cn/s/f459ec712526
提取码：V9US
温馨提示
如果你喜欢本文，请分享到朋友圈，想要获得更多信息，请关注我。