clusters_from_shifts take a vector of shifts edges, and gives the
clustering of the tips induced by them, in a "no homoplasy" model (i.e. no
convergence is allowed).
clusters_from_shifts(tree, edges, part.list = enumerate_tips_under_edges(tree))enumerate_tips_under_edges to compute
the list of tips under each edge.