shifts_to_simmap takes a vector of edges where the shifts occur, and return a
simmap formatted tree, mapped with corresponding regimes.
shifts_to_simmap(tree, shifts_edges)input tree in phylo format
shifts positions on the edges
tree a simmap object
Ancestral state is always 0, and other states are consecutive integers.