sim.coaltree: Simulate a coalescence tree

Description Usage Arguments Details Value Author(s) References See Also Examples

Description

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.

Usage

1
sim.coaltree(nspecies,theta)

Arguments

nspecies

the number of species

theta

the population parameter

Details

theta is the population parameter theta=4N*mu.

Value

The function returns the simulated coalescence tree.

Author(s)

Liang Liu lliu@uga.edu

References

John Wakeley, Coalescent theory: An introduction.

See Also

sim.coaltree.sp

Examples

1
2
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)"

bomeara/phybase documentation built on May 12, 2019, 11:35 p.m.