bd.probability.range: Birth-death probability (Keiding, 1975)
In rachelwarnock/fbdR: Calculate Birth and Death Rates

Description Usage Arguments Value References Examples

Calculate the probabilty of speciaton and extinction rates for a given set of stratigraphic ranges

1	bd.probability.range(frs, b, d, crown = FALSE, tol = NULL)

`frs`	Dataframe of species ranges
`b`	Rate of speciation
`d`	Rate of extinction
`crown`	If TRUE assume the process begins at the first speciation event and not the origin (default = FALSE)
`tol`	rounding-error tolerance for extant taxa

Log likelihood

Keiding, N. 1975. Maximum likelihood estimation in the birth-death process. Annals of Statistics 3: 363-372.

# simulate tree & assume complete sampling
t = TreeSim::sim.bd.taxa(100,1,1,0.1)[[1]]
# simulated mixed speciation
beta = 0.5
lambda.a = 0.1
frs <- FossilSim::sim.taxonomy(tree = t, beta = beta, lambda.a = 0.1)
# calculate birth-death probability
birth = 1
death = 0.1
bd.probability.range(frs, birth, death)