\begin{Schunk}
\begin{Sinput}
 N <- 400
 x <- 20 + rgamma(N, shape=2, scale=5)
 y <- 2 + 0.1 *x + 10 * rnorm(N)
 plot(y ~ x, xlab = "One Mysterious Predictor", ylab="An Output Variable")
\end{Sinput}
\end{Schunk}
