`BAMM`

output`dtRates`

calculates the mean of the marginal posterior
density of the rate of speciation/extinction or trait evolution for
small segments along each branch in a phylogeny.

1 |

`ephy` |
An object of class |

`tau` |
A numeric that specifies the size (as a fraction of tree height) of the segments that each branch will be discretized into. |

`ism` |
An integer vector indexing which posterior samples to include in the calculation. |

`tmat` |
A logical. If |

`dtRates`

bins the phylogeny into windows of time and
calculates average rates of speciation/extinction or phenotypic
evolution along each segment of a branch within a window. The width of
each window is determined by `tau`

. `tau`

is a fraction of
the root to tip distance so a value of `tau = 0.01`

bins the
phylogeny into 100 time windows of equal width.

A `bammdata`

object with a new component named "dtrates", which
is a list with two or three components:

tau The parameter value of

`tau`

used in the calculation.rates If

`ephy$type = "trait"`

: a numeric vector with the phenotypic rates of each segment on each branch. If`ephy$type = "diversification"`

: a list with two components. The first component is a numeric vector of speciation rates. The second component is a numeric vector of extinction rates.tmat A matrix of the starting and ending times of the segments on each branch. Only if

`tmat = TRUE`

.

If there are zero length branches in the input tree `NA`

s will
result.

Mike Grundler

1 2 3 4 5 6 7 8 | ```
data(whales, events.whales)
ed <- getEventData(whales, events.whales, burnin=0.25, nsamples=500)
# use all posterior samples
ed <- dtRates(ed, tau=0.01)
# use specified range of posterior samples
ed <- dtRates(ed, tau=0.01, ism=50:150)
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.