# HG changeset patch
# User testtool
# Date 1488040141 18000
# Node ID be66730c5c3bc22b4a076bb58eda146ac83733db
# Parent  92387cb819622780e5a678f42725302d5d15a2e5
Uploaded

diff -r 92387cb81962 -r be66730c5c3b annotatePeak/annotatePeak.R
--- a/annotatePeak/annotatePeak.R	Sat Feb 25 11:18:27 2017 -0500
+++ b/annotatePeak/annotatePeak.R	Sat Feb 25 11:29:01 2017 -0500
@@ -1,6 +1,7 @@
 require("ChIPseeker", quietly = TRUE)
 require("ChIPpeakAnno", quietly = TRUE)
 require("org.Hs.eg.db", quietly = TRUE)
+require("GenomicRanges", quietly = TRUE)
 options(warn = -1)
 
 args <- commandArgs(trailingOnly = TRUE)
@@ -20,8 +21,7 @@
                     ranges = IRanges
                     (start = DMRInfo[, 2], end = DMRInfo[, 3]))
 
-annotatePeak <-
-  as.data.frame(annotatePeak(DMRPeaks, level = "gene", annoDb = "org.Hs.eg.db"))
-??org.Hs.eg.db
+annotatePeak <- annotatePeak(DMRPeaks, level = "gene", annoDb = "org.Hs.eg.db")
 
-write.csv(annotatePeak, annoPeakTable, row.names = FALSE)
+
+write.csv(as.data.frame(annotatePeak), annoPeakTable, row.names = FALSE)