TreePar
3.3Estimating birth and death rates based on phylogenies
Overview
(i) For a given species phylogeny on present day data which is calibrated to calendar-time, a method for estimating maximum likelihood speciation and extinction processes is provided. The method allows for non-constant rates. Rates may change (1) as a function of time, i.e. rate shifts at specified times or mass extinction events (likelihood implemented as LikShifts, optimization as bd.shifts.optim and visualized as bd.shifts.plot) or (2) as a function of the number of species, i.e. density-dependence (likelihood implemented as LikDD and optimization as bd.densdep.optim) or (3) extinction rate may be a function of species age (likelihood implemented as LikAge and optimization as bd.age.optim.matlab). Note that the methods take into account the whole phylogeny, in particular it accounts for the "pull of the present" effect. (1-3) can take into account incomplete species sampling, as long as each species has the same probability of being sampled. For a given phylogeny on higher taxa (i.e. all but one species per taxa are missing), where the number of species is known within each higher taxa, speciation and extinction rates can be estimated under model (1) (implemented within LikShifts and bd.shifts.optim with groups !=0). (ii) For a given phylogeny with sequentially sampled tips, e.g. a virus phylogeny, rates can be estimated under a model where rates vary across time using bdsky.stt.optim based on likelihood LikShiftsSTT (extending LikShifts and bd.shifts.optim). Furthermore, rates may vary as a function of host types using LikTypesSTT (multitype branching process extending functions in R package diversitree). This function can furthermore calculate the likelihood under an epidemiological model where infected individuals are first exposed and then infectious.
Install
Health
CRAN check results are not tracked yet.
Documentation
- Examples that run
- 100%
- Documented parameters
- 88%
- Return-value docs
- 88%
- References docs
- 94%
Downloads
Dependencies
Nothing depends on this yet.
Code & Tests
People & History
14 releases. Pick two to compare their code metrics. R releases are shown for context.
- RR 4.3.0 released · 2023-04-21
- archivedRemoved from CRAN2022-04-27check issues were not corrected in time length > 1 in && or ||
- RR 4.2.0 released · 2022-04-22
- RR 4.1.0 released · 2021-05-18
- RR 4.0.0 released · 2020-04-24
- RR 3.6.0 released · 2019-04-26
- RR 3.5.0 released · 2018-04-23
- RR 3.4.0 released · 2017-04-21
- RR 3.3.0 released · 2016-05-03
- RR 3.2.0 released · 2015-04-16
- 3.32015-01-02 · diff ↗
- 3.22014-11-01 · diff ↗
- 3.12014-09-18 · diff ↗
- 3.02014-07-28 · diff ↗
- RR 3.1.0 released · 2014-04-10
- 2.52013-05-07 · diff ↗
Show 15 earlier events
- RR 3.0.0 released · 2013-04-03
- 2.32012-08-05 · diff ↗
- RR 2.15.0 released · 2012-03-30
- 2.22011-11-01 · diff ↗
- RR 2.14.0 released · 2011-10-31
- 2.1.22011-08-22 · diff ↗
- 2.12011-07-12 · diff ↗
- 1.42011-05-09 · diff ↗
- RR 2.13.0 released · 2011-04-13
- 1.32011-03-21 · diff ↗
- 1.22011-03-20 · diff ↗
- RR 2.12.0 released · 2010-10-15
- 1.12010-06-17 · diff ↗
- 1.02010-06-17
- RR 2.11.0 released · 2010-04-22
Package metadata
- Total releases
- 14
- License
- GPL-2 OSI
- Download size
- not tracked yet
- Installed size
- not tracked yet
- With dependencies
- not tracked yet
Cite
Cite this package
Run in R for the authors' preferred citation:
citation("TreePar")This is what citation() produces when a package has no citation file of its own. If it prints something else, use that.
Cite the R Observatory
For a number measured here: a download total, a coverage figure, an archival date.
From data release v2026-08-11, which the citation names so these numbers can be found later. More on citing and the projects behind them.