[BioC] multiple Affy normalization results in single object ?
Martin Morgan
mtmorgan at fhcrc.org
Fri Sep 4 15:23:07 CEST 2009
Hi Wolfgang --
Wolfgang Raffelsberger wrote:
> Dear list,
>
> I'm wondering if / how I could store the results from multiple multiple
> Affy summarization/normalization approaches (eg RMA, vsn, etc) in a
> single object.
> Does the Biobase ExpressionSet support this ( I haven't succeded in
> doing so) or do I have to create a now class of objects ?
>
> Of course the simples way would be something like :
> require(affy); require(vsn)
> eset <- ReadAffy()
> combExps <- list(rma=exprs(rma(eset)),vsn= exprs(vsnrma(eset.cel)) #
> and so on for more methods...
>
> I suppose other people have already thought about this and I wonder if
> there are more elegant solutions like the code above...
ExpressionSet assayData must have an 'exprs' matrix, but can also have
others, added by assayDataElement<-. So
abatch <- ReadAffy()
combo = rma(abatch)
assayDataElement(combo, "vsn") = vsnrma(abatch)
and so on. Then things like
combo[,1:10]
do a coordinate subset. Retrieve the data with assayData(combo)[["vsn"]]
Martin
>
> For completeness :
>> sessionInfo()
> R version 2.9.1 (2009-06-26)
> i386-pc-mingw32
>
> locale:
> LC_COLLATE=French_France.1252;LC_CTYPE=French_France.1252;LC_MONETARY=French_France.1252;LC_NUMERIC=C;LC_TIME=French_France.1252
>
>
> attached base packages:
> [1] stats graphics grDevices datasets utils methods base
> other attached packages:
> [1] MASS_7.2-47 fdrtool_1.2.5 limma_2.18.2
> affyPLM_1.20.0 preprocessCore_1.6.0 gcrma_2.16.0
> Biostrings_2.12.7 [8] IRanges_1.2.3 affy_1.22.0
> Biobase_2.4.1
> loaded via a namespace (and not attached):
> [1] affyio_1.12.0 splines_2.9.1 tools_2.9.1
> Thank's in advance,
> Wolfgang
>
>
> . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
> Wolfgang Raffelsberger, PhD
> Laboratoire de BioInformatique et Génomique Intégratives
> CNRS UMR7104, IGBMC, 1 rue Laurent Fries, 67404 Illkirch Strasbourg,
> France
> Tel (+33) 388 65 3300 Fax (+33) 388 65 3276
> wolfgang.raffelsberger (a t) igbmc.fr
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> 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