\begin{Schunk}
\begin{Sinput}
 N <- 500 
 xcat2 <- cut(rnorm(N), breaks=c(-Inf, -2, 0, 0.5, 1, Inf), labels=c("R", "M", "D", "A", "F"))
 xcat2n <- contrasts(xcat2)[xcat2, ]	
 y <- 2 + xcat2n %*% c(3.5, -2.0, 1.3, -0.2) + 1 * rnorm(N)
 plot(y ~ (as.numeric(xcat2) -1), xlab = "One Ordinal Predictor", ylab="An Output Variable", xaxt = "n")
 axis(1, at = 1:5, labels = c("R", "M", "D", "A", "F"))
\end{Sinput}
\end{Schunk}
