\begin{Schunk}
\begin{Sinput}
 set.seed(132)
 b0 <- .1; b1 <- 0.75; b2 <- -0.15
 N <- 100
 x1 <- rnorm(N, m=20, s=5)
 x2 <- rnorm(N, m=10, s=3)
 eta <- b0 + b1 * x1 + b2 * x2
 input <- exp(eta)
 y <- rpois(N, lambda = input)
 dat <- data.frame(x1,x2,y)
\end{Sinput}
\end{Schunk}
