# NOT RUN {
library(MASS)
data(france)
polr.mod <- polr(vote ~ age + male + retnat + lrself, data=france)
typical.france <- data.frame(
age = 35,
retnat = factor(1, levels=1:3, labels=levels(france$retnat)),
stringsAsFactors=TRUE)
oc.res <- ordChange(polr.mod, data=france, typical.dat=typical.france, sim=TRUE)
oc2plot(oc.res)
# }
Run the code above in your browser using DataLab