# NOT RUN {
n <- 5000
set.seed(12345)
#Predictors with normal distribution
set.seed(1235)
scoreNegative <- rnorm(n, mean = 0.25,sd = 0.125)
set.seed(11452)
scorePositive1 <- rnorm(n, mean = 0.55,sd = 0.125)
set.seed(54321)
scorePositive2 <- rnorm(n, mean = 0.65,sd = 0.125)
set.seed(65987)
scorePositive3 <- rnorm(n, mean = 0.75,sd = 0.125)
response = c(rep(c("target"), times = n),rep(c("nontarget"), times = n))
predictor1 = c(scoreNegative,scorePositive1)
predictor2 = c(scoreNegative,scorePositive2)
predictor3 = c(scoreNegative,scorePositive3)
responses <- data.frame(
response = response
)
predictors <- data.frame(
DET1 = predictor1,
DET2 = predictor2,
DET3 = predictor3
)
#We can also plot the ROC curves activating logical attribute 'plotROC'
detcurve <- detc(responses,predictors,
names = names(predictors),
positive="target",
title="Example",
plotROC = TRUE)
#If you want to plot the EER and its CI on the curves:
for (name in names(detcurve)){
points(qnorm(detcurve[[name]]$EER),qnorm(detcurve[[name]]$EER),
pch=19,col = detcurve[[name]]$color,lwd=3)
}
# }
Run the code above in your browser using DataLab