This function fits a birth-death model to a phylogenetic tree with edge lengths.
fit.bd(tree, b=NULL, d=NULL, rho=1, ...)
fit.yule(tree, b=NULL, d=NULL, rho=1, ...)
lik.bd(theta, t, rho=1, N=NULL)
# S3 method for fit.bd
print(x, ...)
object of class "phylo"
.
birth (speciation) rate. Presently doesn't do anything as the rate cannot be fixed.
death (extinction) rate. Presently doesn't do anything as the rate cannot be fixed.
sampling fraction.
vector of b
and d
for likelihood function.
branching times for calculation of the likelihood.
number of tips in the tree.
object of class "fit.bd"
for print
method.
optional arguments.
fit.bd
generates an object of class "fit.bd"
which can be printed.
This function duplicates birthdeath
in ape and make.bd
in the diversitree package.
Nee, S., May, R. M. and Harvey, P. H. (1994) The reconstructed evolutionary process. Philosophical Transactions of the Royal Society of London B, 344, 305-311.
Stadler, T. (2012) How can we improve the accuracy of macroevolutionary rate estimates? Systematic Biology, 62, 321-329.
Revell, L. J. (2012) phytools: An R package for phylogenetic comparative biology (and other things). Methods Ecol. Evol., 3, 217-223.