Description Usage Arguments Value References See Also Examples
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.
1 |
trees |
List of trees to analyse. |
TipVolatility()
returns a named vector listing the volatility
index calculated for each leaf.
Other tip instability functions:
TipInstability()
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.