sarahromanes

library(tidyverse)
library(ggExtra)
library(mvtnorm)

# Generate MVN data 
s <- matrix(c(4,2,2,3), ncol = 2)
s1 <- matrix(c(3, -0.6,-0.6, 1), ncol=2)
s2 <- matrix(c(2,-0.08,-0.08,1), ncol = 2)
m1 <- c(0, 0)
m2 <- c(-3, 4)