view query_microarray_datasets/GetDatasets.R @ 0:cf862e6a6ddb draft

Uploaded
author taoufik
date Thu, 27 Dec 2018 05:54:24 -0500
parents
children
line wrap: on
line source

log_file=file("Query_Geo_Database.log", open = "wt")
sink(log_file)
sink(log_file, type = "out")
suppressWarnings(suppressMessages(library("batch")))
suppressWarnings(suppressMessages(library(GEOquery)))
listArguments = parseCommandArgs(evaluate=FALSE)
GeoCode=listArguments[["GeoCode"]]
GeoCode=toupper(GeoCode)
gethelp.df =tryCatch(eList <- getGEOSuppFiles(GeoCode), error = function(cond)"skip")
if(is.null(gethelp.df)){
  gethelp.df =tryCatch(eList <- getGEOSuppFiles(GeoCode), error = function(cond)"skip")
}
if(is.null(gethelp.df)){
  print('GeoCode with no raw data, choose another code or retry.')
}else{
  if(is.data.frame(gethelp.df)==TRUE){
    rawdata=rownames(eList)[grep('RAW',rownames(eList))][1]
    untar(rawdata, exdir = 'CEL')
  }else{
    if(gethelp.df=="skip"){
      print('Network Trouble Try Again Or Check your GeoCode')
    }else{
      print('GeoCode with no raw data, choose another code or retry.')
    }
  }
}