source("../script/venn_improved.R") source("../UniPeak_1.0/extras/functions.R") regions <- read.unipeak("regions_tf_annot.txt", check.names=F, as.is=T) annot <- data.frame((regions[,"nearest Pol2 peak"] != ""), (regions[,"nearest CAGE peak"] != ""), (rowSums(is.na(regions[,grep("RNA-seq", colnames(regions))])) == 0)) colnames(annot) <- c("Pol II", "CAGE", "TSS") pdf("../figure/venn.pdf", width=(11-2)/2, height=(8.5-2)/2) venn.improved(annot, color=c("yellow", "red", "blue")) # because venn.improved gets them in the wrong order grid.text(sum(rowSums(annot) == 0), 0.8, 0.15) grid.text(paste0("N = ", nrow(annot), " HOT regions"), 0.5, 0.05) dev.off()