# HG changeset patch # User kpbioteam # Date 1529590086 14400 # Node ID 17b814d78ac62ff788da71cc0df8d83686821bf2 # Parent ed66a80e81e85d1613415b0ca64db7f79697b7e5 planemo upload for repository https://github.com/kpbioteam/chipeakanno_annopeaks commit 5ee10ba73cd634a41b84f059f02edeaa04f69b8d-dirty diff -r ed66a80e81e8 -r 17b814d78ac6 chipeakanno_annopeaks.R --- a/chipeakanno_annopeaks.R Wed Jun 20 06:59:04 2018 -0400 +++ b/chipeakanno_annopeaks.R Thu Jun 21 10:08:06 2018 -0400 @@ -1,8 +1,3 @@ -source("https://bioconductor.org/biocLite.R") -biocLite("ChIPpeakAnno") - -biocLite("rtracklayer") - require("ChIPpeakAnno", quietly = TRUE) require("EnsDb.Hsapiens.v75", quietly = TRUE) require("rtracklayer", quietly = TRUE) @@ -21,7 +16,7 @@ DMRInfo <- read.table(input1) peaks <- GRanges(seqnames = DMRInfo[, 1], ranges = IRanges - (start = DMRInfo[, 2], end = DMRInfo[, 3])) + (start = DMRInfo[, 2], end = DMRInfo[, 3],score=DMRInfo[, 9])) annoData <- toGRanges(EnsDb.Hsapiens.v75) seqlevelsStyle(peaks) <- seqlevelsStyle(annoData)