# Multiple Regression Code Supplement
# Michael B. Morrissey and Graeme D. Ruxton

library(mvtnorm)
library(ellipse)

n <- 200
cov_mat <- matrix(c(1, 0.5, 0.5, 1), 2, 2)
x <- rmvnorm(n, c(0, 0), cov_mat)
x1 <- x[,1]
x2 <- x[,2]
y <- rnorm(n, 0.5*x1, sqrt(0.75))