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 <miguel.ve...@ext.uv.es>: > 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/
