mbo1d: Efficient global optimization with iterative point proposals

Description Usage Arguments Details Value Note Author(s) References See Also

Description

Implements the efficient global optimization algorithm based on Kriging. New design points are generated by the expected improvement criterion or expected quantile improvement criterion. The optimizer of the performance criterion can be choosen as generalized simulated annealing or conditional one dimensional optimization.

Usage

1
2
3
4
mbo1d(model, fun, nsteps, lower, upper, parinit, isoInput, 
maxRunsMult=1, repMult=1, tol_input=.Machine$double.eps^0.25, addInfo=TRUE, 
envir=parent.frame(), EIopt="1Dmulti", GenSAmaxCall=100, timeAlloc="constant",
EItype="EI")

Arguments

model

Fitted Kriging model of the experimental design and evaluation measures. See km

fun

Loss function to be minimized.

nsteps

Number of points the efficient global optimizer should compute.

lower

Vector of lower bounds of the tuning parameters. First element is the lower bound of the first tuning parameter, second element the lower bound of the next tuning parameter etc.

upper

Vector of upper bounds of the tuning parameters. First element is the upper bound of the first tuning parameter, second element the upper bound of the next tuning parameter etc.

parinit

Starting value of the one dimensional optimization algorithm. See optimize1dMulti.

isoInput

Force the covariance structure of the km to have only one range parameter. For details see km.

maxRunsMult

Multiplies the base number of iterations of the one dimensional algorithm. Default number of iterations are the number of hyperparameters. Only applied, if the argument EIopt is set to "1Dmulti".

repMult

Multiplies the base number of random starting values of the one dimensional algorithm. Default number of iterations are the number of hyperparameters. Only applied, if the argument EIopt is set to "1Dmulti".

tol_input

Accuracy threshold of the one dimensional optimization algorithm. See optimize1dMulti.

addInfo

Should additional information be displayed during optimization? (logical value). Default is FALSE.

envir

Internal variable to store environments. Default is to look up the one higher level environment. Modification is unnecessary.

EIopt

Specifies which algorithm is used to optimize the expected improvement criterion. Two alternatives are available "1Dmulti" and "GenSA". The former uses the conditional 1D algorithm and the latter generalized, simulated annealing.

GenSAmaxCall

Maximum number of function calls per parameter to estimate in generalized, simulated annealing. Higher values result in more accurate estimates, but the optimization process is slowed.

timeAlloc

Specifies how the new noise variance is influenced by iteration progress. Default is to use "constant"" allocation. The other available option is to specify "zero", which corresponds to the original expected improvement criterion.

EItype

Defines the type of the improvement criterion. The default EI corresponds to the expected improvement. As an alternative EQI the expected quantile improvement is also possible.

Details

In each step iteration the performance criterion is minimized regarding the design input parameters. The loss function is evaluated at the proposed design point. The Kriging model is updated with the new point and refitted with nugget effect to ensure convergence. The starting value for the next step iteration is set to the parameter values with lowest function value.

Value

Retains the parameter values with the lowest function value. The function value (funcVal) and the MBO tuning process history (MBOprogress) are available as attributes.

Note

Function is supplied for model customization and intended for the experienced user. The more user friendly function uses this code as intermediate step.

Author(s)

Thomas Welchowski welchow@imbie.meb.uni-bonn.de

References

Olivier Roustant and David Ginsbourger and Yves Deville, (2012), DiceKriging, DiceOptim: Two R Packages for the Analysis of Computer Experiments by Kriging-Based Metamodeling and Optimization, Journal of Statistical Software, Volume 51, Issue 1

Donald R. Jones and Matthias Schonlau and William J. Welch, (1998), Efficient Global Optimization of Expensive Black-Box Functions, Journal of Global Optimization 13: pages 455-492

Constantino Tsallis and Daniel A. Stariolo, (1996), Generalized Simulated Annealing, Elsevier Physica A: Statistical Mechanics and its Applications, Volume 233, Issues 1-2, 15, Pages 395-406

Victor Picheny, David Ginsbourger, Yann Richet, (2012), Quantile-based optimization of Noisy Computer Experiments with Tunable Precision, HAL-archives-ouvertes.fr, hal-00578550v3

See Also

optimize1dMulti, km, EI, mboAll


kernDeepStackNet documentation built on May 2, 2019, 8:16 a.m.