annotate annotate_peak/annotatePeak.R @ 12:54e19bceff70 draft

Uploaded
author testtool
date Thu, 16 Mar 2017 07:21:37 -0400
parents
children 77c0aee75ac6
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
12
54e19bceff70 Uploaded
testtool
parents:
diff changeset
1 require("ChIPseeker", quietly = TRUE)
54e19bceff70 Uploaded
testtool
parents:
diff changeset
2 require("ChIPpeakAnno", quietly = TRUE)
54e19bceff70 Uploaded
testtool
parents:
diff changeset
3 require("org.Hs.eg.db", quietly = TRUE)
54e19bceff70 Uploaded
testtool
parents:
diff changeset
4 require("GenomicRanges", quietly = TRUE)
54e19bceff70 Uploaded
testtool
parents:
diff changeset
5 options(warn = -1)
54e19bceff70 Uploaded
testtool
parents:
diff changeset
6
54e19bceff70 Uploaded
testtool
parents:
diff changeset
7 args <- commandArgs(trailingOnly = TRUE)
54e19bceff70 Uploaded
testtool
parents:
diff changeset
8 DMR = args[1]
54e19bceff70 Uploaded
testtool
parents:
diff changeset
9 annoPeakTable = args[2]
54e19bceff70 Uploaded
testtool
parents:
diff changeset
10
54e19bceff70 Uploaded
testtool
parents:
diff changeset
11
54e19bceff70 Uploaded
testtool
parents:
diff changeset
12 DMRInfo = read.table(
54e19bceff70 Uploaded
testtool
parents:
diff changeset
13 DMR,
54e19bceff70 Uploaded
testtool
parents:
diff changeset
14 header = FALSE,
54e19bceff70 Uploaded
testtool
parents:
diff changeset
15 sep = "\t",
54e19bceff70 Uploaded
testtool
parents:
diff changeset
16 stringsAsFactors = FALSE,
54e19bceff70 Uploaded
testtool
parents:
diff changeset
17 quote = ""
54e19bceff70 Uploaded
testtool
parents:
diff changeset
18 )
54e19bceff70 Uploaded
testtool
parents:
diff changeset
19
54e19bceff70 Uploaded
testtool
parents:
diff changeset
20 DMRPeaks <- GRanges(seqnames = DMRInfo[, 1],
54e19bceff70 Uploaded
testtool
parents:
diff changeset
21 ranges = IRanges
54e19bceff70 Uploaded
testtool
parents:
diff changeset
22 (start = DMRInfo[, 2], end = DMRInfo[, 3]))
54e19bceff70 Uploaded
testtool
parents:
diff changeset
23
54e19bceff70 Uploaded
testtool
parents:
diff changeset
24 annotatePeak <- as.data.frame(annotatePeak(DMRPeaks, level = "gene", annoDb = "org.Hs.eg.db"))
54e19bceff70 Uploaded
testtool
parents:
diff changeset
25
54e19bceff70 Uploaded
testtool
parents:
diff changeset
26 write.csv(annotatePeak, annoPeakTable, row.names = FALSE)