# NOT RUN {
if (require("MASS")){
.W <- Womenlf
.W$partic <- factor(.W$partic, levels=c("not.work", "parttime", "fulltime"))
poTest(polr(partic ~ hincome + children + region, data=.W))
}
# }
Run the code above in your browser using DataLab