v1 <- factor(rep(1:3, times = 16))
set.seed(1)
x1 <- rnorm(48)
set.seed(2)
y <- x1 + rnorm(48, 0, 3)
m <- lm(y ~ v1 + x1)
out <- DMR(m)
x0 <- out$Best$Hypotheses
y0 <- matrix(c(0, 0, -1, 1, 1, 0, 0, 0), 2, 4)
measures <- roc(y0, x0)
sen <- measures$sensitivity
spe <- measures$specificity
Run the code above in your browser using DataLab