【画图专题】venn图

批量画韦恩图

#**************************
#zhangj
#plot venn pic
#2019-06-05
#**************************

rm(list=ls())
options(stringsAsFactors = F)
library(pacman)
p_load(VennDiagram)

fil <- list.dirs("/home/bmk/zhangj/biowork/mRNA")
fil <- fil[-1]
for(m in 1:length(fil)){
  setwd(paste(fil[m]))
 
  file <- list.files(pattern = "xls")
  data <- list()
  for(i in 1:length(file)){
    temp <- read.table(paste(file[i]),header = T,sep="\t",comment.char = "",check.names = F)
    data[[i]] <- temp[,1]
    
  }
  col <- c('#0099CC','#FF6666','#FFCC99','#0099CC')#圈线颜色
  names(data) <- unlist(strsplit(paste(file),fixed=T,split=".DEG_final.xls"))
  venn.diagram(data,
               filename="venn.png",
               imagetype="png",
               lwd=1,  #圈线粗度
               lty=1,  #圈线类型
               col=col[1:length(file)],
               fill=col[1:length(file)],
               cat.col=col[1:length(file)],
               cat.cex = 0.8,       #A和B的大小
               rotation.degree = 0, #旋转角度
               cex=1.5,             #里面交集字的大小
               alpha = 0.5,         #透明度
               reverse=TRUE,
               width=4000,
               height = 4000,
               resolution =600,
               margin=0.2)
}