Mercurial > repos > testtool > get_geo
view getMETAdata/getMETAdata.R @ 39:5447fc8946ed draft default tip
Uploaded
author | testtool |
---|---|
date | Fri, 13 Oct 2017 11:22:49 -0400 |
parents | 608ab56a90d0 |
children |
line wrap: on
line source
require("data.table", quietly = TRUE) require("GEOquery", quietly = TRUE) options(warn = -1) options("download.file.method"="wget") args <- commandArgs(trailingOnly = TRUE) GSMTable = args[1] MetaTable = args[2] TAB = fread(GSMTable) if (is.null(TAB)) { stop("Must specify input files") } else { options(download.file.method.GEOquery = "wget") GEODataTable <- getGEO(TAB$ID[1], getGPL = FALSE) MetaData <- data.frame(Meta(GEODataTable)) write.table(MetaData, MetaTable, row.names = FALSE, sep = "\t") }