Fits a logistic or probit regression model to an ordered factor response. The default logistic case is <em>proportional odds logistic regression</em>, after which the function is named.