\begin{Schunk}
\begin{Sinput}
 Nsampsize <- 100
 set.seed(12345)
 b0 <- 3; b1 <- 0.25; stde <- 10;
 x <- rnorm(Nsampsize, m=0, sd=10);
 simr <- function(i){
  y <- b0 + b1 * x + stde*rnorm(Nsampsize)
  myr <- lm(y ~ x)
  }
 nsims <- 1000
 mregs <- lapply(1:nsims, simr)
\end{Sinput}
\end{Schunk}
