Dear List Members.
*Does anyone know of any code or function that can help me easily set arbitrary branch lengths according to Nee and Pagel's methods?* I'm trying to analyze trait data on a phylogeny with no branch lengths. Based on archived posts, I've determined that I'd like to experiment with setting arbitrary branch lengths of the following types: 1) All equal (1) R CODE: compute.brlen(tree, 1) 2) Grafen's (1989): branch lengths set to a length equal to the number of descendant tips minus one. R CODE: compute.brlen(tree, method="Grafen") 3) Nee's Method (Purvis 1995): branch lengths set such that the distance from the tips to the current node is equal to the log 10 of the number of tips descending from that node R CODE: Haven't found anything. 4) Pagel's (1992): branch lengths set such that the depth of a node is equal to the maximum of the number of bifurcation levels of the left child and the right child plus one. R CODE: Haven't found anything. Thanks! Brian -- Brian A. Gill VISIT MY WEBSITE: http://gillbriana.wix.com/brian-gill FOLLOW ME ON TWITTER: @CSUBrianGill Colorado State University Biology 1878 Campus Delivery Fort Collins, CO 80523 970-215-7037 [[alternative HTML version deleted]] _______________________________________________ R-sig-phylo mailing list - [email protected] https://stat.ethz.ch/mailman/listinfo/r-sig-phylo Searchable archive at http://www.mail-archive.com/[email protected]/
