annotate GEO/getGEO.R @ 32:2b223a449cf6 draft

Uploaded
author testtool
date Thu, 09 Feb 2017 08:56:42 -0500
parents eb1a8cedf6c5
children a266a7c268e9
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
29
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
1 require("BiocGenerics", quietly = TRUE)
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
2 require("data.table", quietly = TRUE)
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
3 require("GEOquery", quietly = TRUE)
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
4 require("rtracklayer", quietly = TRUE)
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
5 require("FDb.InfiniumMethylation.hg19", quietly = TRUE)
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
6
32
2b223a449cf6 Uploaded
testtool
parents: 29
diff changeset
7 options(warn = -1,"download.file.method"="wget")
2b223a449cf6 Uploaded
testtool
parents: 29
diff changeset
8
29
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
9 args <- commandArgs(trailingOnly = TRUE)
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
10 GSMTable = args[1]
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
11 MetaTable = args[2]
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
12 IlmnTable = args[3]
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
13 gmTable = args[4]
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
14
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
15 TAB = fread(GSMTable)
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
16
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
17 if (is.null(TAB)) {
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
18 stop("Must specify input files")
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
19 } else {
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
20 GEODataTable <- getGEO(TAB$ID[1], getGPL = FALSE)
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
21 IlmnIDTable <- Table(GEODataTable)
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
22 MetaData <- data.frame(Meta(GEODataTable))
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
23
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
24
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
25 write.csv(MetaData, MetaTable, row.names = FALSE, sep = "\t")
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
26 }
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
27
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
28 hm450.hg19 <- getPlatform()
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
29
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
30 IlmnInfo <-
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
31 data.table(
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
32 IlmnID = names(hm450.hg19),
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
33 CHR = as.data.frame(hm450.hg19@seqnames)$value,
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
34 BP = as.numeric(hm450.hg19@elementMetadata$probeStart)
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
35 )
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
36
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
37
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
38 write.csv(IlmnInfo, IlmnTable, row.names = FALSE, sep = "\t")
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
39
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
40
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
41 if (length(TAB$ID) > 1) {
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
42 mysamples <-
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
43 do.call("data.table", lapply(TAB$ID[-1], function(x)
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
44 Table(getGEO(x, getGPL = FALSE))[, -1]))
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
45
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
46 gmSet <- data.table(IlmnIDTable[, -1], mysamples)
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
47
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
48 } else {
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
49 gmSet <- Table(getGEO(TAB$ID, getGPL = FALSE))
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
50 }
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
51
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
52
eb1a8cedf6c5 Uploaded
testtool
parents:
diff changeset
53 write.csv(gmSet, gmTable, row.names = FALSE, sep = "\t")