Use GEOquery to get cell type data via GSM accession ID

r_geo_metadata.R

# requires GEOquery package to be installed
# writes output to stdout
# for more metadata, check the respective classes at http://rileylab.bio.umb.edu/sites/g/files/g1314676/f/201403/GEOquery.pdf
library(GEOquery)

gsm_ids <- c("GSM409307", "GSM1220601", "GSM1120314", "GSM1120315", "GSM1112815", "GSM1112816", "GSM1112817", "GSM1112818", "GSM410808")

for(iter in seq_along(gsm_ids)) {
    gds <- getGEO(gsm_ids[iter])
    cat(gsm_ids[iter], Meta(gds)$title)
}