Abstract: Disclosed are data processing and analysis methods for gene expression data for identifying endogenous reference genes and a composition for the quantitative analysis of gene expression, comprising a pair of primers and/or probes useful in amplifying the identified endogenous reference genes. Introduced with the concepts of “Zero's proportion” and CV, the method allows different datasets to be integrally analyzed, thereby searching for novel reference genes. By the method, 2,087 genes are first found as housekeeping genes which are expressed in most tissues, and the usefulness thereof in the relative quantification of different target genes is determined by analyzing their expression stability. Of the 2,087 genes, 13 genes show higher expression stability with lower expression levels across a wide range of samples than traditional reference genes such as GAPDH and ACTS, and therefore are suitable for the normalization of universal genes having relatively low expression levels.