I have the following data.frame:
dat <- data.frame(bifpar = c(0.5, 0.75, 1.0), sd_distfeed = rep(2,3), sd_nodist = rep(3,3), mean_distfeed = rep(10,3), mea