# NOT RUN {
set.seed(1, kind = "Mersenne-Twister", normal.kind = "Inversion")
tree <- ape::rtree(10)
regimes <- sample(letters[1:3], nrow(tree$edge), replace = TRUE)
PCMTreeSetRegimesForEdges(tree, regimes)
# }
# NOT RUN {
PCMTreePlot(tree) + ggtree::geom_nodelab() + ggtree::geom_tiplab()
# }
# NOT RUN {
regime <- PCMRegimes(tree)[1]
PCMTreeGetTipsInRegime(tree, regime)
print(regime)
# }
Run the code above in your browser using DataCamp Workspace