[BioC] To get expression data from GEO

vegardny vegard.nygaard at medisin.uio.no
Wed Mar 26 14:26:20 CET 2014


Hi Jerry.

The GSE27597 data set seems to consist of data from two different array 
platforms, so getGEO returns a list of two ExpressionSet objects. You 
need to use the list accessor, [[]], to get the individual elements.

Something like this:

library(GEOquery)
data1 = getGEO("GSE27597")
datExpr1 = exprs (data1[[1]])
dim(datExpr1) # the first smaller data set
datExpr2 = exprs (data1[[2]])
dim(datExpr2) # the second larger


And your line
datExpr = exprs (data1)
Just gave an error for me, so You might be using an old version of 
GEOquery.

This is my
sessionInfo()
R version 3.0.2 (2013-09-25)
Platform: x86_64-apple-darwin10.8.0 (64-bit)

locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8

attached base packages:
[1] parallel  stats     graphics  grDevices utils     datasets  methods 
  base

other attached packages:
[1] GEOquery_2.28.0    Biobase_2.22.0     BiocGenerics_0.8.0

loaded via a namespace (and not attached):
[1] RCurl_1.95-4.1 XML_3.95-0.2

Vegard.


On 2014-03-25 23:47, Jerry Cholo wrote:
> Hello,
> 
> 
> I would like to have expression data for GSE27597.  When I ran 
> following
> command line:
> 
> 
> 
> data1 = getGEO("GSE27597")
> 
> 
> 
> I noticed that there are two files related to this GSE (see below):
> 
> 
> 
> 
> 
> Found 2 file(s)
> 
> 
> GSE27597-GPL13243_series_matrix.txt.gz
> 
> Using locally cached version:
> C:\Users\AppData\Local\Temp\RtmpodgO6C/GSE27597-GPL13243_series_matrix.txt.gz
> 
> Using locally cached version of GPL13243 found here:
> 
> C:\Users\AppData\Local\Temp\RtmpodgO6C/GPL13243.soft
> 
> 
> 
> GSE27597-GPL5175_series_matrix.txt.gz
> 
> Using locally cached version:
> C:\Users\AppData\Local\Temp\RtmpodgO6C/GSE27597-GPL5175_series_matrix.txt.gz
> 
> Using locally cached version of GPL5175 found here:
> 
> C:\Users\AppData\Local\Temp\RtmpodgO6C/GPL5175.soft
> 
> 
> 
> When I ran following line
> 
> 
> 
>  datExpr = exprs (data1)
> 
> 
> 
> to get expression data I did not get the main data that includes more
> samples, but I got the smaller dataset that has few samples.
> 
> 
> 
> May someone let me know how I could deal with this problem?
> 
> 
> 
> Thanks,
> 
> Jerry
> 
> 	[[alternative HTML version deleted]]
> 
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives:
> http://news.gmane.org/gmane.science.biology.informatics.conductor



More information about the Bioconductor mailing list