Function for phylogeny resolution

January 13, 2012
By

(This article was first published on Recology - R, and kindly contributed to R-bloggers)

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_tipstrsize_nodesnumpolysnumpolysbytrsize_tipsnumpolysbytrsize_nodesproptipsdescpolypropnodesdich
201340.200.310.70.69
20730.150.430.90.57
201160.300.551.00.45
201340.200.310.70.69
20950.250.561.00.44

To leave a comment for the author, please follow the link and comment on his blog: Recology - R.

R-bloggers.com offers daily e-mail updates about R news and tutorials on topics such as: visualization (ggplot2, Boxplots, maps, animation), programming (RStudio, Sweave, LaTeX, SQL, Eclipse, git, hadoop, Web Scraping) statistics (regression, PCA, time series, trading) and more...



If you got this far, why not subscribe for updates from the site? Choose your flavor: e-mail, twitter, RSS, or facebook...

Comments are closed.