|
71
|
1 require("GEOquery", quietly = TRUE)
|
|
|
2
|
|
|
3
|
|
|
4 options(warn = -1)
|
|
|
5 options("download.file.method"="wget")
|
|
|
6 args <- commandArgs(trailingOnly = TRUE)
|
|
|
7 GSMTable = args[1]
|
|
|
8 MetaTable = args[2]
|
|
|
9
|
|
72
|
10 TAB = read.table(GSMTable)
|
|
71
|
11
|
|
|
12
|
|
|
13 if (is.null(TAB)) {
|
|
|
14 stop("Must specify input files")
|
|
|
15 } else {
|
|
|
16 options(download.file.method.GEOquery = "wget")
|
|
|
17
|
|
|
18
|
|
|
19 GEODataTable <- getGEO(TAB$ID[1], getGPL = FALSE)
|
|
|
20 MetaData <- data.frame(Meta(GEODataTable))
|
|
|
21
|
|
|
22 write.table(MetaData, MetaTable, row.names = FALSE, sep = "\t")
|
|
|
23 }
|
|
|
24
|
