Function for phylogeny resolution

Scott Chamberlain




UPDATE: Yeah, so the treeresstats function had a problem in one of the calculations. I fixed that and added some more calulcations to the function.

I couldn't find any functions to calculate number of polytomies, and related metrics.

Here's a simple function that gives four metrics on a phylo tree object:

Here's output from the gist above:

$trsize_tips
[1] 15

$trsize_nodes
[1] 13

$numpolys
[1] 1

$numpolysbytrsize_tips
[1] 0.06666667

$numpolysbytrsize_nodes
[1] 0.07692308

$proptipsdescpoly
[1] 0.2

$propnodesdich
[1] 0.9230769

And an example with many trees:

trsize_tips trsize_nodes numpolys numpolysbytrsize_tips numpolysbytrsize_nodes proptipsdescpoly propnodesdich
20 13 4 0.20 0.31 0.7 0.69
20 7 3 0.15 0.43 0.9 0.57
20 11 6 0.30 0.55 1.0 0.45
20 13 4 0.20 0.31 0.7 0.69
20 9 5 0.25 0.56 1.0 0.44

Posted in  R ape phylogenetics


blog comments powered by Disqus
Fork me on GitHub