\begin{Schunk}
\begin{Sinput}
 set.seed(234234)
 x1 <- 0.5 * rpois(100, lambda = 7)
 x2 <- rnorm(100, m = 50, sd = 2)
 x3 <- sample(c(0,1), 100, replace = TRUE)
 y <- 2*x1 + 11.45*x2 + 2.8*x1*(x2-45) + 0.3*x1*x3 + rnorm(100, sd = 10)
 y2 <- 2*x1 + 0.15*x2 + 1*x1*x3 + rnorm(100, sd = 10)
 x3 <- factor(x3, labels=c("left", "right"))
 dat <- data.frame(x1, x2, x3, y, y2)
\end{Sinput}
\end{Schunk}
