Description Usage Arguments Author(s) Examples
Gene Ontology enrichment with subsequent clustering of terms
1 2 3 4 | clusteredEnrichmentGO(gomatrix, geneset, universe, geneAlias = NULL,
min.genes = 5, cut_max = 5, dynamic_cut = F, domains = c("BP", "MF",
"CC"), two_sided = F, representative_by_silhouette = F,
sort_by_pval = T, weights = NULL)
|
gomatrix |
Sparse binary matrix, maps gene identifiers (column names) to the GO terms they are annotated with (row names). |
geneset |
Character vector, gene set of interest. |
universe |
Character vector, background set of genes. Note: gene identifiers not included in the gomatrix are dropped! |
geneAlias |
Named character vector, maps gene identifiers (vector names) to alternative identifiers. |
min.genes |
Minimum number of genes of interest that a term needs to contain to be included in the analysis. |
cut_max |
Terms may be combined into a cluster if they contain at most cut_max different genes of interest (i.e. max. edit distance). |
dynamic_cut |
Default: FALSE. If TRUE, evaluate the silhouette index of all clusterings with max. edit distances of [1...cut_max], and choose the optimal cut accordingly. |
domains |
Character vector, select included GO terms by domain: biological process ("BP"), cellular compartment ("CC"), molecular function ("MF"). |
two_sided |
Default: FALSE. If TRUE, performs two-sided Fisher tests to detect depletion as well as enrichment. |
representative_by_silhouette |
Default: FALSE. If TRUE, selects the representative of a cluster based on their silhouette score (i.e. the most specific term for that cluster). Otherwise, the smallest term is chosen, or the term with the shortest description in the presence of ties. |
sort_by_pval |
Default: TRUE. Sort results by p-value. |
Robert Sehlke
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 | ## Not run:
# example for annotation matrix creation
# may be used interchangably with createGoMatrix()
gomatrix = createAnnotationMatrix( list("GO:0000042"=paste0("gene",1:6),
"GO:0000002"=paste0("gene",2:7),
"GO:0000009"=paste0("gene",6:12)) )
# example call for main function
res = clusteredEnrichmentGO(gomatrix,
geneset = paste0("gene",2:6),
universe = paste0("gene",1:12))
# example plot
plotClusteredEnrichment(res$results)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.