BDsim {geiger}R Documentation

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

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.

Value

Population size at each time in TIMES

Author(s)

Richard E. Glor and Luke J. Harmon

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


[Package geiger version 1.2-14 Index]