View source: R/marginalOddsRatioBranches.R
marginalOddsRatioBranches | R Documentation |
Compute marginal posterior-to-prior odds ratio associated with observing one or more rate shift on a given branch.
marginalOddsRatioBranches(ephy, expectedNumberOfShifts)
ephy |
An object of class |
expectedNumberOfShifts |
Expected number of shifts under the prior alone. |
This function returns a copy of a phylogenetic tree where each branch length is equal to the marginal odds ratio in favor of a rate shift on a particular branch. These cannot be interpreted as evidence for a rate shift in an absolute sense. As explained on the website, they are a marginal odds ratio. This function is provided primarily for the purpose of distinguishing core and non-core shifts.
A object of class phylo
but where each branch length is
equal to the marginal shift odds on each branch.
Dan Rabosky
getBranchShiftPriors
,
distinctShiftConfigurations
,
credibleShiftSet
data(whales, events.whales)
ed <- getEventData(whales, events.whales, burnin=0.1, nsamples=500)
marginalOddsRatioBranches(ed, expectedNumberOfShifts = 1)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.