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)
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