Adjust branch lengths with node ages: comparison of two methods
[This article was first published on Recology, and kindly contributed to R-bloggers]. (You can report issue about the content on this page here)
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.
Here is an approach for comparing two methods of adjusting branch lengths on trees: bladj in the program Phylocom and a fxn written by Gene Hunt at the Smithsonian.Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.
Get the code and example files here: http://wp.me/PRT1F-2v
Get phylocom here: http://www.phylodiversity.net/phylocom/
The code uses as input a newick tree file and a text file of node ages.
# This is where the work happens...
# Directory below needs to have at least three items:
# 1. phylocom executable for windows or mac
# 2. tree newick file
# 3. node ages file as required by phylocom, see their manual
# Output: trees_out is a list of three trees, the original, bladj, and Gene Hunt's method
# Also, within the function all three trees are written to file as PDFs
setwd("/Mac/R_stuff/Blog_etc/Bladjing")
trees_out <- AdjBrLens("tree.txt", "nodeages.txt")
# plot trees of three methods together,
# with nodes with age estimates labeled
layout(matrix(1:3, 1, 3))
plot(trees_out[[1]])
nodelabels(trees_out[[1]]$node.label, cex = 0.6)
title("original tree")
plot(trees_out[[2]])
nodelabels(trees_out[[2]]$node.label, cex = 0.6)
title("bladj method")
plot(trees_out[[3]])
nodelabels(trees_out[[3]]$node.label, cex = 0.6)
title("gene hunt method, scalePhylo")It is sort of hard to see the differences in the branch length changes here, but the individual output trees will reveal the differences better.
To leave a comment for the author, please follow the link and comment on their blog: Recology.
R-bloggers.com offers daily e-mail updates about R news and tutorials about learning R and many other topics. Click here if you're looking to post or find an R/data-science job.
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.
