Description Usage Arguments Value Examples
View source: R/KEGG_function.R
The function will plot the correlation distributions for each enriched pathway (result from find_enriched_pathway function), and then Wilcoxon tests would be used to estimate the significance of correlations distribution between genes in each pathway and all genes.
1 2 | plot_pathway_cor(gene_expr, kegg_enriched_pathway, groups = NULL,
side = c("both", "pos", "neg"), alternative = NULL)
|
gene_expr |
the matrix for gene expression, row.names should be NCBI gene ID, such as 67040, 93683 |
kegg_enriched_pathway |
The returned value from find_enriched_pathway function, the enriched pathways. |
groups |
a character used to indicate expression values from different types of samples |
side |
a character string specifying the correlation directions interested, must be one of "both" (default), "pos" or "neg". |
alternative |
a character string specifying the alternative hypothesis, must be one of "two.sided" (default), "greater" or "less". You can specify just the initial letter. |
p values for Wilcoxon tests in each pathway
1 2 3 4 5 | data(pro_pho_expr)
data(pho_sites_count)
genes<-row.names(pho_sites_count)[which(pho_sites_count>=10)]
pho_KEGGresult<-find_enriched_pathway(genes,species='hsa')
result<-plot_pathway_cor(gene_expr=pro_pho_expr,kegg_enriched_pathway=pho_KEGGresult)
|
KEGG.db contains mappings based on older data because the original
resource was removed from the the public domain before the most
recent update was produced. This package should now be considered
deprecated and future versions of Bioconductor may not have it
available. Users who want more current data are encouraged to look
at the KEGGREST or reactome.db packages
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.