[BioC] avereps function

Erika Melissari erika.melissari at bioclinica.unipi.it
Tue Apr 28 14:39:46 CEST 2009


Dear Claus-Jürgen,

on avereps help I read:

A new data object is computed in which each probe is represented by the 
average of its replicate spots. For an MAList object, the components M and A 
are both averaged in this way, as weights and any matrices found in 
object$other.
For an MAList object, ID defaults to MA$genes$ID is that exists, otherwise 
to rownames(MA$M).

If x is of mode "character", then the replicate values are assumed to be 
equal and the first is taken as the average.



Avereps schould works on an MA object, that is after normalization. If his 
object is not an MA list, it works like uniquegenelist function I think, and 
does not perform the average.

Erika


----- Original Message ----- 
From: "Claus-Jürgen Scholz" <scholz at klin-biochem.uni-wuerzburg.de>
To: <erika.melissari at bioclinica.unipi.it>
Cc: <bioconductor at stat.math.ethz.ch>
Sent: Tuesday, April 28, 2009 13:19 PM
Subject: Re: [BioC] avereps function


> Dear Erika,
>
> I think avereps simply calculates a mean value - which might be
> problematic if your expression values are log transformed after
> normalization... Maybe try to apply avereps before normalization and
> compare the results to those without averepsing.
>
> Bests,
> Claus-Jürgen
>
> Erika Melissari schrieb:
>>
>> Dear list,
>>
>> does anyone use avereps function by limma package to average probe
>> copies on microarray?
>> I used this function and obtained a strange result: my differentially
>> expressed gene list sensibly change if I use or not this function (I
>> perform the average before lmFit and after normalization).
>> I know this function only averages probe copies that are on a
>> microarray, also if they are in multiple not fixed number.
>> What is the problem?
>> Thank you very much
>>
>> Erika
>>
>> _______________________________________________
>> 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