Here\'s the data I\'m working with:
data <- data.frame(id = rep(1:3, each = 30), intervention = rep(c("a","b"),each= 2, times=4