16
|
1 require("ChIPseeker", quietly = TRUE)
|
|
2 require("ChIPpeakAnno", quietly = TRUE)
|
|
3 require("EnsDb.Hsapiens.v75", quietly = TRUE)
|
|
4
|
12
|
5 options(warn = -1)
|
17
|
6 options("download.file.method"="wget")
|
|
7
|
12
|
8 args <- commandArgs(trailingOnly = TRUE)
|
|
9 DMR = args[1]
|
|
10 annoPeakTable = args[2]
|
|
11
|
15
|
12 #DMR <- ("test-data/DMR.bed")
|
12
|
13 DMRInfo = read.table(
|
|
14 DMR,
|
|
15 header = FALSE,
|
|
16 sep = "\t",
|
|
17 stringsAsFactors = FALSE,
|
|
18 quote = ""
|
|
19 )
|
|
20
|
14
|
21 peaks <- GRanges(seqnames = DMRInfo[, 1],
|
12
|
22 ranges = IRanges
|
|
23 (start = DMRInfo[, 2], end = DMRInfo[, 3]))
|
17
|
24
|
14
|
25 annoData <- toGRanges(EnsDb.Hsapiens.v75)
|
|
26 seqlevelsStyle(peaks) <- seqlevelsStyle(annoData)
|
17
|
27
|
|
28 anno <- annoPeaks(peaks, annoData=annoData,bindingType ="startSite")
|
|
29
|
15
|
30 #annoPeakTable <- ('test-data/ChIPPeak.csv')
|
17
|
31 write.csv(anno, annoPeakTable, row.names = FALSE) |