# Recent Publications

More Publications

• ggtree: an R package for visualization and annotation of phylogenetic trees with their covariates and other associated data

• ReactomePA: an R/Bioconductor package for reactome pathway analysis and visualization

• ChIPseeker: an R/Bioconductor package for ChIP peak annotation, comparison and visualization

• DOSE: an R/Bioconductor package for disease ontology semantic and enrichment analysis

• Putative cobalt- and nickel-binding proteins and motifs in Streptococcus pneumoniae

• A novel andrographolide derivative AL-1 exerts its cytotoxicity on K562 cells through a ROS-dependent mechanism

• LXtoo: an integrated live Linux distribution for the bioinformatics community

• clusterProfiler: an R Package for Comparing Biological Themes Among Gene Clusters

• Functional similarity analysis of human virus-encoded miRNAs

• A new method for measuring functional similarity of microRNAs

# Posts

### scatterpie for plotting pies on ggplot

Plotting pies on ggplot/ggmap is not an easy task, as ggplot2 doesn’t provide native pie geom. The pie we produced in ggplot2 is actually a barplot transform to polar coordination. This make it difficult if we want to produce a map like the above screenshot, which was posted by Tyler Rinker, the author of R package pacman.

### align genomic features with phylogenetic tree

A question on biostars asking how to generate the following figure:

This can be quite easy to implement in ggtree, I can write a geom layer to layout the alignment. As ggbio already provides many geom for genomic data and I don’t want to re-invent the wheel, I decided to try ggtree+ggbio. This is also the beauty of R that packages complete each others.

### showCategory parameter for visualizing compareCluster output

I am using dotplot() to visualize results from enrichGO(), enrichDO(), enricher() and compareCluster() in clusterProfiler R package. When specifying showCategory, I get the right number of categories except with the results of compareCluser().

In my case, I use compareCluster() on a list of 3 elements:

str(ClusterList)
List of 3
$All : chr [1:1450] "89886" "29923" "100132891" "101410536" ...$ g1  : chr [1:858] "89886" "29923" "100132891" "101410536" ...
\$ g2: chr [1:592] "5325" "170691" "29953" "283392" ...

dotplot(CompareGO_BP, showCategory=10, title="GO - Biological Process")


I ask for 10 categories, but I get 15 categories in All, 8 categories in g1 and 12 categories in g2. None of the categories, neither the sum of the categories are 10…

Is the option showCategory working in the case of comparison? Am I missing something here?

And which categories precisely will it plot? the most significant whatever my 3 cases or the most significant of each case?

The question was posted in Bioconductor support site. It seems quite confusing and I think I need to write a post to clarify it.

### xlim_tree: set x axis limits for only Tree panel

A ggtree user recently asked me the following question in google group:

I try to plot long tip labels in ggtree and usually adjust them using xlim(), however when creating a facet_plot xlim affects all plots and minimizes them.

Is it possible to work around this and only affect the tree and it’s tip labels leaving the other plots in facet_plot unaffected?

This is indeed a desire feature, as ggplot2 can’t automatically adjust xlim for text since the units are in two different spaces (data and pixel).

### facet_plot: a general solution to associate data with phylogenetic tree

ggtree provides gheatmap for visualizing heatmap and msaplot for visualizing multiple sequence alignment with phylogenetic tree.

We may have different data types and want to visualize and align them with the tree. For example, dotplot of SNP site (e.g. using geom_point(shape='|')), barplot of trait values (e.g. using geom_barh(stat='identity')) et al.

To make it easy to associate different types of data with phylogenetic tree, I implemented the facet_plot function which accepts a geom function to draw the input data.frame and display it in an additional panel.

# Teaching

I am a teaching instructor for the following courses at University X: