[BioC] just.gcrma and normalisation
bmb at bmbolstad.com
Tue Jun 27 16:40:31 CEST 2006
I believe you can call just.gcrma() with a normalize=FALSE parameter to
turn off the quantile normalization stage.
If you have affyPLM loaded then you can use the normalize() method on an
exprSet. The methods available are pretty much the same as those in the
affy package for AffyBatch objects. The only thing to note is that you
will have to tell it to transform your data back to natural scale, which
is what most of the underlying normalization routines assume as input,
if your expression values are in log2 scale.
eset <- justGCRMA(normalize=FALSE)
eset.normalized <- normalize(eset,"contrast",transfn="antilog")
On Tue, 2006-06-27 at 13:33 +0100, alex lam (RI) wrote:
> Dear colleagues,
> I am very new to this so please forgive me if this question is too
> basic. I couldn't find an answer in the archive.
> I am trying to pre-process a large number of affy arrays and use the
> expression values in a linkage analysis for expression QTL mapping. I
> ran just.gcrma on the data because running ReadAffy then gcrma didn't
> work on so many arrays. It returned an exprSet object which has been
> normalised. But I have a feeling that the default quantile normalisation
> is not suitable for what I want to do.
> I think I can turn normalisation off, but then how do I apply other
> types of normalisation on an exprSet object? Expresso can only take an
> AffyBatch object, right?
> Alex Lam
> PhD student
> Department of Genetics and Genomics
> Roslin Institute (Edinburgh)
> Midlothian EH25 9PS
> Phone +44 131 5274471
> Web http://www.roslin.ac.uk
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
Ben Bolstad <bmb at bmbolstad.com>
More information about the Bioconductor