\begin{Schunk}
\begin{Sinput}
 incedmod1 <- lm(income~education, data=Prestige)
 plot(income~education, data=Prestige, xlab="Education", ylab="Income")
 abline(incedmod1)
\end{Sinput}
\end{Schunk}
