\begin{Schunk}
\begin{Sinput}
 plot(sei~educ, data=dat, cex=0.5,main="", xlab="Education (years)", ylab="Socio-economic Status",ylim=c(0,100))
 mod1 <- lm(sei~educ, data=dat) 
 abline(mod1) 
 abline(2,3.3,lty=2)
 abline(v=0,lty=2, lwd=0.3) 
 abline(h=0,lty=2, lwd=0.3) 
 legend("topleft",legend=c("OLS","PJ's Guess"),lty=c(1,2),bg="white",lwd=c(1,1))
\end{Sinput}
\end{Schunk}
