Description Usage Arguments Details Value Author(s) See Also Examples
This function prints a compact summary of a phylogenetic tree of class "treeshape"
)
1 2 3 |
object |
an object of class |
... |
further arguments passed to or from other methods. |
summary.treeshape
prints the following information: the number of tips of the tree, its Colless' index, and the expected values and standard deviations of the Colless' index under the PDA and Yule models.The expected value of the Colless' index under the Yule model is given according to the formula: n*log(n)+n*(gamma-1-log(2)) where n is the number of tips of the tree and gamma the Euler's constant. The standard deviation under the Yule model is given by: sqrt((3-pi^2/6-log(2))*n). The expected value of the Colless' index under the PDA model is given according to the formula: sqrt(pi)*n^(3/2). The standard deviation under the PDA model is given by: sqrt(10/3-pi)*n^(3/2).
A NULL value is returned, the results are simply printed.
Michael Blum <michael.blum@imag.fr>
Nicolas Bortolussi <nicolas.bortolussi@imag.fr>
Eric Durand <eric.durand@imag.fr>
Olivier Francois <olivier.francois@imag.fr>
summary
colless
for more informations about the expected values under the Yule and PDA models.
1 2 3 4 5 6 7 8 9 10 11 12 13 | ## Summary of a PDA tree with 100 tips.
summary(rpda(100))
## Note that the standard deviation is very large.
## Summary of a Yule tree with 100 tips.
summary(ryule(100))
## The standard deviation under the Yule model is much smaller than under
## the PDA model.
## Summary of the HIV tree.
data(hivtree.treeshape)
summary(hivtree.treeshape)
## The HIV tree is much closer from the Yule model than from the PDA model.
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.