joyplot for GSEA result

R-bloggers 2017-08-20

(This article was first published on R on Guangchuang YU, and kindly contributed to R-bloggers)

I am very glad to find that someone figure out how to use ggjoy with ggtree.

I really love ggjoy and believe it can be a good tool to visualize gene set enrichment (GSEA) result. DOSE/clusterProfiler support several visualization methods.

running score:

enrichment map:

category-gene-network:

dotplot:


These visualization methods are designed for visualizing/summarizing enrichment results and none of them can incoporate expression values.

In DOSE v>=3.3.2, I defined a joyplot function, which can visualize GSEA result with expression distribution of the enriched categories.

Here is an example:

require(clusterProfiler)data(geneList)x <- gseKEGG(geneList)joyplot(x)

var vglnk = { key: '949efb41171ac6ec1bf7f206d57e90b8' }; (function(d, t) {var s = d.createElement(t); s.type = 'text/javascript'; s.async = true;s.src = '//cdn.viglink.com/api/vglnk.js';var r = d.getElementsByTagName(t)[0]; r.parentNode.insertBefore(s, r); }(document, 'script'));

To leave a comment for the author, please follow the link and comment on their blog: R on Guangchuang YU.

R-bloggers.com offers daily e-mail updates about R news and tutorials on topics such as: Data science, Big Data, R jobs, visualization (ggplot2, Boxplots, maps, animation), programming (RStudio, Sweave, LaTeX, SQL, Eclipse, git, hadoop, Web Scraping) statistics (regression, PCA, time series, trading) and more...