\begin{Schunk}
\begin{Sinput}
 b0 <- 3; b1 <- 0.25; stde <- 10;
 simrhomo <- function(i){
  e <- rnorm(Nsampsize, m=0, s=1)
  yhom <- b0 + b1 * x + stde*e
  mod <- lm(yhom ~ x)
  mod
  }
 nsims <- 1000
 reshomo <- lapply(1:nsims, simrhomo)
\end{Sinput}
\end{Schunk}
