powered by
Simulates a phylogeny arising from the SGD model with exponentially increasing metapopulation size. Notations follow Manceau et al. (2015).
sim_sgd(tau, b, d, nu)
a phylogenetic tree of class "phylo" (see ape documentation)
the simulation time, which corresponds to the length of the phylogeny
the (constant) per-individual birth rate
the (constant) per-individual death rate
the (constant) per-individual mutation rate
M Manceau
Manceau M., Lambert A., Morlon H. (2015) Phylogenies support out-of-equilibrium models of biodiversity Ecology Letters 18: 347-356
tau <- 10 b <- 1e6 d <- b-0.5 nu <- 0.6 tree <- sim_sgd(tau,b,d,nu) plot(tree)
Run the code above in your browser using DataLab