sim.coaltree | R Documentation |
This function can simulate a coalescence tree from a single population with parameter theta. The coalescence times in the tree have exponential distributions. theta
is equal to 4uNe where Ne is the effective population size and u is the mutation rate.
sim.coaltree(nspecies,theta)
nspecies |
the number of species |
theta |
the population parameter |
theta is the population parameter theta=4N*mu.
The function returns the simulated coalescence tree.
Liang Liu lliu@uga.edu
John Wakeley, Coalescent theory: An introduction.
sim.coaltree.sp
sim.coaltree(5,theta=0.2)
##[1] "((5:0.55696,(1:0.34858,3:0.34858):0.20838):2.99874,(2:0.97896,4:0.97896):2.57674)"
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.