Description Usage Arguments Value Author(s) References See Also Examples
The function simulates a gene tree from the species tree using Rannala and Yang's formula. It uses the same population size for all branch lengths. If you want more control, use sim.coaltree.sp().
1 | sim.coaltree.phylo(phy, pop.size=1, nsamples=1)
|
phy |
a tree in phylo format |
pop.size |
the population size in coalescent time units on each branch |
nsamples |
the number of allele samples per species |
the gene tree in phylo format generated from the species tree
Brian O'Meara omeara.brian@gmail.com
Rannala, B. and Z. Yang. 2003. Bayes estimation of species divergence times and ancestral population sizes using DNA sequences from multiple loci. Genetics 164: 1645-1656.
1 2 | species.tree <- rcoal(7)
gene.tree <- sim.coaltree.phylo(species.tree, nsamples=2)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.