The phytools function plot.cophylo
will no plot an object
of class "cophylo"
containing trees without edge lengths
(upon user request).
Here's what I mean:
library(devtools)
install_github("liamrevell/phytools")
Now try it….
library(phytools)
t1<-rtree(n=26,tip.label=LETTERS,br=NULL)
t2<-rtree(n=26,tip.label=LETTERS,br=NULL)
obj<-cophylo(t1,t2)
## Rotating nodes to optimize matching...
## Done.
plot(obj)
This was super easy. I just plugged the line:
if(is.null(tree$edge.length)) tree<-compute.brlen(tree)
into the internally used plotting function, phylogram
, and
that was all there was too it.
No comments:
Post a Comment
Note: due to the very large amount of spam, all comments are now automatically submitted for moderation.