On 7/14/11 2:45 AM, ppi...@uniroma3.it wrote:
Thankyou NIck,
now...I have an error when running it:
Oops. Apologies, this stuff is in-house code, I haven't
organized it. Add these to the text file...
===
get_daughters - function(nodenum, t)
{
daughter_edgenums =
Oops, left that out. Here it is as text and an attached file.
You will have to do
library(ape)
and maybe
library(phangorn)
...ahead of time (hopefully not others...)
chainsaw - function(tr, timepoint=10)
{
# Take a tree and saw it off evenly across a certain
Hi all,
someone knows a smart coding to cut a tree in order to
retain only taxa present in a given time interval bin?
Thanks in advance for any suggestion
best
paolo
___
R-sig-phylo mailing list
R-sig-phylo@r-project.org
I wrote a function called chainsaw to do something like
this, it saws off all the branches at a particular time
point, then you just have to drop.tip on branch tips older
than your time period. Would that be helpful?
On 7/12/11 3:19 PM, ppi...@uniroma3.it wrote:
Hi all,
someone knows a
Do you want a function that gives you a tree of taxa observed in a time bin
(like, in the fossil record, as in Ruta et al., 2007)? Or gives you the tree
of the lineages extant at point in time based on some global phylogeny? If
that is what you want, the following function I wrote for doing
Here's chainsaw(). It also requires sourcing a few other
functions. Please cite me if you use it :-).
=
chainsaw - function(tr, timepoint=10)
{
# Take a tree and saw it off evenly across a certain timepoint.
# This removes any tips above the