enrichment_test {RVenn} | R Documentation |
Perform an enrichment test.
Description
Calculate the p-value of occurrence of an overlap between two sets by chance.
Usage
## S4 method for signature 'Venn'
enrichment_test(venn, set1, set2, univ = "all",
n = 10000, seed = 42)
Arguments
venn |
(Required) A |
set1 |
(Required) The name or the index of the set of interest. |
set2 |
(Required) The name or the index of the set to be checked whether
enriched in |
univ |
(Optional) Population size. Default is "all", implying the union
of all the sets in the |
n |
(Optional) Number of randomly generated sets. Default is 10,000 and minimum is 1,000. |
seed |
(Optional) An integer passed to set.seed function. It is used to fix a seed for reproducibly random number generation. Default is 42. |
Details
This type of analysis can also be performed by hypergeometric test or
Fisher's exact test. Here, the approach is similar to that described in
(Austin et
al., 2016). Briefly, the test is based on randomly generation of sets with
equal size to set1
from the background (universal) set. After creating
n (default is 10,000) random sets, the overlap between these and set2
is calculated to make a null distribution. When this distribution is true,
the probability of seeing an overlap at least as extreme as what was observed
(overlap between set1
and set2
) will be returned as the
p-value.
Value
Returns a list containing the probability (Significance) of occurrence of an overlap between two sets by chance and the number of occurrences (Overlap_Counts) in randomly generated sets.
Examples
set1 = c(1:20, letters[1:10])
set2 = letters[-26]
univ = unique(c(set1, set2, 21:200))
venn = Venn(list(set1, set2, univ))
e = enrichment_test(venn, 1, 2)
e$Significance