getMarginalBranchRateMatrix: Compute mean branch rates for bammdata object

Description

For each sample in the posterior, computes the mean rate for each branch in the focal phylogeny (speciation, extinction, trait evolution). If the bammdata object contains nsamples samples and the target phylogeny has nbranches branches, the function will compute a matrix of nbranches x nsamples.

Usage

1

Arguments

ephy

An object of class bammdata.

verbose

Print progress during processing of bammdata object.

Details

If a type = 'diversification' bammdata object is passed as an argument, the function will return matrices for both speciation and extinction. If type = 'trait' object, the matrix will simply be the corresponding phenotypic rates. Branch-specific rates are the mean rates computed by integrating the relevant rate-through-time function along each branch, then dividing by the length of the branch.

Value

Returns a list with the following components:

Author(s)

Dan Rabosky

Examples

1
2
3
4
5
data(whales)
data(events.whales)
ed <- getEventData(whales, events.whales, nsamples = 10)
mbr <- getMarginalBranchRateMatrix(ed)
dim(mbr$lambda_branch_matrix)

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

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