X-Armed Bandits
Sébastien Bubeck, Rémi Munos, Gilles Stoltz, Csaba Szepesvari
Introduction
In the classical stochastic bandit problem a gambler tries to maximize his revenue by sequentially playing one of a finite number of slot machines that are associated with initially unknown (and potentially different) payoff distributions . Assuming old-fashioned slot machines, the gambler pulls the arms of the machines one by one in a sequential manner, simultaneously learning about the machines’ payoff-distributions and gaining actual monetary reward. Thus, in order to maximize his gain, the gambler must choose the next arm by taking into consideration both the urgency of gaining reward (“exploitation”) and acquiring new information (“exploration”).
Maximizing the total cumulative payoff is equivalent to minimizing the (total) regret, i.e., minimizing the difference between the total cumulative payoff of the gambler and the one of another clairvoyant gambler who chooses the arm with the best mean-payoff in every round. The quality of the gambler’s strategy can be characterized as the rate of growth of his expected regret with time. In particular, if this rate of growth is sublinear, the gambler in the long run plays as well as the clairvoyant gambler. In this case the gambler’s strategy is called Hannan consistent.
Bandit problems have been studied in the Bayesian framework , as well as in the frequentist parametric [25; 2] and non-parametric settings , and even in non-stochastic scenarios [5; 10]. While in the Bayesian case the question is whether the optimal actions can be computed efficiently, in the frequentist case the question is how to achieve low rate of growth of the regret in the lack of prior information, i.e., it is a statistical question. In this paper we consider the stochastic, frequentist, non-parametric setting.
Although the first papers studied bandits with a finite number of arms, researchers have soon realized that bandits with infinitely many arms are also interesting, as well as practically significant. One particularly important case is when the arms are identified by a finite number of continuous-valued parameters, resulting in online optimization problems over continuous finite-dimensional spaces. Such problems are ubiquitous to operations research and control. Examples are “pricing a new product with uncertain demand in order to maximize revenue, controlling the transmission power of a wireless communication system in a noisy channel to maximize the number of bits transmitted per unit of power, and calibrating the temperature or levels of other inputs to a reaction so as to maximize the yield of a chemical process” . Other examples are optimizing parameters of schedules, rotational systems, traffic networks or online parameter tuning of numerical methods. During the last decades numerous authors have investigated such “continuum-armed” bandit problems [3; 21; 6; 22; 12]. A special case of interest, which forms a bridge between the case of a finite number of arms and the continuum-armed setting, is formed by bandit linear optimization, see and the references therein.
In many of the above-mentioned problems, however, the natural domain of some of the optimization parameters is a discrete set, while other parameters are still continuous-valued. For example, in the pricing problem different product lines could also be tested while tuning the price, or in the case of transmission power control different protocols could be tested while optimizing the power. In other problems, such as in online sequential search, the parameter-vector to be optimized is an infinite sequence over a finite alphabet [13; 7].
The motivation for this paper is to handle all these various cases in a unified framework. More precisely, we consider a general setting that allows us to study bandits with almost no restriction on the set of arms. In particular, we allow the set of arms to be an arbitrary measurable space. Since we allow non-denumerable sets, we shall assume that the gambler has some knowledge about the behavior of the mean-payoff function (in terms of its local regularity around its maxima, roughly speaking). This is because when the set of arms is uncountably infinite and absolutely no assumptions are made on the payoff function, it is impossible to construct a strategy that simultaneously achieves sublinear regret for all bandits problems (see, e.g., [9, Corollary 4]). When the set of arms is a metric space (possibly with the power of the continuum) previous works have assumed either the global smoothness of the payoff function [3; 21; 22; 12] or local smoothness in the vicinity of the maxima . Here, smoothness means that the payoff function is either Lipschitz or Hölder continuous (locally or globally). These smoothness assumptions are indeed reasonable in many practical problems of interest.
In this paper, we assume that there exists a dissimilarity function that constrains the behavior of the mean-payoff function, where a dissimilarity function is a measure of the discrepancy between two arms that is neither symmetric, nor reflexive, nor satisfies the triangle inequality. (The same notion was introduced simultaneously and independently of us by [23, Section 4.4] under the name “quasi-distance.”) In particular, the dissimilarity function is assumed to locally set a bound on the decrease of the mean-payoff function at each of its global maxima. We also assume that the decision maker can construct a recursive covering of the space of arms in such a way that the diameters of the sets in the covering shrink at a known geometric rate when measured with this dissimilarity.
Our work generalizes and improves previous works on continuum-armed bandits.
In particular, Kleinberg and Auer et al. focused on one-dimensional problems, while we allow general spaces. In this sense, the closest work to the present contribution is that of Kleinberg et al. , who considered generic metric spaces assuming that the mean-payoff function is Lipschitz with respect to the (known) metric of the space; its full version relaxed this condition and only requires that the mean-payoff function is Lipschitz at some maximum with respect to some (known) dissimilarity. The present paper paper is a concurrent and independent work with respect to the paper of Kleinberg, Slivkins, and Upfal . An extended abstract of the latter was published in May 2008 at STOC’08, while the NIPS’08 version of the present paper was submitted at the beginning of June 2008. At that time, we were not aware of the existence of the full version , which was released in September 2008. Kleinberg et al. proposed a novel algorithm that achieves essentially the best possible regret bound in a minimax sense with respect to the environments studied, as well as a much better regret bound if the mean-payoff function has a small “zooming dimension”.
Our contribution furthers these works in two ways:
our algorithms, motivated by the recent successful tree-based optimization algorithms [24; 18; 13], are easy to implement;
we show that a version of our main algorithm is able to exploit the local properties of the mean-payoff function at its maxima only, which, as far as we know, was not investigated in the approach of Kleinberg et al. .
The precise discussion of the improvements (and drawbacks) with respect to the papers by Kleinberg et al. requires the introduction of somewhat extensive notations and is therefore deferred to Section 5. However, in a nutshell, the following can be said.
First, by resorting to a hierarchical approach, we are able to avoid the use of the doubling trick, as well as the need for the (covering) oracle, both of which the so-called zooming algorithm of Kleinberg et al. relies on. This comes at the cost of slightly more restrictive assumptions on the mean-payoff function, as well as a more involved analysis. Moreover, the oracle is replaced by an a priori choice of a covering tree. In standard metric spaces, such as the Euclidean spaces, such trees are trivial to construct, though, in full generality they may be difficult to obtain when their construction must start from (say) a distance function only. We also propose a variant of our algorithm that has smaller computational complexity of order compared to the quadratic complexity of our basic algorithm. However, the cheaper algorithm requires the doubling trick to achieve an anytime guarantee (just like the zooming algorithm).
Second, we are also able to weaken our assumptions and to consider only properties of the mean-payoff function in the neighborhoods of its maxima; this leads to regret bounds scaling as \widetilde{O}\bigl{(}\sqrt{n}\bigr{)} We write when up to a logarithmic factor. when, e.g., the space is the unit hypercube and the mean-payoff function has a finite number of global maxima around which it is locally equivalent to a function with some known degree . Thus, in this case, we get the desirable property that the rate of growth of the regret is independent of the dimensionality of the input space. (Comparable dimensionality-free rates are obtained under different assumptions in .)
Finally, in addition to the strong theoretical guarantees, we expect our algorithm to work well in practice since the algorithm is very close to the recent, empirically very successful tree-search methods from the games and planning literature [16; 17; 27; 11; 15].
In Section 2 we formalize the –armed bandit problem.
In Section 3 we describe the basic strategy proposed, called HOO (hierarchical optimistic optimization).
We present the main results in Section 4. We start by specifying and explaining our assumptions (Section 4.1) under which various regret bounds are proved. Then we prove a distribution-dependent bound for the basic version of HOO (Section 4.2). A problem with the basic algorithm is that its computational cost increases quadratically with the number of time steps. Assuming the knowledge of the horizon, we thus propose a computationally more efficient variant of the basic algorithm, called truncated HOO and prove that it enjoys a regret bound identical to the one of the basic version (Section 4.3) while its computational complexity is only log-linear in the number of time steps. The first set of assumptions constrains the mean-payoff function everywhere. A second set of assumptions is therefore presented that puts constraints on the mean-payoff function only in a small vicinity of its global maxima; we then propose another algorithm, called local-HOO, which is proven to enjoy a regret again essentially similar to the one of the basic version (Section 4.4). Finally, we prove the minimax optimality of HOO in metric spaces (Section 4.5).
In Section 5 we compare the results of this paper with previous works.
Problem setup
The mean-payoff function thus defined is assumed to be measurable. For simplicity, we shall also assume that all have bounded supports, included in some fixed bounded intervalMore generally, our results would also hold when the tails of the reward distributions are uniformly sub-Gaussian. , say, the unit interval $f$.
A decision maker (the gambler of the introduction) that interacts with a stochastic bandit problem plays a game at discrete time steps according to the following rules. In the first round the decision maker can select an arm and receives a reward drawn at random from . In round the decision maker can select an arm based on the information available up to time , i.e., , and receives a reward drawn from , independently of given . Note that a decision maker may randomize his choice, but can only use information available up to the point in time when the choice is made.
Formally, a strategy of the decision maker in this game (“bandit strategy”) can be described by an infinite sequence of measurable mappings, , where maps the space of past observations,
to the space of probability measures over . By convention, does not take any argument. A strategy is called deterministic if for every , is a Dirac distribution.
The goal of the decision maker is to maximize his expected cumulative reward. Equivalently, the goal can be expressed as minimizing the expected cumulative regret, which is defined as follows. Let
be the best expected payoff in a single round. At round , the cumulative regret of a decision maker playing is
Finally, we define the cumulative pseudo-regret as
that is, the actual rewards used in the definition of the regret are replaced by the mean-payoffs of the arms pulled. Since (by the tower rule)
As it is argued in , in many real-world problems, the decision maker is not interested in his cumulative regret but rather in its simple regret. The latter can be defined as follows. After rounds of play in a stochastic bandit problem , the decision maker is asked to make a recommendation based on the obtained rewards . The simple regret of this recommendation equals
In this paper we focus on the cumulative regret , but all the results can be readily extended to the simple regret by considering the recommendation , where is drawn uniformly at random in . Indeed, in this case,
The Hierarchical Optimistic Optimization (HOO) strategy
The HOO strategy (cf. Algorithm 1) incrementally builds an estimate of the mean-payoff function over . The core idea (as in previous works) is to estimate precisely around its maxima, while estimating it loosely in other parts of the space . To implement this idea, HOO maintains a binary tree whose nodes are associated with measurable regions of the arm-space such that the regions associated with nodes deeper in the tree (further away from the root) represent increasingly smaller subsets of . The tree is built in an incremental manner. At each node of the tree, HOO stores some statistics based on the information received in previous rounds. In particular, HOO keeps track of the number of times a node was traversed up to round and the corresponding empirical average of the rewards received so far. Based on these, HOO assigns an optimistic estimate (denoted by ) to the maximum mean-payoff associated with each node. These estimates are then used to select the next node to “play”. This is done by traversing the tree, beginning from the root, and always following the node with the highest –value (cf. lines 4–14 of Algorithm 1). Once a node is selected, a point in the region associated with it is chosen (line 16) and is sent to the environment. Based on the point selected and the received reward, the tree is updated (lines 18–33).
The tree of coverings which HOO needs to receive as an input is an infinite binary tree whose nodes are associated with subsets of . The nodes in this tree are indexed by pairs of integers ; node is located at depth from the root. The range of the second index, , associated with nodes at depth is restricted by . Thus, the root node is denoted by . By convention, and are used to refer to the two children of the node . Let be the region associated with node . By assumption, these regions are measurable and must satisfy the constraints
As a corollary, the regions at any level cover the space ,
In the algorithm listing the recursive computation of the –values (lines 28–33) makes a local copy of the tree; of course, this part of the algorithm could be implemented in various other ways. Other arbitrary choices in the algorithm as shown here are how tie breaking in the node selection part is done (lines 9–12), or how a point in the region associated with the selected node is chosen (line 16). We note in passing that implementing these differently would not change our theoretical results.
To facilitate the formal study of the algorithm, we shall need some more notation. In particular, we shall introduce time-indexed versions (, , , , , etc.) of the quantities used by the algorithm. The convention used is that the indexation by is used to indicate the value taken at the end of the round.
In particular, is used to denote the finite subtree stored by the algorithm at the end of round . Thus, the initial tree is and it is expanded round after round as
where is the node selected in line 15. We call the node played in round . We use to denote the point selected by HOO in the region associated with the node played in round , while denotes the received reward.
Node selection works by comparing –values and always choosing the node with the highest –value. The –value, , at node by the end of round is an estimated upper bound on the mean-payoff function at node . To define it we first need to introduce the average of the rewards received in rounds when some descendant of node was chosen (by convention, each node is a descendant of itself):
Here, denotes the set of all descendants of a node in the infinite tree,
and is the number of times a descendant of is played up to and including round , that is,
A key quantity determining is , an initial estimate of the maximum of the mean-payoff function in the region associated with node :
In the expression corresponding to the case , the first term added to the average of rewards accounts for the uncertainty arising from the randomness of the rewards that the average is based on, while the second term, , accounts for the maximum possible variation of the mean-payoff function over the region . The actual bound on the maxima used in HOO is defined recursively by
The role of is to put a tight, optimistic, high-probability upper bound on the best mean-payoff that can be achieved in the region . By assumption, . Thus, assuming that (resp., ) is a valid upper bound for region (resp., ), we see that \max\bigl{\{}B_{h+1,2i-1}(n),B_{h+1,2i}(n)\bigr{\}} must be a valid upper bound for region . Since is another valid upper bound for region , we get a tighter (less overoptimistic) upper bound by taking the minimum of these bounds.
Obviously, for leafs of the tree , one has , while close to the root one may expect that ; that is, the upper bounds close to the root are expected to be less biased than the ones associated with nodes farther away from the root.
Note that at the beginning of round , the algorithm uses to select the node to be played (since will only be available at the end of round ). It does so by following a path from the root node to an inner node with only one child or a leaf and finally considering a child of the latter; at each node of the path, the child with highest –value is chosen, till the node with infinite –value is reached.
Figure 1 illustrates the computation done by HOO in round , as well as the correspondence between the nodes of the tree constructed by the algorithm and their associated regions. Figure 2 shows trees built by running HOO for a specific environment.
At the end of round , the size of the active tree is at most , making the storage requirements of HOO linear in . In addition, the statistics and –values of all nodes in the active tree need to be updated, which thus takes time . HOO runs in time at each round , making the algorithm’s total running time up to round quadratic in . In Section 4.3 we modify HOO so that if the time horizon is known in advance, the total running time is , while the modified algorithm will be shown to enjoy essentially the same regret bound as the original version.
Main results
We start by describing and commenting on the assumptions that we need to analyze the regret of HOO. This is followed by stating the first upper bound, followed by some improvements on the basic algorithm. The section is finished by the statement of our results on the minimax optimality of HOO.
The main assumption will concern the “smoothness” of the mean-payoff function. However, somewhat unconventionally, we shall use a notion of smoothness that is built around dissimilarity functions rather than distances, allowing us to deal with function classes of highly different smoothness degrees in a unified manner. Before stating our smoothness assumptions, we define the notion of a dissimilarity function and some associated concepts.
However, it is also natural to wonder what is the class of functions for which the algorithm (given a fixed tree) can achieve non-trivial regret bounds; a similar question for regression was investigated e.g., by Yang . We shall indicate below how to construct a subset of such a class, right after stating our assumptions connecting the tree, the dissimilarity, and the environment (the mean-payoff function). Of these, Assumption A2. will be interpreted, discussed, and equivalently reformulated below into (4), a form that might be more intuitive. The form (3) stated below will turn out to be the most useful one in the proofs.
There exists such that for all integers ,
for all , there exists such that
for all .
The mean-payoff function satisfies that for all ,
We show next how a tree induces in a natural way first a dissimilarity and then a class of environments. For this, we need to assume that the tree of coverings –in addition to (1a) and (1b)– is such that the subsets and are disjoint whenever and that none of them is empty. Then, each corresponds to a unique path in the tree, which can be represented as an infinite binary sequence , where
For points with respective representations and , we let
It is not hard to see that this dissimilarity satisfies A1.. Thus, the associated class of environments is formed by those with mean-payoff functions satisfying A2. with the so-defined dissimilarity. This is a “natural class” underlying the tree for which our tree-based algorithm can achieve non-trivial regret. (However, we do not know if this is the largest such class.)
In general, Assumption A1. ensures that the regions in the tree of coverings shrink exactly at a geometric rate. The following example shows how to satisfy A1. when the domain is a –dimensional hyper-rectangle and the dissimilarity is some positive power of the Euclidean (or supremum) norm.
We now argue that Assumption A1. is satisfied. With no loss of generality we take . Then, for all integers and ,
It is now easy to see that Assumption A1. is satisfied for the indicated dissimilarity, e.g., with the choice of the parameters and \nu_{1}=b\,\bigl{(}2\sqrt{D}\bigr{)}^{a} for HOO, and the value .
This time, Assumption A1. is satisfied, e.g., with the choice of the parameters and for HOO, and the value .
Here is another interpretation of these two facts; it will be useful when considering local assumptions in Section 4.4 (a weaker set of assumptions). First, concerning the behavior around global maxima, Assumption A2. implies that for any set with ,
denotes the set of –optimal arms. This second property essentially states that there is no sudden and large drop in the mean-payoff function around the global maxima (note that this property can be satisfied even for discontinuous functions).
Figure 3 presents an illustration of the two properties discussed above.
Before stating our main results, we provide a straightforward, though useful consequence of Assumptions A1. and A2., which should be seen as an intuitive justification for the third term in (2).
is called the suboptimality factor of node . Depending whether it is positive or not, a node is called suboptimal () or optimal ().
Under Assumptions A1. and A2., if the suboptimality factor of a region is bounded by for some , then all arms in are –optimal, that is,
Proof For all , we denote by an element of such that
By the weak Lipschitz property (Assumption A2.), it then follows that for all ,
Letting and substituting the bounds on the suboptimality and on the diameter of (Assumption A1) concludes the proof.
2 Upper bound for the regret of HOO
Auer et al. [6, Assumption 2] observed that the regret of a continuum-armed bandit algorithm should depend on how fast the volumes of the sets of –optimal arms shrink as . Here, we capture this by defining a new notion, the near-optimality dimension of the mean-payoff function. The connection between these concepts, as well as with the zooming dimension defined by Kleinberg et al. , will be further discussed in Section 5. We start by recalling the definition of packing numbers.
The following example shows that using a dissimilarity (rather than a metric, for instance) may sometimes allow for a significant reduction of the near-optimality dimension.
Note that in all these cases the –near-optimality dimension of is independent of the value of .
We can now state our first main result. The proof is presented in Section A.1.
Note that if is infinite, then the bound is vacuous. The constant in the theorem depends on and on all other parameters of HOO and of the assumptions, as well as on the bandit environment . (The value of is determined in the analysis; it is in particular proportional to .) The next section will exhibit a refined upper bound with a more explicit value of in terms of all these parameters.
The tuning of the parameters of HOO is critical for the assumptions to be satisfied, thus to achieve a good regret; given some environment, one should select the parameters of HOO such that the near-optimality dimension of the mean-payoff function is minimized. Since the mean-payoff function is unknown to the user, this might be difficult to achieve. Thus, ideally, these parameters should be selected adaptively based on the observation of some preliminary sample. For now, the investigation of this possibility is left for future work.
3 Improving the running time when the time horizon is known
A deficiency of the basic HOO algorithm is that its computational complexity scales quadratically with the number of time steps. In this section we propose a simple modification to HOO that achieves essentially the same regret as HOO and whose computational complexity scales only log-linearly with the number of time steps. The needed amount of memory is still linear. We work out the case when the time horizon, , is known in advance. The case of unknown horizon can be dealt with by resorting to the so-called doubling trick, see, e.g., [10, Section 2.3], which consists of periodically restarting the algorithm for regimes of lengths that double at each such fresh start, so that the instance of the algorithm runs for rounds.
We consider two modifications to the algorithm described in Section 3. First, the quantities of (2) are redefined by replacing the factor by , that is, now
(This results in a policy which explores the arms with a slightly increased frequency.) The definition of the –values in terms of the is unchanged. A pleasant consequence of the above modification is that the –value of a given node changes only when this node is part of a path selected by the algorithm. Thus at each round , only the nodes along the chosen path need to be updated according to the obtained reward.
However, and this is the reason for the second modification, in the basic algorithm, a path at round may be of length linear in (because the tree could have a depth linear in ). This is why we also truncate the trees at a depth of the order of . More precisely, the algorithm now selects the node to pull at round by following a path in the tree , starting from the root and choosing at each node the child with the highest –value (with the new definition above using ), and stopping either when it encounters a node which has not been expanded before or a node at depth equal to
(It is assumed that so that .) Note that since no child of a node located at depth will ever be explored, its –value at round simply equals .
We call this modified version of HOO the truncated HOO algorithm. The computational complexity of updating all –values at each round is of the order of and thus of the order of . The total computational complexity up to round is therefore of the order of , as claimed in the introduction of this section.
As the next theorem indicates this new procedure enjoys almost the same cumulative regret bound as the basic HOO algorithm.
Fix a horizon such that . Then, the regret bound of Theorem 6 still holds true at round for truncated HOO up to an additional additive factor.
4 Local assumptions
In this section we further relax the weak Lipschitz assumption and require it only to hold locally around the maxima. Doing so, we will be able to deal with an even larger class of functions and in fact we will show that the algorithm studied in this section achieves a bound on the regret regret when it is used for functions that are smooth around their maxima (e.g., equivalent to for some known smoothness degree ).
For the sake of simplicity and to derive exact constants we also state in a more explicit way the assumption on the near-optimality dimension. We then propose a simple and efficient adaptation of the HOO algorithm suited for this context.
Further, there exists such that for all and ,
There exist and such that for all ,
Assumption A3. is not a real assumption but merely a reformulation of the definition of near optimality (with the small added ingredient that the limit can be achieved, see the second step of the proof of Theorem 6 in Section A.1).
Proof Note that is such that . Therefore, for all ,
which concludes the proof of the second part of A2’..
4.2 Modified HOO algorithm
We now describe the proposed modifications to the basic HOO algorithm.
We first consider, as a building block, the algorithm called –HOO, which takes an integer as an additional parameter to those of HOO. Algorithm –HOO works as follows: it never plays any node with depth smaller or equal to and starts directly the selection of a new node at depth . To do so, it first picks the node at depth with the best –value, chooses a path and then proceeds as the basic HOO algorithm. Note in particular that the initialization of this algorithm consists (in the first rounds) in playing once each of the nodes located at depth in the tree (since by definition a node that has not been played yet has a –value equal to ). We note in passing that when , algorithm –HOO coincides with the basic HOO algorithm.
Algorithm local-HOO employs the doubling trick in conjunction with consecutive instances of –HOO. It works as follows. The integers will index different regimes. The regime starts at round and ends when the next regime starts; it thus lasts for rounds. At the beginning of regime , a fresh copy of –HOO, where , is initialized and is then used throughout the regime.
Note that each fresh start needs to pull each of the nodes located at depth at least once (the number of these nodes is ). However, since round lasts for time steps (which is exponentially larger than the number of nodes to explore), the time spent on the initialization of –HOO in any regime is greatly outnumbered by the time spent in the rest of the regime.
In the rest of this section, we propose first an upper bound on the regret of –HOO (with exact and explicit constants). This result will play a key role in proving a bound on the performance of local-HOO.
4.3 Adaptation of the regret bound
In the following we write for the smallest integer such that
and consider the algorithm –HOO, where . In particular, when is chosen, the obtained bound is the same as the one of Theorem 6, except that the constants are given in analytic forms.
then the following bound holds for the expected regret of –HOO:
The proof, which is a modification of the proof to Theorem 6, can be found in Section A.3 of the Appendix. The main complication arises because the weakened assumptions do not allow one to reason about the smoothness at an arbitrary scale; this is essentially due to the threshold used in the formulation of the assumptions. This is why in the proposed variant of HOO we discard nodes located too close to the root (at depth smaller than ). Note that in the bound the second term arises from playing in regions corresponding to the descendants of “poor” nodes located at level . In particular, this term disappears when , in which case we get a bound on the regret of HOO provided that holds.
Under the prescribed assumptions, the rate of convergence is of order no matter the ambient dimension . Although the rate is independent of , the latter impacts the performance through the multiplicative factor in front of the rate, which is exponential in . This is, however, not an artifact of our analysis, since it is natural that exploration in a –dimensional space comes at a cost exponential in . (The exploration performed by HOO naturally combines an initial global search, which is bound to be exponential in , and a local optimization, whose regret is of the order of .)
The following theorem is an almost straightforward consequence of Theorem 9 (the detailed proof can be found in Section A.4 of the Appendix). Note that local-HOO does not require the knowledge of the parameter in A2’..
5 Minimax optimality in metric spaces
In this section we provide two theorems showing the minimax optimality of HOO in metric spaces. The notion of packing dimension is key.
The reader interested in the explicit expressions of and is referred to the last lines of the proof of the theorem in the Appendix.
Discussion
In this section we would like to shed some light on the results of the previous sections. In particular we generalize the situation of Example 5, discuss the regret that we can obtain, and compare it with what could be obtained by previous works.
We equip with a norm . We assume that the mean-payoff function has a finite number of global maxima and that it is locally equivalent to the function –with degree – around each such global maximum of ; that is,
This means that there exist such that for all satisfying ,
Smoothness overestimated: Now, if the true smoothness is overestimated by choosing or and , then the assumption of weak Lipschitzness is violated and we are unable to provide any guarantee on the behavior of HOO. The latter, when used with an overestimated smoothness parameter, may lack exploration and exploit too heavily from the beginning. As a consequence, it may get stuck in some local optimum of , missing the global one(s) for a very long time (possibly indefinitely). Such a behavior is illustrated in the example provided in and showing the possible problematic behavior of the closely related algorithm UCT of . UCT is an example of an algorithm overestimating the smoothness of the function; this is because the –values of UCT are defined similarly to the ones of the HOO algorithm but without the third term in the definition (2) of the –values. This corresponds to an assumed infinite degree of smoothness (that is, to a locally constant mean-payoff function).
2 Relation to previous works
Several works [3; 21; 12; 6; 22] have considered continuum-armed bandits in Euclidean or, more generally, normed or metric spaces and provided upper and lower bounds on the regret for given classes of environments.
Kleinberg considered mean-payoff functions on the real line that are Hölder continuous with degree . The derived regret bound is \Theta\bigl{(}n^{(\alpha+1)/(\alpha+2)}\bigr{)}.
Auer et al. extended the analysis to classes of functions that are equivalent to around their maxima , where the allowed smoothness degree is also larger: . They derived the regret bound
where the parameter is such that the Lebesgue measure of –optimal arm is .
For an illustration, consider again the example of Section 5.1. The result of Auer et al. shows that for , the regret is (since here , with the notation above). Our result extends the rate of the regret bound to any dimension .
On the other hand the analysis of Kleinberg et al. does not apply because in this example is controlled only when is close in some sense to (i.e., when ), while (9) requires such a control over the whole set . However, note that the local weak-Lipschitz assumption A2’. requires an extra condition in the vicinity of compared to (9) as it is based on the notion of weak Lipschitzness. Thus, A2’. and (9) are in general incomparable (both capture a different phenomenon at the maxima).
We now compare our results to those of and under Assumption A2. (which does not cover the example of Section 5.1 unless is large). Under this assumption, our algorithms enjoy essentially the same theoretical guarantees as the zooming algorithm of [22; 23]. Further, the following hold.
Our algorithms do not require the oracle needed by the zooming algorithm.
Our truncated HOO algorithm achieves a computational complexity of order , whereas the complexity of a naive implementation of the zooming algorithm is likely to be much larger.The zooming algorithm requires a covering oracle that is able to return a point which is not covered by the set of active strategies, if there exists one. Thus a straightforward implementation of this covering oracle might be computationally expensive in (general) continuous spaces and would require a ‘global’ search over the whole space.
Both truncated HOO and the zooming algorithms use the doubling trick. The basic HOO algorithm, however, avoids the doubling trick, while meeting the computational complexity of the zooming algorithm.
The fact that the doubling trick can be avoided is good news since an algorithm that uses the doubling trick must start from tabula rasa time to time, which results in predictable, yet inevitable, sharp performance drops –a quite unpleasant property. In particular, for this reason algorithms that rely on the doubling trick are often neglected by practitioners. In addition, the fact that we avoid the oracle needed by the zooming algorithm is attractive as this oracle might be difficult to implement for general (non-metric) dissimilarities.
Acknowledgements
We thank one of the anonymous referee for his valuable comments, which helped us to provide a fair and detailed comparison of our work to prior contributions.
This work was supported in part by French National Research Agency (ANR, project EXPLO-RA, ANR-08-COSI-004), the Alberta Ingenuity Centre of Machine Learning, Alberta Innovates Technology Futures (formerly iCore and AIF), NSERC and the PASCAL2 Network of Excellence under EC grant no. 216886.
Appendix A Proofs
We begin with three lemmas. The proofs of Lemmas 15 and 16 rely on concentration-of-measure techniques, while the one of Lemma 14 follows from a simple case study. Let us fix some path , , , of optimal nodes, starting from the root. That is, denoting , we mean that for all , the suboptimality of equals and is a child of .
Let be a suboptimal node. Let be the largest depth such that is on the path from the root to . Then for all integers , we have
Proof Consider a given round . If , then this is because the child of on the path to had a better –value than its brother . Since by definition, –values can only increase on a chosen path, this entails that . This is turns implies, again by definition of the –values, that . Thus,
But, once again by definition of –values,
and the argument can be iterated. Since up to round no more than nodes have been played (including the suboptimal node ), we know that has not been played so far and thus has a –value equal to . (Some of the previous optimal nodes could also have had an infinite –value, if not played so far.) We thus have proved the inclusion
Now, for any integer it holds that
where we used for the inequality the fact that the quantities are constant from to , except when , in which case, they increase by 1; therefore, on the one hand, at most of the can be smaller than and on the other hand, can only happen if . Using (11) and then taking expectations yields the result.
Let Assumptions A1. and A2. hold. Then, for all optimal nodes and for all integers ,
Proof On the event that was not played during the first rounds, one has, by convention, . In the sequel, we therefore restrict our attention to the event \bigl{\{}T_{h,i}(n)\geqslant 1\bigr{\}}.
Lemma 3 with ensures that for all arms . Hence,
We take care of the last term with a union bound and the Hoeffding-Azuma inequality for martingale differences.
To do this in a rigorous manner, we need to define a sequence of (random) stopping times when arms in were pulled:
Note that , hence it holds that . We denote by the arm pulled in the region corresponding to . Its associated corresponding reward equals and
where we used a union bound to get the last inequality.
is a martingale w.r.t. the filtration \mathcal{G}_{t}=\sigma\bigl{(}\widetilde{X}_{1},Z_{1},\ldots,\widetilde{X}_{t},Z_{t},\widetilde{X}_{t+1}\bigr{)}. This follows, via optional skipping (see [14, Chapter VII, adaptation of Theorem 2.3]), from the facts that
is a martingale w.r.t. the filtration and that the events are –measurable for all .
Applying the Hoeffding-Azuma inequality for martingale differences (see ), using the boundedness of the ranges of the induced martingale difference sequence, we then get, for each ,
For all integers , for all suboptimal nodes such that , and for all integers such that
Proof The mentioned in the statement of the lemma are such that
Now it follows from the same arguments as in the proof of Lemma 15 (optional skipping, the Hoeffding-Azuma inequality, and a union bound) that
where we used the stated bound on to obtain the last inequality.
Combining the results of Lemmas 14, 15, and 16 leads to the following key result bounding the expected number of visits to descendants of a “poor” node.
Under Assumptions A1. and A2., for all suboptimal nodes with , we have, for all ,
Proof We take as the upper integer part of and use union bounds to get from Lemma 14 the bound
Lemmas 15 and 16 further bound the quantity of interest as
Proof (of Theorem 6) First, let us fix . The statement will be proven in four steps.
First step. For all , denote by the set of those nodes at depth that are –optimal, i.e., the nodes such that . (Of course, .) Then, let be the union of these sets when varies. Further, let be the set of nodes that are not in but whose parent is in . Finally, for we denote by the nodes in that are located at depth in the tree (i.e., whose parent is in ).
Lemma 17 bounds in particular the expected number of times each node is visited. Since for these nodes , we get
Second step. We bound the cardinality of . We start with the case . By definition, when , one has , so that by Lemma 3 the inclusion holds. Since by Assumption A1., the sets contain disjoint balls of radius , we have that
We prove below that there exists a constant such that for all ,
Thus we obtain the bound |\mathcal{I}_{h}|\leqslant C\,\bigl{(}\nu_{2}\rho^{h}\bigr{)}^{-d^{\prime}} for all . We note that the obtained bound |\mathcal{I}_{h}|\leqslant C\,\bigl{(}\nu_{2}\rho^{h}\bigr{)}^{-d^{\prime}} is still valid for , since .
It only remains to prove (13). Since , where is the near-optimality of , we have, by definition, that
and thus, there exists such that for all ,
which in turn implies that for all ,
The result is proved with if . Now, consider the case . Given the definition of packing numbers, it is straightforward that for all \varepsilon\in\bigl{[}\varepsilon_{d^{\prime}},\,\nu_{2}\bigr{]},
therefore, for all \varepsilon\in\bigl{[}\varepsilon_{d^{\prime}},\,\nu_{2}\bigr{]},
for the choice C=\max\bigl{\{}1,\,\,u_{d^{\prime}}\,\nu_{2}^{d^{\prime}}\bigr{\}}. Because we take the maximum with 1, the stated inequality also holds for , which concludes the proof of (13).
Third step. Let be an integer to be chosen later. We partition the nodes of the infinite tree into three subsets, , as follows. Let the set contain the descendants of the nodes in (by convention, a node is considered its own descendant, hence the nodes of are included in ); let ; and let contain the descendants of the nodes in . Thus, and are potentially infinite, while is finite.
We recall that we denote by the node that was chosen by HOO in round . From the definition of the algorithm, each node is played at most once, thus no two such random variables are equal when varies. We decompose the regret according to which of the sets the nodes belong to:
The contribution from is easy to bound. By definition any node in is –optimal. Hence, by Lemma 3, the corresponding domain is included in . By definition of a tree of coverings, the domains of the descendants of these nodes are still included in . Therefore,
For , consider a node . It belongs to and is therefore –optimal. By Lemma 3, the corresponding domain is included in . By the result of the second step of this proof and using that each node is played at most once, one gets
We finish by bounding the contribution from . We first remark that since the parent of any element is in , by Lemma 3 again, we have that . We now use the first step of this proof to get
Now, it follows from the fact that the parent of is in that when . Substituting this and the bound on obtained in the second step of this proof, we get
Fourth step. Putting the obtained bounds together, we get
(recall that ). Note that all constants hidden in the symbol only depend on , , and .
Now, by choosing such that is of the order of , that is, is of the order of , we get the desired result, namely,
A.2 Proof of Theorem 8 (regret bound for truncated HOO)
The proof follows from an adaptation of the proof of Theorem 6 and of its associated lemmas; for the sake of clarity and precision, we explicitly state the adaptations of the latter.
Adaptations of the lemmas. Remember that denotes the maximum depth of the tree, given horizon . The adaptation of Lemma 14 is done as follows. Let be a suboptimal node with and let be the largest depth such that is on the path from the root to . Then, for all integers , one has
As for Lemma 15, its straightforward adaptation states that under Assumptions A1. and A2., for all optimal nodes with and for all integers ,
Similarly, the same changes yield from Lemma 16 the following result for truncated HOO. For all integers , for all suboptimal nodes such that and , and for all integers such that
Combining these three results (using the same methodology as in the proof of Lemma 17) shows that under Assumptions A1. and A2., for all suboptimal nodes such that and , one has
(We thus even improve slightly the bound of Lemma 17.)
Adaptation of the proof of Theorem 6. The main change here comes from the fact that trees are cut at the depth . As a consequence, the sets , , , and are defined only by referring to nodes of depth smaller than . All steps of the proof can then be repeated, except the third step; there, while the bounds on the regret resulting from nodes of and go through without any changes (as these sets were constructed by considering all descendants of some base nodes), the bound on the regret associated with the nodes calls for a modified proof since at this stage we used the property that each node is played at most once. But this is not true anymore for nodes located at depth , which can be played several times. Therefore the proof is modified as follows.
Consider a node at depth . Then, by definition of ,
Since the considered nodes are –optimal, the corresponding domains are –optimal by Lemma 3, thus also –optimal. The instantaneous regret incurred when playing any of these nodes is therefore bounded by ; and the associated cumulative regret (over rounds) can be bounded by . In conclusion, with the notations of Theorem 6, we get the new bound
The rest of the proof goes through and only this additional additive factor of is suffered in the final regret bound. (The additional factor can be included in the notation.)
A.3 Proof of Theorem 9 (regret bound for z𝑧z–HOO)
We start with the following equivalent of Lemma 3 in this new local context. Remember that is the smallest integer such that
Under Assumptions A1. and A2’., for all , if the suboptimality factor of a region is bounded by for some , then all arms in are –optimal, that is,
When , i.e., the node is optimal, the bound improves to
Since and , the second part of Assumption A2’. then yields
It follows from the definition of that f^{*}-f(x)+\varepsilon=\max\bigl{\{}f^{*}-f(x),\,\nu_{1}\rho^{h}\bigr{\}}, and this implies
But , i.e., , we thus have proved
In conclusion, for all sufficiently small . Letting concludes the proof.
that is, .
We now provide an adaptation of Lemma 17 (actually based on adaptations of Lemmas 14 and 15), providing the same bound under local conditions that relax the assumptions of Lemma 17 to some extent.
Consider a depth . Under Assumptions A1. and A2.’, the algorithm –HOO satisfies that for all and all suboptimal nodes with and ,
Proof We consider some path of optimal nodes, starting at depth . We distinguish two cases, depending on whether there exists such that or not.
In the first case, we denote the largest such . The argument of Lemma 14 can be used without any change and shows that for all integers ,
In the second case, we denote by the ancestor of located at depth . By definition of –HOO, at some round only if and since –values can only increase on a chosen path, can only happen if . Repeating again the argument of Lemma 14, we get that for all integers ,
Now, notice that Lemma 16 is valid without any assumption. On the other hand, with the modified assumptions, Lemma 15 is still true but only for optimal nodes with . Indeed, the only point in its proof where the assumptions were used was in the fourth line, when applying Lemma 3; here, Lemma 18 with provides the needed guarantee.
The proof is concluded with the same computations as in the proof of Lemma 17.
Proof (of Theorem 9) We follow the four steps in the proof of Theorem 6 with some slight adjustments. In particular, for , we use the sets of nodes and defined therein.
First step. Lemma 19 bounds the expected number of times each node is visited. Since for these nodes , we get
Second step. We bound here the cardinality . By Lemma 18 with , when and , one has .
Now, by Assumption A1. and by using the same argument as in the second step of the proof of Theorem 6,
Assumption A3. can be applied since and yields the inequality |\mathcal{I}_{h}|\leqslant C\bigl{(}\nu_{2}\rho^{h}\bigr{)}^{-d}.
Third step. We consider some integer to be defined by the analysis in the fourth step. We define a partition of the nodes located at a depth equal to or larger than ; more precisely,
contains the nodes of and their descendants,
,
contains the nodes and their descendants,
is formed by the nodes located at depth not belonging to , i.e., such that , and their descendants.
As in the proof of Theorem 6 we denote by the regret resulting from the selection of nodes in , for .
As for , we also use here that nodes in belong to some , with ; in particular, they are the child of some element of and as such, firstly, they are –optimal (by Lemma 18) and secondly, their number is bounded by |\mathcal{J}_{h}|\leqslant 2|\mathcal{I}_{h-1}|\leqslant 2C\bigl{(}\nu_{2}\rho^{h-1}\bigr{)}^{-d}. Thus,
where we used the bound of Lemma 19. Finally, for , we use that it contains at most nodes, each of them being associated with a regret controlled by Lemma 19; therefore,
Fourth step. Putting things together, we have proved that
where (using that in the second inequality)
It remains to define the parameter . In particular, we propose to choose it such that the terms
are balanced. To this end, let be the smallest integer such that ; in particular,
Note from the inequality that this is such that
and thus this satisfies in view of the assumption of the theorem indicating that is large enough. The final bound on the regret is then
A.4 Proof of Theorem 10 (regret bound for local-HOO)
Proof We use the notation of the proof of Theorem 9. Let be a positive integer such that for , one has
we can therefore apply the result of Theorem 9 in regimes indexed by . For previous regimes, we simply upper bound the regret by the number of rounds, that is, . For round , we denote by the index of the regime where lies in (regime ). Since regime terminates at round , we have
where denote some constants depending only on the parameters but not on . Note that for the last equality we used that the first term in the sum of the two terms that depend on dominates the second term.
A.5 Proof of Theorem 12 (uniform upper bound on the regret of HOO against the class of all weak Lipschitz environments)
A.6 Proof of Theorem 13 (minimax lower bound in metric spaces)
The set contains “difficult” environments for the –armed bandit problem.
For any strategy suited to the –armed bandit problem, one can construct a strategy for the –armed bandit problem such that
Proof We only deal with the case of deterministic strategies. The extension to randomized strategies can be done using Fubini’s theorem (by integrating also w.r.t. the auxiliary randomizations used).
First step. Let be a real number and be an integer, both to be defined during the course of the analysis. The set only contains elements, denoted by and given by product distributions. For , the distribution is obtained as the product of the when and where
One can extract the following result from the proof of the lower bound of [10, Section 6.9].
For all strategies for the –armed bandit (where ), one has
Second step. We now need to construct such that item (ii) is satisfied. We assume that is such that contains disjoint balls with radius . (We shall quantify later in this proof a suitable value of .) Denoting by the corresponding centers, these disjoint balls are then .
With each of these balls we now associate a bandit environment over , in the following way. For all , we introduce a mapping on defined by
for all . This mapping is used to define an environment over , as follows. For all ,
Let be the corresponding mean-payoff function; its values equal
Third step. We describe how to associate with each (deterministic) strategy on a (random) strategy on the finite set of arms . Each of these strategies is indeed given by a sequence of mappings,
where for , the mappings and should only depend on the past up to the beginning of round . Since the strategy is deterministic, the mapping takes only into account the past rewards and is therefore a mapping . (In particular, equals a constant.)
We use the notations and for, respectively, the arms pulled and the rewards obtained by the strategy at each round . The arms are drawn at random according to the distributions
which we now define. (Actually, they will depend on the obtained payoffs only.) To do that, we need yet another mapping that links elements in to probability distributions over . Denoting by the Dirac probability on , the mapping is defined as
for all . Note that this definition is legitimate because the balls are disjoint when varies between and .
Finally, is defined as follows. For all ,
Before we proceed, we study the distribution of the reward obtained under (for ) by the choice of a random arm drawn according to , for some . Since can only take the values 0 or 1, its distribution is a Bernoulli distribution whose parameter we compute now. The computation is based on the fact that under , the Bernoulli distribution corresponding to arm has as an expectation, except if , in which case it is . Thus, for all ,
That is, on .
Fourth step. We now prove that the distributions of the regrets of under and of under are equal for all . On the one hand, the expectations of rewards associated with the best arms equal under the two environments. On the other hand, one can prove by induction that the sequences and have the same distribution. (In the argument below, conditioning by empty sequences means no conditioning. This will be the case only for .)
Under and given , the distribution of is obtained by definition as the two-step random draw of and then, conditionally on this first draw, . By the above results, the distribution of is thus a Bernoulli distribution with parameter .
At the same time, under and given , the choice of
yields a reward distributed according to , that is, by definition and with the notations above, a Bernoulli distribution with parameter .
The argument is concluded by induction and by using the fact that rewards are drawn independently in each round.
Fifth step. We summarize what we proved so far. For , provided that there exist disjoint balls in , we could construct, for all strategies for the –armed bandit problem, a strategy for the –armed bandit problem such that, for all and all ,
But by the assumption on the packing dimension, there exists such that for all , the choice of guarantees the existence of such disjoint balls. Substituting this value, and using the results of the first and fourth steps of the proof, we get
the left-hand side is smaller than the maximal regret w.r.t. all weak Lipschitz environments;
the right-hand side can be lower bounded and then optimized over in the following way.
By definition of and the fact that it is larger than 2, one has
where C=\sqrt{\bigl{(}4\ln(4/3)\bigr{)}\,\big{/}\,c}. We can optimize the final lower bound over .
To that end, we choose, for instance, such that , that is,
To ensure that this choice of is valid we need to show that . Since the latter requirement is equivalent to
it suffices to choose the right-hand side to be ; we then get that indeed holds for all , thus concluding the proof of the theorem.