I am currently trying to estimate a ordinal logistic regression model with hierarchical data in R (my dependent variable has three categories and the data stems from the 50