bactdate | R Documentation |
Main function of BactDating package to date the nodes of a phylogenetic tree
bactdate(
tree,
date,
initMu = NA,
initAlpha = NA,
initSigma = NA,
updateMu = T,
updateAlpha = T,
updateSigma = T,
updateRoot = T,
nbIts = 10000,
thin = ceiling(nbIts/1000),
useCoalPrior = T,
model = "arc",
useRec = F,
minbralen = 0.1,
showProgress = F
)
tree |
Tree wih branches measured in unit of substitutions |
date |
Sampling dates for the leaves of the tree |
initMu |
Initial rate of substitutions per genome (not per site), or zero to use root-to-tip estimate |
initAlpha |
Initial coalescent time unit |
initSigma |
Initial std on per-branch substitution rate (only used in relaxed or mixed-relaxed models) |
updateMu |
Whether or not to update the substitution rate |
updateAlpha |
Whether or not to update the coalescent time unit |
updateSigma |
Whether or not to per-branch substitution rate (only used in relaxed or mixed-relaxed models) |
updateRoot |
Whether or not to use the root finding algorithm |
nbIts |
Number of MCMC iterations to perform |
thin |
Thining interval between recorded MCMC samples |
useCoalPrior |
Whether or not to use a coalescent prior on the tree |
model |
Which model to use (poisson or negbin or strictgamma or relaxedgamma or mixedgamma or arc or carc or mixedcarc) |
useRec |
Whether or not to use results from previous recombination analysis |
minbralen |
Minimum branch length of the phylogenetic tree (in number of substitutions) |
showProgress |
Whether or not to show a progress bar |
Dating results
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.