data("Laurasiatherian")
set.seed(123)
bs <- bootstrap.phyDat(Laurasiatherian,
FUN=function(x)upgma(dist.ml(x)), bs=100)
tree_compat <- allCompat(bs, rooted=TRUE) |>
add_edge_length(bs)
plot(tree_compat)
add_boxplot(tree_compat, bs)
Run the code above in your browser using DataLab