udata <- unique(sort(userValue)) delta <- min(diff(udata))/2 cutpts <- c(udata - delta, udata[length(udata)] + delta)
for (i in 1:length(cutpts)) { pred <- rule(data, cutpts[i]) sens[i] <- mean(pred[group==1]) spec[i] <- mean(1 - pred[group==2]) }
h <- diff(x) lx <- length(x) 0.5 * sum(h * (y[-1] + y[-lx]))