I searched a lot and could not find a package that had functions to predict the outcome of ordinal multilevel models. So I coded the functions in mrp/functions.R with what I learned about the equations behind ordinal regressions. Would love to get some feedback on it.