BayesianOptimization
Bayesian optimization results
Description
A BayesianOptimization
object contains the
results of a Bayesian optimization. It is the output of bayesopt
or a fit function that accepts the
OptimizeHyperparameters
namevalue pair such as fitcdiscr
. In addition, a BayesianOptimization
object contains data for each iteration of
bayesopt
that can be accessed by a plot function or an output
function.
Creation
Create a BayesianOptimization
object by using the
bayesopt
function or one of the following
fit functions with the OptimizeHyperparameters
namevalue
argument.
Classification fit functions:
fitcdiscr
,fitcecoc
,fitcensemble
,fitcgam
,fitckernel
,fitcknn
,fitclinear
,fitcnb
,fitcnet
,fitcsvm
,fitctree
Regression fit functions:
fitrensemble
,fitrgam
,fitrgp
,fitrkernel
,fitrlinear
,fitrnet
,fitrsvm
,fitrtree
Properties
Problem Definition Properties
ObjectiveFcn
— ObjectiveFcn
argument used by
bayesopt
function handle
This property is readonly.
ObjectiveFcn
argument used by
bayesopt
, specified as a function handle.
If you call
bayesopt
directly,ObjectiveFcn
is thebayesopt
objective function argument.If you call a fit function containing the
'OptimizeHyperparameters'
namevalue pair argument,ObjectiveFcn
is a function handle that returns the misclassification rate for classification or returns the logarithm of one plus the crossvalidation loss for regression, measured by fivefold crossvalidation.
Data Types: function_handle
VariableDescriptions
— VariableDescriptions
argument that bayesopt
used
vector of optimizableVariable
objects
This property is readonly.
VariableDescriptions
argument that
bayesopt
used, specified as a vector of optimizableVariable
objects.
If you called
bayesopt
directly,VariableDescriptions
is thebayesopt
variable description argument.If you called a fit function with the
OptimizeHyperparameters
namevalue pair,VariableDescriptions
is the vector of hyperparameters.
Options
— Options that bayesopt
used
structure
This property is readonly.
Options that bayesopt
used, specified as a
structure.
If you called
bayesopt
directly,Options
is the options used inbayesopt
, which are the namevalue pairs Seebayesopt
Input Arguments.If you called a fit function with the
OptimizeHyperparameters
namevalue pair,Options
are the defaultbayesopt
options, modified by theHyperparameterOptimizationOptions
namevalue pair.
Options
is a readonly structure containing the
following fields.
Option Name  Meaning 

AcquisitionFunctionName  Acquisition function name. See Acquisition Function Types. 
IsObjectiveDeterministic  true means the objective function
is deterministic, false
otherwise. 
ExplorationRatio  Used only when
AcquisitionFunctionName is
'expectedimprovementplus' or
'expectedimprovementpersecondplus' .
See Plus. 
MaxObjectiveEvaluations  Objective function evaluation limit. 
MaxTime  Time limit. 
XConstraintFcn  Deterministic constraints on variables. See Deterministic Constraints — XConstraintFcn. 
ConditionalVariableFcn  Conditional variable constraints. See Conditional Constraints — ConditionalVariableFcn. 
NumCoupledConstraints  Number of coupled constraints. See Coupled Constraints. 
CoupledConstraintTolerances  Coupled constraint tolerances. See Coupled Constraints. 
AreCoupledConstraintsDeterministic  Logical vector specifying whether each coupled constraint is deterministic. 
Verbose  Commandline display level. 
OutputFcn  Function called after each iteration. See Bayesian Optimization Output Functions. 
SaveVariableName  Variable name for the
@assignInBase output function.

SaveFileName  File name for the @saveToFile
output function. 
PlotFcn  Plot function called after each iteration. See Bayesian Optimization Plot Functions 
InitialX  Points where bayesopt evaluated
the objective function. 
InitialObjective  Objective function values at
InitialX . 
InitialConstraintViolations  Coupled constraint function values at
InitialX . 
InitialErrorValues  Error values at InitialX . 
InitialObjectiveEvaluationTimes  Objective function evaluation times at
InitialX . 
InitialIterationTimes  Time for each iteration, including objective function evaluation and other computations. 
Data Types: struct
Solution Properties
MinObjective
— Minimum observed value of objective function
real scalar
This property is readonly.
Minimum observed value of objective function, specified as a real scalar. When there are coupled constraints or evaluation errors, this value is the minimum over all observed points that are feasible according to the final constraint and Error models.
Data Types: double
XAtMinObjective
— Observed point with minimum objective function value
1
byD
table
This property is readonly.
Observed point with minimum objective function value, specified as a
1
byD
table, where
D
is the number of variables.
Data Types: table
MinEstimatedObjective
— Estimated objective function value
real scalar
This property is readonly.
Estimated objective function value at
XAtMinEstimatedObjective
, specified as a real
scalar.
MinEstimatedObjective
is the mean value of the
posterior distribution of the final objective model. The software
estimates the MinEstimatedObjective
value by
passing XAtMinEstimatedObjective
to the object
function predictObjective
.
Data Types: double
XAtMinEstimatedObjective
— Point with minimum upper confidence bound of objective function value
1
byD
table
This property is readonly.
Point with the minimum upper confidence bound of the objective
function value among the visited points, specified as a
1
byD
table, where
D
is the number of variables. The software uses
the final objective model to find the upper confidence bounds of the
visited points.
XAtMinEstimatedObjective
is the same as the best
point returned by the bestPoint
function with the
default criterion
('minvisitedupperconfidenceinterval'
).
Data Types: table
NumObjectiveEvaluations
— Number of objective function evaluations
positive integer
This property is readonly.
Number of objective function evaluations, specified as a positive integer. This includes the initial evaluations to form a posterior model as well as evaluation during the optimization iterations.
Data Types: double
TotalElapsedTime
— Total elapsed time of optimization in seconds
positive scalar
This property is readonly.
Total elapsed time of optimization in seconds, specified as a positive scalar.
Data Types: double
NextPoint
— Next point to evaluate if optimization continues
1
byD
table
This property is readonly.
Next point to evaluate if optimization continues, specified as a
1
byD
table, where
D
is the number of variables.
Data Types: table
Trace Properties
XTrace
— Points where the objective function was evaluated
T
byD
table
This property is readonly.
Points where the objective function was evaluated, specified as a
T
byD
table, where
T
is the number of evaluation points and
D
is the number of variables.
Data Types: table
ObjectiveTrace
— Objective function values
column vector of length T
This property is readonly.
Objective function values, specified as a column vector of length
T
, where T
is the number of
evaluation points. ObjectiveTrace
contains the
history of objective function evaluations.
Data Types: double
ObjectiveEvaluationTimeTrace
— Objective function evaluation times
column vector of length T
This property is readonly.
Objective function evaluation times, specified as a column vector of
length T
, where T
is the number of
evaluation points. ObjectiveEvaluationTimeTrace
includes the time in evaluating coupled constraints, because the
objective function computes these constraints.
Data Types: double
IterationTimeTrace
— Iteration times
column vector of length T
This property is readonly.
Iteration times, specified as a column vector of length
T
, where T
is the number of
evaluation points. IterationTimeTrace
includes both
objective function evaluation time and other overhead.
Data Types: double
ConstraintsTrace
— Coupled constraint values
T
byK
array
This property is readonly.
Coupled constraint values, specified as a
T
byK
array, where
T
is the number of evaluation points and
K
is the number of coupled constraints.
Data Types: double
ErrorTrace
— Error indications
column vector of length T
of 1
or 1
entries
This property is readonly.
Error indications, specified as a column vector of length
T
of 1
or
1
entries, where T
is the
number of evaluation points. Each 1
entry indicates
that the objective function errored or returned NaN
on the corresponding point in XTrace
. Each
1
entry indicates that the objective function
value was computed.
Data Types: double
FeasibilityTrace
— Feasibility indications
logical column vector of length T
This property is readonly.
Feasibility indications, specified as a logical column vector of
length T
, where T
is the number of
evaluation points. Each 1
entry indicates that the
final constraint model predicts feasibility at the corresponding point
in XTrace
.
Data Types: logical
FeasibilityProbabilityTrace
— Probability that evaluation point is feasible
column vector of length T
This property is readonly.
Probability that evaluation point is feasible, specified as a column
vector of length T
, where T
is the
number of evaluation points. The probabilities come from the final
constraint model, including the error constraint model, on the
corresponding points in XTrace
.
Data Types: double
IndexOfMinimumTrace
— Which evaluation gave minimum feasible objective
column vector of integer indices of length
T
This property is readonly.
Which evaluation gave minimum feasible objective, specified as a
column vector of integer indices of length T
, where
T
is the number of evaluation points. Feasibility
is determined with respect to the constraint models that existed at each
iteration, including the error constraint model.
Data Types: double
ObjectiveMinimumTrace
— Minimum observed objective
column vector of length T
This property is readonly.
Minimum observed objective, specified as a column vector of length
T
, where T
is the number of
evaluation points.
Data Types: double
EstimatedObjectiveMinimumTrace
— Estimated objective
column vector of length T
This property is readonly.
Estimated objective, specified as a column vector of length
T
, where T
is the number of
evaluation points. The estimated objective at each iteration is
determined with respect to the objective model at that iteration. At
each iteration, the software uses the object function predictObjective
to
estimate the objective function value at the point with the minimum
upper confidence bound of the objective function among the visited
points.
Data Types: double
UserDataTrace
— Auxiliary data from the objective function
cell array of length T
This property is readonly.
Auxiliary data from the objective function, specified as a cell array
of length T
, where T
is the number
of evaluation points. Each entry in the cell array is the
UserData
returned in the third output of the
objective function.
Data Types: cell
Object Functions
bestPoint  Best point in a Bayesian optimization according to a criterion 
plot  Plot Bayesian optimization results 
predictConstraints  Predict coupled constraint violations at a set of points 
predictError  Predict error value at a set of points 
predictObjective  Predict objective function at a set of points 
predictObjectiveEvaluationTime  Predict objective function run times at a set of points 
resume  Resume a Bayesian optimization 
Examples
Create a BayesianOptimization
Object Using bayesopt
This example shows how to create a BayesianOptimization
object by using bayesopt
to minimize crossvalidation loss.
Optimize hyperparameters of a KNN classifier for the ionosphere
data, that is, find KNN hyperparameters that minimize the crossvalidation loss. Have bayesopt
minimize over the following hyperparameters:
Nearestneighborhood sizes from 1 to 30
Distance functions
'chebychev'
,'euclidean'
, and'minkowski'
.
For reproducibility, set the random seed, set the partition, and set the AcquisitionFunctionName
option to 'expectedimprovementplus'
. To suppress iterative display, set 'Verbose'
to 0
. Pass the partition c
and fitting data X
and Y
to the objective function fun
by creating fun
as an anonymous function that incorporates this data. See Parameterizing Functions.
load ionosphere rng default num = optimizableVariable('n',[1,30],'Type','integer'); dst = optimizableVariable('dst',{'chebychev','euclidean','minkowski'},'Type','categorical'); c = cvpartition(351,'Kfold',5); fun = @(x)kfoldLoss(fitcknn(X,Y,'CVPartition',c,'NumNeighbors',x.n,... 'Distance',char(x.dst),'NSMethod','exhaustive')); results = bayesopt(fun,[num,dst],'Verbose',0,... 'AcquisitionFunctionName','expectedimprovementplus')
results = BayesianOptimization with properties: ObjectiveFcn: [function_handle] VariableDescriptions: [1x2 optimizableVariable] Options: [1x1 struct] MinObjective: 0.1197 XAtMinObjective: [1x2 table] MinEstimatedObjective: 0.1213 XAtMinEstimatedObjective: [1x2 table] NumObjectiveEvaluations: 30 TotalElapsedTime: 29.7824 NextPoint: [1x2 table] XTrace: [30x2 table] ObjectiveTrace: [30x1 double] ConstraintsTrace: [] UserDataTrace: {30x1 cell} ObjectiveEvaluationTimeTrace: [30x1 double] IterationTimeTrace: [30x1 double] ErrorTrace: [30x1 double] FeasibilityTrace: [30x1 logical] FeasibilityProbabilityTrace: [30x1 double] IndexOfMinimumTrace: [30x1 double] ObjectiveMinimumTrace: [30x1 double] EstimatedObjectiveMinimumTrace: [30x1 double]
Create a BayesianOptimization
Object Using a Fit Function
This example shows how to minimize the crossvalidation loss in the ionosphere
data using Bayesian optimization of an SVM classifier.
Load the data.
load ionosphere
Optimize the classification using the 'auto'
parameters.
rng default % For reproducibility Mdl = fitcsvm(X,Y,'OptimizeHyperparameters','auto')
=====================================================================================================  Iter  Eval  Objective  Objective  BestSoFar  BestSoFar  BoxConstraint KernelScale    result   runtime  (observed)  (estim.)    =====================================================================================================  1  Best  0.22222  9.7944  0.22222  0.22222  64.836  0.0015729   2  Accept  0.35897  0.073885  0.22222  0.23074  0.036335  5.5755   3  Best  0.13105  3.9262  0.13105  0.14145  0.0022147  0.0023957   4  Accept  0.35897  0.089783  0.13105  0.13108  5.1259  98.62   5  Best  0.12536  0.21521  0.12536  0.12538  0.0010294  0.037239   6  Accept  0.12821  0.15008  0.12536  0.12502  0.0010316  0.014375   7  Accept  0.1339  0.68533  0.12536  0.12424  0.025081  0.016111   8  Accept  0.12536  0.14609  0.12536  0.12631  0.0010618  0.016299   9  Best  0.11681  0.22373  0.11681  0.12174  0.001133  0.02714   10  Accept  0.12251  0.29918  0.11681  0.12148  0.001008  0.026832   11  Accept  0.12536  0.18019  0.11681  0.12261  0.0010061  0.028162   12  Accept  0.14245  0.58455  0.11681  0.11764  733.87  5.6856   13  Accept  0.1396  0.22235  0.11681  0.12257  0.004017  0.016529   14  Accept  0.1339  8.3508  0.11681  0.12287  974.87  0.49411   15  Accept  0.1396  0.17699  0.11681  0.12324  985.69  119.34   16  Accept  0.35897  0.29536  0.11681  0.12301  0.0010027  8.6693   17  Accept  0.35897  0.15066  0.11681  0.12197  984.97  967.71   18  Accept  0.12536  0.40566  0.11681  0.12202  996.7  37.438   19  Accept  0.13105  0.1635  0.11681  0.12219  0.034548  0.07008   20  Accept  0.13105  0.75845  0.11681  0.12241  0.59807  0.073309  =====================================================================================================  Iter  Eval  Objective  Objective  BestSoFar  BestSoFar  BoxConstraint KernelScale    result   runtime  (observed)  (estim.)    =====================================================================================================  21  Accept  0.14245  0.1707  0.11681  0.1227  0.0011223  0.13671   22  Accept  0.12251  0.17267  0.11681  0.12248  0.0011076  0.043248   23  Accept  0.12251  0.12665  0.11681  0.12232  0.0010439  0.042774   24  Accept  0.13105  0.83651  0.11681  0.12235  0.15094  0.035846   25  Accept  0.12251  0.52615  0.11681  0.12234  997.65  56.317   26  Accept  0.35897  0.20878  0.11681  0.12217  0.55688  823.33   27  Accept  0.1339  0.34271  0.11681  0.12214  1.4203  0.34163   28  Accept  0.1339  0.77884  0.11681  0.12215  8.2318  0.28112   29  Accept  0.12251  0.2529  0.11681  0.12206  0.10082  0.18816   30  Accept  0.13105  0.154  0.11681  0.1221  0.23533  0.16241 
__________________________________________________________ Optimization completed. MaxObjectiveEvaluations of 30 reached. Total function evaluations: 30 Total elapsed time: 55.4351 seconds Total objective function evaluation time: 30.4624 Best observed feasible point: BoxConstraint KernelScale _____________ ___________ 0.001133 0.02714 Observed objective function value = 0.11681 Estimated objective function value = 0.12305 Function evaluation time = 0.22373 Best estimated feasible point (according to models): BoxConstraint KernelScale _____________ ___________ 0.0010061 0.028162 Estimated objective function value = 0.1221 Estimated function evaluation time = 0.17879
Mdl = ClassificationSVM ResponseName: 'Y' CategoricalPredictors: [] ClassNames: {'b' 'g'} ScoreTransform: 'none' NumObservations: 351 HyperparameterOptimizationResults: [1x1 BayesianOptimization] Alpha: [99x1 double] Bias: 4.2973 KernelParameters: [1x1 struct] BoxConstraints: [351x1 double] ConvergenceInfo: [1x1 struct] IsSupportVector: [351x1 logical] Solver: 'SMO' Properties, Methods
The fit achieved about 12% loss for the default 5fold cross validation.
Examine the BayesianOptimization
object that is returned in the HyperparameterOptimizationResults
property of the returned model.
disp(Mdl.HyperparameterOptimizationResults)
BayesianOptimization with properties: ObjectiveFcn: @createObjFcn/inMemoryObjFcn VariableDescriptions: [5x1 optimizableVariable] Options: [1x1 struct] MinObjective: 0.1168 XAtMinObjective: [1x2 table] MinEstimatedObjective: 0.1221 XAtMinEstimatedObjective: [1x2 table] NumObjectiveEvaluations: 30 TotalElapsedTime: 55.4351 NextPoint: [1x2 table] XTrace: [30x2 table] ObjectiveTrace: [30x1 double] ConstraintsTrace: [] UserDataTrace: {30x1 cell} ObjectiveEvaluationTimeTrace: [30x1 double] IterationTimeTrace: [30x1 double] ErrorTrace: [30x1 double] FeasibilityTrace: [30x1 logical] FeasibilityProbabilityTrace: [30x1 double] IndexOfMinimumTrace: [30x1 double] ObjectiveMinimumTrace: [30x1 double] EstimatedObjectiveMinimumTrace: [30x1 double]
Version History
Introduced in R2016b
Open Example
You have a modified version of this example. Do you want to open this example with your edits?
MATLAB Command
You clicked a link that corresponds to this MATLAB command:
Run the command by entering it in the MATLAB Command Window. Web browsers do not support MATLAB commands.
Select a Web Site
Choose a web site to get translated content where available and see local events and offers. Based on your location, we recommend that you select: .
You can also select a web site from the following list:
How to Get Best Site Performance
Select the China site (in Chinese or English) for best site performance. Other MathWorks country sites are not optimized for visits from your location.
Americas
 América Latina (Español)
 Canada (English)
 United States (English)
Europe
 Belgium (English)
 Denmark (English)
 Deutschland (Deutsch)
 España (Español)
 Finland (English)
 France (Français)
 Ireland (English)
 Italia (Italiano)
 Luxembourg (English)
 Netherlands (English)
 Norway (English)
 Österreich (Deutsch)
 Portugal (English)
 Sweden (English)
 Switzerland
 United Kingdom (English)