\begin{Schunk}
\begin{Sinput}
 gradratehat <- predict(mod1, interval="confidence")
 plot(gradratehat ~ model.matrix(mod1)[,"private"], data=dat, xlab="Private or Public?", ylab="Graduation Rate")