Mercurial > repos > taoufik > microarrays
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.') } } }