treeTrace
trees : | TimeTree[] (pass by value) |
Vector of trees. | |
burnin : | Probability (pass by value) |
The fraction/number of samples to discard as burnin. | |
Default : 0.25 |
# Read in a vector of trees
trees = readTrees("trees.nex")
# Create a tree trace
tree_trace = treeTrace(trees, burnin=0.25)
# Create a distribution of trees from the tree trace
tree ~ dnEmpiricalTree(tree_trace)
# Add an MCMC move
moves[1] = mvEmpiricalTree(tree)