sim.coaltree: Simulate a coalescence tree

In lliu1871/phybase: Basic functions for phylogenetic analysis

Simulate a coalescence tree

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

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 [e-mail]

References

John Wakeley, Coalescent theory: An introduction.

See Also

sim.coaltree.sp

Examples

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

lliu1871/phybase documentation built on April 21, 2024, 3:16 a.m.