The rank product is a biologically motivated rank test for the detection of differentially expressed genes in replicated microarray experiments.It is a simple non-parametric statistical method based on ranks of fold changes. In addition to its use in expression profiling, it can be used to combine ranked lists in various application domains, including proteomics, metabolomics, statistical meta-analysis, and general feature selection.
Given n genes and k replicates, let
rg,i
Compute the rank product via the geometric mean:
kr | |
RP(g)=(\Pi | |
g,i |
)1/k
See main article: article and Permutation test. Simple permutation-based estimation is used to determine how likely a given RP value or better is observed in a random experiment.
ERP(g)=c/p
pfp(g)=ERP(g)/rank(g)
rank(g)
RP
Permutation re-sampling requires a computationally demanding number of permutations to get reliable estimates of the p-values for the most differentially expressed genes, if n is large. Eisinga, Breitling and Heskes (2013) provide the exact probability mass distribution of the rank product statistic. Calculation of the exact p-values offers a substantial improvement over permutation approximation, most significantly for that part of the distribution rank product analysis is most interested in, i.e., the thin right tail. However, exact statistical significance of large rank products may take unacceptable long amounts of time to compute. Heskes, Eisinga and Breitling (2014) provide a method to determine accurate approximate p-values of the rank product statistic in a computationally fast manner.