Change function equivalent_shifts_values, to actualize on incidence matrix rather than on shifts values (handle non-parametric trees).