韦恩图绘制
2019-08-29 本文已影响0人
落寞的橙子
直接上代码
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"))