Dear R users, I would like to fit an ordinal (proportional odds) multilevel model. Is it possible with lme4? Or another R package? Or should I go to SAS? (Glimmix and nlmixed can do it). Thanks for any advice. Hester Lingsma