Rev Language Reference
mapTree
Finds the maximum a posteriori (MAP) topology from a trace of trees and summarizes branch lengths.
Usage
mapTree(TraceTree trace, String file, Bool ccAges, Bool ccp, Bool conditionalAges, Probability hpd, Bool mean, Bool sampledAncestors, Bool positiveBranchLengths)
Arguments
| trace : | TraceTree (pass by reference) |
| The samples of trees from the posterior. | |
| file : | String (pass by value) |
| The name of the file where to store the tree. | |
| ccAges : | Bool (pass by value) |
| Annotate conditional clade ages? | |
| Default : FALSE | |
| ccp : | Bool (pass by value) |
| Annotate conditional clade probabilities? | |
| Default : FALSE | |
| conditionalAges : | Bool (pass by value) |
| Annotate node ages conditional on the topology? | |
| Default : FALSE | |
| hpd : | Probability (pass by value) |
| The probability mass of the highest posterior density node age interval. | |
| Default : 0.95 | |
| mean : | Bool (pass by value) |
| Annotate node ages using the mean age instead of the median? | |
| Default : TRUE | |
| sampledAncestors : | Bool (pass by value) |
| Annotate sampled ancestor probs? | |
| Default : TRUE | |
| positiveBranchLengths : | Bool (pass by value) |
| Force negative branch lengths to be short but positive? | |
| Default : FALSE |
Return Type
Example
# Read in tree trace
tree_trace = readTreeTrace("output/my.trees", burnin=0.25)
# Generate the MAP tree
map_tree = mapTree(trace=tree_trace, file="map.tree")
Edit this page