A better multiRFfunction

multiRF.R

multiRF <- function(X, Y, mtry_vector,...) {
  
  stopifnot(require(doRNG))
  stopifnot(require(randomForest))

  stopifnot(is.numeric(mtry_vector))
  stopifnot(all(mtry_vector>1))
  
  foreach(i=mtry_vector,.combine=randomForest::combine,.packages='randomForest',.export=c('X','Y'),.inorder=FALSE) %dorng% {
    randomForest(X,Y,mtry=i,...)
  }
}

library(doParallel)
data(iris)
a <- c(rep(3,10),rep(4,10),rep(5,10))

registerDoParallel(cores=1)
system.time(output<-multiRF(iris[,-5], iris[,5], a, ntree=500))

registerDoParallel(cores=8)
system.time(output<-multiRF(iris[,-5], iris[,5], a, ntree=500))