## No-PASt-BO: Normalized Portfolio Allocation Strategy for Bayesian Optimization

Thiago de P. Vasconcelos, Daniel A. R. M. A. de Souza, C´ esar L. C. Mattos and Jo˜ ao P. P. Gomes

Department of Computer Science

Federal University of Cear´ a

Fortaleza-CE, Brazil thiagodpvasconcelos@hotmail.com, danielramos@lia.ufc.br, cesarlincoln@dc.ufc.br, jpaulo@dc.ufc.br

Abstract -Bayesian Optimization (BO) is a framework for black-box optimization that is especially suitable for expensive cost functions. Among the main parts of a BO algorithm, the acquisition function is of fundamental importance, since it guides the optimization algorithm by translating the uncertainty of the regression model in a utility measure for each point to be evaluated. Considering such aspect, selection and design of acquisition functions are one of the most popular research topics in BO. Since no single acquisition function was proved to have better performance in all tasks, a well-established approach consists of selecting different acquisition functions along the iterations of a BO execution. In such an approach, the GP-Hedge algorithm is a widely used option given its simplicity and good performance. Despite its success in various applications, GP-Hedge shows an undesirable characteristic of accounting on all past performance measures of each acquisition function to select the next function to be used. In this case, good or bad values obtained in an initial iteration may impact the choice of the acquisition function for the rest of the algorithm. This fact may induce a dominant behavior of an acquisition function and impact the final performance of the method. Aiming to overcome such limitation, in this work we propose a variant of GP-Hedge, named No-PASt-BO, that reduce the influence of far past evaluations. Moreover, our method presents a built-in normalization that avoids the functions in the portfolio to have similar probabilities, thus improving the exploration. The obtained results on both synthetic and realworld optimization tasks indicate that No-PASt-BO presents competitive performance and always outperforms GP-Hedge.

Index Terms -Bayesian Optimization, Acquisition Functions, Portfolio Allocation

## I. INTRODUCTION

Bayesian optimization (BO, [1]) is a sequential design strategy based on the Bayesian framework for global optimization of unknown functions, i.e., possibly noisy functions without known closed-form expression and gradient information. It has been widely used in the most diverse tasks, such as selection of hyperparameters for machine learning algorithms [2]-[6], control policies in robotics [7], [8], automated circuit design [9], [10], etc.

The BO approach is especially useful in scenarios where the objective function is costly to evaluate. The inherent uncertainty considered by the Bayesian methodology allows

This article was part of a project financed through the R&amp;D program of the Brazilian Electric Energy Agency (ANEEL), with number PD-000633045 and funded by CPFL Energy and Companhia Energetica Rio das Antas (CERAN). The authors also acknowledge the support from Delfos Intelligent Maintenance.

for a more efficient exploration of the optimization domain with respect to the number of queried points. This is critical in applications where each individual evaluation usually involves substantial financial and/or computational effort. The uncertainty with respect to the function to be optimized usually is modelled using the Gaussian Process (GP) framework [11].

One of the main ingredients of the BO methodology is the so-called acquisition function, which translates the uncertainty in the task domain to a simple to evaluate function that quantifies the expected utility of each point. Such function is then optimized to select the next point to be evaluated, i.e., the next candidate solution. However, although several acquisition functions have been proposed in the literature, there is not a single choice that is always better for any task [12].

Hoffman et al. [12] tackled the aforementioned issue by proposing a hierarchical hedging approach for managing an adaptive portfolio of acquisition functions based on their past performances. Similar approaches have been pursued by other authors. Shahriari et al. [13] expands the original hedge by proposing a choice criterion based on information theoretic considerations. Although in a multi-armed bandit learning context, Shen et al. [14] also consider a portfoliobased strategy for balancing exploration and exploitation during the sequential decision procedure. Recently, Lyu et al. [9] proposed an alternative strategy that considers the multiobjective optimization of multiple acquisition functions to obtain a Pareto front from where candidate points can be sampled in a batch fashion.

Despite the above compelling recent work on the acquisition function choice problem, the resulting solutions stray from the simplicity and applicability of the original hedging strategy presented in [12], named GP-Hedge. Furthermore, GP-Hedge presents some undesirable properties. For instance, since it accounts for the historical performance of the individual acquisition functions to select the next candidate solution, initial discrepant values of either good or bad performance can compromise the quality of the selection strategy.

In this work we aim to propose a modified portfolio-based BO methodology that overcomes the GP-Hedge limitations while maintaining its ease of use. Our approach reduces the influence of far past evaluations to enable the recovery of initially bad performing acquisition functions and to avoid the dominance of initially good performing functions. Moreover, our proposal, named No-PASt-BO (Normalized Portfolio Allocation Strategy for Bayesian Optimization), presents a built-in normalization mechanism to avoid the functions within the portfolio to have similar probabilities of being chosen, which promotes exploration.

The new No-PASt-BO approach is empirically evaluated in the task of optimizing synthetic benchmark functions. We also consider the task of optimizing the hyperparameters of machine learning models in real world applications. The obtained results indicate that No-PASt-BO presents competitive performance and always outperform GP-Hedge.

The remaining of the paper is organized as follows. Section II summarizes the required theoretical background; Section III illustrates the GP-Hedge limitations and details the proposed No-PASt-BO methodology; Section IV presents and discuss the performed computational experiments; Section V concludes the paper with pointers for further investigations.

## II. THEORETICAL BACKGROUND

In this section we summarize the main theoretical aspects of the BO framework, including GP basics, the role of the acquisition function and the original GP-Hedge algorithm.

## A. Gaussian Process Basics

In a standard single dimension output regression setting, we aim to obtain a mapping f : R D → R from a set of N inputs x i ∈ R D , organized in a matrix X ∈ R N × D , to a set of N correspondent outputs f i ∈ R . However, we usually only observe y ∈ R N , a noisy version of the vector f . Considering a Gaussian observation noise glyph[epsilon1] ∼ N (0 , σ 2 ) we have:

$$y _ { i } = f _ { i } + \epsilon, \text{ where } f _ { i } = f ( x _ { i } ), \text{ } 1 \leq i \leq N. \text{ } \text{ } \text{ } ( 1 ) \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{  \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{} \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{
}$$

In the GP modeling framework we choose a multivariate Gaussian prior for the latent (non-observed) vector f . If a zero mean prior is chosen, we get [11]:

$$p ( f | X ) & = \mathcal { N } ( f | 0, K _ { f } ), \quad \text{$\text{$an eas}$} \\ p ( y | f, X ) & = \mathcal { N } ( y | f, \sigma ^ { 2 } I ) \mathcal { N } ( f | 0, K _ { f } ), \quad \text{$\text{$quisiti}$} \\ p ( y | X ) & = \mathcal { N } ( y | 0, K _ { f } + \sigma ^ { 2 } I ), \quad \text{$(2)} \quad \text{means} \\ \quad \text{me$\mathcal{ }A\sim} }$$

where in Eq. (2) we were able to analytically integrate out f . The elements of the covariance matrix K f ∈ R N × N are calculated by [ K f ] ij = k ( x i , x j ) , ∀ i, j ∈ { 1 , · · · , N } , where k ( · , · ) is the so-called covariance (or kernel ) function.

The kernel hyperparameters are usually optimized following the gradients of the log-marginal likelihood, i.e., the logarithm of Eq. (2), also called the evidence of the model.

Given a new input x ∗ ∈ R D , the posterior predictive distribution of the related output f ∗ ∈ R is calculated analytically using standard Gaussian distribution conditioning properties:

$$\mu _ { \ast } & p ( f _ { * } | y, X, x _ { * } ) & = \mathcal { N } \left ( f _ { * } | \mu _ { * }, \sigma _ { * } ^ { 2 } \right ), & ( 3 ) \quad \text{literatu} \\ \mu _ { * } & = k _ { * f } ( K _ { f } + \sigma _ { y } ^ { 2 } I ) ^ { - 1 } y, & \quad \text{detail t} \\ \sigma _ { * } ^ { 2 } & = K _ { * } - k _ { * f } ( K _ { f } + \sigma _ { y } ^ { 2 } I ) ^ { - 1 } k _ { f * }, & \quad \text{detail t} \\ \mu _ { * } & \quad \text{propos} \quad$$

where k f ∗ = [ k ( x ∗ , x 1 ) , · · · , k ( x ∗ , x N )] glyph[latticetop] ∈ R N , k ∗ f = k glyph[latticetop] f ∗ and K ∗ = k ( x ∗ , x ∗ ) ∈ R . Importantly, each prediction is a fully defined distribution, instead of a point estimate, which reflects the inherent uncertainty of the regression problem.

## B. The Bayesian Optimization Framework

BO is a general framework for black-box optimization. Mathematically, the problem consists in finding a global minimizer (or maximizer) of an unknown loss function [15]

$$x ^ { * } = \underset { x \in \mathcal { X } } { \arg \min } l ( x ),$$

where l ( ) · : X → R denotes the loss function (or objective, when maximizing), X is the search space, usually given by a compact subset of R D , and x ∗ ∈ R D denotes the optimal solution. The loss function is usually assumed to be either hard to compute or have no simple closed form, but it can be evaluated at an arbitrary query point x . Similar to Eq. (1), the BO framework also considers cases in which we do not observe l ( ) · directly, but rather noisy observations.

BO solves the problem in Eq. (4) by sequentially querying the loss function as we keep the best-so-far candidate solution x + . In doing so, at iteration t , we select a new location x t +1 at which the loss function l ( x t +1 ) is evaluated. As the method iterates, the querying points and its corresponding loss values D t = { x i , l ( x i ) } t i =1 are used for modeling the loss function. When a stopping criterion has been achieved, we return x + as an approximation to the actual minimizer x ∗ . A fundamental aspect is how to provide location guesses along iterations. For that, a probabilistic model is necessary, since l ( ) · is unknown.

The most common probabilistic approach used in the BO framework is the GP model, summarized in Section II-A. The GP model is able to quantify the uncertainty with respect the function l ( ) · , especially at locations at which the function was not yet evaluated and knowledge is scarce.

The uncertainty about the loss is used to select the most promising candidate point for evaluation. It is achieved by an easy-to-compute acquisition function . Typically, such acquisition functions are model-derived and are used to tradeoff exploration and exploitation; in the sense that exploration means investigate non-explored areas (with high uncertainty), and exploitation refers to considering regions where the model prediction is high. The goal is to maximize the acquisition function, and it can be achieved because such functions are assumed to be cheap to evaluate, usually with gradient information available.

We refer the reader to the comprehensive survey in [16] for more details and challenges in the general BO framework.

## C. Acquisition Functions

Several acquisition functions have been proposed in the literature. Each proposal aims at measuring the quality of the candidates to be queried following a specific strategy. Next we detail three of the most used acquisition functions in practice.

1) Probability of improvement (PI): The PI function, firstly proposed in [17], focus on choosing the domain point with the highest probability of being lower than µ -= min i µ ( x i ) , where µ ( x i ) indicates the GP predicted mean at the input x i . This formulation focus on exploitation [12], which can be balanced with the trade-off hyperparameter ξ ≥ 0 as follows:

$$\text{PI} ( x ) = P ( f ( x ) \leq \mu ^ { - } - \xi ) = \Phi \left ( \frac { \mu ^ { - } - \xi - \mu ( x ) } { \sigma ( x ) } \right ), \ \ ( 5 ) \ \overline { \text{Algorit} }$$

where Φ is the CDF of a standard Gaussian distribution.

2) Expected improvement (EI): The EI function, introduced in [18], considers the probability of an evaluation being lower than the current best known evaluation, but it also takes into account the magnitude of the improvement. Let µ -= min i µ ( x i ) , the EI function is zero if σ ( x ) = 0 , otherwise it is given by

$$\text{EI} ( x ) & = \tau ( x ) \Phi \left ( \frac { \tau ( x ) } { \sigma ( x ) } \right ) + \sigma ( x ) \phi \left ( \frac { \tau ( x ) } { \sigma ( x ) } \right ) \quad ( 6 ) \quad \text{A} \\ \text{where } \tau ( x ) & = \mu ^ { - } - \xi - \mu ( x )$$

In Eq. (6), Φ and φ are respectively the CDF and the PDF of a standard Gaussian distribution.

3) GP Lower confidence bound (LCB): In [19] it is introduced the 'Sequential Design for Optimization' (SDO), which selects a point to evaluate based on the posterior mean and variance, minimizing µ ( x ) -κσ ( x ) . In the original paper, κ is a hyperparameter, but no clues are given in how to select it. In [20], the SDO algorithm is revisited and distinct approaches to select a value for κ are discussed. Thus, the so-called GP-LCB formulation is presented below:

$$\text{GP-LCB} ( x ) = \mu ( x ) - \sqrt { \nu \beta _ { t } } \sigma ( x ), \quad \ \ ( 7 ) \, \begin{array} { c c } \text{or one} \\ \text{Thus}, \end{array}$$

where we have considered κ = √ νβ t , β t = 2 log( t D/ 2+2 π / δ 2 3 ) varies with the sequential iteration t and ν, δ &gt; 0 [15].

## D. GP-Hedge

As previously mentioned, there is not a single acquisition function which is the best choice for all possible optimization tasks. Therefore, a strategy which can choose among a set of acquisition functions may be a good direction to handle this issue. The GP-Hedge, introduced in [12], follows such approach.

In the GP-Hedge framework, a set of predefined acquisition functions is considered, comprising a portfolio . Each function nominates a candidate for the next point x of the domain to be evaluated. The candidates are then selected with a probability proportional to how good the posterior mean of the previous points the corresponding acquisition function has suggested before.

The aforementioned approach follows the hedge strategy. According to [21], such method consists in choosing the action j among J options with probability p j ∝ exp( ηG j ( )) t , where η is a hyperparameter and G j ( ) t = ∑ t t ′ =1 score ( j x j ( t ′ )) is the total score of the action j up to the time t .

The original GP-Hedge algorithm proposed in [12] was defined to solve a maximization problem, so the reward of each acquisition function is equal to the sum of the previous posterior means, i.e., G j ( ) = t ∑ t t ′ =1 µ j ( x j ( t ′ )) . Since we define our tasks as minimization problems, we multiply each score by -1 before adding it, i.e., G j ( ) = t -∑ t t ′ =1 µ j ( x j ( t ′ )) . The full algorithm is detailed in Algorithm 1.

## Algorithm 1 GP-Hedge

Select hyperparameter η ∈ R +

Set

G

j

(0) = 0

for

j

= 1 2

,

, . . . , J

for t = 1,2,. . . do

Nominate points from each acquisition function h j :

$$\begin{array} { c } \text{$num_poms_normal_coverage_common_name.n_{j}$}. \\ \text{$x_{j}(t)=argmax_{x}h_{j}(x)$.} \\ \text{$Select a nominee $x(t)=x_{j}(t)$ with probability} \\ p_{j}(t)=\frac{\exp(\eta G_{j}(t-1))} \\ \sum _ { j ^ { \prime } = 1 } ^ { j } \exp(\eta G _ { j ^ { \prime } } ( t - 1 ) ) }. \\ \text{$Commute $u(t)$} \text{hv evaluating the objective on noint $x(t)$.} \end{array}$$

Compute y t ( ) by evaluating the objective on point x ( ) t .

Augment the data

D

t

with the new pair

(

x

( )

t

, y

( ))

t

.

Update the surrogate GP model.

Update the rewards

G

j

( ) =

t

G

j

(

the updated GP posterior.

end for

## III. PROPOSED METHODOLOGY

## A. GP-Hedge limitations

The original GP-Hedge algorithm relies on the cumulative performance of each acquisition function in the portfolio during previous iterations to favor the choice of a function over the others. However, for large enough horizons, the influence of the first iterations may not be relevant or even desirable. Thus, we want to reduce the importance of an iteration in the reward function as the former becomes more distant from the current iteration. Importantly, it is emphasized in [13] that the reward function is critical for the GP-Hedge performance, which encourages our argument.

In Fig. 1 it is illustrated an example of how the scores of each acquisition function evolve with the iterations. We consider a 3-function portfolio comprised by the PI, EI and GP-LCB functions and the well known Hartmann 6 benchmark. The corresponding probabilities of being chosen are also presented. It is possible to note that for the GP-Hedge the GPLCB function finishes with a large score lead over the other 2 functions. In comparison with the versions with a memory factor, we note that this lead is smaller and can still be lost. Both score graphs for the versions with memory factors are very similar. However, the probability graphs are different, since in the normalized version it is more uncommon to obtain situations of equal probabilities, which would result in a completely random choice between the 3 available functions. Both the memory factor and normalization mechanisms will be detailed in the next sections.

## B. Memory Factor

In order to tackle the presented GP-Hedge issues, we propose to change the reward function update by including a memory factor. The main goal of the memory factor is to decrease the influence of previous iterations in the reward evaluation as new iterations are completed, while still considering past experiences. Thus, we aim to avoid that discrepant

t

-

1)

-

µ

(

x

j

( ))

t

from

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00361v1-figure-c7b0bce5-de81-4a64-bb95-fa82f30c7c3f.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00361v1-figure-ec563b90-2429-4355-aa3b-2e3b599ab4d9.jpg)

## 0.8 memory factor normalized

Fig. 1. In the first row, it is presented the scores of each acquisition function along the 100 iterations of a typical run for the Hartmann 6 benchmark function. The first column corresponds to the GP-Hedge, the second column corresponds to the GP-Hedge with a memory factor of 0.8, and, finally the third column corresponds to the memory factor of 0.8 using the normalization strategy, which consists in our approach. In the second row, the correspondent probabilities of each acquisition function be chosen is shown.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00361v1-figure-f4c1706e-804f-4ba7-bb47-3ba4ea5b093f.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00361v1-figure-8ebcbef6-e3df-43ad-93bd-b4b355be66e7.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00361v1-figure-3b047de3-a756-48bf-a862-227ef162de83.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00361v1-figure-a4619089-2331-411c-af92-13b6619dfafb.jpg)

values in the first few iterations of being determinant in the acquisition function selection during all the later iterations. Our approach also enables more probability for acquisition functions that are better in the recent past, according to the memory factor value, enabling 'recovery' from far past mistakes.

iterations, as illustrated in Fig. 1. In such scenario, all the acquisition functions have about the same probability of being chosen. In the extreme case where all the functions have equal probability we would get an undesired completely random portfolio behavior.

More specifically, we propose to change the reward update computation as follows:

$$G _ { j } ( t ) = m G _ { j } ( t - 1 ) - \mu ( \boldsymbol x _ { j } ( t ) ), \text{ \quad \ \ } ( 8 )$$

where 0 ≤ m ≤ 1 is the memory factor hyperparameter, µ ( x j ( )) t is the GP posterior mean for the input x j ( ) t suggested by acquisition function h j at iteration t . Typical values for the memory factor hyperparameter are between 0 7 . and 1 , where the latter recovers the original GP-Hedge.

Eq. (8) imposes a decrease in relevance to past rewards. For m &lt; 1 we achieve two behaviors that are difficult to observe in the original GP-Hedge: ( ) initially bad acquisition i functions may receive some probability in the later iterations if they improve along the optimization steps; ( ii ) acquisition functions that go very well in the beginning may lose the preference if they are not able to keep the good performance. Those behaviors can be observed in Fig. 1, where we can note, for instance, that the EI function is able to recover some probability after the initial iterations in the variants that include the memory factor.

## C. Rewards Normalization

The use of the previously presented memory factor may cause the rewards of the portfolio to be very close at some

In order to solve this problem, we propose to normalize the reward functions values before computing the choice probabilities as follows:

$$r _ { j } ( t ) & = \frac { G _ { j } ( t ) - r _ { \max } ( t ) } { r _ { \max } ( t ) - r _ { \min } ( t ) }, \quad \text{(9)} \\ & \text{where } r _ { \min } ( t ) & = \min _ { j } ( G _ { j } ( t ) ), \\ & \text{$r_{\max} ( t ) = \max _ { j } ( G _ { j } ( t ) ). \\ \ell \quad \cdot \quad \cdot \quad \ell \cdot \quad \ell \cdot \quad \cdot \quad \cdot.$$

In the former expressions, the term r j ( ) t indicates the normalized reward for the acquisition function j after the iteration t . Those values, computed for each acquisition function h j | J j =1 , are then used to obtain the probabilities of each acquisition function being chosen as follows:

$$p _ { j } ( t ) = \frac { \exp ( \eta r _ { j } ( t - 1 ) } { \sum _ { j ^ { \prime } = 1 } ^ { J } \exp ( \eta r _ { j ^ { \prime } } ( t - 1 ) ) }. \quad \quad ( 1 0 )$$

Note that the normalized r j ( ) t values are considered only to compute Eq. (10). The original values of G j ( ) t are not overwritten.

The proposed normalization step constrains the terms r j ( ) t to be between -1 and 0 , with the highest value always being 0 and the lowest always being -1 . To solve a possible division by 0 , if the highest and lowest rewards are equal, all probabilities are set equally.

The proposed changes in the rewards function of the base GP-Hedge approach result in the proposal of the No-PASt-BO (Normalized Portfolio Allocation Strategy for Bayesian Optimization) algorithm, which is detailed in Algorithm 2. Note that the proposed method presents the same computational cost of the original GP-Hedge.

## Algorithm 2 No-PASt-BO

Select hyperparameter η ∈ R + Select hyperparameter m ∈ [0 , 1] Set G j (0) = 0 for j = 1 2 , , . . . , J for t = 1,2,. . . do Nominate points from each acquisition function h j : x j ( ) = arg max t x h j ( x ) Compute r min ( t -1) = min ( j G j ( t -1)) Compute r max ( t -1) = max ( j G j ( t -1)) Compute the normalized rewards: r j ( t -1) = G j ( t -1) -r max ( t -1) r max ( t -1) -r min ( t -1) Select a nominee x ( ) = t x j ( ) t with probability p j ( ) = t exp( ηr j ( t -1) ∑ J j ′ =1 exp( ηr j ′ ( t -1)) . Compute y t ( ) by evaluating the objective on point x ( ) t . Augment the data D t with the new pair ( x ( ) t , y ( )) t . Update the surrogate GP model. Update the rewards G j ( ) = t mG t j ( -1) -µ ( x j ( )) t from the updated GP posterior. end for

## IV. EXPERIMENTS

In order to evaluate the proposed No-PASt-BO performance a battery of experiments were made in synthetic benchmark function optimization. We also considered the real world application of BO in the task of optimizing the hyperparameters of machine learning models. For each test 25 runs were executed and the mean logarithmic error is reported along with the correspondent obtained confidence interval, i.e., the standard deviation scaled by the square root of the number of runs.

As baselines, we also evaluate a random portfolio (RP) approach, which chooses randomly among a set of predefined acquisition functions, the original GP-Hedge and standard BO with a single acquisition function. For all of the experiments, we have used the GpyOpt package, a general BO framework introduced in [22]. After preliminary experiments, the hyperparameter η was set to 4 for all of the No-PASt-BO experiments, while for the GP-Hedge the η was defined using the strategy suggested in [12].

## A. Memory Factor Sensibility

Before the comparison results, we wish to explore the impact of the memory factor m in the No-PASt-BO method. Thus, we perform an initial experiment with 7 different values for the memory factor: { 0 7 . , 0 75 0 8 . , . , 0 85 0 9 . , . , 0 95 1 . , } . Note that when the memory factor is equal to 1 , it means that previous executions do not lose importance. With the exception of the normalizing step (see Section III-C), the latter is similar to the GP-Hedge. As the memory factor decreases, the importance of previous evaluations also decreases.

The impact of different the memory factor values is presented in Fig. 2 for the Hartmann 6 benchmark function. It is possible to verify how the rewards change over the iterations, and how it affects the probability of the correspondent acquisition function being chosen. Higher values of memory factor make it more difficult to select an acquisition function as the most probable one. Moreover, lower values for the memory factor imply in more quickly varying probabilities, which may improve the diversity in the acquisition function selection step.

## B. Synthetic Benchmark Functions

In this section 3 standard benchmark functions were used: Branin, Hartmann 3 and Hartmann 6 1 . Their domains are respectively 2, 3 and 6 dimensional. These functions were the same chosen in the original GP-Hedge paper [12].

In all experiments, the portfolio-based strategies used the same set of acquisition functions: PI, EI and GP-LCB. We have set ξ = 0 01 . for PI and EI, and δ = 0 1 . and ν = 0 2 . for GP-LCB. Furthermore, we evaluate the proposed NoPASt-BO with 7 different memory factors. For each of the experiments, 5 initial points were selected using the Latin Hypercube Sampling (LHS) approach [23].

The results of these experiments can be found in Fig. 3. For all of scenarios, a portfolio strategy obtained the best performance. A version of the No-PASt-BO achieved the best result in two of the three case, and was very close to the best version in the case it loses. Also, the No-PASt-BO with memory factor of 0 7 . outperforms GP-Hedge in all of the executed experiments.

A second battery of tests was run for the same synthetic functions, but using 9 acquisition functions in the portfolio variants. These functions come from the same set of PI, EI and GP-LCB adding ξ = 0 1 . , ξ = 1 0 . , ν = 0 1 . and ν = 1 to the already studied values. The results are show in Fig. 4. Overall, the Random Portfolio performance was worse than the previous experiment with only 3 functions, with the remarkable exception of the Hartmann 3 experiment, due to the fact that as more functions were inserted in the portfolios, some of them are worse to the task at hand. In the same way, we can see that GP-Hedge got more competitive in general, but the No-PAST-BO still reaches lower optimized values faster.

## C. Real World Problems

In order to evaluate the No-PASt-BO and compare it with the available baselines in a real world problem, we first consider a regression setting with the standard Boston Housing dataset 2 . The Support Vector Regression (SVR) model [24] is used to predict houses prices from the available attributes. The BO strategies have the task of optimizing the 3 SVR

1 Functions definitions available at https://www.sfu.ca/ ∼ ssurjano/ optimization.html.

2 Available at https://www.cs.toronto.edu/ ∼ delve/data/boston/bostonDetail. html.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00361v1-figure-08c2ea3f-0d2e-4012-93ac-3b2f362bd196.jpg)

100

1

1

1

0.75

0 50

0 25

100

0 50

000

100

0.75

0 50

0 25

000

1

1

1

0.75

0 25

000

100

0 75

0 50

0.25

000

100

0 50

0 25

0.00

100

50

GP-LCB

75

100

Iterations

No-PASt 0.8

GP-LCB

WW

50

75

Iterations

No-PASt 0.9

GP-LCB

25

50

Iterations

1

0.6

0 4

0 2

00

75

No-PASt 1

25

50

50

GP-LCB

75

100

Iterations

No-PASt 0.85

GP-LCB

50

75

100

Iterations

No-PASt 0.95

GP-LCB

25

GP-LCB

75

100

Iterations

Fig. 2. The effect of different memory factors in the NO-PASt-BO execution. As the memory factor increases, it takes more time to select an acquisition function as the one with higher probability.

50

75

Iterations

100

100

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00361v1-figure-66423373-2f7f-41fc-a947-bf66cce106b0.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00361v1-figure-3c1c913e-a0cf-42c8-abe3-74bce2cda5cd.jpg)

Hartmann 6

Fig. 3. Results of the minimization of the synthetic benchmark using 3 acquisition functions. Note that only the 2 best choices for the NO-PAStBO memory factor are shown. Also, only the best non-portfolio acquisition function is presented.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00361v1-figure-af357730-bca3-48f6-bb5b-1b01c7ae4019.jpg)

hyperparameters (gamma, C and epsilon), as implemented in the scikit learn toolbox [25].

In this experiment we have used a 10 -fold cross-validation procedure, where the objective of each BO step is to minimize the average root mean square error (RMSE). For each experiment 100 optimization iterations were performed and each experiment was repeated 25 times. The obtained average RMSE values and their confidence intervals are shown in Fig. 5. We can see that the NO-PASt-BO variants achieved the best results both in terms of best final solution and faster lower

Fig. 4. Results of the minimization of the synthetic benchmark using 9 acquisition functions. Note that only the 2 best choices for the NO-PAStBO memory factor are shown. Also, only the best non-portfolio acquisition function is presented.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00361v1-figure-6f6e2cce-bf81-411a-9085-d0845e5f5fc0.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00361v1-figure-7b640262-cbd8-4897-b5df-95b629eb0fe5.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00361v1-figure-bd251e9b-ec42-497b-b0b5-294545e0f8d4.jpg)

values. Moreover, the GP-Hedge performed comparable to the the random portfolio strategy.

We tackled a second real world problem to compare the evaluated BO methods in a regression setting which consists in predicting the average gearbox high-speed shaft temperature in a Wind Turbine Generator (WTG) from a given set of measures. Following the same methodology presented in [26], [27], the input variables were: average active power, average rotor speed, average nacelle (the turbine outer casing) temperature and average outdoor temperature. A limit filter

## Boston Housing SVM

Fig. 5. Results for automatic tuning of SVM and MLP hyperparameters. For both cases a No-PASt-BO version has achieved the best performance.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00361v1-figure-05338db5-a74b-421f-99e5-c45d0a2601b1.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00361v1-figure-f99302e1-c15c-479b-910b-ad4b6cf25db9.jpg)

was applied to the available measured data, removing values that are physically absurd. After that, we applied a clustering filter to remove some outliers from the dataset and, finally, we applied a continuity filter in order to remove isolated points. Details of these preprocessing steps can be seen in [26], [27].

We considered a nonlinear model comprised of a Multilayer Perceptron (MLP) network with a single hidden layer and hyperbolic tangent hidden activations. The model was trained via stochastic gradient descent optimizer. The task of the BO was to minimize the RMSE of the predicted temperature in a hold-out validation set, varying over the scikit-learn MLPRegressor hyperparameters [25]: neurons in the hidden layer, learning rate, alpha (weight decay hyperparameter) and momentum. The data used in this experiment was obtained from a WTG located in Brazil. Results are presented in Fig. 5. Although all the methods obtained somehow comparable results, the detailed ispection indicate that the NO-PASt-BO variants obtained the best RMSE values at the end.

## V. CONCLUSIONS AND FURTHER WORK

In this work we have proposed No-PASt-BO, a new methodology for performing black-box Bayesian optimization of unknown functions. Our approach builds on the popular GPHedge method [12] to create an adaptive portfolio of acquisition functions. We aimed to tackle some of the GP-Hedge's limitations, such as the sometimes undesirable high influence of far past evaluations, by incorporating a limited memory and a normalization mechanisms.

We evaluate the proposed approach on both synthetic and real-world optimization tasks, where No-PASt-BO obtained competitive performance with respect to the other evaluated strategies. Importantly, our approach always outperforms GPHedge, while maintaining its simplicity and general applicability. The latter features may enable No-PASt-BO to become the default off-the-shelf method for portfolio-based BO . 3

Further investigations shall verify how to incorporate nonmyopic concepts to our portfolio BO methodology, such as the ones explored in [28]. That would enable more clever decision making when we have a known limited budget in terms of number of objective evaluations.

Another interesting subject for investigation is the task of function optimization with constraints known a priori . Strategies such as the ones presented in [29] may turn our portfolio framework even more applicable in real world problems.

## REFERENCES

- [1] J. Mockus and L. Mockus, 'Bayesian approach to global optimization and application to multiobjective and constrained problems,' Journal of Optimization Theory and Applications , vol. 70, no. 1, pp. 157-172, 1991.
- [2] J. Snoek, H. Larochelle, and R. P. Adams, 'Practical Bayesian optimization of machine learning algorithms,' in Advances in Neural Information Processing Systems 25 (NIPS) . Lake Tahoe, Nevada, USA: NIPS Foundation, 2012, pp. 2951-2959.
- [3] A. Klein, S. Falkner, S. Bartels, P. Hennig, and F. Hutter, 'Fast bayesian optimization of machine learning hyperparameters on large datasets,' arXiv preprint arXiv:1605.07079 , 2016.
- [4] L. Kotthoff, C. Thornton, H. H. Hoos, F. Hutter, and K. LeytonBrown, 'Auto-weka 2.0: Automatic model selection and hyperparameter optimization in weka,' The Journal of Machine Learning Research , vol. 18, no. 1, pp. 826-830, 2017.
- [5] S. Falkner, A. Klein, and F. Hutter, 'Bohb: Robust and efficient hyperparameter optimization at scale,' in International Conference on Machine Learning , 2018, pp. 1436-1445.
- [6] M. Feurer and F. Hutter, 'Hyperparameter optimization,' in Automated Machine Learning . Springer, 2019, pp. 3-33.
- [7] R. Calandra, A. Seyfarth, J. Peters, and M. P. Deisenroth, 'Bayesian optimization for learning gaits under uncertainty,' Annals of Mathematics and Artificial Intelligence , vol. 76, no. 1-2, pp. 5-23, 2016.
- [8] K. Chatzilygeroudis, R. Rama, R. Kaushik, D. Goepp, V. Vassiliades, and J.-B. Mouret, 'Black-box data-efficient policy search for robotics,' in Intelligent Robots and Systems (IROS), 2017 IEEE/RSJ International Conference on . IEEE, 2017, pp. 51-58.
- [9] W. Lyu, F. Yang, C. Yan, D. Zhou, and X. Zeng, 'Batch Bayesian optimization via multi-objective acquisition ensemble for automated analog circuit design,' in International Conference on Machine Learning , 2018, pp. 3312-3320.
- [10] H. M. Torun, M. Swaminathan, A. K. Davis, and M. L. F. Bellaredj, 'A global bayesian optimization algorithm and its application to integrated system design,' IEEE Transactions on Very Large Scale Integration (VLSI) Systems , vol. 26, no. 4, pp. 792-802, 2018.

3 The No-PASt-BO source code is available at https://github.com/ thiago-vasconcelos/no-past-bo.

- [11] C. Rasmussen and C. Williams, Gaussian Processes for Machine Learning , 1st ed. Cambridge, MA, USA: MIT Press, 2006.
- [12] M. D. Hoffman, E. Brochu, and N. de Freitas, 'Portfolio allocation for bayesian optimization.' in UAI . Citeseer, 2011, pp. 327-336.
- [13] B. Shahriari, Z. Wang, M. W. Hoffman, A. Bouchard-Cˆ ot´ e, and N. de Freitas, 'An entropy search portfolio for bayesian optimization,' Proc. NIPS Workshop Bayesian Optim. , 2014.
- [14] W. Shen, J. Wang, Y.-G. Jiang, and H. Zha, 'Portfolio choices with orthogonal bandit learning,' in Twenty-Fourth International Joint Conference on Artificial Intelligence , 2015.
- [15] E. Brochu, V. M. Cora, and N. De Freitas, 'A tutorial on bayesian optimization of expensive cost functions, with application to active user modeling and hierarchical reinforcement learning,' arXiv preprint arXiv:1012.2599 , 2010.
- [16] B. Shahriari, K. Swersky, Z. Wang, R. P. Adams, and N. De Freitas, 'Taking the human out of the loop: A review of bayesian optimization,' Proceedings of the IEEE , vol. 104, no. 1, pp. 148-175, 2015.
- [17] H. J. Kushner, 'A new method of locating the maximum point of an arbitrary multipeak curve in the presence of noise,' Journal of Basic Engineering , vol. 86, no. 1, pp. 97-106, 1964.
- [18] J. Mo´ ckus, V. Tiesis, and A. Zilinskas, ´ 'The application of bayesian methods for seeking the extremum. vol. 2,' 1978.
- [19] D. D. Cox and S. John, 'Sdo: A statistical method for global optimization,' in in Multidisciplinary Design Optimization: State-of-the-Art , 1997, pp. 315-329.
- [20] N. Srinivas, A. Krause, S. Kakade, and M. Seeger, 'Gaussian process optimization in the bandit setting: No regret and experimental design,' in Proceedings of the 27th International Conference on Machine Learning , no. CONF. Omnipress, 2010.
- [21] P. Auer, N. Cesa-Bianchi, Y. Freund, and R. E. Schapire, 'Gambling in a rigged casino: The adversarial multi-armed bandit problem,' in Proceedings of IEEE 36th Annual Foundations of Computer Science . IEEE, 1995, pp. 322-331.
- [22] The GPyOpt authors, 'GPyOpt: A bayesian optimization framework in python,' http://github.com/SheffieldML/GPyOpt, 2016.
- [23] M. D. McKay, R. J. Beckman, and W. J. Conover, 'Comparison of three methods for selecting values of input variables in the analysis of output from a computer code,' Technometrics , vol. 21, no. 2, pp. 239-245, 1979.
- [24] H. Drucker, C. J. Burges, L. Kaufman, A. J. Smola, and V. Vapnik, 'Support vector regression machines,' in Advances in neural information processing systems , 1997, pp. 155-161.
- [25] F. Pedregosa, G. Varoquaux, A. Gramfort, V. Michel, B. Thirion, O. Grisel, M. Blondel, P. Prettenhofer, R. Weiss, V. Dubourg, J. Vanderplas, A. Passos, D. Cournapeau, M. Brucher, M. Perrot, and E. Duchesnay, 'Scikit-learn: Machine learning in Python,' Journal of Machine Learning Research , vol. 12, pp. 2825-2830, 2011.
- [26] P. Bangalore, S. Letzgus, D. Karlsson, and M. Patriksson, 'An artificial neural network-based condition monitoring method for wind turbines, with application to the monitoring of the gearbox,' Wind Energy , vol. 20, no. 8, pp. 1421-1438, 2017. [Online]. Available: https://onlinelibrary.wiley.com/doi/abs/10.1002/we.2102
- [27] D. A. R. M. A. Souza, T. P. Vasconcelos, C. L. C. Mattos, and J. P. P. Gomes, 'Evaluation of data based normal behavior models for fault detection in wind turbines,' in Brazilian Conference on Intelligent Systems , 2019, p. to appear.
- [28] J. Gonz´ alez, M. Osborne, and N. D. Lawrence, 'Glasses: Relieving the myopia of bayesian optimisation,' 2016.
- [29] J. R. Gardner, M. J. Kusner, Z. E. Xu, K. Q. Weinberger, and J. P. Cunningham, 'Bayesian optimization with inequality constraints.' in ICML , 2014, pp. 937-945.
Thiago de P. Vasconcelos Daniel A. R. M. A. de Souza César L. C. Mattos João P. P. Gomes
No-PASt-BO: Normalized Portfolio Allocation Strategy for Bayesian Optimization
Bayesian Optimization (BO) is a framework for black-box optimization that is
especially suitable for expensive cost functions. Among the main parts of a BO
algorithm, the acquisition function is of fundamental importance, since it
guides the optimization algorithm by translating the uncertainty of the
regression model in a utility measure for each point to be evaluated.
Considering such aspect, selection and design of acquisition functions are one
of the most popular research topics in BO. Since no single acquisition function
was proved to have better performance in all tasks, a well-established approach
consists of selecting different acquisition functions along the iterations of a
BO execution. In such an approach, the GP-Hedge algorithm is a widely used
option given its simplicity and good performance. Despite its success in
various applications, GP-Hedge shows an undesirable characteristic of
accounting on all past performance measures of each acquisition function to
select the next function to be used. In this case, good or bad values obtained
in an initial iteration may impact the choice of the acquisition function for
the rest of the algorithm. This fact may induce a dominant behavior of an
acquisition function and impact the final performance of the method. Aiming to
overcome such limitation, in this work we propose a variant of GP-Hedge, named
No-PASt-BO, that reduce the influence of far past evaluations. Moreover, our
method presents a built-in normalization that avoids the functions in the
portfolio to have similar probabilities, thus improving the exploration. The
obtained results on both synthetic and real-world optimization tasks indicate
that No-PASt-BO presents competitive performance and always outperforms
GP-Hedge.