\begin{Schunk}
\begin{Sinput}
 dat <- dat[dat$instcstP1000 < 25, ]
 datpublic  <- dat[dat$schtype %in% levels(dat$schtype)[1], ]
 datprivate <- dat[dat$schtype %in% levels(dat$schtype)[2], ]
 par(mfrow=c(1,2))
 plot(gradrate ~ instcstP1000, data = datpublic, main = "Public", xlab = "Instructional Cost/1000", ylab = "Graduation Rate")
 plot(gradrate ~ instcstP1000, data = datprivate, main = "Private", xlab = "Instructional Cost/1000", ylab = "Graduation Rate")
\end{Sinput}
\end{Schunk}
