Title: | An R/shiny package for microbiome functional enrichment analysis |
---|---|
Description: | This is an R/shiny package to perform functional enrichment analysis for microbiome data. This package was based on clusterProfiler. Moreover, MicrobiomeProfiler support KEGG enrichment analysis, COG enrichment analysis, Microbe-Disease association enrichment analysis, Metabo-Pathway analysis. |
Authors: | Guangchuang Yu [cre, aut] , Meijun Chen [aut] |
Maintainer: | Guangchuang Yu <[email protected]> |
License: | GPL-2 |
Version: | 1.13.0 |
Built: | 2024-10-30 08:24:18 UTC |
Source: | https://github.com/bioc/MicrobiomeProfiler |
COG enrichment analysis for microbiome data
enrichCOG( gene, dtype = "category", pvalueCutoff = 0.05, pAdjustMethod = "BH", universe, minGSSize = 10, maxGSSize = 500, qvalueCutoff = 0.2 )
enrichCOG( gene, dtype = "category", pvalueCutoff = 0.05, pAdjustMethod = "BH", universe, minGSSize = 10, maxGSSize = 500, qvalueCutoff = 0.2 )
gene |
a vector of COG ids. |
dtype |
one of "category", "pathway" |
pvalueCutoff |
adjusted pvalue cutoff on enrichment tests to report. |
pAdjustMethod |
one of "holm","hochberg","hommel","bonferroni","BH", "BY","fdr","none". |
universe |
universe background genes. If missing,use the all COGs. |
minGSSize |
minimal size of genes annotated by KEGG term for testing. |
maxGSSize |
maximal size of genes annotated for testing. |
qvalueCutoff |
qvalue cutoff on enrichment tests to report. |
A enrichResult
instance.
data(Psoriasis_data) cog <- enrichCOG(Psoriasis_data,dtype="category")
data(Psoriasis_data) cog <- enrichCOG(Psoriasis_data,dtype="category")
Metabolism enrichment analysis for microbiome data
enrichHMDB( metabo_list, pvalueCutoff = 0.05, pAdjustMethod = "BH", universe, minGSSize = 10, maxGSSize = 500, qvalueCutoff = 0.2 )
enrichHMDB( metabo_list, pvalueCutoff = 0.05, pAdjustMethod = "BH", universe, minGSSize = 10, maxGSSize = 500, qvalueCutoff = 0.2 )
metabo_list |
a vector of metabolites in HMDB.ID |
pvalueCutoff |
adjusted pvalue cutoff on enrichment tests to report. |
pAdjustMethod |
one of "holm", "hochberg", "hommel", "bonferroni", "BH", "BY", "fdr", "none". |
universe |
universe background genes. If missing, use HMDB db. |
minGSSize |
minimal size of genes annotated by KEGG term for testing. |
maxGSSize |
maximal size of genes annotated for testing. |
qvalueCutoff |
qvalue cutoff on enrichment tests to report. |
A enrichResult
instance.
x1 <- c("HMDB0000001","HMDB0000005","HMDB0000008") x2 <- enrichHMDB(x1)
x1 <- c("HMDB0000001","HMDB0000005","HMDB0000008") x2 <- enrichHMDB(x1)
KO enrichment for microbiome data
enrichKO( gene, pvalueCutoff = 0.05, pAdjustMethod = "BH", universe, minGSSize = 10, maxGSSize = 500, qvalueCutoff = 0.2 )
enrichKO( gene, pvalueCutoff = 0.05, pAdjustMethod = "BH", universe, minGSSize = 10, maxGSSize = 500, qvalueCutoff = 0.2 )
gene |
a vector of K gene id (e.g. K00001) or EC id (e.g. 1.1.1.27). |
pvalueCutoff |
adjusted pvalue cutoff on enrichment tests to report. |
pAdjustMethod |
one of "holm","hochberg","hommel","bonferroni","BH", "BY","fdr","none". |
universe |
universe background genes. If missing, use all K genes. |
minGSSize |
minimal size of genes annotated by KEGG term for testing. |
maxGSSize |
maximal size of genes annotated for testing. |
qvalueCutoff |
qvalue cutoff on enrichment tests to report. |
A enrichResult
instance.
data(Rat_data) ko <- enrichKO(Rat_data) head(ko)
data(Rat_data) ko <- enrichKO(Rat_data) head(ko)
Metabolism enrichment analysis for microbiome data
enrichMBKEGG( metabo_list, pvalueCutoff = 0.05, pAdjustMethod = "BH", universe, minGSSize = 10, maxGSSize = 500, qvalueCutoff = 0.2 )
enrichMBKEGG( metabo_list, pvalueCutoff = 0.05, pAdjustMethod = "BH", universe, minGSSize = 10, maxGSSize = 500, qvalueCutoff = 0.2 )
metabo_list |
a vector of metabolites in KEGG.ID |
pvalueCutoff |
adjusted pvalue cutoff on enrichment tests to report. |
pAdjustMethod |
one of "holm", "hochberg", "hommel", "bonferroni", "BH", "BY", "fdr", "none". |
universe |
universe background genes. If missing, use KEGG as default. |
minGSSize |
minimal size of genes annotated by KEGG term for testing. |
maxGSSize |
maximal size of genes annotated for testing. |
qvalueCutoff |
qvalue cutoff on enrichment tests to report. |
A enrichResult
instance.
mblist3 <- c("C00019","C00020","C00022") mb3 <- enrichMBKEGG(mblist3) head(mb3)
mblist3 <- c("C00019","C00020","C00022") mb3 <- enrichMBKEGG(mblist3) head(mb3)
Microbe-Disease associations enrichment analysis
enrichMDA( microbe_list, pvalueCutoff = 0.05, pAdjustMethod = "BH", universe, minGSSize = 10, maxGSSize = 500, qvalueCutoff = 0.2 )
enrichMDA( microbe_list, pvalueCutoff = 0.05, pAdjustMethod = "BH", universe, minGSSize = 10, maxGSSize = 500, qvalueCutoff = 0.2 )
microbe_list |
a vector of microbe ncbi tax ids. |
pvalueCutoff |
adjusted pvalue cutoff on enrichment tests to report. |
pAdjustMethod |
one of "holm", "hochberg", "hommel", "bonferroni", "BH", "BY", "fdr", "none". |
universe |
universe background genes. If missing, use disbiome as default. |
minGSSize |
minimal size of genes annotated by KEGG term for testing. |
maxGSSize |
maximal size of genes annotated for testing. |
qvalueCutoff |
qvalue cutoff on enrichment tests to report. |
A enrichResult
instance.
data(microbiota_taxlist) mda <- enrichMDA(microbiota_taxlist) head(mda)
data(microbiota_taxlist) mda <- enrichMDA(microbiota_taxlist) head(mda)
Module enrichment for microbiome data
enrichModule( gene, pvalueCutoff = 0.05, pAdjustMethod = "BH", universe, minGSSize = 10, maxGSSize = 500, qvalueCutoff = 0.2 )
enrichModule( gene, pvalueCutoff = 0.05, pAdjustMethod = "BH", universe, minGSSize = 10, maxGSSize = 500, qvalueCutoff = 0.2 )
gene |
a vector of K gene id (e.g. K00001). |
pvalueCutoff |
adjusted pvalue cutoff on enrichment tests to report. |
pAdjustMethod |
one of "holm","hochberg","hommel","bonferroni","BH", "BY","fdr","none". |
universe |
universe background genes. If missing, use all K genes. |
minGSSize |
minimal size of genes annotated by KEGG term for testing. |
maxGSSize |
maximal size of genes annotated for testing. |
qvalueCutoff |
qvalue cutoff on enrichment tests to report. |
A enrichResult
instance.
data(Rat_data) ko <- enrichModule(Rat_data) head(ko)
data(Rat_data) ko <- enrichModule(Rat_data) head(ko)
Metabolism enrichment analysis for microbiome data
enrichSMPDB( metabo_list, pvalueCutoff = 0.05, pAdjustMethod = "BH", universe, minGSSize = 10, maxGSSize = 500, qvalueCutoff = 0.2 )
enrichSMPDB( metabo_list, pvalueCutoff = 0.05, pAdjustMethod = "BH", universe, minGSSize = 10, maxGSSize = 500, qvalueCutoff = 0.2 )
metabo_list |
a vector of metabolites in smpdb Metabolite.ID |
pvalueCutoff |
adjusted pvalue cutoff on enrichment tests to report. |
pAdjustMethod |
one of "holm", "hochberg", "hommel", "bonferroni", "BH", "BY", "fdr", "none". |
universe |
universe background genes. If missing, use SMPDB db. |
minGSSize |
minimal size of genes annotated by KEGG term for testing. |
maxGSSize |
maximal size of genes annotated for testing. |
qvalueCutoff |
qvalue cutoff on enrichment tests to report. |
A enrichResult
instance.
smp <- enrichSMPDB(c("PW_C000164","PW_C000078","PW_C000040")) head(smp)
smp <- enrichSMPDB(c("PW_C000164","PW_C000078","PW_C000040")) head(smp)
download compound annotation of the latest version of KEGG pathway and stored in a 'GSON' object
gson_cpd()
gson_cpd()
a 'GSON' object
download compound annotation of the latest version of KEGG pathway to enzyme and stored in a 'GSON' object
gson_enzyme()
gson_enzyme()
a 'GSON' object
download KO annotation of the latest version of KEGG pathway and stored in a 'GSON' object
gson_KO()
gson_KO()
a 'GSON' object
download compound annotation of the latest version of KEGG Module and stored in a 'GSON' object
gson_module(db = "ko")
gson_module(db = "ko")
db |
ko or enzyme |
a 'GSON' object
This example data was reported on Forslund K,et al.2016 (doi: 10.1038/nature15766) and used for Microbe-Disease Association analysis.
a vector with 54 genera tax ids
https://www.nature.com/articles/nature15766
data(microbiota_taxlist)
data(microbiota_taxlist)
This example data was reported on Xiao S,et al.2021 (doi:10.3389/fcimb.2021.605825) and used for COG enrichment analysis.
a vector with 134 COGs
https://www.frontiersin.org/articles/10.3389/fcimb.2021.605825/full
data(Psoriasis_data)
data(Psoriasis_data)
This example data was reported on Fang S,et al.2019 (doi: 10.1111/1751-7915.13485) and used for KEGG enrichment analysis.
a vector with 91 KEGG Orthologies (KOs)
https://sfamjournals.onlinelibrary.wiley.com/doi/10.1111/1751-7915.13485
data(Rat_data)
data(Rat_data)
Run the Shiny Application
run_MicrobiomeProfiler( onStart = NULL, options = list(), enableBookmarking = NULL, uiPattern = "/", ... )
run_MicrobiomeProfiler( onStart = NULL, options = list(), enableBookmarking = NULL, uiPattern = "/", ... )
onStart |
A function that will be called before the app is actually run.
This is only needed for |
options |
Named options that should be passed to the |
enableBookmarking |
Can be one of |
uiPattern |
A regular expression that will be applied to each |
... |
arguments to pass to golem_opts. See '?golem::get_golem_options' for more details. |
Shiny application object.
if (interactive()) {run_MicrobiomeProfiler()}
if (interactive()) {run_MicrobiomeProfiler()}