Description Usage Arguments Details Value Author(s) References Examples
For a given set of candidate genes, reference genes and a list object of gene sets (e.g. Gene Ontology terms or gene sets from pathways) a one-sided Fisher's exact test ("greater") is performed for each gene set in the collection.
1 2 |
enrichment(genes, reference, genesets, adj = "fdr", verbose = FALSE)
|
genes |
A character vector of gene identifiers used as a candidate gene list that is assessed by the functional enrichment analysis. The candidate gene list is a subset of the reference gene list. |
reference |
A character vector of gene identifiers used as the reference gene list. Note all candidate genes are included in the reference gene list. |
genesets |
A named list object of a collection of gene sets. The identifiers used for the candidate and reference genes need to match the identifier types used for the gene sets. An example list of gene sets is given in data(exgensets) showing an example list object of gene sets from pathways with gene symbols. $'Reactome:REACT_115566:Cell Cycle' [1] "APITD1" "TAOK1" "CDC23" [4] "RB1" "PRKCA" "HIST1H4J" [7] "MCM10" "PPP1CC" "NUP153" ... $'Reactome:REACT_152:Cell Cycle, Mitotic' [1] "APITD1" "TAOK1" "CDKN2C" [4] "RB1" "PRKCA" "MCM10" [7] "HIST1H2BH" "NUP153" "TUBGCP3" [10] "APEX1" "RPA2" "PRKACA" ... |
adj |
The default value is <fdr> (False discovery rate using the Benjamini-Hochberg approach). Multiple testing correction based on the function stats::p.adjust() with available options for "holm", "hochberg", "hommel", "bonferroni", "BH", "BY", "fdr" and "none" |
verbose |
The default value is <FALSE>. If this option is set <TRUE> the number and name of the gene sets during their processing is reported. |
The enrichment analysis is based on a one-sided Fisher's exact test.
The function returns a data.frame object with the columns
"TermID" the given name of a gene set i from the named gene set collection list object S. "genes" the number of candidate genes present in the gene set i "all" the number of all genes present in the gene set i "pval" the nominal p-value from a one-sided fisher's exact test "padj" the adjusted p-value to consider for multiple testing
Ricardo de Matos Simoes
de Matos Simoes R, Tripathi S, Emmert-Streib F. Organizational structure and the periphery of the gene regulatory network in B-cell lymphoma. BMC Syst Biol. 2012; 6:38.
Inference and Analysis of Gene Regulatory Networks in R: Applications in Biology, Medicine, and Chemistry, DOI: 10.1002/9783527694365.ch10 In book: Computational Network Analysis with R, 2016, pp.289-306
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 | # In the following candidate genes are defined from a
# giant connected component of an example igraph network
# where we use the remaining genes of a given network
# as a reference list.
data(exanet)
data(exgensets)
candidate=V(getgcc(exanet))$name
reference=V(exanet)$name
# hypergeometric test is performed for
# each defined set of genes in the
# list object exgensets
tab.hypg=enrichment(candidate, reference, exgensets, verbose=TRUE)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.