\begin{Schunk}
\begin{Sinput}
 dat$schnum <- as.numeric(dat$schtype) - 1
 plot(gradrate ~ schnum, data = dat, xlab = "Private or Public?", ylab = "Graduation Rate")
 abline(mod1)
\end{Sinput}
\end{Schunk}
