I am running an Ordinal Logistic Regression, as my Dependent Variable has 4 possible levels. I am using the function polr from the MASS package in R. However, now I should include instrumental variables, and the only R packages I have found so far are AER (function ivreg ) and ivtools (function ivglm ). Both do not seem to support the polr function. Does anyone know a solution for this?