\begin{Schunk}
\begin{Sinput}
 incedmod1 <- lm(income~education, data=Prestige)
 summary(incedmod1)
 anova(incedmod1, test="F")
 vcov(incedmod1)
 confint(incedmod1)
 plot(incedmod1)
 termplot(incedmod1, se=T, partial=T)
\end{Sinput}
\end{Schunk}
