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