library(MASS)
data(iris)
x1 <- iris[iris[,5]=="setosa",1:2]
x2 <- iris[iris[,5]=="versicolor",1:2]
Rhat <- kroc(x1=x1, x2=x2, approx.ref="kcde")
## approx.ref="kcde" is faster but less accurate than "pmvnorm"
plot(Rhat)
abline(0,1, lty=2, col="grey")Run the code above in your browser using DataLab