Search

Equivalence testing is of emerging importance in genomics studies but has hitherto been little studied in this content. In this paper, we define the notion of equivalence of gene expression and determine a `strength of evidence' measure for gene equivalence. It is common practice in genome-wide studies to rank genes according to observed gene-specific P-values or adjusted P-values, which are assumed to measure the strength of evidence against the null hypothesis of no differential gene expression. We show here, both empirically and formally, that the equivalence P-value does not satisfy the basic consistency requirements for a valid strength of evidence measure for equivalence. This means that the widely-used q-value (Storey, 2002) defined for each gene to be the minimum positive false discovery rate that would result in the inclusion of the corresponding P-value in the discovery set, cannot be translated to the equivalence testing framework. However, when represented as a posterior probability, we find that the q-value does satisfy some basic consistency requirements needed to be a credible measure of evidence for equivalence. We propose a simple estimate for the q-value from posterior probabilities of equivalence, and analyse data from a mouse stem cell microarray experiment which demonstrate the theory and methods presented here., Comment: 26 pages, 5 figures

In microarray experiments, it is often of interest to identify genes which have a pre-specified gene expression profile with respect to time. Methods available in the literature are, however, typically not stringent enough in identifying such genes, particularly when the profile requires equivalence of gene expression levels at certain time points. In this paper, the authors introduce a new methodology, called gene profiling, that uses simultaneous differential and equivalent gene expression level testing to rank genes according to a pre-specified gene expression profile. Gene profiling treats the vector of true gene expression levels as a linear combination of appropriate vectors, i.e., vectors that give the required criteria for the profile. This gene-profile model is fitted to the data and the resultant parameter estimates are summarized in a single test statistic that is then used to rank the genes. The theoretical underpinnings of gene profiling (equivalence testing, intersection-union tests) are discussed in this paper, and the gene profiling methodology is applied to our motivating stem cell experiment., Comment: 20 pages, 6 figures, final version of the paper is now published online in Biostatistics