# Iris data
data(iris)
x <- subset(iris, select = -Species)
y <- iris$Species
method <- c("randomForest","svm","pcalda","knn")
pars <- valipars(sampling="boot", niter = 10, nreps=4)
res <- maccest(x, y, method=method, pars=pars, comp="anova")
res
htest <- mc.anova(res$acc.iter)
oldpar <- par(mar = c(5,10,4,2) + 0.1)
plot(htest$tukey,las=1) ## plot the tukey results
par(oldpar)
Run the code above in your browser using DataLab