Description Usage Arguments Value References Examples
View source: R/range.calculations.tree.wide.R
By default the function returns tree-wide estimates of speciation and extinction rates, which are calculated by summing taxon counts across all intervals. This approach relies on having equal length intervals.
1 2 | boundary.crosser.rates(fossils, max.age, strata, continuous = TRUE,
return.intervals = FALSE)
|
fossils |
Dataframe of sampled fossils (sp = unique species ID. h = ages.) |
max.age |
Maximum age of the oldest stratigraphic interval |
strata |
Number of stratigraphic intervals |
continuous |
If TRUE calculate continuous rates (i.e. account for interval length) |
return.intervals |
If TRUE return per interval estimates |
named list with the overall speciation rate, overall extinction rate and a dataframe of per interval estimtes if return.intervals = TRUE. Note this approach does not return rates for the first interval.
Foote, M. 2000. Origination and extinction components of taxonomic diversity: General problems. Paleobiology 26: 74-102.
1 2 3 4 5 6 7 8 9 10 11 | # simulate tree & fossils
t = TreeSim::sim.bd.taxa(100,1,1,0.5)[[1]]
# budding speciation
sp = FossilSim::sim.taxonomy(t)
# simulate fossils
max = FossilSim::tree.max(t)
f = FossilSim::sim.fossils.intervals(taxonomy = sp, max.age = max, strata = 10, probabilities = rep(0.5,10), use.exact.times = FALSE)
# add extant occurrences
f = FossilSim::sim.extant.samples(f, taxonomy = sp)
# calculate speciation and extinction rates
boundary.crosser.rates(f, max, 10)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.