Calculating ordinal regression models in SAS and S-plus