data(iris)
rf <- function(formula, ...)
ml_model$new(formula, info="grf::probability_forest",
fit=function(x,y, ...) grf::probability_forest(X=x, Y=y, ...),
pred=function(fit, newdata) predict(fit, newdata)$predictions, ...)
args <- expand.list(num.trees=c(100,200), mtry=1:3,
formula=c(Species ~ ., Species ~ Sepal.Length + Sepal.Width))
models <- lapply(args, function(par) do.call(rf, par))
x <- models[[1]]$clone()
x$estimate(iris)
predict(x, newdata=head(iris))
# Reduce Ex. timing
a <- targeted::cv(models, data=iris)
cbind(coef(a), attr(args, "table"))
Run the code above in your browser using DataLab