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