View source: R/simulatePhylogenies.R
createPhyloBiSSE | R Documentation |
The BiSSE diversification model has 6 parameters:
λ_0 speciation rate of state 0
λ_1 speciation rate of state 1
μ_0 extinction rate of state 0
μ_1 extinction rate of state 1
q_{01} transition rate from state 0 to 1
q_{10} transition rate from state 1 to 0
createPhyloBiSSE(n.tips, params)
n.tips |
number of tips |
params |
list
|
We assume 4 constraints:
λ_1 = 2λ_0
μ_0 = 0
μ_1 = 0
q_{10} = q_{01}
phylogeny (ape format)
n.tips <- 10 # 10 extant species lambda0 <- 1 # speciation rate q <- 0.1 # transition rate createPhyloBiSSE(n.tips, list(lambda0 = lambda0, q = q))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.