geiger (version 0.2-6)

BDsim: Birth-death population simulator

Description

Simulates species richness (or population growth) under a uniform 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

References

Geiger

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