Dear List
I am trying to plot two trees facing each other with pie charts at the nodes
indicating marginal likelihoods of an ancestral state. I have broken the
plotting frame into 3 sections, with one in the centre for the tip labels. I am
using text(rep()) to plot the tip labels into this
The way I've handled this is probably inefficient, but I get the labels
to match up by writing the tree as I have plotted it, then reading that
back in: then, the tip labels are in order with the plotted tips, bottom
to top:
a = rtree(20)
plot(a)
a$tip.label ## labels match, bottom to top
Dear list.
I am advertising for a postdoctoral position in my lab at the University
of Massachusetts Boston. Details below. Please contact me for more
information.
- Liam
--
Postdoctoral research associate in phylogenetic comparative methods
A postdoctoral
Hi Henry,
I suppose the two trees do not have the same topology (otherwise there
would be no difficulty). Have you considered using rotateConstr() on the
second tree? That'd be something like:
tree2 - rotateConstr(tree2, tree1$tip.label)
The node numbering won't be changed, so you can still
I tried Brownian rate parameter estimation by tranforming the branch
lengths according to the lambda values of trait in each group. It gives
me lower varaince for the group that has higher lambda value and that is
confusing. In fact with the same trait values, if I transform the branch