runGOEnrich {rliger} | R Documentation |
Run Gene Ontology enrichment analysis on differentially expressed genes.
Description
This function forms genesets basing on the differential expression result, and calls gene ontology (GO) analysis method provided by gprofiler2.
Usage
runGOEnrich(
result,
group = NULL,
useBg = TRUE,
orderBy = "padj",
logFCThresh = 1,
padjThresh = 0.05,
splitReg = FALSE,
...
)
Arguments
result |
Data frame of unfiltered output from |
group |
Selection of one group available from |
useBg |
Logical, whether to set all genes involved in DE analysis
(before threshold filtering) as a domain background of GO analysis. Default
|
orderBy |
Name of DE statistics metric to order the gene list for each
group. Choose from |
logFCThresh |
The log2FC threshold above which the genes will be used.
Default |
padjThresh |
The adjusted p-value threshold less than which the genes
will be used. Default |
splitReg |
Whether to have queries of both up-regulated and
down-regulated genes for each group. Default |
... |
Additional arguments passed to |
Value
A list object where each element is a result list for a group. Each result list contains two elements:
result |
data.frame of main GO analysis result. |
meta |
Meta information for the query. |
See gprofiler2::gost()
. for detailed explanation.
References
Kolberg, L. et al, 2020 and Raudvere, U. et al, 2019
Examples
res <- runMarkerDEG(pbmcPlot)
# Setting `significant = FALSE` because it's hard for a gene list obtained
# from small test dataset to represent real-life biology.
if (requireNamespace("gprofiler2", quietly = TRUE)) {
go <- runGOEnrich(res, group = 0, significant = FALSE)
}