Mercurial > repos > testtool > annotate_peak

diff annotate_peak/annotatePeak.R @ 15:53df7871db21 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author testtool
date Thu, 16 Mar 2017 11:03:24 -0400
parents 4a489b0d247f
children ee3fc24b2108
line wrap: on
line diff
--- a/annotate_peak/annotatePeak.R	Thu Mar 16 10:59:08 2017 -0400
+++ b/annotate_peak/annotatePeak.R	Thu Mar 16 11:03:24 2017 -0400
@@ -3,7 +3,7 @@
 DMR = args[1]
 annoPeakTable = args[2]
 
-DMR <- ("test-data/DMR.bed")
+#DMR <- ("test-data/DMR.bed")
 DMRInfo = read.table(
   DMR,
   header = FALSE,
@@ -28,7 +28,7 @@
 ## do annotation by nearest TSS
 anno <- annotatePeakInBatch(peaks, AnnotationData=annoData)
 anno[1:2]
-annoPeakTable <- ('test-data/ChIPPeak.csv')
+#annoPeakTable <- ('test-data/ChIPPeak.csv')
 write.csv(anno, annoPeakTable, row.names = FALSE)