See mbo_parallel for all parallelization options.

mbo(fun, design = NULL, learner = NULL, control,
  show.info = getOption("mlrMBO.show.info", TRUE), more.args = list())

Arguments

fun

[smoof_function]
Fitness function to optimize. For one dimensional target functions you can obtain a smoof_function by using makeSingleObjectiveFunction. For multi dimensional functions use makeMultiObjectiveFunction. It is possible to return even more information which will be stored in the optimization path. To achieve this, simply append the attribute “extras” to the return value of the target function. This has to be a named list of scalar values. Each of these values will be stored additionally in the optimization path.

design

[data.frame]
Initial design as data frame. If the y-values are not already present in design, mbo will evaluate the points. If the parameters have corresponding trafo functions, the design must not be transformed before it is passed! Functions to generate designs are available in ParamHelpers: generateDesign, generateGridDesign, generateRandomDesign. Default is NULL, which means generateDesign is called and a design of size 4 times number of all parameters is created The points are drawn via maximinLHS to maximize the minimal distance between design points.

learner

[Learner]
Regression learner from mlr, which is used as a surrogate to model our fitness function. If NULL (default), the default learner is determined as described here: mbo_default_learner.

control

[MBOControl]
Control object for mbo.

show.info

[logical(1)]
Verbose output on console? Default is TRUE.

more.args

[list]
Further arguments passed to fitness function.

Value

[MBOSingleObjResult | MBOMultiObjResult]

Examples

# simple 2d objective function obj.fun = makeSingleObjectiveFunction( fn = function(x) x[1]^2 + sin(x[2]), par.set = makeNumericParamSet(id = "x", lower = -1, upper = 1, len = 2) ) # create base control object ctrl = makeMBOControl() # do three MBO iterations ctrl = setMBOControlTermination(ctrl, iters = 3L) # use 500 points in the focussearch (should be sufficient for 2d) ctrl = setMBOControlInfill(ctrl, opt.focussearch.points = 500) # create initial design des = generateDesign(n = 5L, getParamSet(obj.fun), fun = lhs::maximinLHS) # start mbo res = mbo(obj.fun, design = des, control = ctrl)
#> Computing y column(s) for design. Not provided.
#> [mbo] 0: x=0.494,0.412 : y = 0.645 : 0.0 secs : initdesign
#> [mbo] 0: x=-0.122,0.0783 : y = 0.093 : 0.0 secs : initdesign
#> [mbo] 0: x=-0.208,-0.325 : y = -0.276 : 0.0 secs : initdesign
#> [mbo] 0: x=0.897,-0.988 : y = -0.0307 : 0.0 secs : initdesign
#> [mbo] 0: x=-0.979,0.69 : y = 1.6 : 0.0 secs : initdesign
#> [mbo] 1: x=1,-0.322 : y = 0.683 : 0.0 secs : infill_cb
#> [mbo] 2: x=-0.181,-0.691 : y = -0.605 : 0.0 secs : infill_cb
#> [mbo] 3: x=0.103,-1 : y = -0.831 : 0.0 secs : infill_cb
print(res)
#> Recommended parameters: #> x=0.103,-1 #> Objective: y = -0.831 #> #> Optimization path #> 5 + 3 entries in total, displaying last 10 (or less): #> x1 x2 y dob eol error.message exec.time cb #> 1 0.4941278 0.41208499 0.64468297 0 NA <NA> 0 NA #> 2 -0.1216173 0.07832955 0.09304025 0 NA <NA> 0 NA #> 3 -0.2077841 -0.32457760 -0.27573422 0 NA <NA> 0 NA #> 4 0.8966086 -0.98750787 -0.03074897 0 NA <NA> 0 NA #> 5 -0.9794215 0.69022501 1.59597716 0 NA <NA> 0 NA #> 6 0.9999967 -0.32226298 0.68327963 1 NA <NA> 0 -0.6290274 #> 7 -0.1807007 -0.69122233 -0.60482669 2 NA <NA> 0 -0.4985452 #> 8 0.1029667 -0.99996360 -0.83084918 3 NA <NA> 0 -1.0001956 #> error.model train.time prop.type propose.time se mean lambda #> 1 <NA> NA initdesign NA NA NA NA #> 2 <NA> NA initdesign NA NA NA NA #> 3 <NA> NA initdesign NA NA NA NA #> 4 <NA> NA initdesign NA NA NA NA #> 5 <NA> NA initdesign NA NA NA NA #> 6 <NA> 0.086 infill_cb 0.179 0.3197429 -0.3092844 1 #> 7 <NA> 0.047 infill_cb 0.189 0.3779856 -0.1205596 1 #> 8 <NA> 0.038 infill_cb 0.164 0.1795866 -0.8206090 1
# NOT RUN { plot(res) # }