R语言韦恩图绘制

R语言直接上代码
参考文章一
参考文章二
绘制韦恩图(Venn plot)和花瓣图(Flower plot)的几个常用R包及在线网站

我的代码一

library (VennDiagram)
venn_dir="~/Venn"
venn.diagram(x= list(Group1=na.omit(pc_list),Group2=na.omit(mRNA_list,Group3=na.omit(lnc_list)), 
             filename = paste0(venn_dir,"/results.png"), 
             height = 1000, width = 1000,resolution =500, imagetype="png", col="transparent",
             fill=c("cornflowerblue","darkorchid1","grey"),alpha = 0.50, cex=0.45, cat.cex=0.45)
overlap_genes<-as.character(Reduce(intersect,list(na.omit(mRNA_list), 
                                                     na.omit(lnc_list),na.omit(pc_list))))
write.csv(overlap_genes,paste0(venn_dir,"/overlap_genes.csv"))

我的代码二

suppressMessages(library(VennDiagram))
venn.diagram(x= list(as.character(human_nc[,1]),as.character(mouse_nc[,1])),
             category.names = c("Human" , "Mouse" ),
             filename = paste0(your_dir,"/Human_mouse_overlap.png"),
             height = 600, width = 600,resolution =500, imagetype="png", 
             col=c( "#0A5EB9","#DF3D8C"),
             fill=c(alpha("#0A5EB9",0.3), alpha('#DF3D8C',0.3)),
             cex=0.3, cat.cex=0.3,
             compression = "lzw",
             lwd = 1,
             cat.default.pos = "outer",
             output = TRUE,
             cat.pos = c(27, -27),
             cat.dist = c(0.055, 0.055),
             inverted=T
)

保存为pdf

suppressMessages(library(gridSVG))
suppressMessages(library(VennDiagram))
dev.off()
my_plot <-venn.diagram(x= list(human_com,mouse_com_human),
                       category.names =  c("Human" , "Mouse" ),
                       filename = NULL,height = 600, width = 600,
                       col=c( "#E64B35FF","#4DBBD5FF"),
                       fill=c(alpha("#E64B35FF",0.3), alpha('#4DBBD5FF',0.3)),
                       cex=1, cat.cex=1,
                       compression = "lzw",
                       lwd = 2,
                       cat.default.pos = "outer",
                       output = TRUE,
                       cat.pos = c(27, -27),
                       cat.dist = c(0.055, 0.055),
                       inverted=T
)
grid.draw(my_plot)
#overlap.pdf 4/4
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容