Mercurial > repos > kpbioteam > chipeakanno_annopeaks

diff chipeakanno_annopeaks.R @ 1:a90b6d85463c draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/kpbioteam/chipeakanno_annopeaks commit b7298370ebae64340e792262c4e4f35568757871-dirty
author kpbioteam
date Wed, 02 May 2018 10:23:51 -0400
parents efd9e1315ed5
children 9ec5e9e2e7b2
line wrap: on
line diff
--- a/chipeakanno_annopeaks.R	Wed May 02 09:24:03 2018 -0400
+++ b/chipeakanno_annopeaks.R	Wed May 02 10:23:51 2018 -0400
@@ -13,7 +13,10 @@
 input4 = as.numeric(args[4])
 output1  = args[5]
 
-peaks <- get(load(input1))
+DMRInfo <-  read.table(input1)
+peaks <- GRanges(seqnames = DMRInfo[, 1],
+                    ranges = IRanges
+                    (start = DMRInfo[, 2], end = DMRInfo[, 3]))
 
 annoData <- toGRanges(EnsDb.Hsapiens.v75)
 seqlevelsStyle(peaks) <- seqlevelsStyle(annoData)