BDsim {geiger} | R Documentation |
Simulates species richness (or population growth) under a uniform birth-death process
BDsim(nStart, b, d, times)
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 |
Population size at each time in TIMES
Richard E. Glor and Luke J. Harmon
Geiger
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")