16
|
1 require("ChIPseeker", quietly = TRUE)
|
|
2 require("ChIPpeakAnno", quietly = TRUE)
|
|
3 require("EnsDb.Hsapiens.v75", quietly = TRUE)
|
18
|
4 require("data.table", quietly = TRUE)
|
16
|
5
|
12
|
6 options(warn = -1)
|
17
|
7 options("download.file.method"="wget")
|
|
8
|
12
|
9 args <- commandArgs(trailingOnly = TRUE)
|
|
10 DMR = args[1]
|
|
11 annoPeakTable = args[2]
|
|
12
|
18
|
13 DMRInfo = fread(
|
12
|
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
|
18
|
30 write.table(anno, annoPeakTable, row.names = FALSE) |