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.
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")
Fitted Kriging model of the experimental design and evaluation measures. See km
Loss function to be minimized.
Number of points the efficient global optimizer should compute.
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.
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.
Starting value of the one dimensional optimization algorithm. See optimize1dMulti
.
Force the covariance structure of the km to have only one range parameter. For details see km
.
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"
.
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"
.
Accuracy threshold of the one dimensional optimization algorithm. See optimize1dMulti
.
Should additional information be displayed during optimization? (logical value). Default is FALSE.
Internal variable to store environments. Default is to look up the one higher level environment. Modification is unnecessary.
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.
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.
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.
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.
Retains the parameter values with the lowest function value. The function value (funcVal
) and the MBO tuning process history (MBOprogress
) are available as attributes.
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.
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