Major new fitting function orm = ordinal regression model