b

DiscoverSearch
About
My stuff
K-NN active learning under local smoothness assumption
2020·arXiv
Abstract
Abstract

There is a large body of work on convergence rates either in passive or in active learning. Here we first outline some of the main results that have been obtained, more specifically in a nonparametric setting under assumptions about the smoothness of the regression function (or the boundary between classes) and the margin noise. We discuss the relative merits of these underlying assumptions by putting active learning in perspective with recent work on passive learning. In particular, the regression function is assumed to be smooth with respect to the marginal probability of the data. This is more general than what was previously used in the literature. Our contribution is an active nearest neighbor learning algorithm, that is able to deal with this setting, and that outperforms its passive counterpart. Our algorithm works for a larger class of probability distributions than those previously used, especially for distributions of probability for which the density function is not necessarily bounded below, and also for discrete distributions.

Keywords Nonparametric learning  ·active learning  ·nearest-neighbors  ·smoothness condition.

Active learning is a machine learning approach for reducing the data labeling effort. Given an instance space X or a pool of unlabeled data  {X1, . . . , Xw}provided by a distribution  PX, the learner focuses its labeling effort only on the most “informative” points so that a model built from them can achieve the best possible guarantees (Dasgupta, 2011). Such guarantees are particularly interesting when they are significantly better than those obtained in passive learning (Hanneke and Yang, 2015). In the context of this work, we consider binary classification (where the label Y of X takes its value in {0, 1}) in a nonparametric setting. Extensions to multiclass classification and adaptive algorithms are discussed at the end of this paper (Section 6).

The nonparametric setting has the advantage of providing guarantees with many informations such as the dependence on the dimensional and distributional parameters by using some hypotheses on the regularity of the decision boundary (Castro and Nowak, 2008), on the regression function (Minsker, 2012; Locatelli et al., 2017), and on the geometry of instance space (called strong density assumption) (Audibert and Tsybakov, 2007; Locatelli et al., 2017; Minsker, 2012). One of the initial works on nonparametric active learning (Castro and Nowak, 2008) assumed that the decision boundary is the graph of a smooth function, that a margin assumption very similar to Tsybakov’s noise assumption (Mammen and Tsybakov, 1999) holds, and that distribution  PXis uniform. This led to a better guarantee than in passive learning. Instead of the assumption on the decision boundary, other works (Minsker, 2012; Locatelli et al., 2017) supposed rather that the regression function is smooth (in some sense). This assumption, along with Tsybakov’s noise assumption and the strong density assumption also gave a better guarantee than in passive learning. Moreover, unlike in (Castro and Nowak, 2008), they provided algorithms that are adaptive with respect to the margin’s noise and to the smoothness parameters.

However, recent work (Chaudhuri and Dasgupta, 2014) pointed out some disadvantages of the preceding smoothness assumption, and extended it in the context of passive learning with k-nearest neighbors (k-NN) by using a more general smoothness assumption that is able to sharply characterize the rate of convergence for all probability distributions that satisfy it.

In this paper, we thus extend the work of (Chaudhuri and Dasgupta, 2014) to the active learning setting, and provide a novel algorithm that outputs a clas-sifier with the same rate of convergence as other recent algorithms with more restrictive hypotheses, as for example (Minsker, 2012; Locatelli et al., 2017). Section 2 introduces general definitions, Section 3 presents previous work on convergence rates in active and passive non-parametric learning, with a special emphasis on the assumptions related to our work. Section 4 provides an outline of our algorithm while Section 5 describes its theoretical motivations and Section 6 contains the conclusion and some perspectives for future work.

We begin with some general definitions and notations about active learning in binary classification, then recall the concept of k-NN classifiers. Finally, the main assumptions that are used in nonparametric active learning are explained.

2.1 Active learning setting

Let (X, ρ) be a metric space. In this paper we set  X ⊂ Rdand refer to it as the instance space, and take  ρas the Euclidean metric. Let Y = {0, 1} the label space. We assume that the pairs (X, Y ) are random variables distributed according to an unknown probability P over  X × Y. Let us denote  PXthe marginal distribution of P over X.

Given  w ∈ Nand an i.i.d. sample (X1, Y1), . . . , (Xw, Yw) drawn according to probability P, the learning problem consists in minimizing the risk R(f) = P(Y ̸= f(X)) over all measurable functions, called classifiers,  f : X → Y.

In active learning, the labels are not available from the beginning but we can request iteratively at a certain cost (to a so-called oracle) a given number n of samples, called the budget (n ≤ w). In passive learning, all labels are available and n = w. At any time, we choose to request the label of a point X according to the previous observations. The point X is chosen to be most “informative”, which amounts to belonging to a region where classification is difficult and requires more labeled data to be collected. Therefore, the goal of active learning is to design a sampling strategy that outputs a classifier �fn,wwhose excess risk (see below) is as small as possible with high probability over the requested samples, as reviewed in (Dasgupta, 2011; Hanneke and Yang, 2015; Dasgupta, 2017).

Given x in X, let us introduce the regression function  η(x) = E(Y |X = x) = P(Y = 1| X = x). It is easy to show (Lugosi, 2002) that the function  f ∗(x) = 1η(x)≥1/2achieves the minimum risk and that  R(f ∗) =  EX(min(η(X), 1 −η(X))). Because P is unknown, the function  f ∗is unreachable and thus the aim of a learning algorithm is to return a classifier �fn,wwith minimum excess risk R( �fn,w)−R(f ∗) with high probability over the sample (X1, Y1), . . . , (Xw, Yw).

2.2 k-Nearest Neighbors (k-NN) classifier

Given two integers k, n such that k < n, and a test point  X ∈ X, the k-NN classifier predicts the label of X by giving the majority vote of its k nearest neighbors amongst the sample  X1, . . . , Xn. For k = 1, the k-NN classifier returns the label of the nearest neighbor of X amongst the sample  X1, . . . , Xn. If k is allowed to grow with n, the method is called  kn-NN. For a complete discussion of nearest neighbors classification, see for example (Biau and Devroye, 2015; Shalev-Shwartz and Ben-David, 2014; Chaudhuri and Dasgupta, 2014).

2.3 Regularity, noise and strong density assumptions

Let B(x, r) =  {x′ ∈ X, ρ(x, x′) < r}and ¯B(x, r) =  {x′ ∈ X, ρ(x, x′) ≤ r}the open and closed balls (with respect to the Euclidean metric  ρ), respectively, centered at  x ∈ Xwith radius r > 0. Let supp(PX) =  {x ∈ X, ∀r >0, PX(B(x, r)) > 0}the support of the marginal distribution  PX.

image

Let  η : X →[0, 1] be the regression function defined as  η(x) = P(Y = 1|X = x). We say that  ηis (α, L)-H¨older continuous(0  < α ≤ 1,and  L ≥1) if  ∀x, x′ ∈ X,

image

The notion of H¨older continuity ensures that the proximity between two closest (according to the metric  ρ) points is reflected in a similar value of the conditional probability  η.

This definition remains true for a general metric space, but when  ρis the Euclidean metric, we should always have 0  < α ≤1, otherwise  ηbecomes constant(Pugh, 2002).

In most of the previous works (for example (Audibert and Tsybakov, 2007; Minsker, 2012; Ga¨ıffas, 2007)), the definition H1a is used along with the following notion (H1b) for technical reasons.

image

Let P be the probability distribution defined over  X × Yand  PXthe marginal distribution of P over X. We say that P satisfies the strong density assumption if there exists some constants  r0 >0,  c0 >0,  pmin >0 such that for all x ∈supp(PX):

image

where  pXis the density function of the marginal distribution  PXand  λis the Lebesgue measure.

The strong density assumption ensures that, given a realisation X = x according to  PX, there exists an infinite number of realisations  X1 = x1, . . . , Xm= xm, . . .in a neighborhood of x.

Sometimes, the notion of strong density is used to geometrically characterize the set where the classification is difficult (Locatelli et al., 2017), and then combined with the following definition of Margin noise, allows to nicely control the error of classification beyond a given number of label requests.

image

The probability distribution P satisfies the margin noise assumption (sometimes called  Tsybakov’s noiseassumption (Audibert and Tsybakov, 2007)) with parameter  β ≥0 if for all 0  < ǫ ≤1, there is  C = C(β) ∈[1, +∞[ such that

image

The margin noise assumption gives a bound on the probability that the label of the points in the neigborhood of a point x differs from the label of x given by the conditional probability  η(x). It also describes the behavior of the regression function in the vicinity of the decision boundary  η(x) = 12. When  βgoes to infinity, we observe a “jump” of  ηaround the decision boundary, and then we obtain Massart’s noise condition (Massart and N´ed´elec, 2006). Small values of βallow for  ηto “cuddle” 12when we approach the decision boundary.

image

Let 0  < α ≤1 and L > 1. The regression function is (α, L)-smoothif for all x, z ∈supp(PX) we have:

image

where d is the dimension of the instance space.

Equivalently, (H3) can be rewritten as:

image

It is important to note that the (α, L)-smooth assumption (H3) is more general than the H¨older continuity assumption (H1a), as stated in Theorem 1 below.

image

Suppose that  X ⊂ Rd, that the regression function  ηis (αh, Lh)-H¨older continuous, and that  PXsatisfies H1b. Then there is a constant L > 1 such that for any  x, z ∈supp(PX), we have:

image

The marginal distribution  PXis a doubling-probability if there exists a constant  Cdb >0 such that for any  x ∈ X, and r > 0, we have:

image

This notion was initially introduced for geometric purposes in the setting of measure theory (Heinonen, 2012; Federer, 2014). It helps for constructing a subcover of a metric space by also minimizing the overlap between the elements of the subcover. Doubling-probability has been used in a machine learning context, particularly k-NN classification (or regression), where the constant  Cdbis interpreted as the intrinsic dimension of the region where the data belong (Kpotufe, 2011). This allows to reduce considerably the complexity of the clas-sification problem and to bypass the so-called curse of dimension. Also, it is also proved (Kpotufe, 2011) that the notion of doubling-probability generalizes the strong density assumption H1b. It is thus more universal, and does not require a probability density.

In this paper, doubling-probability is used only for geometrical purposes. It is later relaxed, so that it becomes sufficient to consider only balls B(x, r) with  PX(B(x, r)) sufficiently large to satisfy the doubling-probability condition (H4).

3.1 Previous work

Active learning theory has been mostly studied during the last decades in a parametric setting, see for example (Balcan et al., 2010; Hanneke, 2011; Dasgupta, 2017) and references therein. One of the pioneering works studying the achievable limits in active learning in a nonparametric setting (Castro and Nowak, 2008) required that the decision boundary is the graph of a H¨older continuous function with parameter  α(H1a). Using a notion of margin noise (with parameter  β) very similar to (H2), the following minimax rate was obtained:

image

where  γ= d−1αand d is the dimension of instance space (X = Rd).

Note that this result assumes the knowledge of the smoothness and margin noise parameters, whereas an algorithm that achieves the same rate, but that adapts to these parameters was proposed recently in (Locatelli et al., 2018).

In passive learning, by assuming that the regression function is H¨older continuous (H1a), along with (H1b) and (H2), the following minimax rate was established (Audibert and Tsybakov, 2007):

image

In active learning, using the same assumptions (H1a), (H1b) and (H2), with the additional condition  αβ < d, the following minimax rate was obtained (Locatelli et al., 2017)

image

where ˜O indicates that there may be additional logarithmic factors. This active learning rate given by (3) thus represents an improvement over the passive learning rate (2) that uses the same hypotheses.

With another assumption on the regression function relating the  L2and  L∞approximation losses of certain piecewise constant or polynomial approximations of  ηin the vicinity of the decision boundary, the same rate (3) was also obtained (Minsker, 2012).

3.2 Link with k-NN classifiers

For practicals applications, an interesting question is whether k-NN classifiers attain the rate given by (2) in passive learning and by (3) in active learning. In passive learning, under assumptions (H1a), (H1b) and (H2), and for suitable  kn, it was shown in (Chaudhuri and Dasgupta, 2014) that  kn-NN indeed achieves the rate (2).

In active learning a pool-based algorithm that outputs a k-NN classifier has been proposed in (Kontorovich et al., 2016), but its assumptions differ from ours in terms of smoothness and noise. Similarly, the algorithm proposed in (Hanneke) outputs a 1-NN classifier based on a subsample of a given pool of data, such that the label of each instance of this subsample is determined with high probability by the labels of its neighbors within the pool. The number of neighbors is adaptively chosen for each instance in the subsample, leading to the minimax rate (3) under the same assumptions as in (Locatelli et al., 2017).

To obtain more general results on the rate of convergence for k-NN classifiers in metric spaces under minimal assumptions, the more general smoothness assumption given by (H3) was used in (Chaudhuri and Dasgupta, 2014). By using a k-NN algorithm, and under assumptions (H2) and (H3), the rate of convergence obtained in (Chaudhuri and Dasgupta, 2014) is also of the order of (2). Additionally, using assumption (H3) instead of (H1a) removes the need for the strong density assumption (H1b), which therefore allows for more probability classes.

3.3 Contributions of the current work

In this work, we will use the assumptions that were used in the context of passive learning in (Chaudhuri and Dasgupta, 2014), and show that is is possible to use them in active learning as well.

For the sake of clarity, let us restate here these assumptions that will be used throughout this paper. We assume that the assumptions (H3), (H2), (H4) simultaneously hold respectively with parameters (α, L), (β, C), Cdb.

In this paper, we provide an active learning algorithm under assumptions (H3), (H2) that were used in passive learning in (Chaudhuri and Dasgupta, 2014). We additionally assume that the underlying marginal probability  PXsatisfies (H4) mostly for geometrical convenience. Our algorithm has several advantages:

image

According to the assumption (H2), as we will see, our algorithm also (as in (Minsker, 2012)) covers the most interesting case where the regression function is allowed to cross the boundary decision  {x, η(x) = 1

image

As explained in Section 2.1, we consider an active learning setting with a pool of i.i.d. unlabeled examples  K = {X1, X2, . . . , Xw}. Let  n ≤ wthe budget, that is the maximum number of points whose label we are allowed to query to the oracle. The objective of the algorithm is to build a 1-NN classifier, based on a labelled set  Sacof carefully chosen points. This set contains a subset of most informative points in K and is called the active set. More precisely, a point Xtis considered informative if its label cannot be inferred (see below) from the previous observations  Xt′(with  t′ < t). The set  Sacstarts with  Xt1 = X1chosen arbitrarily in K and stops when the budget n is reached or when  Xwis attained.

When a point  Xtis informative, instead of requesting directly its label to the (noisy) oracle, we infer it by requesting the labels of its nearest neighbors in K, as was done in (Hanneke). This is reasonable for practical situations where the uncertainty about the label of  Xthas to be overcome, and it is related to the assumption (H3). Note that it differs from the setting of (Locatelli et al., 2018), where the label of  Xtis requested several times. The number of neighbors  ktused for inferring that label of  Xtis determined such that, while respecting the budget, we can predict with high confidence the true label as  f ∗(Xt) of Xtby the empirical mean of the labels of its  ktnearest neighbors.

The labelled active set  Sacoutput by the algorithm will comprise only the informative points on which we have sufficient guarantees when considering the inferred label as the right label. Finally, we show that the labelled active set Sacis sufficient to predict the label of any new point with a 1-NN classification rule �fn,w.

4.2 Algorithm

The KALLS algorithm (Algorithm 1) aims at determining the active set defined in Section 4.1 and the related 1-NN classifier �fn,wunder the assumption (H3) and (H2).

Before beginning the description of KALLS, let us introduce some variables and notations, whose precise form will be justified in Section 5. The latter contains the proof sketch of the convergence of KALLS, while the complete proofs are in Appendix A.

For  ǫ, δ ∈(0, 1),  k ≥1, set:

image

where

image

Let

image

For  Xs ∈ K = {X1, . . . , Xw}, we denote henceforth by  X(k)sits k-th nearest neighbor in K, and  Y (k)sthe corresponding label. For an integer  k ≥1, let

image

The inputs of KALLS are a pool K of unlabelled data of size w, the budget n, the smoothness parameters (α, L) from (H3), the margin noise parameters (β, C) from (H2), a confidence parameter  δ ∈(0, 1) and an accuracy parameter ǫ ∈(0, 1). For the moment, these parameters are fixed from the beginning but adaptive algorithms such as (Locatelli et al., 2017) could be exploited, in particular for the  αand  βparameters.

At any given stage, the current version of the labelled active set  Sacis denoted by �S. Based on  Sac, with high confidence, the 1-NN classifier �fn,wagrees with the Bayes classifier at points that lie beyond some margin  ∆o >0 of the decision boundary. Formally, given  x ∈ Xsuch that  |η(x) − 1/2| > ∆0, we have �fn,w(x) =  1η(x)≥1/2with high confidence. We will show in Section 5 that, with a suitable choice of  ∆o, the assumption (H2) leads to the desired rate of convergence (3).

KALLS uses two main subroutines : Reliable and ConfidentLabel, which are detailed below in Sections 4.3 and 4.4, respectively.

4.3 Reliable subroutine

The Reliable subroutine is a binary test that checks if the label of a current point X can be inferred with high confidence from some previously informative points before reaching X. These points are obtained via a set �Scalled current active set. Each element of �Scan be seen as a triplet (X′, �Y ′, c) where  X′is an informative point, �Y ′its inferred label, and c > 0 can be thought as

image

a guarantee for predicting the right label Y of  X′as �Y ′. Formally, we have O(c) ≤ |η(X′) − 12|when (X′, �Y ′, c) ∈ �Sand  X′is relatively far from the decision boundary. If  Reliable(X, δ, α, L, �S) outputs True, the point X is not considered to be informative, and �Swill not be updated. By convention, Reliable(X, δ, α, L, ∅) always returns False.

The inputs are the current point X, a confidence parameter  δ, the smoothness parameters (α, L) from (H3), and the set �Sbefore examining the point X.

If  |η(X) − 12|entails the same confidence lower bound O(c) as that of some previous informative point  X′(with (X′, �Y ′, c) ∈ �S), there is a low degree of uncertainty on the label of X, and X is considered to be uninformative.

Using the assumption (H3), it suffices to have

image

Because the  PXappearing in (9) are unknown, it has to be replaced by an estimate. We will show that it can be estimated with arbitrary precision and confidence using only unlabelled data from K. The Reliable subroutine uses  EstProb(X, r, ǫo,50, δ) (inspired from (Kontorovich et al., 2016)) as follows:

image

1. Call the subroutine  BerEst(ǫo, δ,50). 2. To draw a single  piin  BerEst(ǫo, δ,50), sample randomly an example  Xifrom K, and set  pi = 1Xi∈B(X,r).

The subroutine BerEst consists in estimating adaptively with high probability the expectation of a Bernoulli variable  Z ∼ p. In our setting, we estimate a probability-ball, so that a realisation of Z can be set as  pi= 1Xi∈B(x,r). The variables  p1, . . . , p4are sampled at the beginning for theoretical analysis where we want a concentration inequality to hold for a number of samples greater than 4 (see (Kontorovich et al., 2016; Maurer and Pontil, 2009) for more details).

However, it is not dramatic if  PXis supposed to be known by the learner. This is not a limitation, since it can be assumed that the pool K of data is large enough such that  PXcan be estimated to any desired accuracy.

image

If a point X is considered informative, it is introduced in the ConfidentLabel (Algorithm(5)), along with an integer  k′, a budget parameter t and a confi-dence parameter  δ. This subroutine infers with high confidence (at least 1−δ) the label of X, by using the labels of its  k′nearest neighbors, knowing that we can request at most t labels. The parameter  k′is chosen such that, with high probability, the empirical majority of the  k′-NN labels differs from the majority in expectation by less than some margin, and all the  k′-NN are at most at some distance from X. The ConfidentLabel subroutine outputs �Y, Q where Q represents the set of labeled nearest neighbors in the subroutine, and �Yrepresents the majority label in Q.

image

This Section provides the main results and theoretical motivations behind the KALLS algorithm. Let us recall  K = {X1, . . . , Xw}is the pool of unlabeled data and n is the budget.

Let us denote by  Aa,wthe set of active learning algorithms on K, and  P(α, β) the set of probabilities that satisfy assumption (H3) and (H2).

Additionally, let us introduce the set of probabilities  P′(α, β) on  X × Y. A probability  P ∈ P′(α, β) if  P ∈ P(α, β) and its marginal probability  PXis a doubling-probability. For  A ∈ Aa,w, we denote by �fA,n,w:= �fn,wthe classifier that is provided by A.

Theorem 2 and its equivalent form in Theorem 3 are the main results of this paper. They provide bounds on the excess risk for the KALLS algorithm in terms of the set  P′(α, β). The main idea of the proof is sketched in Section 5.2, while a detailed proof can be found in Appendix A.

5.1 Main results

image

Let the set  P′(α, β) such that  αβ < dwhere d is the dimension of the input space  X ⊂ Rd. For  P ∈ P′(α, β), if �fn,wis the 1  − NNclassifier provided by KALLS, then we have:

image

The result (10) is also stated below (Theorem 3) in a more practical form using label complexity. This latter form will be used in the proof.

image

Before proving this theorem, a couple of important remarks should be made:

1. The rate of convergence (10) obtained in Theorem 2 is an improvement over

image

2. The rate (10) is also minimax. Indeed, let us introduce the set of probabilities ¯P(α, β) that satisfy the H¨older continuous assumption (H1a) (with parameter  α), strong density assumption (H1b), and margin noise assumption (H2).

Let us assume that  αβ < d. It was proven in (Minsker, 2012) that if supp(PX) ⊂[0, 1]d, there exists a constant  γ >0 such that for all n large enough and for any active classifier �fn, we have:

image

Moreover, the strong density assumption implies the doubling-probability assumption (Kpotufe, 2011), and according to Theorem 1, the lower bound obtained in (17) is also valid for the family of probabilities  P′(α, β).

image

For a classifier �fn,w, it is well known(Lugosi, 2002) that the excess of risk is:

image

We thus aim to prove that (11) is a sufficient condition to guarantee (with probability  ≥ 1−δ), that �fn,wagrees with  f ∗on the set  {x, |η(x)−1/2| > ∆o}, for a suitable choice of  ∆o >0.

Introducing  ∆oin (18) leads to:

image

Therefore, if  ∆o ≤ ǫ2then we have immediately,  R( �fn,w) − R(f ∗) ≤ ǫ. On the other hand, if  ∆o > ǫ2, by hypothesis (H2), we have  R( �fn,w) − R(f ∗) ≤2C∆β+1o .In the latter case, setting  ∆o = � ǫ2C� 1β+1guarantees  R( �fn,w) −R(f ∗) ≤ ǫ. Altogether, using for  ∆othe value  ∆= max( ǫ2,� ǫ2C� 1β+1) guarantees R( �fn,w) − R(f ∗) ≤ ǫ. This explains the expression (6).

We present the proof sketch of Theorem 3 in three main steps, and refer to the corresponding Lemmas and Theorems in the Appendix A for more details.

image

will be violated after at most ˜k(ǫ, δs) requests, with ˜k(ǫ, δs) ≤ k(ǫ, δs). Also, the label inferred after ˜k(ǫ, δs) label requests corresponds to the true label  f ∗(Xs). The intuition behind is to adapt the number of labels requested with respect to the noise; i.e., fewer label requests on a less noisy point (i.e.,  |η(Xs)− 12| ≥ 12∆), and more label requests on a noisy point. This provides significant savings in the number of requests needed to predict with high probability the correct label.

In the event  A1 ∩ A2, any informative point  Xsfalls in a high density region such that all the  k(ǫ, δs) nearest neighbors of  Xsare within at most some distance to  Xs, and the condition (12) is sufficient to have k(ǫ, δs) ≤ w.

image

(21) In this case, let  Xs′be such a point that satisfies (20) and (21), we can easily prove that when  Xs′is relatively far from the boundary, i.e.,  |η(Xs′)− 12| ≥∆

image

and easily deduce by using the smoothness assumption, (21) and (22), that the points  Xsand  Xs′have the same label, then we do not need to use  Xsin the subroutine ConfidentLabel (Algorithm 5) and  Xsis an uninformative point. In addition, (13) is a sufficient condition such that the number of

image

The set I is introduced in KALLS (Algorithm(1)) as the set of informative points indexes. Let  sI= max I, the index of the last informative point.

For  ǫ, δ ∈(0, 12), ∆defined in (6), and (α, L) the smoothness parameters, let us introduce

image

where  X(1)xis the nearest neighbor of x in �Sac, and �fn,wthe 1-NN classifier on �Sac. Additionally, on  A1 ∩ A2 ∩ A3 ∩ A4 ∩ A5we prove that (11) is sufficient to obtain (23).

image

In this paper we first reviewed the main results for convergence rates in a non-parametric setting for active learning, with a special emphasis on the relative merits of the assumptions about the smoothness and the margin noise. Then, by putting active learning in perspective with recent work on passive learning, we provided a novel active learning algorithm with a particular smoothness assumption customized for k-NN.

We showed that our algorithm has a convergence rate comparable to state-of-the art active learning algorithms, but using less restrictive assumptions.

This removes unnecessary restrictions on the distribution that would exclude important densities (e.g., Gaussian). Additionally, our algorithm can readily be extended to multi-class classifi-cation, and then compared to recent results obtained in passive learning by (Reeve and Brown, 2017) which extended the work of (Chaudhuri and Dasgupta, 2014) to multi-class classification. Finally, an important direction for further work is to extend our results to the case where the key parameters of the problem (defining smoothness and noise) are unknown. Ongoing work in this direction builds upon previous results in an adaptive setting (Locatelli et al., 2017), (Minsker, 2012), (Balcan and Hanneke, 2012), (Hanneke, 2011).

Jean-Yves Audibert and Alexandre B Tsybakov. Fast learning rates for plug-in classifiers. The Annals of statistics, 35(2):608–633, 2007.

Maria Florina Balcan and Steve Hanneke. Robust interactive learning. In Conference on Learning Theory, pages 20–1, 2012.

Maria-Florina Balcan, Steve Hanneke, and Jennifer Wortman Vaughan. The true sample complexity of active learning. Machine learning, 80(2-3):111– 139, 2010.

G´erard Biau and Luc Devroye. Lectures on the nearest neighbor method. Springer, 2015.

Rui M Castro and Robert D Nowak. Minimax bounds for active learning. IEEE Transactions on Information Theory, 54(5):2339–2353, 2008.

Kamalika Chaudhuri and Sanjoy Dasgupta. Rates of convergence for near- est neighbor classification. In Advances in Neural Information Processing Systems, pages 3437–3445, 2014.

Sanjoy Dasgupta. Two faces of active learning. Theoretical computer science, 412(19):1767–1781, 2011.

Sanjoy Dasgupta. Active learning theory. Encyclopedia of Machine Learning and Data Mining, pages 14–19, 2017.

Maik D¨oring, L´aszl´o Gy¨orfi, and Harro Walk. Rate of convergence of k-nearest- neighbor classification rule. The Journal of Machine Learning Research, 18 (1):8485–8500, 2017.

Gerald A Edgar. Packing measure in general metric space. Real Analysis Exchange, 26(2):831–852, 2000.

Herbert Federer. Geometric measure theory. Springer, 2014.

St´ephane Ga¨ıffas. Sharp estimation in sup norm with random design. Statistics & probability letters, 77(8):782–794, 2007.

Steve Hanneke. Nonparametric active learning, part 1: Smooth regression functions.

Steve Hanneke. Rates of convergence in active learning. The Annals of Statistics, 39(1):333–361, 2011.

Steve Hanneke and Liu Yang. Minimax analysis of active learning. The Journal of Machine Learning Research, 16(1):3487–3602, 2015.

Juha Heinonen. Lectures on analysis on metric spaces. Springer Science & Business Media, 2012.

Wassily Hoeffding. Probability inequalities for sums of bounded random vari- ables. Journal of the American Statistical Association, 58(301):13–30, 1963.

Emilie Kaufmann, Olivier Capp´e, and Aur´elien Garivier. On the complexity of best-arm identification in multi-armed bandit models. The Journal of Machine Learning Research, 17(1):1–42, 2016.

Aryeh Kontorovich, Sivan Sabato, and Ruth Urner. Active nearest-neighbor learning in metric spaces. In Advances in Neural Information Processing Systems, pages 856–864, 2016.

Samory Kpotufe. k-nn regression adapts to local intrinsic dimension. In Advances in Neural Information Processing Systems, pages 729–737, 2011.

Andrea Locatelli, Alexandra Carpentier, and Samory Kpotufe. Adaptivity to noise parameters in nonparametric active learning. Proceedings of Machine Learning Research vol, 65:1–34, 2017.

Andrea Locatelli, Alexandra Carpentier, and Samory Kpotufe. An adaptive strategy for active learning with smooth decision boundary. In Algorithmic Learning Theory, pages 547–571, 2018.

G´abor Lugosi. Pattern classification and learning theory. In Principles of nonparametric learning, pages 1–56. Springer, 2002.

Enno Mammen and Alexandre B Tsybakov. Smooth discrimination analysis. The Annals of Statistics, 27(6):1808–1829, 1999.

Pascal Massart and ´Elodie N´ed´elec. Risk bounds for statistical learning. The Annals of Statistics, 34(5):2326–2366, 2006.

Andreas Maurer and Massimiliano Pontil. Empirical bernstein bounds and sample variance penalization. arXiv preprint arXiv:0907.3740, 2009.

Stanislav Minsker. Plug-in approach to active learning. Journal of Machine Learning Research, 13(Jan):67–90, 2012.

Wolfgang Mulzer. Five proofs of Chernoff’s bound with applications. arXiv preprint arXiv:1801.03365, 2018.

Charles Chapman Pugh. Real mathematical analysis. Springer, 2002.

Henry WJ Reeve and Gavin Brown. Minimax rates for cost-sensitive learn- ing on manifolds with approximate nearest neighbours. In International Conference on Algorithmic Learning Theory, pages 11–56, 2017.

Shai Shalev-Shwartz and Shai Ben-David. Understanding machine learning: From theory to algorithms. Cambridge university press, 2014.

Mathukumalli Vidyasagar. Learning and generalisation: with applications to neural networks. Springer Science & Business Media, 2013.

A Detailed proof of Theorem 3

This Appendix is organized as follows: in Section A.1, we introduce some additional no-

tations. In Section A.2 we adaptively determine the number of label requests needed to

accurately predict the label of an informative point that is relatively far from the bound-

ary decision. In Section A.4, we provide some lemmas that give a sufficient condition for

a point to be informative. In Section A.5, we give theorems that allow us to classify each

instance relatively far from the decision boundary. Finally in Section A.6, we provide the

label complexity and establish Theorem 3.

A.1 Notations

Some notations that will be used throughout the proofs are listed here for convenience.

As defined in Section 2.3, let  B(x, r) = {x′ ∈ X, ρ(x, x′) < r} and ¯B(x, r) = {x′ ∈

X, ρ(x, x′) ≤ r}the open and closed balls with respect to the Euclidean metric  ρ, respec-

tively, centered at  x ∈ Xwith radius r > 0. Let supp(PX) = {x ∈ X, ∀r > 0, PX(B(x, r)) >

0} the support of the marginal distribution  PX.

For  p ∈ (0, 1], and x ∈ supp(PX), let us define

image

Let us recall for  Xs ∈ K = {X1, . . . , Xw}, we denote by  X(k)s its k-th nearest neighbor in

K, and Y (k)sthe corresponding label.

For an integer  k ≥ 1, let

image

A.2 Adaptive label requests on informative points

Lemma 1 (Chernoff bounds, (Mulzer, 2018))

Suppose  X1, . . . , Xmare independent random variables taking value in {0, 1}. Let X denote

their sum and  µ = E(X)its expected value. Then,

image

Lemma 2 (Logarithmic relationship, (Vidyasagar, 2013))

Suppose a, b, c > 0, abec/a > 4 log2(e), and u ≥ 1. Then:

image

Lemma 3 (Chaudhuri and Dasgupta, 2014)

For  p ∈ (0, 1], and x ∈ supp(PX), let us define  rp(x) = inf{r > 0, PX(B(x, r)) ≥ p}. For

all  p ∈ (0, 1], and x ∈ supp(PX), we have:

image

Theorem 4

Let  ǫ, δ ∈ (0, 1). Set ∆ = max(ǫ,� ǫ2C� 1β+1 ), and pǫ =�31∆1024L�d/α, where α, L, β, C are

parameters used in (H2) and (H3).

For  p ∈ (0, 1], and x ∈ supp(PX), let us introduce  rp(x) = inf{r > 0, PX(B(x, r)) ≥ p} and

ks := k(ǫ, δs) defined in (5) (where δs = δ32s2 ).

image

1 −δ16 , such that on  A1, for all 1 ≤ s ≤ w, if

image

then the  ksnearest neighbors of  Xs(in the pool K) belong to the ball  B(Xs, rpǫ(Xs)).

Additionally, the condition

image

is sufficient to have (29).

Proof

Fix  x ∈ supp(PX). For k ∈ N, let us denote  X(k)x , the kth nearest neighbor of x in the pool.

we have,

image

Then, by using Lemma 1 and Lemma 3, and if  kssatisfies (29), we have:

image

Fix  x = Xs. Given Xs, there exists an event  A1,s, such that  P (A1,s) ≥ 1 − δ/(32s2), and

on  A1,s, if

image

we have  B(Xs, rpǫ(Xs))∩{X1, . . . , Xw} ≥ ks. By setting  A1 = ∩s≥1A1,s, we have P (A1) ≥

1  − δ/16, and on A1, for all 1  ≤ s ≤ w, if ks ≤ (1 − τks,s)pǫ(w − 1), then B(Xs, rpǫ(Xs)) ∩

{X1, . . . , Xw} ≥ ks.

Now, let us proof that the condition (30) is sufficient to guarantee (29).

The relation (29) implies

We can see by a bit of calculus, that  τks,s ≤ 12, and then

image

where  b = 4c�1024L31

image

Now, we are going to apply the Lemma 2. If we set in Lemma 2

image

we can easily see that  c ≥ a, a ≥4 and then

image

Then, the relation

image

is sufficient to guarantee (32).

Let us note that the guarantee obtained in the preceding theorem corresponds to that

obtained in passive setting (w = n).

A.3 Motivation for choosing  ksfor  Xs

Lemma 4 (Hoeffding’s inequality,(Hoeffding, 1963))

image

Second version: Let  X1, . . . , Xmbe independent random variables such that  −1 ≤ Xi ≤ 1, (i = 0, . . . , m).We define the empirical mean of these variables by

image

Lemma 5 (Kaufmann et al., 2016)

Let  ζ(u) = �

image

tributed, such that, for all  v > 0, E(evX1) ≤ ev2σ2/2. For every positive integer t, let

St = X1 + . . . + Xt. Then, for all  γ > 1 and r ≥

image

Lemma 6

Let  m ≥ 1 and u ≥ 20. Then we have:

image

Proof

Define  φ(m) = m − u log(log(m)), and let  m0 = 2u log(log(u)). We have:

image

It can be shown numerically that  φ(m0) ≥ 0 for u ≥ 20.

Also, we have:  φ′(m) = m log(m)−um log(m) ≥ 0 for all m ≥ m0(notice that  m0 ≥ u for u ≥ 20).

Then it is easy to see that  φ(m) ≥ φ(m0) for all m ≥ m0. This establishes the lemma.

Theorem 5

Let  δ ∈ (0, 1), and ǫ ∈ (0, 1). Let us assume that  w satisfies (12). For Xs, set ˜k(ǫ, δs) (with

δs = δ32s2 ) as

image

where  c ≥ 7.106. For k ≥ 1, s ≤ w, let ∆ = max( ǫ2, � ǫ2C� 1β+1 ) and bδs,k defined in (4).

Then, there exists an event  A2, such that  P (A2) ≥ 1 − δ/8, and on A1 ∩ A2, we have:

1. For  k ≥ 1, �ηk(Xs) and ¯ηk(Xs) defined in (26), for all s ∈ {1, . . . , w},

image

Where ¯ksis the number of requests made in  ConfidentLabel(Xs).

Proof

1. Let us begin with the proof of the first part of Theorem 5. Here, we follow the proof of Theorem 8 in (Kaufmann et al., 2016), with few additional modifications.

image

η(X(k)s)) = 0, and the random variables�Y (i)s − η(X(i)s ), i = 1, . . . , k�are independent. Then by Lemma 4, given  {X1, . . . , Xw}, as Y (1)s − η(X(1)s) takes values in [−1, 1], wehave  E(ev(Y (1)s −η(X(1)s ))) ≤ ev2/2 for all v >0. Furthermore, set  z = log( 32s2δ ), andr = z + 3 log(z). We have r ≥ 8(e−1)2, and by Lemma 5, with  γ = 3/2, we have: P

image

Then, we have, given  s ∈ {1, . . . , w}, there exists an event  A′2,s such that P (A′2,s) ≥1  − δ/32s2, and simultaneously for all  k ≥1, we have:

image

image

(39) On the other hand, the right-hand side is smaller than:

image

image

We can apply the Lemma 6 in (42) by taking:  m = ek and u = 73728e|η(Xs)− 12 |2 . We havem ≥ 1 and u ≥20 and then, a sufficient condition to have (42) is:

image

As  |η(Xs) − 12| ≥ 12∆, we can easily see that ˜k(ǫ, δs) ≤ k(ǫ, δs). By taking the minimum value ¯ks = ¯k(ǫ, δs) that satisfies (44), we can see that when the budget allows us, the subroutine  ConfidentLabel requests ¯kslabels, and we have:

image

By setting  A′′2 = ∩s≥1A′′2,s, we have P (A′′2 ) ≥ 1 − δ/16, and we can deduce (34). We have on  A′2, for all s ≤ w, k ≤ k(ǫ, δs),

image

image

Thus we can easily deduce (35). By setting  A2 = A′2 ∩ A′′2, we have P (A2) ≥ 1 − δ/8 and on A1 ∩ A2, the item 1 and item 2 hold simultaneously.

A.4 Sufficient condition to be an informative point

As noticed in Section 4.3, a sufficient condition for a point  Xt (with t ≤ w) to be considered

as not informative is:

image

for some previous informative point  Xs (with (Xs, �Ys, �LBs) ∈ �S the current active set just

before attaining  Xt in KALLS(Algorithm(1))). Because  PXis unknown, we provide a com-

putational scheme sufficient to obtain (49).

Firstly we follow the general procedure used in (Kontorovich et al., 2016) to estimate adap-

tively the expectation of a Bernoulli random variable. And secondly, we apply it to the

Bernoulli variable  1A where A = {x, x ∈ B(X, r)} for r > 0 and X ∈ X.

Lemma 7 (Kontorovich et al., 2016)

Let δ′ ∈ (0, 1), ǫo > 0, t ≥ 7 and set g(t) = 1 + 83t +�

Bernoulli random variables with expectation  p. Let �pbe the output of  BerEst(ǫo, δ′, t). There

exists an event  A′, such that  P (A′) ≥ 1 − δ′, and on A′, we have:

1. If  �p ≤ ǫog(t) then p ≤ ǫo, otherwise, we have  p ≥ 2−g(t)g(t) ǫ0.

2. The number of random draws in the BerEst subroutine (Algorithm 4) is at most

where  ψ := max(ǫo, pg(t) ).

Lemma 8

Let  ǫ, δ ∈ (0, 1), r > 0. Let us assume that w satisfies (13).

There exists an event  A3, such that  P (A3) ≥ 1 − δ/16, we have, on  A3, for all s ≤ w:

If there exists  1 ≤ s′ < s, such that  Xs′is an informative point, and  (Xs′, �Ys′, �LBs′) ∈ �S

(the current active set just before attaining  Xs defined in KALLS(Algorithm(1))), and that

satisfies:

image

where

image

and

image

then

image

Otherwise, if (50) does not holds, i.e:

image

then min(PX(B(Xs, ρ(Xs′, Xs))), PX(B(Xs′, ρ(Xs′, Xs)))) ≥ 2847

Proof

By following the scheme of subroutine Estprob, this Lemma is a direct application of Lemma

7 by taking for all  s ≤ w, t = 50, ǫo =� 164L �LBs′�d/α, δ′ = δs, r = ρ(Xs, Xs′), A3,s := A′.

And then, if we set  A3 = ∩s≥1A3,s, we have P (A3) ≥ 1 − δ/16, and on the event  A3, we

can easily deduce (51) and (52) in each cases.

On the other hand, for all  s ≤ w, the number of draws in  Estprob(Xs, ρ(Xs, Xs′),� 164L �LBs′�d/α, 50, δs)

(respectively  Estprob(Xs′ , ρ(Xs, Xs′),� 164L �LBs′�d/α, 50, δs)) is always lower than w. In-

deed, by Lemma 7, the number of draws is at most:

image

Then we have:

image

In equation (53),  bδs′ ,|Qs′ |is defined by (4), and  |Qs′|represents the number of label requests

used in the subroutine ConfidentLabel(Algorithm (5)) at the stage  s′.

A.5 Label the instance space

image

Let I the set of indexes of informative points used in KALLS (Algorithm 1). Let us consider

its last update in KALLS (Algorithm 1) and also denoted it by I.

Then, set  sI = max Ithe index of the last informative point. Let �Sacbe the active set

obtained in KALLS (Algorithm 1) and denote by �fn,wthe output  1NN(�Sac). There exists an

event  A4 such that P (A4) ≥ 1 − δ/8, and on A1 ∩ A2 ∩ A3 ∩ A4, we have

image

2. If w satisfies (12) and (13) and the following condition holds

image

then, for all  x ∈ supp(PX) such that |η(x) − 12 | > ∆, there exists  s := s(x) ∈ I suchthat:

image

Proof This proof is based on results from (Hanneke) with few additional modifications. 1. Let us begin by proving the first part of Theorem 6. For  x ∈ supp(PX), let us introduce

image

By Lemma 3, we have  PX(B(x, r˜pǫ(x)) ≥ ˜pǫ. Then each ¯X ∈ {X1, . . . , XTǫ,δ} belongsto  B(x, r˜pǫ(x)) with probability at least ˜pǫ. If we denote �Pthe probability over the data, we have:

image

Then, there exists an event  A4, such that  P (A4) ≥ 1 − δ/8 and (55) holds on  A4. Andthen, we can easily conclude the first part.

2. For the second part of Theorem 6, let  x ∈ supp(PX). By (55), on  A4there exists  Xx ∈{X1, . . . , XTǫ,δ }such that:

image

As  sI ≥ Tǫ,δ, then there exists  s′ such that Xx := Xs′ and Xs′passes through the subroutine Reliable. We have two cases: a)  Xs′ is uninformative. Then there exists  s < s′, such that  Xsis an informative point, and

image

image

that contradicts (62), then we have  |η(Xs) − 12 | ≥ 132 ∆. Therefore, by (63), (64), we have:

image

b)  Xs′ is informative. In this case,  s = s′ and then we always obtains the equation (69), which becomes

image

On  A1 ∩ A2, by Theorem 5, the subroutine  ConfidentLabel(Xs) uses at most ˜k(ǫ, δs)request labels, and returns the correct label (with respect to the Bayes classifier) of  Xs.Let us proof that  f∗(x) = f∗(Xs). Let us assume without loss of generality that  η(Xs)−12 ≥0. We will show that  η(x) − 12 ≥0. We have:

image

Then  f∗(x) = f∗(Xs).As  |η(Xs) − 12| ≥ 12∆, by using Theorem 5 (the second part), we can easily see that (Xs, �Ys) ∈ �Sac (where �Ysis the inferred label of  Xsprovided by the subroutine ConfidentLabel in KALLS(Algorithm (1))). Let  X(1)xthe nearest neighbor of  x in �Sac. We have:

image

Then,  |η(X(1)x ) − 12 | ≥ (1 − 791985 )|η(x) − 12| ≥ 12∆and by Theorem 5, the subroutine ConfidentLabel(X(1)x ) outputs

image

A.6 Label complexity

Lemma 9

Let us assume that  w satisfies (12), (13), and w ≥ Tǫ,δ. Then, there exists an event  A5such that  P (A5) ≥ 1 − δ/8, and on A1 ∩ A2 ∩ A3 ∩ A5. The condition (11) is sufficient to guarantee (56).

Before beginning the proof, let us define a notion that will be used through the proof.

Definition 5

Let a set  F ⊂ supp(PX). Let {x1, . . . , xm} ⊂ F and p (0,1]. We say that the set  {x1, . . . , xm} ⊂

F is a p-probability-packing set of F if:

image

where  rpis defined by (25), and  a ∨ b = max(a, b) for a, b ∈ R

This notion of p-probability-packing comes from the Definition 1.4 in (Edgar, 2000). It will

be used on a particular set of the form  {x ∈ supp(PX), γ ≤ |η(x) − 12| ≤ γ′} (where

0  < γ < γ′). This allows us to upper bound the number of informative points where we have

a very high confidence for inferring their labels.

Proof

Let us consider the last update of I, the set of indexes of informative points used in

KALLS(Algorithm 1).

Set  sI = max I, the index of the last informative point. We consider two cases:

1.  First case: sI = w: we can easily see that (56) is satisfied, and we have trivially that the condition (11) is sufficient to guarantee (56). 2.  Second case: sI < w: then the total number of label requests up to  sI is:

image

ConfidentLabel(Xs, t) implicitly assumes that the process of label request do not takes into account the constraint related to the budget n (very large budget with respect to k(ǫ, δs)) such that  ConfidentLabel(Xs, t)=ConfidentLabel(Xs, t = ∞) . Then we have:

image

On the other hand, we want to guarantee the condition (56). For this, necessary for all s ∈ I, such that  s ≤ Tǫ,δ, and s < sI, at the end of the subroutine  ConfidentLabel(Xs, t),the budget n is not yet reached and then we can replace the relation (77) by

image

Then, necessarily, (56) holds when (78) holds. Also, for  s ∈ I, by theorem5, if we assume that  |η(Xs) − 12 | ≥ 12 ∆, we have that  |Qs| ≤˜k(ǫ, δs), and the subroutine  ConfidentLabel(Xs, t), (with t = n − �si∈I,si<s |Qsi|)terminates when the cut-off condition (34) is satisfied. The right hand side of (78) is equal to:

image

image

Let  Xs, Xs′any two informative points that belong to  Ij. Let us assume that  s < s′.As  Xs ∈ Ij, we have |η(Xs)− 12 | ≥ ∆2 and by Theorem 5, the number of label requests ¯ks used in ConfidentLabel(Xs) satisfies:

image

where  �η¯ks := �η¯ks(Xs) and bδs,¯ksare respectively defined by (26) and (4). Then

image

Additionally, as  Xs and Xs′are both informative points, by Lemma 8, we necessary have on event  A3(see Lemma 8), that

image

(89) Let us introduction the quantity ¯η¯ks := ¯η¯ks(Xs) defined by (26). We have, by Theorem 5, on the event  A2(see Theorem 5),

image

image

(93) As the same way, we also obtain (93) if  s′ < s.Then, if we set

image

Therefore, the set of informative points that belong to  Ij = {x, γj−1 ≤ |η(x)− 12| ≤γj} forms an pj-probability-packing set.

(b) As second step, let us determine an upper bound of the cardinal of any  pj-probability-packing set of  Ij. Let Λj = {x1, . . . , xMj } any pj-probability-packing set of  Ij.For all  s, s′ ≤ Mj, we obviously have:

image

On the other hand, if  z ∈ B(xs, rpj (xs)/2) for some  s ≤ Mj, by the assumption (H3), equation (25), and the fact that  |η(xs) − 12 | ≤ γj, we have:

image

max�ǫ2, � ǫ2C� 1β+1�.Now, it remains to upper bound the second term in (79). Denote it by  T2. By Lemma1, (equation(28)), there exists an event  A5 such that P (A5) ≥ 1 − δ/8, and on A5, wehave:

image

Because  |Qs| ≤ k(ǫ, δ) (according to the subroutine  ConfidentLabel) for all s ≤ Tǫ,δand where ˜B = {x, |η(x) − 12 | ≤ ∆2 } and k(ǫ, δ) is defined in (5). Consequently, we have:

image

As  αβ ≤ d, ∆ ≤ 1, C ≥1, the term (128)d/α−β 64βC � 1∆�d/α−βin (103) is greater than 1. Thus, (103) becomes:

image

Thus, if the label budget n satisfies n  ≥ b0

image

We have that n satisfies (78), and (56) is necessary satisfied.

Proof of Theorem 3

Finally, we are able to prove Theorem 3.

Previously, we have designed five events  A1, A2, A3, A4, A5, such that if the event A :=

A1 ∩ A2 ∩ A3 ∩ A4 ∩ A5, and if the label budget  n ≥ ˜O�� 1ǫ� 2α+d−αβα(β+1) �(105) (where ˜O

assumes the existence of polylogarithmic factor in 1ǫ and 1δ ), then (56) is necessary satisfied,

therefore , by Theorem 6, if w satisfies (12) and (13), and the event A holds, the classifier

provided by the algorithm KALLS agrees with the Bayes classifier on  {x, |η(x) − 12| > ∆}.

Consequently, we have:

image

Thus with probability at least 1  − P (Ac) ≥ 1 − δ, (14) holds.


Designed for Accessibility and to further Open Science