I am performing a comparative study (across species) and I would like to test how a single continuous X-variable affects a single Y-variable which is ordinal (0 < 1 < 2), while controlling for the phylogeny. I would like to be able to test differences in both intercepts and slopes.