\begin{Schunk}
\begin{Sinput}
 dat3 <- genCorrelatedData(N=150, beta = c(0, 0.15, 0.25, 0.1), stde = 150)
 dat3$x3 <- rpois(150, lambda = 7)
 dat3$x4 <- rgamma(150, 2, 1)
 m1 <- lm(y ~ x1 + x2 + x3 + x4, data=dat3)
\end{Sinput}
\end{Schunk}
