Is there a function to make prediction in ordinal regression model with bootstrap in R?

Question

I have an ordinal dependent variable and several independent variables I am using the polr function from the MASS package to build my model.

First I have the datafram