Mercurial > repos > testtool > anno_peak_figure
changeset 5:1de0d1429a3c draft default tip
Uploaded
author | testtool |
---|---|
date | Mon, 20 Mar 2017 06:49:22 -0400 |
parents | 57082662d2e5 |
children | |
files | annoPeakFigure/annoPeakFigure.R |
diffstat | 1 files changed, 3 insertions(+), 3 deletions(-) [+] |
line wrap: on
line diff
--- a/annoPeakFigure/annoPeakFigure.R Fri Mar 17 07:20:20 2017 -0400 +++ b/annoPeakFigure/annoPeakFigure.R Mon Mar 20 06:49:22 2017 -0400 @@ -1,6 +1,7 @@ require("ChIPseeker", quietly = TRUE) require("ChIPpeakAnno", quietly = TRUE) require("EnsDb.Hsapiens.v75", quietly = TRUE) +require("data.table", quietly = TRUE) options(warn = -1) options("download.file.method"="wget") @@ -9,8 +10,7 @@ DMR = args[1] annoPeakFigure = args[2] -#DMR <- ("test-data/DMR.bed") -DMRInfo = read.table( +DMRInfo = fread( DMR, header = FALSE, sep = "\t", @@ -25,9 +25,9 @@ annoData <- toGRanges(EnsDb.Hsapiens.v75) seqlevelsStyle(peaks) <- seqlevelsStyle(annoData) + anno <- annoPeaks(peaks, annoData=annoData,bindingType ="startSite") -#annoPeakFigure <- ('test-data/annoPie.png') png(file = annoPeakFigure, width = 1200,