This article has not been peer reviewed.Deposited research article
A new estimate of the proportion unchanged genes in a microarray experiment
Biological Sciences, AstraZeneca R&D Lund, S-221 87 Lund, Sweden
Genome Biology 2004, 5:P10 doi:10.1186/gb-2004-5-5-p10
This was the first version of this article to be made available publicly.Published: 1 April 2004
In the analysis of microarray data one generally produces a vector of p-values that for each gene give the likelihood of obtaining equally strong evidence of change by pure chance. The distribution of these p-values is a mixture of two components corresponding to the changed genes and the unchanged ones. The basic question 'What proportion of genes is changed' is a non-trivial one, with implications for the way that such experiments are analysed. An estimate not requiring any assumptions on the distributions is proposed and evaluated. The approach relies on the concept of a moment generating function.
A simulation model of real microarray data was used to assess the proposed method. The method fared very well, and gave evidence of low bias and very low variance.
The approach opens up a new possibility of sharpening the inference concerning microarray experiments, including more stable estimates of the false discovery rate.