simulateBirthDeathRich {evolved} | R Documentation |
Simulating richness through birth-death processes
Description
simulateBirthDeathRich
calculates the number of species at a certain
point in time, following a birth-death process.
Usage
simulateBirthDeathRich(t, S = NULL, E = NULL, K = NULL, R = NULL)
Arguments
t |
Point in time which richness will be simulated. |
S |
A numeric representing the per-capita speciation rate (in number
of events per lineage per million years). Must be larger than |
E |
A numeric representing the per-capita extinction rate (in number
of events per lineage per million years). Must be smaller than |
K |
A numeric representing the extinction fraction (i.e.,
|
R |
A numeric representing the per-capita Net Diversification
rate (i.e., |
Details
The function only accepts as inputs S
and E
, or
K
and R
.
Value
The number of simulated species (i.e., the richness).
Author(s)
Matheus Januario, Daniel Rabosky, Jennifer Auler
References
Raup, D. M. (1985). Mathematical models of cladogenesis. Paleobiology, 11(1), 42-52.
Examples
# running a single simulation:
SS <- 0.40
EE <- 0.09
tt <- 10 #in Mya
simulateBirthDeathRich(t = tt, S = SS, E = EE)
#running many simulations and graphing results:
nSim <- 1000
res <- vector()
for(i in 1:nSim){
res <- c(res,
simulateBirthDeathRich(t = tt, S = SS, E = EE))
}
plot(table(res)/length(res),
xlab="Richness", ylab="Probability")