First column is the element (possibly gene) identifier, and the second is its value on which to sort. For differential gene expression, values are often -log10(P-value) * sign(effect).
data.frame. Same as list1.
Controls the resolution of the test: how many items between any two overlap tests.
Character vector with two elements: the labels of the two lists.
Either "enrichment" for a one sided test, or "two.sided" for a two sided test. See Details section.
Logical. Should output plots be returned?
Path name where plots ae returned.
Logical. Should Benjamini-Yekutieli FDR corrected pvalues be computed?
Logical. Should pvalues be reported and plotted in -log10 scale and not -log scale?
Following the method in the reference, the function computes the number of overlapping elements in the first i*stepsize and j*stepsize elements of each list, and return the observed significance of this overlap using a hypergeometric test (see fisher.test).
The output is returned as a list of matrices including: the overlap in the first i*stepsize,j*stepsize elements and the significance of this overlap.
If plots=TRUE then plots of these matrices are stored in .jpg format.
In the case of alternative='two.sided' the pvalue plots are signed, just like in , thus distinguishing between over and under enrichment.
Matrix of -log(pvals) of the test for the first i,j elements of the lists.
Counts of the number of agreements in the first i,j elements of the lists.
An optional output of the B-Y corrected p-values of hypermat
Matrix of the type of deviation from the null. Negative for underenrichment and positive for overenrichment.
By default, pvalues are reported in (minus) the natural log scale and not in (minus) log 10 scale.
This behaviour is governed by log10.ind.
The p-values of the two-sided hypothesis test differ from those in reference . This is because the two-sided p-values suggested in , are based on taking either the upper or lower tail of the distribution without appropriately using both tails. This method does not correctly control the type I error rate. In the implementation here, for a two-sided test we sum the probabilities from both tails of the hypergeometric distribution.
See the package vignette for a small simulation.
Jonathan Rosenblatt and Jason Stein
 Plaisier, Seema B., Richard Taschereau, Justin A. Wong, and Thomas G. Graeber. "Rank-rank Hypergeometric Overlap: Identification of Statistically Significant Overlap Between Gene-expression Signatures." Nucleic Acids Research 38, no. 17(September 1, 2010)
 Benjamini, Y., and D. Yekutieli. 2001. "The Control of the False Discovery Rate in Multiple Testing Under Dependency." ANNALS OF STATISTICS 29 (4): 1165-1188.
 Stein JL(*), de la Torre-Ubieta L(*), Tian Y, Parikshak NN, Hernandez IA, Marchetto MC, Baker DK, Lu D, Lowe JK, Wexler EM, Muotri AR, Gage FH, Kosik KS, Geschwind DH. "A quantitative framework to evaluate modeling of cortical development by neural stem cells." Manuscript in press at Neuron.
(*) Authors contributed equally to this work.
R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
Loading required package: grid
> png(filename="/home/ddbj/snapshot/RGM3/R_BC/result/RRHO/RRHO.Rd_%03d_medium.png", width=480, height=480)
> ### Name: RRHO
> ### Title: Rank-Rank Hypergeometric Overlap Test
> ### Aliases: RRHO
> ### Keywords: htest
> ### ** Examples
> list.length <- 100
> list.names <- paste('Gene',1:list.length, sep='')
> gene.list1<- data.frame(list.names, sample(100))
> gene.list2<- data.frame(list.names, sample(100))
> # Enrichment alternative
> RRHO.example <- RRHO(gene.list1, gene.list2, alternative='enrichment')
> # Two sided alternative
> RRHO.example <- RRHO(gene.list1, gene.list2, alternative='two.sided')
Created & Maintained by Osamu Ogasawara (email@example.com) and