bd.ms: estimate net diversification rate

Description Usage Arguments Details Value Author(s) References Examples

View source: R/diversification.R

Description

estimating net diversification rate with confidence limits and testing diversities

Usage

1
2
3
4
5
6
bd.ms(phy=NULL, time, n, missing = 0, crown=TRUE, epsilon = 0)
bd.km(phy=NULL, time, n, missing = 0, crown=TRUE)
crown.p(phy=NULL, time, n, r, epsilon)
stem.p(phy=NULL, time, n, r, epsilon)
crown.limits(time, r, epsilon, CI=0.95)
stem.limits(time, r, epsilon, CI=0.95)

Arguments

phy

a phylogenetic tree of class 'phylo' (see Details)

time

time interval (can be a vector)

n

number of extant species

r

net diversification rate, birth - death

epsilon

extinction rate as a fraction of speciation rate

missing

number of taxa missing from tree

crown

whether time is treated as crown age (or otherwise as stem age)

CI

confidence level for estimated parameters

Details

bd.ms uses the Magallon and Sanderson (2000) method to calculate net diversification rate for a clade given extant diversity and age. bd.km computes the Kendall-Moran estimate of speciation rate, which assumes a complete phylogenetic tree.

Associated functions crown.p and stem.p also calculate the probability of obtaining a clade with at least n species given a net diversification rate (r), extinction fraction (epsilon), and time interval. Associated functions stem.limits and crown.limits generate confidence limits on extant diversity given a net diversification rate (r), extinction fraction (epsilon), and time interval.

Where a function calls for a time and an n element, a tree may be given instead (as argument phy). The argument n is taken from the number of tips in the tree. The method will attempt to discern whether the model should be fitted assuming crown or stem. If the tree has a non-NULL phy$root.edge element, the length will be assumed for the stem and crown is assumed to be FALSE (see also read.tree).

Value

Author(s)

LJ Harmon and C Brock

References

Magallon S and MJ Sanderson. 2000. Absolute diversification rates in angiosperm clades. Evolution 55:1762-1780.

Examples

1
2
3
4
5
6
7
geo=get(data(geospiza))

# Assuming no extinction
bd.ms(phy=geo$phy, missing=1, epsilon=0)

# Assuming high extinction
bd.ms(phy=geo$phy, missing=1, epsilon=0.9)

geiger documentation built on Jan. 25, 2020, 5:06 p.m.