a funtion group deg data frame and make enrichment analysis in one step
enrich_kegg(
deg_data,
x,
y,
cut_FC = 1,
cut_P = 0.05,
organism = "hsa",
keyType = "kegg",
OrgDb = "org.Hs.eg.db",
pvalueCutoff = 0.05,
pAdjustMethod = "BH",
minGSSize = 10,
maxGSSize = 500,
qvalueCutoff = 0.2,
use_internal_data = FALSE,
label = c("Down", "Stable", "Up"),
label_ns = "Stable",
mc.cores = 1L
)
a DEG data frame contains logFC and p value
which column is log FC
which column is P value
a single number character or numeric vector in threshold value of log FC
threshold value of P value, can set for every cut_FC number in numeric vector format
supported organism listed in 'http://www.genome.jp/kegg/catalog/org_list.html'
one of "kegg", 'ncbi-geneid', 'ncib-proteinid' and 'uniprot'
OrgDb
adjusted pvalue cutoff on enrichment tests to report
one of "holm", "hochberg", "hommel", "bonferroni", "BH", "BY", "fdr", "none"
minimal size of genes annotated by Ontology term for testing.
maximal size of genes annotated for testing
qvalue cutoff on enrichment tests to report as significant. Tests must pass i) pvalueCutoff on unadjusted pvalues, ii) pvalueCutoff on adjusted pvalues and iii) qvalueCutoff on qvalues to be reported.
logical, use KEGG.db or latest online KEGG data
symbol word for groups
which group is the stable group
param for mclapply, choose cores to use in non-Windows machine
enrichKEGG result
if (FALSE) {
enrich_kegg(deg_data = DEG_df, x = "log2FoldChange", y = "pvalue")
}