[BioC] reading many affy U133 Plus_2.0
holger.schw at gmx.de
Fri Oct 5 15:11:55 CEST 2007
since we also had a similar problem (actually, we wanted to apply fitPLM to more than 400 Plus2 chips), I wrote a just-Version of fitPLM. This just-version works for up to about 100 Plus2 chips on an AMD Athlon XP 3000+ with 1 GB RAM. I have also implemented a function that allows to fit PLM models to several hundreds of Plus2 arrays by splitting up the problem in several chunks. (The background correction and the quantile normalization is done chipwise, afterwards the chips are split in several chunks such that the resulting PLM signals are the same as the ones that would have been returned by fitPLM if fitPLM could have been used.) For details, see
Chapter 5 of
If you are interested in any of these functions then let me know.
But please note that these functions were written more than a year ago using the BioC 1.8 (??) version of affyPLM. Haven't used these function since then. So no guarantee that they also work using the BioC 2.1 version of affyPLM. (But should work if Ben Bolstad did not recode the C code and the internal affyPLM functions in fitPLM. One problem might be that the resulting object is an object of class exprSet and not ExpressionSet.)
I was thinking about publishing these functions. But didn't do it because of the R packages BufferedMatrix and oligo, in which the idea of splitting the problem into several chunks is implemented in a more elegant way.
-------- Original-Nachricht --------
> Datum: Fri, 5 Oct 2007 13:57:53 +0200 (CEST)
> Von: darteta001 at ikasle.ehu.es
> An: bioconductor at stat.math.ethz.ch
> Betreff: [BioC] reading many affy U133 Plus_2.0
> Dear list,
> In order to normalise my data using fitPLM(), I need to create an
> affybatch using ReadAffy(). This time I think justRMA() is not an
> option I can use.
> Has anyone been successful at reading in a large group of Plus chips
> (say 100) using ReadAffy() and the normalise them using fitPLM() on
> a "normal" PC on Windows? What amount of memory is needed for it?
> Any examples on how many chips you might have loaded with memory
> details would be very appreciated.
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> Search the archives:
Psssst! Schon vom neuen GMX MultiMessenger gehört?
Der kanns mit allen: http://www.gmx.net/de/go/multimessenger
More information about the Bioconductor