Description Usage Arguments Details Value References See Also
This function calculates the mean and variance of the posterior mean of labeled substitution counts over a prespecified set of branches on a phylogeny using a simulation-free dynamic programming algorithm and Monte Carlo sampling.
1 2 3 | postmean.moments.phylojumps(tree, rate.mat, label.mat,
edge.set, root.dist,
scale = FALSE, N = 100000)
|
tree |
A tree object of class |
rate.mat |
A reversible continuous-time Markov chain (CTMC) rate matrix. |
label.mat |
A 0-1 matrix of the same size as |
edge.set |
An integer vector of branch indices specifying the branch set over which the moments will be calculated. |
root.dist |
A numeric vector defining the root distribution of the evolutionary process on the tree. If this process starts at stationarity (as is commonly assumed), then the root distribution is equal to the CTMC stationary distribution. |
scale |
A logical indicating whether to scale the time dimension. If
|
N |
An integer specifying the number of Monte Carlo samples to draw. |
The mean obtained in this function is computed exactly, while the variance is partially approximated using Monte Carlo simulations.
More information on phylogenetic tree objects of class "phylo"
can be
found at http://ape-package.ird.fr/misc/FormatTreeR_24Oct2012.pdf.
A numeric vector holding the mean and variance of the posterior mean of labeled substitution counts.
Dhar A and Minin VN (2017) “Calculating Higher-Order Moments of Phylogenetic Stochastic Mapping Summaries in Linear Time”, Journal of Computational Biology, 24(5):377-399.
prior.moments.phylojumps
,
post.moments.phylojumps
,
joint.postmean.moments.phylojumps
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.