# simulate data
x <- matrix(rnorm(100), 50, 2)
y <- x %*% c(1, -2)
dat <- data.frame(y = y, x = x, z = as.factor(rbinom(50, 2, 0.5)))
# fit vine regression model
fit <- vinereg(y ~ ., dat)
plot_effects(fit)
Run the code above in your browser using DataLab