Learn R Programming

geiger (version 0.2-5)

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

Run this code
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")

Run the code above in your browser using DataLab