# 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)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab