geiger (version 1.0-91)

BDsim: Birth-death population simulator

Description

Simulates species richness (or population growth) under a uniform, time-homogeneous birth-death process.

Usage

BDsim(nStart, b, d, times)

Arguments

nStart

Number of taxa at starting time zero

b

Per-lineage birth (speciation) rate

d

Per-lineage death (extinction) rate

times

Vector of times where extant species are counted

Value

• Population size at each time in TIMES

Details

This function simulates species diversification under a uniform birth-death process. This differs from birthdeath.tree in that only the number of species, and not their phylogenetic affinities, are stored. This function relates to GEIGER's rate.estimate (and associated functions), which are also non-phylogenetic.

References

Yule, G. U. 1924. A mathematical theory of evolution based on the conclusions of Dr. J. C. Willis, FRS. Philos. Trans. R. Soc. London Ser. B 213: 21-87

See Also

rate.estimate, stem.p, crown.p, stem.limits, crown.limits

Examples

pop1=BDsim(nStart=10, b=0.1, d=0, times=1:10)
pop2=BDsim(nStart=10, b=0, d=0.1, times=1:10)
pop3=BDsim(nStart=10, b=0.1, d=0.1, times=1:10)

plot(pop1, type="l", ylim=c(0,20))
lines(pop2, col="red")
lines(pop3, col="blue")