\begin{Schunk}
\begin{Sinput}
 mod3 <- lm(gradrate ~ instcst, data = dat, subset = instcst < 25000)
 plotSlopes(mod3, plotx = "instcst", interval = "confidence")
\end{Sinput}
\end{Schunk}
