AbstractBecause quantitative reverse transcription PCR (RT-qPCR) gene expression data are compositional, amounts of quantified RNAs must be normalized using reference genes. However, the two most used methods to select reference genes (NormFinder and geNorm) ignore the compositional nature of RT-qPCR data, and often lead to different results making reliable reference genes selection difficult. We propose a method, based on all pairwise equivalence tests on ratio of gene expressions, to select genes that are stable enough to be used as reference genes among a set a candidate genes. This statistical procedure controls the error of selecting an inappropriate gene. Application to 30 candidate reference genes commonly used in human studies, assessed by RT-qPCR in RNA samples from lymphoblastoid cell lines of 14 control subjects and 26 patients with bipolar disorder, allowed to select 7 reference genes. This selection was consistent with geNorm’s ranking, less with NormFinder’s ranking. Our results provide an important fundamental basis for reference genes identification using sound statistics taking into account the compositional nature of RT-qPCR data. The method, implemented in the SARP.compo package for R (available on the CRAN), can be used more generally to prove that a set of genes shares a common expression pattern.