Mercurial > repos > testtool > get_gsm

view getGSMdata/test-data/getGSMdata.R @ 38:af9cb4b9f01e draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author testtool
date Mon, 07 Aug 2017 08:33:28 -0400
parents 31fe37b30836
children
line wrap: on
line source

require("rentrez", quietly = TRUE)

args <- commandArgs(trailingOnly = TRUE)

GSE = args[1]
GSMTable = args[2]



search1 <- entrez_search(db = "gds", term = GSE, retmax = 1)

rec2 <- entrez_summary(db = "gds",
                       id = search1$ids,
                       version = "2.0")

extract <- extract_from_esummary(rec2,rownames(summary(rec2)))

samples <- extract$samples
colnames(samples)[which(names(samples) == "accession")] <- "ID"
colnames(samples)[which(names(samples) == "title")] <-  "Phenotype"

samples$Phenotype<-gsub("\\s", "",samples$Phenotype) 

write.table(samples, GSMTable, row.names = FALSE,quote=FALSE,sep="\t")