--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/getGSMdata/.Rhistory	Mon Feb 06 05:17:56 2017 -0500
@@ -0,0 +1,17 @@
+GSMTable <- ("test-data/GSMTable.txt")
+write.table(samples, GSMTable, row.names = FALSE, sep = "\t")
+?write.table
+GSE <- "GSE51547"
+search1 <- entrez_search(db = "gds", term = GSE, retmax = 1)
+require("rentrez", quietly = TRUE)
+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"
+write.table(samples, GSMTable, row.names = FALSE, sep = "\t",quote = FALSE)
+write.table(samples, file=GSMTable.txt, row.names = FALSE, sep = "\t",quote = FALSE)
+write.table(samples, file="GSMTable.txt", row.names = FALSE, sep = "\t",quote = FALSE)