# Pure-birth tree
p1<-birthdeath.tree(b=0.1, d=0, time.stop=20)
plot(p1)
# Birth-death tree with extinct taxa
p2<-birthdeath.tree(b=0.2, d=0.05, time.stop=20)
plot(p2)
# Previous tree with extinct taxa removed
prune.extinct.taxa(p2)->p3
if(!is.null(p3)) plot(p3) else cat("No survivors
")
#Sometimes there are no extant taxa
Run the code above in your browser using DataLab