Observation <- 1:14
Dose <- c(3.4,5.2,7,8.5,10.5,13,18,21,28,6.5,10,14,21.5,29)
Prep <- factor(c(rep(0,9),rep(1,5)))
Convulse <- c(0,5,11,14,18,21,23,30,27,2,10,18,21,27)
Total <- c(33,32,38,37,40,37,31,37,30,40,30,40,35,37)
NoConvulse <- Total-Convulse
mice <- data.frame(Observation,Dose,Prep,Convulse,Total,NoConvulse)
brMask(data=mice,blinded="Prep")
Run the code above in your browser using DataLab