[BioC] R: How to use GEOquery to extract more than the default information from a GSE

[Manca Marco (PATH)]

Dear James,

thank you for your prompt and kind reply.

I was doing like the following and I wasn't able to see my annotation associated to the filesL
library("GEOquery")
gse <- getGEO("GSE9820")
gse

...following your suggestion I get exactly the same output as you.

Nevertheless I would love to be able to build my own ExprSet from a GSE using GEOquery with the option GSEMatrix=FALSE and then selecting the variables I want to import/include. In GEOquery's vignette there is an example of this but I am not able to find a document listing the options and the language/naming I should use to personalize the final file (the vignette only mentions that personalizing everything is quite difficult, but possible anyway).

Thank you.

Best regards,
Marco
________________________________________
Da: James F. Reid [james.reid at ifom-ieo-campus.it]
Inviato: venerdì 24 luglio 2009 15.38
A: Manca Marco (PATH)
Cc: sdavis2 at mail.nih.gov; bioconductor mailing list
Oggetto: Re: [BioC] How to use GEOquery to extract more than the default information from a GSE

Hi Marco,

I'm not sure what you mean by 'more than default information'.

Using GEOquery can be a bit complicated if the GEO series (GSE) contains
multiple platforms, but in your case you're fine because there is only one.

If you can get a complete ExpressionSet which stores samples annotation,
platform annotation and expression values by doing:

library("GEOquery")
gse <- getGEO("GSE9820")
names(gse)
##[1] "GSE9820_series_matrix.txt.gz"
gse[[1]]

which prints out:
ExpressionSet (storageMode: lockedEnvironment)
assayData: 20589 features, 153 samples
   element names: exprs
phenoData
   sampleNames: GSM247703, GSM247704, ..., GSM247855  (153 total)
   varLabels and varMetadata description:
     title: NA
     geo_accession: NA
     ...: ...
     data_row_count: NA
     (33 total)
featureData
   featureNames: ILMN_10000, ILMN_10001, ..., ILMN_9999  (20589 total)
   fvarLabels and fvarMetadata description:
     ID: NA
     GB_ACC: NA
     ...: ...
     SYNONYM: NA
     (6 total)
   additional fvarMetadata: Column, Description
experimentData: use 'experimentData(object)'
Annotation: GPL6255

fvarLabels(gse[[1]])
[1] "ID"         "GB_ACC"     "SYMBOL"     "DEFINITION" "ONTOLOGY"
[6] "SYNONYM"

contains all the information for the platform, varLabels will give you
the labels of the sample information and you can get to the expression
values by means of exprs(gse[[1]]).

HTH,
J.


Manca Marco (PATH) wrote:
> Dear Sean and dear bioconductors,
>
> I am writing you to ask a source of inspiration (code pieces, notes, references, whatever you might think appropriate) to import array annotation and other data from the GSE I am trying to work with (namely the GSE9820) into my eset.
>
> I have read on GEOquery's vignette that this is actually possible, despite being a bit tricky:
>
> "So, using a combination of lapply on the GSMList, one can extract as many columns of interest as necessary to build the data structure of choice. Because the GSM data from the GEO website are fully downloaded and included in the GSE object, one can extract foreground and background as well as quality for two-channel arrays, for example. Getting array annotation is also a bit more complicated, but by replacing \platform" in the lapply call to get platform information for each array, one can get other information associated with each array. Future work with this package will likely focus on better tools for manipulating GSE data" From http://www.bioconductor.org/packages/2.4/bioc/vignettes/GEOquery/inst/doc/GEOquery.pdf Page 22 of 22
>
> ...but I can't find anywhere any hint.
>
> Thank you in advance for your patience and support.
>
> My best regards,
> Marco
>