how about this trick:

library(ggtree)
set.seed(1234)
tr=rtree(50)

geom_strip = function(...) geom_cladelabel(label="", barsize=5, align=T, offset=.5, ...)

ggtree(tr) + geom_strip(58, color='red') + geom_strip(84, color='green') + geom_strip(52, color='blue')

This sort of addresses the problem. However I am still not sure how to deal with this as my data relates tip labels to taxonomic names (at the phyla level). Thus I first have to find the MRCA of a given phyla however often its not mono-phyletic. Therefore the results are incorrect.

It would be easier if I could use the %<+% opperator to associate the phyla names to the tips and then just call a geom_strip function.

Do you have any suggestions?

So I found a workaround - a package called MonoPhy that will estimate the MRCA despite non monophyletic clades. So I can use that... Its a pretty big work around though, it would be much easier if you could just associate a color with the tips as a bar.

it's supported in ggtree >= 1.5.1

it's similar to geom_cladelabel, except that:

• the input is not internal node number but taxa1 and taxa2.
• label is optional.

It works for Monophyletic, Polyphyletic or Paraphyletc Taxa, so we can think of it as a super-geom_cladelabel.

wow, just realize you request a pull. Thanks, I will review your code tomorrow.

I think I like your solution with geom_strip much better for my research! That said the pull request may be helpful as it fixes what i would consider a small bug with gheatmap function.

Thank you again for the wonderful work you are doing!

• Justin

On May 9, 2016, at 11:58 AM, Guangchuang Yu notifications@github.com wrote:

wow, just realize you request a pull. Thanks, I will review your code tomorrow.

—
You are receiving this because you authored the thread.
Reply to this email directly or view it on GitHub #52 (comment)

merged & thanks