aldous.test {apTreeshape} | R Documentation |
A graphical test to decide if tree data fit the Yule or the PDA models.
aldous.test(tree, xmin=20, ...)
tree |
An object of class |
xmin |
An object of class |
... |
further arguments to be passed to |
A binary tree contains a set of splits
(m,i)=(size of parent clade, size of smaller daughter clade)
which can be plotted as a scatter diagram. Aldous' proposal for studying tree balance is that, given a large phylogenetic tree, one should estimate the median size of the smaller daughter clade as a function of the parent clade and use this function as a descriptor of balance or imbalance of the tree. It is convenient to make a log-log plot and to ignore small parent clades. The scatter diagram shows lines giving the approximate median values of the size of smaller daughter clade predicted by the beta-splitting model for two values of beta, the value for the Yule (beta=0) and PDA (beta=-1.5) models. In other words, if the null model were true, then the scatter diagram for a typical tree would have about half the points above the line and half below the line, throughout the range.
The green line represents the median regression estimated from the tree data.
The function provides a graphical display of results.
Michael Blum michael.blum@imag.fr
Nicolas Bortolussi nicolas.bortolussi@imag.fr
Eric Durand eric.durand@imag.fr
Olivier Francois olivier.francois@imag.fr
Aldous, D. J. (1996) Probability Distributions on Cladograms. pp.1-18 of Random Discrete Structures eds D. Aldous and R. Pemantle, IMA Volumes Math. Appl. 76.
Aldous, D. J. (2001) Stochastic Models and Descriptive Statistics for Phylogenetic Trees, from Yule to Today. Statistical Science, 16, 23 – 24.
library(quantreg) aldous.test(rbiased(2000, p=.5)) ## Test with a huge balanced tree: aldous.test(rbiased(2000, p=.5))