其他
R绘图模板——Venn图绘制！
点击上方
“科研后花园”
关注我们
代码如下：
1、设置工作环境并加载R包：
rm(list=ls())#clear Global Environmentsetwd('D:\\Venn')#设置工作路径
#安装包#install.packages("venn")# install.packages("VennDiagram")#加载R包library(VennDiagram) # Generate High-Resolution Venn and Euler Plots library (venn)library( RColorBrewer) # not installed on this machinelibrary(ggplot2) # Create Elegant Data Visualisations Using the Grammar of Graphicslibrary(ggprism) # A 'ggplot2' Extension Inspired by 'GraphPad Prism'2、加载作图数据并对其进行预处理：
#读取数据,为OTU水平的丰度表data <- read.table(file="otu.txt",sep="\t",header=T,check.names=FALSE,row.names = 1)#查看前6行head(data)#组内合并df <- data.frame(A=rowSums(data[,c(1:3)]), B=rowSums(data[,c(4:6)]), C=rowSums(data[,c(7:9)]), D=rowSums(data[,c(10:12)]))head(df)#创建空列表df1 <- list()#获取每个样本（组）所有的OTUfor (i in 1:length(colnames(df))){ group<- colnames(df)[i] df1[[group]] <- rownames(df)[which(df[,i]!= 0)]}3、绘制Venn图：
#Venn包绘制venn(df1, zcolor=c('red','yellow','blue','green'),#颜色设置 opacity = 0.6,#颜色透明度 box=F,#边框去除 sncs=1.5,#组名字体大小 ilcs=0.8)#图片中数字大小4、计算各组总OTU个数并绘制柱状图进行展示：
df2<-dfdf2[df2>0]=1df3<-rbind(df2,Total=colSums(df2))#提取作图数据df4<-as.data.frame(t(df3[3010,])) df4$group<-rownames(df4)#绘图ggplot(df4,aes(x =group, y = Total)) + geom_col(aes(fill=group),width = 0.8,alpha=0.6)+ geom_text(aes(label=Total, y=Total+20), position=position_dodge(0.9), vjust=0)+ labs(x = NULL, y = NULL)+ theme_prism(palette = "candy_bright", base_fontface = "plain", base_family = "serif", base_size = 16, base_line_size = 0.8, axis_text_angle = 45)+ scale_y_continuous(expand = c(0,0),limits = c(0,1500))+ theme(legend.position = "none")+ scale_fill_manual(values = c('red','yellow','blue','green'))5、查看并导出交集结果：
df_inter <- get.venn.partitions(df1)for (i in 1:nrow(df_inter)) df_inter[i,'values'] <- paste(df_inter[[i,'..values..']], collapse = ', ')df_inter<-df_inter[-c(5, 6)]write.table(df_inter, 'df_Venn.txt', row.names = FALSE, sep = '\t', quote = FALSE)6、利用AI对图片进行拼接及美化：
温馨提示
如果你喜欢本文，请分享到朋友圈，想要获得更多信息，请关注我。