This package provides the visualization of bayesian network inferred from gene expression data. The networks are based on enrichment analysis results inferred from packages including clusterProfiler and ReactomePA. The networks between pathways and genes inside the pathways can be inferred and visualized.
インストールする前に、エージェントに与えられる指示の中身を確認できます。
Infer and plot a Bayesian network of genes within a specific enriched pathway using expression data.
library(CBNplot)
library(bnlearn)
library(org.Hs.eg.db)
# 1. Load User Inputs & Prepare Data
data(gaussian.test)
kegg <- org.Hs.egPATH2EG
mapped <- mappedkeys(kegg)
genes <- as.list(kegg[mapped])[["00532"]]
counts <- head(gaussian.test, length(genes))
row.names(counts) <- genes
pway <- clusterProfiler::enrichKEGG(gene = genes)
pway <- clusterProfiler::setReadable(pway, org.Hs.eg.db, keyType="ENTREZID")
# 2. Visualize the gene-level Bayesian network for the first pathway
bngeneplot(results = pway, exp = counts, pathNum = 1, expRow="ENTREZID")
# 3. Extract the underlying network structure, strength, and direction
ret <- bngeneplot(results = pway, exp = counts, pathNum = 1, returnNet=TRUE, expRow="ENTREZID")
head(ret$str)
# 4. Convert the averaged network to an igraph object and perform centrality analysis
g <- bnlearn::as.igraph(ret$av)
igraph::evcent(g)$vector
# 5. Customize the gene network visualization
bngeneplotCustom(results = pway, exp = counts, expRow="ENTREZID", pathNum=1, fontFamily="sans", glowEdgeNum=NULL, hub=1)
Inputs/Outputs: Takes an enrichment analysis result object and a normalized expression data frame as inputs, and outputs Bayesian network plots and network structure data frames.
bngeneplot.bnpathplot.bngeneplotCustom or bnpathplotCustom.igraph object using bnlearn::as.igraph for centrality analysis.minet or GENIE3 instead).pathview instead).clusterProfiler::enrichKEGG or ReactomePA).counts) where row names match the gene identifier type (e.g., ENTREZID) specified in the enrichment results.clusterProfiler::enrichKEGG).bnpathplot.keyType of the enrichment results (e.g., convert using clusterProfiler::setReadable).mappedkeys to map pathway identifiers to genes when preparing test or custom pathway gene lists.returnNet = TRUE in bngeneplot to inspect the underlying network strength and direction before relying solely on the visual plot.clusterProfiler::setReadable to align key types (e.g., to "ENTREZID").pathview for overlaying expression data onto static KEGG pathway diagrams.enrichplot for standard enrichment visualizations (e.g., cnetplot, emapplot) without Bayesian network inference.bnlearn for general Bayesian network inference without direct integration with enrichment results.まだレビューはありません。使ってみた感想をお寄せください。
概要と使いどころ
The package contains functions for exploratory oligonucleotide array analysis. The dependence on tkWidgets only concerns few convenience functions. 'affy' is fully functional without it.
日本語の概要は準備中です。原文の説明を表示しています。
This tool uses the label-tree function from HyPhy to annotate a phylogenetic tree. It allows users to select a subset of leaves using either a regular expression or a list of sequence names, and then apply a specified label to these selecte
日本語の概要は準備中です。原文の説明を表示しています。
This package provides a client for the Bioconductor AnnotationHub web resource. The AnnotationHub web resource provides a central location where genomic files (e.g., VCF, bed, wig) and other resources from standard locations (e.g., UCSC, En
日本語の概要は準備中です。原文の説明を表示しています。
apeglm provides Bayesian shrinkage estimators for effect sizes for a variety of GLM models, using approximation of the posterior for individual coefficients.
日本語の概要は準備中です。原文の説明を表示しています。
ASURAT is a software for single-cell data analysis. Using ASURAT, one can simultaneously perform unsupervised clustering and biological interpretation in terms of cell type, disease, biological process, and signaling pathway activity. Inputting a single-cell RNA-seq data and knowledge-based databases, such as Cell Ontology, Gene Ontology, KEGG, etc., ASURAT transforms gene expression tables into original multivariate tables, termed sign-by-sample matrices (SSMs).
日本語の概要は準備中です。原文の説明を表示しています。
AUCell allows to identify cells with active gene sets (e.g. signatures, gene modules...) in single-cell RNA-seq data. AUCell uses the "Area Under the Curve" (AUC) to calculate whether a critical subset of the input gene set is enriched with
日本語の概要は準備中です。原文の説明を表示しています。