shifts.list_to_matrix
takes the list description of the shifts
to give the matricial representation of the shifts : the b th element of
the lth line has the value of the shift on character l occuring on that branch b
shifts.list_to_matrix(phy, shifts, p = nrow(shifts$values))
shifts.matrix_to_list