clusterEnrichment {transcriptogramer} | R Documentation |
If species
is a character, this method uses the biomaRt package
to build a gene2GO list, if species
is a data.frame, it will be used
instead.
The gene2GO list will be used with the
topGO package to detect the most significant terms of each cluster
present in the DE slot of the object
.
clusterEnrichment(object, universe = NULL, species, ontology = "biological process", algorithm = "classic", statistic = "fisher", pValue = 0.05, adjustMethod = "BH", nCores = 1L) ## S4 method for signature 'Transcriptogram' clusterEnrichment(object, universe = NULL, species, ontology = "biological process", algorithm = "classic", statistic = "fisher", pValue = 0.05, adjustMethod = "BH", nCores = 1L)
object |
An object of class Transcriptogram. |
universe |
A character vector containing ENSEMBL Peptide IDs, or NULL,
if the universe
is composed by all the proteins present in the transcriptogramS2 slot of
|
species |
A character string specifying the species; or a data.frame containing two columns, the first one with ENSEMBL Peptide IDs (character), which may, or not, to contain the taxonomy ID of the species as prefix, and the second containing its respective Gene Ontology term (character). |
ontology |
A character string specifying the Gene Ontology domain, ignoring case sensitivity, the possible values are 'biological process', 'cellular component' and 'molecular function', the default value of this argument is 'biological process'. |
algorithm |
Character string specifying which algorithm to use, the possible values are 'classic', 'elim', 'weight', 'weight01', 'lea' and 'parentchild', the default value of this argument is 'classic'. |
statistic |
Character string specifying which test to use, the possible values are 'fisher', 'ks', 't', 'sum' and 'globaltest', the default value of this argument is 'fisher'. |
pValue |
A numeric value between 0 and 1 giving the required family-wise error rate or false discovery rate, the default value is 0.05. |
adjustMethod |
Character string specifying p-value adjustment method, the possible values are 'none', 'BH', 'fdr' (equivalent to 'BH'), 'BY', 'hochberg', 'hommel', 'bonferroni', and 'holm', the default value of this argument is 'BH'. |
nCores |
An integer number, referring to the number of processing cores to be used; or a logical value, TRUE indicating that all processing cores should be used, and FALSE indicating the use of just one processing core. The default value of this argument is 1. |
A data.frame containing the most significant terms of each cluster.
Diego Morais
differentiallyExpressed, transcriptogramPreprocess, GSE9988, GPL570, Hs900, HsBPTerms, association, transcriptogramStep1, transcriptogramStep2
transcriptogram <- transcriptogramPreprocess(association, Hs900, 50) ## Not run: transcriptogram <- transcriptogramStep1(transcriptogram, GSE9988, GPL570) transcriptogram <- transcriptogramStep2(transcriptogram) levels <- c(rep(FALSE, 3), rep(TRUE, 3)) transcriptogram <- differentiallyExpressed(transcriptogram, levels, 0.01) terms <- clusterEnrichment(transcriptogram, species = "Homo sapiens", pValue = 0.005) ## this call also works terms <- clusterEnrichment(transcriptogram, species = HsBPTerms, pValue = 0.005) ## End(Not run)