Ordinal logistic regression

lrm.R

library(rms)  # This is the library that does all the orginal logistic regression magic
EMC <- latest.EMC.2011.clean.v2  # Make a copy of the data frame with a shorter name
model <- lrm(lifeoverall ~ fire + med + streetrepair + streetclean, data=EMC)  # Fit an ordinal logistic regression model
model  # General model summary
exp(coefficients(model))  # Factor change coefficients?
anova(model)  # Wald statistics
Predict(model, fire=0, med=0, streetrepair=0, streetclean=0)  # Predict yhat based on given values