[BioC] No replicates and differential analysis !!
magnus.rattray at manchester.ac.uk
Fri Jan 27 14:22:23 CET 2006
I tried to send this message to the mailing list yesterday but it seems
not to have got through.
Like others I'd caution against using no replicates at all, but I think
we have some methods that are useful with no or very few replicates.
Like Ben Bolstad said in a previous message, I would suggest that you
use the probe-level analysis to estimate the technical variance. We have
a method, multi-mgMOS, which estimates the log-signal and associates
this estimate with credibility intervals (either variances or
percentiles of the posterior log-signal). A paper describing this method
is available from,
and the most recent code and other publications are on the project web-site,
An older version of mmgMOS is available in bioconductor.
More recently we have developed a Bayesian t-test, PPLR, which allows
this technical variance to be included in determining differential
expression from replicated conditions. The PPLR code also works with no
replicates, in which case it just uses the technical variance. R code is
available from the above website (the paper is submitted and we will get
PPLR into bioconductor shortly).
I suggest you try mmgMOS, followed by mean or median centering of the
log-signal estimates, followed by PPLR. Let me know how you get on -
these are new methods and feedback is very welcome.
Dr Magnus Rattray
School of Computer Science,
University of Manchester,
Manchester M13 9PL, UK.
> Does anybody know a R package or function to compare expression level
> (affy data) of two groups with no replicates in each group ? In fact,
> just compare one array to an other.
> The purpose is to find differentially expressed genes.
> We cannot used statistical test (not enougth replicates), but we can
> used graphical approach based on scatter plot, and outliers detection
> Thanks for your help,
More information about the Bioconductor