Learn R Programming

BioGeoBEARS (version 0.2.1)

average_tr_tips: Average the heights of (non-fossil) tips to make ultrametric-ish.

Description

When you have a digitized tree, or other slightly uneven source tree, average the tips to get them all to line up at 0 my before present. This makes an ultrametric tree if and only if there are no fossil tips in the tree.

Usage

average_tr_tips(tr, fossils_older_than = 0.6)

Arguments

tr
An ape phylo object
fossils_older_than
Tips that are older than fossils_older_than will be excluded from the tips that are going to be averaged. This is not currently set to 0, because Newick files can have slight precision issues etc. that mean not all tips quite come to zero (which is why you need average_tr_tips in the first place!). Obviously you should be cautious about the value of , depending on the absolute timescale of your tree. Make sure you do not inappropriately average in fossils!!

Value

edge_times_bp A 2-column matrix with the age (from the present) of the top and bottom of each edge.

Details

If the user includes fossils accidentally, this function can easily lead to pathological results (negative branch lengths etc.), so use with care!!

References

http://phylo.wikidot.com/matzke-2013-international-biogeography-society-poster

Matzke_2012_IBS

See Also

prt, chainsaw2, extend_tips_to_ultrametricize

Examples

Run this code
test=1

Run the code above in your browser using DataLab