zoomPathway | R Documentation |
Plots the gene to gene interactions for a given pathway in igraph.
zoomPathway(x, pway, graph_layout = NULL)
x |
A NetGSA object returned from calling |
pway |
Name of pathway to plot |
graph_layout |
(Optional) Layout function to pass to igraph plots. This function should only take one parameter (an igraph object). For example one might create a custom layout by setting the spring.length and spring.constant with: |
Generates igraph plot for gene to gene interactions for a given pathway
No return value, called for side effects
Michael Hellstern
Ma, J., Shojaie, A. & Michailidis, G. (2016) Network-based pathway enrichment analysis with incomplete network information. Bioinformatics 32(20):165–3174.
plot.NetGSA
## Not run:
## load the data
data("breastcancer2012_subset")
## consider genes from just 2 pathways
genenames <- unique(c(pathways[["Adipocytokine signaling pathway"]],
pathways[["Adrenergic signaling in cardiomyocytes"]]))
sx <- x[match(rownames(x), genenames, nomatch = 0L) > 0L,]
db_edges <- obtainEdgeList(rownames(sx), databases = c("kegg", "reactome"))
adj_cluster <- prepareAdjMat(sx, group, databases = db_edges, cluster = TRUE)
out_cluster <- NetGSA(adj_cluster[["Adj"]], sx, group,
pathways_mat[c(1,2), rownames(sx)], lklMethod = "REHE", sampling = FALSE)
plot(out_cluster)
my_layout <- function(graph) layout_with_graphopt(graph = graph,
spring.length = 1000,
spring.constant = 0.00004)
zoomPathway(out_cluster, "Adipocytokine signaling pathway", my_layout)
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.