data(politicalInformation)
table(politicalInformation$y,exclude=NULL)
op <- MASS::polr(y ~ collegeDegree + female + log(age) + homeOwn + govt + log(length),
data=politicalInformation,
Hess=TRUE,
method="probit")
Run the code above in your browser using DataLab