Navigation Menu

Skip to content

Instantly share code, notes, and snippets.

@dritoshi
Created January 16, 2012 12:25
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save dritoshi/1620637 to your computer and use it in GitHub Desktop.
Save dritoshi/1620637 to your computer and use it in GitHub Desktop.
get series_id from ExpressionSet (GEOquery)
library(GEOquery)
gse <- getGEO("GSE20861")
a <- t(sapply(sampleNames(gse[[1]]), function(x) { gsm <- getGEO(x); gsm@header$series_id} ))
a
[,1] [,2]
GSM521700 "GSE20859" "GSE20861"
GSM521701 "GSE20859" "GSE20861"
GSM521702 "GSE20859" "GSE20861"
GSM521703 "GSE20859" "GSE20861"
GSM521704 "GSE20859" "GSE20861"
GSM521705 "GSE20859" "GSE20861"
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment