view getGEO/getGEO.R @ 15:c035531adb1b draft

Uploaded
author testtool
date Thu, 19 Jan 2017 06:44:40 -0500
parents fafbd5bb1f05
children f49d9ccacc73
line wrap: on
line source

require("BiocGenerics", quietly = TRUE)
require("data.table", quietly = TRUE)
require("GEOquery", quietly = TRUE)
require("rtracklayer", quietly = TRUE)
require("FDb.InfiniumMethylation.hg19", quietly = TRUE)

options(warn = -1)
args <- commandArgs(trailingOnly = TRUE)
GSMTable = args[1]
MetaTable = args[2]
IlmnTable = args[3]
gmTable = args[4]

TAB = fread(GSMTable)

if (is.null(TAB)) {
  stop("Must specify input files")
} else {
  GEODataTable <- getGEO(TAB$ID[1], getGPL = FALSE)
  IlmnIDTable <- Table(GEODataTable)
  MetaData <- data.frame(Meta(GEODataTable))

  
  write.csv(MetaData, MetaTable, row.names = FALSE, sep = "\t")
}

hm450.hg19 <- getPlatform()

IlmnInfo <-
  data.table(
    IlmnID = names(hm450.hg19),
    CHR = as.data.frame(hm450.hg19@seqnames)$value,
    BP = as.numeric(hm450.hg19@elementMetadata$probeStart)
  )


write.csv(IlmnInfo, IlmnTable, row.names = FALSE, sep = "\t")


if (length(TAB$ID) > 1) {
  mysamples <-
    do.call("data.table", lapply(TAB$ID[-1], function(x)
      Table(getGEO(x, getGPL = FALSE))[, -1]))

  gmSet <- data.table(IlmnIDTable[, -1], mysamples)

} else {
  gmSet <- Table(getGEO(TAB$ID, getGPL = FALSE))
}


write.csv(gmSet, gmTable, row.names = FALSE, sep = "\t")