Simulated data to fit example models against
sim_data_interactive_cat
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) X3 <- rbinom(n, 5, .23) Z_interactive_cat <- b0 + b1 * X1 + b2 * X3 + b3 * (X1 * X3) Y_interactive_cat <- rnorm(n, Z_interactive_cat, 1) sim_data_interactive_cat <- data.frame(cbind(X1, X3, Y = Y_interactive_cat))