\begin{Schunk}
\begin{Sinput}
 stde <- 2
 dat$y2 <- 3 + 0.08 * dat$x + stde * rnorm(500)
 m2 <- lm(y2 ~ x, data=dat); summary(m2)
 plot(y2 ~ x, data=dat); abline(m2)
\end{Sinput}
\end{Schunk}
