\begin{Schunk}
\begin{Sinput}
 m1 <- lm(ycatn ~ xcat2n)
 plot.default(ycatn ~ xcat2n, xlim=c(-0.5, 1.5),ylim = c(-0.5, 1.5), xlab = "One Categorical Predictor", ylab="An Output Variable", axes = F, bty="L")
 axis(1, at = c(0, 1), labels = c("L=0","R=1"))
 axis(2, at = c(0, 1), labels = c("L=0", "H=1"))
\end{Sinput}
\end{Schunk}
