\begin{Schunk}
\begin{Sinput}
 stde <- 1
 dat <- data.frame(x=rpois(500, lambda=200))
 dat$y <- 3 + 0.08 * dat$x + stde * rnorm(500)
 m1 <- lm(y ~ x, data=dat); summary(m1)
 plot(y ~ x, data=dat); abline(m1)
\end{Sinput}
\end{Schunk}
