Dear Miguel & Eduardo.
This is correct, except that you should keep in mind that pbtree (like
most, but not all, such simulators) generates trees from the first
speciation event. Thus, the correct procedure must first involve drawing
a wait time from an exponential distribution, deciding if the next event
was a speciation or an extinction, attaching a stem to the base tree,
and then (if the event was a speciation) simulating a birth-death tree
to attach to this stem.
The algorithm would thus be something like:
1. Decide on a desired total height for the tree, h. (If, say, the tree
has a current total height/depth of 100, this could be any value >100.)
2. Go to a tip in the base tree. Compute it's height, t1.
3. Draw a random value from an exponential distribution, t2, with rate
b+d for birth rate b & death rate d.
4. If t1+t2Otherwise, add h-t1 to the current tip and go to 2 for the next tip.
5. Decide if the current event is a speciation or extinction event. Draw
a random value from a uniform distribution on (0,1). If that value is
greater than b/(b+d) you are done with this tip: go back to step 2 for
the next tip. Otherwise, advance to step 6.
6. Compute t3=h-(t1+t2). Simulate a birth-death tree for t3 with birth &
death rates b & d. Attach the b/d tree to the current tip. Advance to
the next tip & go back to step 2.
All the best, Liam
Liam J. Revell, Associate Professor of Biology
University of Massachusetts Boston
web: http://faculty.umb.edu/liam.revell/
On 7/30/2018 8:58 AM, Eduardo Ascarrunz wrote:
Maybe you could use a birth-death model. You can estimate the speciation
rate and extinction rate from your original tree, and use those parameters
to run simulations of trees for 5 My with the pbtree function from
phytools. Run as many simulations as you have living tips in your original
tree, and bind each simulated tree to one of the living tips. It shouldn't
be too difficult to write a wrapper function to perform all those steps.
Cheers,
Eduardo
2018-07-29 11:06 GMT+02:00 :
Hi list:
Is there any function to simulate a tree from a given tree?. For example,
I have a tree encompassing 10 Million years and I would like to simulate
the "growth" of such tree for the following 5 MY.
Thanks in advance
Miguel Verdú
--
Centro de Investigaciones sobre Desertificacion
(CSIC-UV-GV)
Carretera Moncada - Náquera, Km. 4,5
Apartado Oficial
46113 Moncada (Valencia)
Spain
Tel +34 96 3424204
Fax +34 96 3424160
www.uv.es/verducam
___
R-sig-phylo mailing list - R-sig-phylo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
Searchable archive at http://www.mail-archive.com/r-
sig-ph...@r-project.org/
[[alternative HTML version deleted]]
___
R-sig-phylo mailing list - R-sig-phylo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
Searchable archive at http://www.mail-archive.com/r-sig-phylo@r-project.org/
___
R-sig-phylo mailing list - R-sig-phylo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
Searchable archive at http://www.mail-archive.com/r-sig-phylo@r-project.org/
