Stemwardness | R Documentation |
Functions to describe the position of a leaf relative to the root. "Stemmier" leaves ought to exhibit a smaller root-node distance and a larger sister size.
SisterSize(tree, tip)
## S3 method for class 'numeric'
SisterSize(tree, tip)
## S3 method for class 'character'
SisterSize(tree, tip)
RootNodeDistance(tree, tip)
## S3 method for class 'numeric'
RootNodeDistance(tree, tip)
## S3 method for class 'character'
RootNodeDistance(tree, tip)
RootNodeDist(tree, tip)
tree |
A tree of class |
tip |
Either a numeric specifying the index of a single tip, or a character specifying its label. |
RootNodeDistance()
calculates the number of nodes between the chosen leaf
and the root of tree
.
This is an unsatisfactory measure, as the range of possible
distances is a function of the shape of the tree
\insertCiteAsher2020TreeTools.
As an example, leaf X1 in the tree (.,(.,(.,(.,(X1,(a,b))))))
falls outside the clade (a, b) and has a root-node distance of 4,
whereas leaf X2 in the tree (.,((.,(.,.)),(b,(X2,a))))
falls within the clade (a, b), so should be considered more "crownwards",
yet has a smaller root-node distance (3).
Stemwardness.pngTreeTools par(mfrow = c(1, 2), mar = rep(0.3, 4)) plot(ape::read.tree(text="(.,(.,(.,(.,(X1,(a,b))))));")) ape::nodelabels(1:4, 9:12) ape::edgelabels(1:2, 11:12)
plot(ape::read.tree(text="(.,((.,(.,.)),(b,(X2,a))));")) ape::nodelabels(1:3, c(9, 12, 13)) ape::edgelabels(1, 12)
SisterSize()
measures the number of leaves in the clade that is sister to
the chosen leaf, as proposed by \insertCiteAsher2020;textualTreeTools.
In the examples above, X1 has a sister size of 2 leaves, whereas X2,
which is "more crownwards", has a smaller sister size (1 leaf), as desired.
SisterSize()
returns an integer specifying the number of leaves
in the clade that is sister to tip
.
RootNodeDist()
returns an integer specifying the number of nodes between
tip
and the root node of tree
.
Martin R. Smith (martin.smith@durham.ac.uk)
Other tree characterization functions:
CladisticInfo()
,
Consensus()
,
J1Index()
,
TotalCopheneticIndex()
bal8 <- BalancedTree(8)
pec8 <- PectinateTree(8)
SisterSize(bal8, 3)
SisterSize(pec8, "t3")
SisterSize(RootTree(pec8, "t3"), "t3")
RootNodeDist(bal8, 3)
RootNodeDist(pec8, "t3")
RootNodeDist(RootTree(pec8, "t3"), "t3")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.