diff annoPeakFigure/annoPeakFigure.R @ 1:6278004aa9b2 draft

Uploaded

--- a/annoPeakFigure/annoPeakFigure.R	Mon Feb 06 10:58:35 2017 -0500
+++ b/annoPeakFigure/annoPeakFigure.R	Tue Feb 28 05:42:05 2017 -0500
@@ -5,7 +5,8 @@
 
 args <- commandArgs(trailingOnly = TRUE)
 DMR = args[1]
-annoPeakFigure = args[2]
+annoPeakFigure1 = args[2]
+annoPeakFigure2 = args[3]
 
 #DMR <- ("test-data/DMR.bed")
 DMRInfo = read.table(
@@ -22,12 +23,22 @@
 
 annotatePeak <- (annotatePeak(DMRPeaks, level = "gene", annoDb = "org.Hs.eg.db"))
 
-#annoPeakFigure <- ('test-data/annoPeakFigure.png')
+#annoPeakFigure1 <- ('test-data/annoPeakUpset.png')
+#annoPeakFigure2 <- ('test-data/annoPeakAnnoBar.png')
 
-png(file = annoPeakFigure,
+png(file = annoPeakFigure1,
+    width = 1200,
+    height = 600)
+
+plotAnnoBar(annotatePeak)
+
+dev.off()
+
+png(file = annoPeakFigure2,
     width = 1200,
     height = 600)
 
 upsetplot(annotatePeak, vennpie = TRUE)
 
 dev.off()
+