Sunday, June 16, 2013

To plot a nifty plot with the rates:

ave.rates(tree.trim,min_split,extract.clade(tree.trim, node=min_split$node)$tip.label,colMeans(mcmc)["sig1"], colMeans(mcmc)["sig2"],min_split) 

No comments:

Post a Comment