Calculate the volatility of each tip: namely, the impact on the mean phylogenetic information distance \insertCiteSmith2020Rogue between trees when that tip is removed. Effective when the number of trees is small.
List of trees to analyse.
TipVolatility() returns a named vector listing the volatility
index calculated for each leaf.
Other tip instability functions:
1 2 3 4 5 6 7 8 9
library("TreeTools", quietly = TRUE) trees <- AddTipEverywhere(BalancedTree(8), 'Rogue') trees <- lapply(trees, AddTip, 'Rogue', 'Rogue2') sb <- TipVolatility(trees) sbNorm <- 1 + (99 * (sb - min(sb)) / (max(sb - min(sb)))) col <- hcl.colors(128, 'inferno')[sbNorm] plot(consensus(trees), tip.color = col) plot(ConsensusWithout(trees, names(sb[sb == max(sb)])))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.