data(plethodon)
# suppose plethodon basal divergence occurred 23 mya:
svec <- scaleBranchingtimes(plethodon, basal = 23)
plotLtt(svec)
# plots lineages through time in units of time, rather than genetic
# distance
pureBirth(svec)
# returns speciation rate in units of lineages/million years
Run the code above in your browser using DataLab