Hi Rafael. When I load your workspace & run the following code I can
produce the attached plot. Does this not work for you?
library(phytools)
dat<-dat[!is.na(dat)]
tips<-setdiff(tr$tip.label,names(dat))
tr<-drop.tip(tr,tips)
pdf(file="tree.pdf")
plotTree.wBars(tree=tr, x=dat, lwd=1, type='fan',
Hi Liam,
That did solve that problem and now a tree is plotted. But it is just a
plain tree without bars at the tips. Here's what I'm doing and the error I
get:
> plotTree.wBars(tree=tr, x=dat, type='fan', scale=0.05, width=1, col='red')
Error in if (min(x) < 0) h <- max(nodeHeights(tree)) :