\begin{Schunk}
\begin{Sinput}
 mod1 <- lm(y ~ x1, data = dat)
 plot(y ~ x1, data = dat)
 abline(mod1)
\end{Sinput}
\end{Schunk}
