Parameter tuning with mlrHyperopt

[This article was first published on r-bloggers on Machine Learning in R, and kindly contributed to R-bloggers]. (You can report issue about the content on this page here)
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.

Hyperparameter tuning with mlr is rich in options as they are multiple tuning methods:

  • Simple Random Search
  • Grid Search
  • Iterated F-Racing (via irace)
  • Sequential Model-Based Optimization (via mlrMBO)

Also the search space is easily definable and customizable for each of the 60+ learners of mlr using the ParamSets from the ParamHelpers Package.

The only drawback and shortcoming of mlr in comparison to caret in this regard is that mlr itself does not have defaults for the search spaces. This is where mlrHyperopt comes into play.

mlrHyperopt offers

  • default search spaces for the most important learners in mlr,
  • parameter tuning in one line of code,
  • and an API to add and access custom search spaces from the mlrHyperopt Database.

Installation

# version >= 1.11 needed.
devtools::install_github("berndbischl/ParamHelpers") 
devtools::install_github("jakob-r/mlrHyperopt", dependencies = TRUE)

Tuning in one line

Tuning can be done in one line relying on the defaults. The default will automatically minimize the missclassification rate.

library(mlrHyperopt)
res = hyperopt(iris.task, learner = "classif.svm")
res
## Tune result:
## Op. pars: cost=1.42e+04; gamma=7.81e-05
## mmce.test.mean=0.0333333

We can find out what hyperopt did by inspecting the res object.

Depending on the parameter space mlrHyperopt will automatically decide for a suitable tuning method:

res$opt.path$par.set
##          Type len Def    Constr Req Tunable Trafo
## cost  numeric   -   0 -15 to 15   -    TRUE     Y
## gamma numeric   -  -2 -15 to 15   -    TRUE     Y
res$control
## Tune control: TuneControlMBO
## Same resampling instance: TRUE
## Imputation value: 1
## Start: <NULL>
## 
## Tune threshold: FALSE
## Further arguments:

As the search space defined in the ParamSet is only numeric, sequential Bayesian optimization was chosen. We can look into the evaluated parameter configurations and we can visualize the optimization run.

tail(as.data.frame(res$opt.path))
##          cost       gamma mmce.test.mean dob eol error.message exec.time
## 20  0.1382035  -0.2092684     0.05333333  20  NA          <NA>     0.116
## 21 -1.5577127   1.5836927     0.07333333  21  NA          <NA>     0.122
## 22 14.9974853 -12.1007129     0.04000000  22  NA          <NA>     0.115
## 23 14.9958728 -14.9974155     0.03333333  23  NA          <NA>     0.105
## 24 10.4035727 -13.9625830     0.04666667  24  NA          <NA>     0.101
## 25 13.7909439 -13.6440308     0.03333333  25  NA          <NA>     0.105
plotOptPath(res$opt.path)

The upper left plot shows the distribution of the tried settings in the search space and contour lines indicate where regions of good configurations are located. The lower right plot shows the value of the objective (the miss-classification rate) and how it decreases over the time. This also shows nicely that wrong settings can lead to bad results.

Using the mlrHyperopt API with mlr

If you just want to use mlrHyperopt to access the default parameter search spaces from the Often you don’t want to rely on the default procedures of mlrHyperopt and just incorporate it into your mlr-workflow. Here is one example how you can use the default search spaces for an easy benchmark:

lrns = c("classif.xgboost", "classif.nnet")
lrns = makeLearners(lrns)
tsk = pid.task
rr = makeResampleDesc('CV', stratify = TRUE, iters = 10)
lrns.tuned = lapply(lrns, function(lrn) {
  if (getLearnerName(lrn) == "xgboost") {
    # for xgboost we download a custom ParConfig from the Database
    pcs = downloadParConfigs(learner.name = getLearnerName(lrn))
    pc = pcs[[1]]
  } else {
    pc = getDefaultParConfig(learner = lrn)
  }
  ps = getParConfigParSet(pc)
  # some parameters are dependend on the data (eg. the number of columns)
  ps = evaluateParamExpressions(ps, 
    dict = mlrHyperopt::getTaskDictionary(task = tsk))
  lrn = setHyperPars(lrn, par.vals = getParConfigParVals(pc))
  ctrl = makeTuneControlRandom(maxit = 20)
  makeTuneWrapper(learner = lrn, resampling = rr, par.set = ps, 
                  control = ctrl)
})
res = benchmark(learners = c(lrns, lrns.tuned), tasks = tsk, 
                resamplings = cv10)
plotBMRBoxplots(res) 

As we can see we were able to improve the performance of xgboost and the nnet without any additional knowledge on what parameters we should tune. Especially for nnet improved performance is noticable.

Additional Information

Some recommended additional reads

To leave a comment for the author, please follow the link and comment on their blog: r-bloggers on Machine Learning in R.

R-bloggers.com offers daily e-mail updates about R news and tutorials about learning R and many other topics. Click here if you're looking to post or find an R/data-science job.
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.

Never miss an update!
Subscribe to R-bloggers to receive
e-mails with the latest R posts.
(You will not see this message again.)

Click here to close (This popup will not appear again)