Simulated data to fit example models against
sim_data_interactive
a data.frame
set.seed(123456) b0 <- 0.2 # true value for the intercept b1 <- 0.5 # true value for first beta b2 <- 0.7 # true value for second beta b3 <- -0.3 # true value for second beta n <- 500 # sample size X1 <- runif(n, -1, 1) X2 <- runif(n, -1, 1) Z_interactive <- b0 + b1 * X1 + b2 * X2 + b3 * (X1 * X2) Y_interactive <- rnorm(n, Z_interactive, 1) sim_data_interactive <- data.frame(cbind(X1, X2, Y = Y_interactive))