Nuclear discrepancy for single-shot batch active learning

Active learning algorithms propose what data should be labeled given a pool of unlabeled data. Instead of selecting randomly what data to annotate, active learning strategies aim to select data so as to get a good predictive model with as little labeled samples as possible. Single-shot batch active learners select all samples to be labeled in a single step, before any labels are observed. We study single-shot active learners that minimize generalization bounds to select a representative sample, such as the maximum mean discrepancy (MMD) active learner. We prove that a related bound, the discrepancy, provides a tighter worst-case bound. We study these bounds probabilistically, which inspires us to introduce a novel bound, the nuclear discrepancy (ND). The ND bound is tighter for the expected loss under optimistic probabilistic assumptions. Our experiments show that the MMD active learner performs better than the discrepancy in terms of the mean squared error, indicating that tighter worst case bounds do not imply better active learning performance. The proposed active learner improves significantly upon the MMD and discrepancy in the realizable setting and a similar trend is observed in the agnostic setting, showing the benefits of a probabilistic approach to active learning. Our study highlights that assumptions underlying generalization bounds can be equally important as bound-tightness, when it comes to active learning performance. Code for reproducing our experimental results can be found at https://github.com/tomviering/NuclearDiscrepancy.


Introduction
Supervised machine learning models require enough labeled data to obtain good generalization performance.For many practical applications such as medical diagnosis or video topic prediction labeling data can be expensive or time consuming (Settles 2012).Often in these settings unlabeled data is abundant.In active learning an algorithm chooses unlabeled samples for labeling (Cohn et al. 1994).The idea is that models can perform better with less labeled data if the labeled data is chosen carefully instead of randomly.Active learning makes the most of a small labeling budget and can reduce labeling costs.
Several works use upperbounds on the expected loss to motivate particular active learning strategies (Gu and Han 2012;Ganti and Gray 2012;Gu et al. 2012Gu et al. , 2014;;Wang and Ye 2013).We study pool-based active learners that choose queries that explicitly minimize generalization bounds and investigate the relation between bounds and active learning performance.We evaluate generalization with respect to the surrogate loss in the classification setting and use the kernel regularized least squares model (Rifkin et al. 2003), a popular model in active learning (Huang et al. 2010;Wang and Ye 2013).Our focus is on active learners that select a batch of queries in a single shot (Contardo et al. 2017).This means that there is no label information available at the time the batch of queries is determined.Since the active learners have only have unlabeled data at their disposal they aim to select the most representative subset of the unlabeled pool.This is different from batch mode or sequential active learning, where after requesting labels from the oracle the algorithm has to determine new queries, creating a feedback loop.The advantage of zero-shot active learning is that all queries can be computed ahead of time, and collected labels do not have to be fed into the active learner.
For applications this can be very convenient: it simplifies the annotation setup.Furthermore, active learning algorithm may require substantial amounts of time to compute the next query.In situations where annotation have to be done by domain experts whose time is costly this can be impractical.For example, if we were to apply active learning to to the problem of Esteva et al. (2017), who build a deep learning model to classify skin cancer, sequential or batch mode active learning strategies usually train a model as intermediate step before being able to determine the next query.For deep models this could take several hours.With zero-shot active learning the dermatologist can annotate all queries without waiting once.
Another example where requesting labels is costly is personalized machine learning models such as for movie recommendation.Here applications may ask feedback from end-users to improve their service.This problem can also be studied using the active learning framework (Harpale and Yang 2008).Asking end-users for feedback usually interrupts their activity in the application.Therefore, we may only interrupt the user a limited amount of times.Using zero-shot active learning users only have to be interrupted once and can answer multiple queries without waiting for new queries to be determined.
The Maximum Mean Discrepancy (MMD) is used for batch-mode active learning by Chattopadhyay et al. (2012) to match the marginal distribution of the selected samples to the marginal distribution of all unlabeled samples.This active learner has been shown to minimize a generalization bound (Wang and Ye 2013).The MMD is a divergence measure (Gretton et al. 2012) which is closely related to the Discrepancy divergence measure of Mansour et al. (2009), both have been used in domain adaptation (Huang et al. 2007;Cortes and Mohri 2014).
Using the Discrepancy, we show that we can get a tighter worst case generalization bound than the MMD in the realizable setting.Tighter bounds are generally considered better as they estimate the expected loss more accurately.One might therefore expect the Discrepancy to lead to better queries in active learning.
We show, however, that the Discrepancy and MMD generalization bounds can be derived, using a probabilistic analysis, from pessimistic assumptions.We subsequently apply the principle of maximum entropy to derive probabilistic assumptions that are more optimistic, inspiring us to introduce the Nuclear Discrepancy (ND) bound.Under these optimistic assumptions the ND provides a tighter bound on the expected loss than the MMD, while the Discrepancy bound is the loosest.
We compare the active learning performance of the proposed ND bound to the existing MMD and Discrepancy bounds.Our hypothesis is that we often find ourselves in a more optimistic average-case scenario than a worst-case scenarios.To this end we empirically study the behavior of the active learners on 13 datasets, and we investigate whether probabilistic assumptions or worst-case assumptions better model observed behavior in our experiments.
In the realizeable setting a model from the model class can perfectly predict the groundtruth labels, as in this setting there is no model mismatch or model misspecification.For this we show that the tightness relations between the generalization bounds is strict.As such, for the realizeable case, our theory gives the strongest predictions for the ranking of the active learners in terms of performance.In the agnostic case, where no such model may exist, the tightness relations can change, which renders our theory less applicable.We perform experiments in both settings to see the effect of the theoretical assumptions not being fulfilled.
We study the realizable setting since it is more amendable to theoretical analysis.This setting is often studied in active learning and is still a topic of active investigation (Tosh and Dasgupta 2017).The general case of the agnostic case is much harder to analyze.To illustrate this, we remark that it has been observed that if a model class is sufficiently wrongly chosen, active learning can even decrease model performance (Settles 2011;Attenberg and Provost 2011;Loog and Yang 2016;Yang and Loog 2018).
These counter-intuitive behaviors further underline the need for further theoretical studies.We believe that by improving our understanding of simpler active learning settings (realizeable case) will contribute to improved understanding of more difficult active learning settings (agnostic case).
To this end, our study provides new quantitative tightness relations between the MMD, Discrepancy and ND bound under different probabilistic assumptions.We investigates the connection between bound tightness and active learning performance.Our most important conclusion is that not only bound tightness is important for performance, but that appropriate assumptions are equally important.

Overview and contributions
First we discuss related work in Sect. 2. In Sect. 3 we describe the considered active learning setting and notation.We present our theoretical results regarding the MMD and Discrepancy in Sect. 4. In Sect. 5 we motivate our novel Nuclear Discrepancy bound.We evaluate the proposed active learners experimentally in Sect.6.In Sect.7 we give a discussion and in Sect.8 we give the conclusions of this work.All proofs, additional background theory and experimental results are given in the Appendix.The main contributions of this work are: 1.An improved MMD bound for active learning and a more informed way to choose the kernel of the MMD in the context of learning.2. A proof that the Discrepancy bound on the worst case loss is tighter than the MMD bound.3. A probabilistic interpretation of the MMD bound.This table gives an overview of the newly proven tightness relations between the generalization bounds and the experimental results.Observe that the tightness relations under the 'Average-Case' correlate well with the experimental performance of the active learners.Therefore, we stipulate that the 'Average-Case' is the most accurate assumption for our considered active learning setting.Note that the tightness relations only hold under the conditions of Theorem 2, and that the experimental performance shown here best reflect the performance in the realizable setting.In the agnostic setting the ranking of the active learning methods is less clear, but the same trend is observed 4. The Nuclear Discrepancy (ND) bound that provides the tightest bound on the expected loss under probabilistic assumptions that follow from the principle of maximum entropy. 5.A probabilistic analysis that explains the differences in empirical performance (in terms of the mean squared error) achieved by the active learners.
In Table 1 we give a visual summary of our work.It shows all formal results and shows in which sections to find them.It also shows the relation between the theory and experiments, and the main findings of the experiments.

Related work
Many active learning methods have been proposed, Settles (2012) provides an excellent introduction and overview.Our work is related to active learning methods that select representative samples (Xu et al. 2003).Most active learning strategies of this kind are combined with an uncertainty criteria (Xu et al. 2003;Chattopadhyay et al. 2012;Wang and Ye 2013;Huang et al. 2010), and often the representative component is used to diversify queries when chosen in batches in order to avoid redundancy (Xu et al. 2003;Wang and Ye 2013).This is different from our considered setting: since there is no labeled data and we have to choose all queries in one shot, our only option is to select representative samples, since uncertainty criteria can only be computed if some labels are known.
A closely related well-known concept to our work is that of (Transductive or) Optimal Experimental Design (Yu et al. 2006).Here also no labeled data is required to select queries for the case of the linear regression model.These methods aim to minimize some form of posterior variance of the model.A closely related statistical approach relies on maximization of the Fisher Information to reduce model uncertainty (Hoi et al. 2006).However, for these approaches it is often required to explicitly specify a noise model (such as Gaussian i.i.d.noise), while in this work we consider deterministic labeling functions.
Our work is motivated by several active learners that minimize generalization bounds.Gu and Han (2012) uses the Transductive Rademacher Complexity generalization bound to perform active learning on graphs.Gu et al. (2012) show that the strategy of Yu et al. (2006) also minimizes a generalization bound, and extend the method to work with a semi-supervised model.Ganti and Gray (2012) introduce an active learning strategy that uses importance weighting to ensure asymptotic consistency of the actively learned model.Their strategy minimizes a generalization bound for the squared loss under some conditions on the data distribution.Gu et al. (2014) introduce an strategy that minimizes a generalization bound on the risk for logistic regression.Wang and Ye (2013) also use a generalization bound based on the MMD to perform active learning, but we will describe this work later in more detail when discussing all methods that use the MMD.
Many theoretical active learning works motivate algorithms by generalization bounds, for example one of the first active learning algorithms 'CAL' (Cohn et al. 1994) and its agnostic generalization A 2 (Balcan et al. 2009) have been thoroughly analyzed using generalization bounds by making use of the Disagreement Coefficient (Hanneke 2007).Most of these theoretical works consider worst-case performance guarantees, where the distribution is chosen by an adversary subject to constraints.Balcan and Urner (2016) provides a short and concise overview of these and other recent theoretical active learning works.In contrast with our work, these algorithms consider generalization in terms of zero-one loss instead of squared loss and do not apply to one shot active learning.
A straightforward approach to one shot active learning is through clustering: cluster the data and request the labels of the cluster centers (BodÃ et al. 2011;Hu et al. 2010;Zhu et al. 2008;Nguyen and Smeulders 2004).However, unlike our work, these methods are not motivated by generalization bounds.Obtaining bounds for such approaches may be difficult because the clustering algorithm and machine learning model may rely on different assumptions.To still get bounds one can use the clustering algorithm instead to also provide predictions for new samples (Urner et al. 2013).Instead, we stick to the regularized least squares model and use the MMD and Discrepancy to get bounds for this model.Our approach can be used to derive bounds and corresponding active learning strategies for any kernelized L 2 regularized model, however, in this work we only focus on the squared loss.
Our work is closely related to that of Chattopadhyay et al. (2012): we use a greedy version of their proposed active learning algorithm.Chattopadhyay et al. (2012) are the first to use the MMD for active learning in a batch-mode setting.An in-depth empirical analysis shows that the MMD outperforms other active learning criteria as judged by the zero-one error when used with kernelized SVMs.They show that the MMD easily can be combined with uncertainty-based active learning approaches and transfer learning.Since we consider oneshot active learning we don't consider the uncertainty-based component of their algorithm.In follow up work active learning and transfer learning is solved jointly using the MMD (Chattopadhyay et al. 2013).
Our theoretical analysis of the MMD bound extends the analysis of Wang and Ye (2013).Wang and Ye (2013) show that active learning by minimization of the MMD and the empirical risk can be seen as minimizing a generalization bound on the true risk.They introduce an active learner that balances exploration (distribution matching using MMD) with exploitation (a form of uncertainty sampling).They show empirically that their proposed algorithm is competitive with several other active learning strategies as evaluated by the zero-one error using kernelized SVMs.
We build upon the generalization bound of Wang and Ye (2013) and improve it.Their bound considers the underlying distribution of the unlabeled pool and labeled (queried) sample, however, this is problematic because the labeled sample is non-i.i.d.due to dependence of the queries of the active learner.We resolve this issue and introduce an additional term η that measures the error of approximating the worst-case loss function.Mansour et al. (2009) introduce the Discrepancy generalization bound for domain adaptation with general loss functions.In a follow up work, Cortes and Mohri (2014) contrast the Discrepancy with the MMD generalization bound: they argue that the Discrepancy is favorable from a theoretical point of view because it takes the loss function and hypothesis set of the model into account, while the MMD does not.This means that the MMD bound for an SVM and regularized least squares model would be exactly the same, while the Discrepancy bound specializes to the chosen model and surrogate loss.They derive an efficient domain adaptation algorithm and empirically show that the Discrepancy improves upon the MMD in several regression adaptation tasks.
Prior to our work, the Discrepancy measure (Cortes and Mohri 2014) has not yet been used to perform active learning.We show that by choosing the kernel for the MMD carefully, we can adapt the MMD to take the hypothesis set and loss into account, addressing one of the theoretical limitations of the MMD identified by Cortes and Mohri (2014).Under these conditions we find that we can compare the MMD and Discrepancy bounds in terms of tightness.This quantitative comparison of these bounds is novel and was not considered before.Germain et al. (2013) adapt the Discrepancy for the zero-one loss to a PAC-Bayes setting in order to do domain adaptation.Their analysis is specifically for the zero-one loss, while we consider the squared loss.Their PAC-Bayes framework is significantly different from our analysis: instead of minimizing a surrogate loss, they use a Gibbs classifier, and they minimize bounds on the expected risk directly.This involves a non-convex optimization problem.Instead, we simply minimize the empirical risk and consider deterministic models, similar to most PAC style analysis.This makes our analysis is simpler.Furthermore, they propose a framework to jointly minimize the empirical risk and domain divergence.To this end, their algorithm requires labeled data which is unavailable in zero-shot active learning, making it unsuitable for our zero-shot setting.
In Cortes et al. (2019) a new domain adaptation algorithm based on a new divergence measure, the Generalized Discrepancy, is introduced.The algorithm consists of two stages: first it minimizes the Discrepancy, afterward it minimizes the empirical risk and the Generalized Discrepancy jointly.The strategy of Cortes et al. (2019) is difficult to apply to active learning for two reasons.First of all, their algorithm requires labeled data to minimize the empirical risk and the General Discrepancy jointly, which is impossible in our zero-shot active learning setting.Second, their algorithm requires i.i.d.samples from the unlabeled pool to estimate the hyperparameter r .This would require costly random queries in the active learning setting.Because of these reasons, we believe their algorithm is more suitable to a joint active and domain adaptation setting (such as considered by Chattopadhyay et al. (2013)) where more labeled data is available.
Our theoretical analysis is substantially different from the analysis of Cortes et al. (2019).Because Cortes et al. (2019) use labeled data, they can make a more accurate characterization of possible worst case scenario's, refining the worst-case scenario of the Discrepancy to obtain tighter bounds.We take an orthogonal approach: we consider probabilistic generalization bounds that hold in expectation.Instead of considering a worst-case, we make probabilistic assumptions to get to a plausible average-case.Cortes et al. (2019) compare the Generalized Discrepancy and Discrepancy bounds in terms of tightness.We compare the tightness of the bounds of the MMD, Discrepancy and Nuclear Discrepancy.We show several orderings of the tightness of the bounds under different probabilistic assumptions, while Cortes et al. (2019) only takes a worst-case approach.
In summary, our work differs from previous works by considering instead of worst-case analysis (Cortes et al. 2019;Cortes and Mohri 2014), a probabilistic analysis of generalization bounds.Unlike most other works that use generalization bounds for domain adaptation (Cortes et al. 2019;Cortes and Mohri 2014;Germain et al. 2013), we use bounds to perform active learning.For the MMD active learner, studied by Chattopadhyay et al. (2012); Wang and Ye (2013), we give new theoretical results: an improved bound for active learning and we provide a principled way to choose the kernel for the MMD.We give new quantitative comparisons of bound tightness for the MMD and Discrepancy in multiple settings, while before these bounds were compared only qualitatively (Cortes and Mohri 2014).Furthermore, we study the novel question: how does bound tightness relate to active learning performance?

Setting and notation
Let X = R d denote the input space and Y the output space.Like Cortes and Mohri (2014) we assume there is a function f : X → Y that determines the outputs and there is an unknown distribution with density P over X from which we get an independent and identically distributed (i.i.d.) unlabeled sample P = (x 1 , . . ., x n P ) ∈ X n P .We study single-shot batch active learners that given the unlabeled pool P selects a batch Qn ⊂ P of n samples before observing any labels.The active learner submits the batch to the labeling oracle that provides the labels of the batch.A kernel regularized least squares (KRLS) model is trained on Qlab n , where lab indicates a labeled dataset.
We take the kernel of the model K to be positive semi-definite (PSD), and denote the reproducing kernel Hilbert space (RKHS) as H where ||h|| K denotes the norm in H.A model corresponds to h ∈ H and is obtained by minimizing for h ∈ H when trained on Qlab , where we follow the convention of Cortes and Mohri (2014).L Q (h, f ) is the average empirical loss of h on Q with outputs given by f : where l : R × R → R is a loss function.For KRLS l is the squared loss: ) is the mean squared error (MSE) on Q. Model complexity is controlled by the regularization parameter μ > 0. We choose as our hypothesis set where Training KRLS always leads to a solution h ∈ H (Mohri et al. 2012, Lemma 11.1).
In classification typically we are interested in the zero-one error (accuracy), however, our study focuses on the squared loss (the surrogate loss).We use the squared loss because we can relate the bounds of the MMD, Nuclear Discrepancy and Discrepancy in closed form and compare them quantitatively.Since our goal is to investigate the correlation between bound tightness and performance, this is essential to our study.
We have made the standard assumption that the data comes from an unknown distribution P. The goal of the active learner is to choose a batch of queries in such a way as to minimize the expected loss of the model under this distribution P: (1) Ideally we would want to train our model on Plab , since small L P (h, f ) will lead to small L P (h, f ) if the model complexity is appropriate, as illustrated by the following theorem (Mohri et al. 2012, p. 240).
Theorem 1 (Generalization bound Squared Loss (Mohri et al. 2012)) Let l be the squared loss.For any δ > 0, with probability at least 1 − δ over an i.i.d.sample P of size n P from P, the following inequality holds for all h ∈ H: Here R m (H ) is the Rademacher complexity of the hypothesis set H , and M is a constant such that |h(x) − f (x)| ≤ M for all x ∈ X and all h ∈ H.
If the model complexity is appropriate R m (H ) will be small.The third term is small when the pool P is large.If both of these criteria are met, it is unlikely that we overfit as reflected by a tight bound.Then training on Plab will likely minimize L P (h, f ).
Ideally we would train on Plab , however, since we only have access to the unlabeled sample P this is impossible.Therefore we upperbound L P (h, f ) instead.This upperbound is minimized by the active learners.The studied bounds are of the form Due to training L Q (h, f ) will be relatively small.The term η is a constant that cannot be minimized during active learning since it depends on Plab .However, if the model misspecification is small, η will be small.Therefore we ignore this term during active learning, this is also (sometimes implicitly) done in other works (Huang et al. 2007;Chattopadhyay et al. 2012;Cortes and Mohri 2014).Thus the active learners choose the batch Q to minimize obj( P, Q).This objective can be the MMD, disc or disc N which will be introduced in the next sections.This term measures the similarity between the unlabeled pool P and the batch Q. Minimizing it leads to selecting a representative sample.
We consider two settings.In the agnostic setting binary labels are used, i.e., Y = {−1, +1}, and generally we have f / ∈ H .In the realizable setting f ∈ H , so a model of our hypothesis set can perfectly reproduce the labels as there is no model misspecification.In this case Y is a subset of R. In the realizeable setting η can become zero under some conditions, which allows us to compare the tightness of the bounds and enables our probabilistic analysis.
K (x, x ) indicates the kernel function between x and x .We mainly use the Gaussian where σ , the bandwidth, is a hyperparameter of the kernel.For the MMD we require a second PSD kernel, K L .We indicate its RKHS and bandwidth (for a Gaussian kernel) by H L and σ L , respectively.All vectors are column vectors.X P and X Q are the n P × d and n Q × d matrices of the sets P and Q.

Analysis of existing bounds
First we provide an improved MMD generalization bound for active learning which is inspired by Cortes et al. (2019).Then we review a bound in terms of the Discrepancy of Cortes et al. (2019) and we review how to compute the Discrepancy quantity (Mansour et al. 2009).We show that the MMD can be computed using a novel eigenvalue analysis, and thereby making the MMD and Discrepancy bounds comparable.We wrap up the section with a probabilistic interpretation of both bounds.As a roadmap for the reader we give an overview of the tightness relations in Table 1 which will be proven in this section and the next section.

Improved MMD bound for active learning
The MMD measures the similarity between the two unlabeled samples Q and P. Using this criterion we give a generalization bound similar to the one given by Wang and Ye (2013) suitable for active learning.The empirical MMD quantity is given by here l is the worst-case function from a set of functions H L .We take the standard choice In Appendix A.1 we revisit how to compute the MMD quantity.We extend the technique of Cortes et al. (2019) to give a generalization bound in terms of the MMD.To get a bound for the MMD we approximate the loss function Proposition 1 (Agnostic MMD worst case bound) Let l be any loss function l : R × R → R.
Then for all h ∈ H and any labeling function f : X → Y we have Here η MMD measures the approximation error since we may have that g(h, f ) / ∈ H L .Our MMD bound above differs in two aspects from the bound of Wang and Ye (2013).Wang and Ye (2013) estimate the MMD between the distributions P and Q.However, to estimate the MMD between distributions i.i.d.samples are required (Gretton et al. 2012, Appendix A.2).The sample Q is not i.i.d.since it is chosen by an active learner.
Our bound allows for non-i.i.d.samples since it estimates the MMD between empirical samples and is therefore better suited for active learning.The second novelty is that we measure the error of approximating the loss function g(h, f ) using the term η MMD .This allows us to adjust the MMD to the hypothesis set H and loss l similar to the Discrepancy measure of Cortes and Mohri (2014).We give the theorem below with a small proof sketch for the simplified case of the linear kernel.See the Appendix for the full proof.

Theorem 2 (Adjusted MMD) Let l be the squared loss and assume f ∈ H (realizable setting). If K
Proof sketch Here we give a proof sketch for the case where K is the linear kernel: are all monomials of degree 2 (Shawe-Taylor and Cristianini 2004, chap.9.1).Therefore H L can be used to model any quadratic function such as g(h, f ).Therefore if Λ L is chosen appropriately we have g(h, f ) ∈ H L .

Corollary 1 Let l be the squared loss and f ∈ H and let K be a Gaussian kernel with bandwidth
Compared to other works Theorem 2 gives a more informed way to choose the MMD kernel in the context of learning.1Typically, a Gaussian kernel is used for the MMD with σ L = σ .However, Corollary 1 shows that if σ L = σ , we may have that η MMD = 0 even in the realizable setting, since σ L is too large-the true loss function g(h, f ) is less smooth than the functions in H L .This is undesirable since η MMD cannot be minimized during active learning.Our choice for σ L is preferable, as it ensures η MMD = 0 in the realizable setting.

Discrepancy bound
The Discrepancy is defined as (4) Observe it depends on H and l and therefore automatically adjusts to the loss and hypothesis set.We give a bound of Cortes et al. (2019) in terms of the Discrepancy.
Theorem 3 (Agnostic Discrepancy worst case bound (Cortes et al. 2019)) Assume that for all x ∈ X and for all h ∈ H that l(h(x), f (x)) ≤ C and let l be the squared loss.Then for all h ∈ H and any labeling function f : X → Y we have where Here η disc measures the model misspecification.In the realizable setting, f ∈ H , and η disc = 0.

Eigenvalue analysis
We show the relation between the Discrepancy and MMD using a novel eigenvalue analysis.To this end we introduce the matrix M P, Q to compute the Discrepancy.
For notational convenience we will often write M instead of M P, Q .The matrix M measures the difference between two sets of samples using their second-order moment.Considering its kernelized version such comparison can implicitly take higher-order moments into account as well.In particular, for a Gaussian kernel all moments of the samples are compared and we have that M = 0 only if P = Q.
In the following we will look at the eigendecomposition of M. Since M is the difference between two covariance matrices, it can have positive and negative eigenvalues.A positive (negative) eigenvalue means that in the direction of the corresponding eigenvector P has more (less) variance than Q.Recall that in active learning, our aim is to approximate P using representative samples Q, and thus small absolute eigenvalues are desirable, because this would indicate that in the direction of the corresponding eigenvector P is well approximated by Q.
Theorem 4 (Discrepancy computation (Mansour et al. 2009)) Assume K is the linear kernel, K (x i , x j ) = x T i x j , and l is the squared loss, then where λ i are the eigenvalues of M, and λ is the vector of eigenvalues of M.
Note that h will later play the role of f , the true labeling function.The theorem shows that in the worst case, the h and h that maximize the Discrepancy in Eq. 4 are chosen exactly in the direction where Q and P differ most, i.e., the direction of the largest absolute eigenvalue.Cortes and Mohri (2014) show that we can replace M by M K to compute the Discrepancy for any PSD kernel. 2efore we can give our main result we require some additional notation.Assume that the eigenvalues λ i of M are ordered by absolute value where |λ 1 | is the largest absolute eigenvalue.λ indicates the vector of eigenvalues, with r = rank(M) non-zero eigenvalues.e i is the normalized (unit-length) eigenvector corresponding to λ i .By careful analysis we can realize the relationship between M and the featuremap of the squared kernel to show that the MMD can be computed as follows.
Theorem 5 This theorem shows that the MMD measures differences between the samples Q and P differently.The Discrepancy only measures similarity along one dimension, namely the direction where the samples differ the most.The MMD considers all dimensions to compare the samples Q and P. Due to the square in the Euclidean norm, the MMD gives directions that differ more more weight in the comparison.
Under these conditions the Discrepancy bound (Theorem 3) is tighter than the MMD bound (Proposition 1), since η MMD = η disc = 0. Since the Discrepancy bound is tighter, one may expect that active learning by minimization of the Discrepancy may result in better active learning queries than minimization of the MMD, in particular if η MMD and η disc are small or zero.

Probabilistic analysis
We show the MMD can provide a tighter bound on the expected loss under certain probabilistic assumptions.From this point on we assume the conditions of Theorem 2 and take h to be the model trained on the set Q, and f to be the true labeling function.In addition, define u = h − f and U = {u ∈ H : ||u|| K ≤ 2Λ} and let ūi = u T e i , where e i is the eigenvector of M. Then It is more convenient to work with ū, since then the matrix M diagonalizes: The difference u is the unknown error our trained model h makes compared with the true model f .By making different probabilistic assumptions about the distribution of u we can arrive at different bounds.We now provide the building block for our probabilistic bounds.

By noting that L
and by making use of the triangle inequality, we find the following.
Lemma 1 (Probabilistic bound) Assume3 u is distributed according to a pdf p(u) over U .Then where we defined G(u, M) ) is generally small due to the training procedure of the model, thus generally E u L Q (h, f ) will be small as well.Therefore we focus our probabilistic analysis on the term E u G(u, M).By giving bounds on this quantity, we derive several probabilistic bounds that hold in expectation w.r.t.u.
The Discrepancy can be interpreted to put all probability mass on u = 2Λe 1 .
Proposition 2 (Worst case: Probabilistic Discrepancy) Given the pdf p(u) = δ(u − 2Λe 1 ) where δ(x) is the Dirac delta distribution.Then Only one u ∈ U can be observed under this pdf.This is a worst case distribution because this p(u) maximizes E u G(u, M).The Discrepancy assumes that the model error u points exactly in the direction that causes us to make the biggest error on P.Under this distribution the Discrepancy gives a tighter bound on the expected loss than the MMD because of Corollary 2. Under a different p(u) the MMD bound is tighter.
Theorem 6 (Pessimistic case: Probabilistic MMD) Let p(u) be a pdf on U s such that4 then Unlike for the distribution of the Discrepancy, for the above p(u) it is possible to observe different model errors u.However, the model error u in this case is biased: Equation 9suggests that u is more likely to point in the direction of eigenvectors with large absolute eigenvalues.This assumption is pessimistic since large absolute eigenvalues can contribute more to E u G(u, M).Another way to interpret this is that model errors are more likely to occur in directions where Q and P differ more.Because Q and P differ more in those directions, these model errors can count more towards the MSE on P.
For this p(u) the MMD bound is tighter.If the probabilistic assumption of the MMD is more accurate, we can expect that the MMD active learner will yield better active learning queries than the Discrepancy.

Nuclear discrepancy
In this section we motivate the optimistic probabilistic assumption that leads to the Nuclear Discrepancy (ND) bound.First, let us introduce the Nuclear Discrepancy quantity In the absence of any prior knowledge, we choose the pdf p(u) according to the well established principle of maximum entropy.This principle dictates that in case nothing is known about a distribution, the distribution with the largest entropy should be chosen (Jaynes 1957).Accordingly, we choose p(u) uniform over U , which leads to the following.
Theorem 7 (Optimistic case: Probabilistic ND) Let p(u) be uniform over all u ∈ U s , then 4 In addition we have that disc Under the uniform distribution, u is unbiased: each direction for the model error is equally likely.This is more optimistic than the assumption of the MMD, where u was biased towards directions that could larger errors on P. Because now u is not biased, E u G(u, M) is smaller under this p(u) than in Theorems 2 and 6 and so this p(u) is more optimistic.The Nuclear Discrepancy (ND) owns its name to the fact that it is proportional to the nuclear matrix norm of M.
An appealing property of this choice of p(u) is that, given a fixed P, any choice of Q does not influence p(u).For the Discrepancy and the MMD, choosing different Q leads to different p(u).Thus choosing queries changes the distribution of p(u) and thus also implicitly the distribution of h and f .Instead, for the ND, our queries don't influence the distribution of h and f .This assumption seems reasonable, since f is usually assumed to be fixed and independent of our actions.
Under the uniform distribution the ND provides the tightest bound on the expected loss, while the MMD bound is looser, and the Discrepancy bound is the loosest.Therefore, if this probabilistic assumption is the most accurate, minimization of the Nuclear Discrepancy may lead to the best queries for active learning, followed by the MMD and Discrepancy, in that order.5

Experiments
We explain the setup and baselines, afterward we review our main results: the realizable setting.We discuss the results and examine the probabilistic assumptions empirically.Somewhat similar results are observed in the agnostic setting which we will briefly discuss.An additional experiment investigates the influence of subsampling of datasets on our results.This subsampling experiment and all results of the agnostic case are discussed in detail in the Appendix.

Experimental setup and baselines
An overview of the experimental procedure is given in Algorithm 1.A training set (65%) and test set (35%) are used-the training set corresponds to P and we indicate the testset by T .We use the active learners to select batches of size n = 1, 2, . . ., 50.For computational reasons we select batches in a sequential greedy fashion.Initially at t = 0 the batch is empty: Q0 = ∅.In iteration 1 ≤ t ≤ n the active learner selects a sample x t from the unlabeled pool Ût−1 = P \ Qt−1 according to x t = arg min s∈ Ût−1 obj( P, Qt−1 ∪ s).We perform experiments multiple times to ensure significance of the results.We call each repetition a run, and for each run a new training and test split is used.During one run, we evaluate each active learner using the described procedure of Algorithm 1.As baseline we use random sampling and a greedy version of the state-of-the-art MMD active learner (Chattopadhyay et al. 2012;Wang and Ye 2013).We compare the baselines with our novel active learners: the Discrepancy active learner and the Nuclear Discrepancy active learner.
The methods are evaluated on 13 datasets that originate either from the UCI Machine Learning repository (Lichman 2013) or were provided by Cawley and Talbot (2004).See Appendix E for the dataset names and characteristics.Furthermore, we perform an experiment on the image dataset MNIST.The MNIST dataset (LeCun et al. 1998) consists of images of handwritten digits of size 28×28 pixels.By treating each pixel as a feature, the dimensionality of this dataset is 784 which is relatively high dimensional.Like Yang and Loog (2018) we construct 3 difficult binary classification problems: 3vs5, 7vs9 and 5vs8.
To make datasets conform to the realizable setting we use the approach of Cortes and Mohri (2014): we fit a model of our hypothesis set to the whole dataset and use its outputs as labels.
To set reasonable hyperparameters we use a similar procedure as Gu et al. (2012).We use labeled data before any experiments are performed to perform model selection to determine hyperparameters (σ and μ of the KRLS model).This can be motivated by the fact that in practice a related task or dataset may be available in order to obtain a rough estimate of the hyperparameter settings.This procedure makes sure η MMD and η disc are small in the agnostic setting.Recall that the active learners minimize bounds on L P (h, f ).Therefore active learners then implicitly also minimizes a bound on L P (h, f ), see Theorem 1.By choosing hyperparameters in the described way above, we ensure that the Rademacher complexity term R m (H ) is not too large and we don't overfit.We measure performance on an independent test set in order to get an unbiased estimate of L P (h, f ).
To aid reproducibility we give all hyperparameters and additional details in Appendix E. We set σ L according to our analysis in Corollary 1.

Realizable setting
First we benchmark the active learners in the realizable setting.In this setting we are assured that η = 0 in all bounds and therefore we eliminate unexpected effects that can arise due to model misspecification.We study this scenario to validate our theoretical results and gain more insight, furthermore, note that this scenario is also studied in adaptation (Cortes and Mohri 2014).
Several learning curves are shown in Fig. 1, all curves can be found in Appendix H.1.The MSE of the active learner minus the mean performance (per query) of random sampling is displayed on the y-axis (lower is better).The curve is averaged over 100 runs.Error bars represent the 95% confidence interval of the mean computed using the standard error.
We summarize results on all datasets using the Area Under the (mean squared error) Learning Curve (AULC) in Table 2.The AULC is a different metric than the well known AUROC or AUPRC measures.The AUROC measure summarize the performance of a model for different misclassification costs (type I and type II costs) and the AUPRC is useful when one class is more important than the other, such as in object detection.
By contrast, AULC is specifically suited to active learning, and summarizes the performance of an active learning algorithm for different number of labeling budgets (O'Neill et al. 2017;Huijser and van Gemert 2017;Settles and Craven 2008).Low AULC is obtained when Significance improvement is judged by a paired two tailed t-test (significance level p = 0.05).We may use a paired test since during one run all active learners are evaluated using the same training and test split.
In the majority of the cases the MMD improves upon the Discrepancy (see Table 2).The results on the ringnorm dataset are remarkable, here the Discrepancy sometimes performs worse than random sampling, see Fig. 1.We observe that generally the Discrepancy performs the worst.These results illustrates that tighter worst case bounds do not guarantee improved performance.The proposed ND active learner significantly improves upon the MMD in 9 out of the 13 datasets tested.Here we counted MNIST once, while we remark that on all subproblems the ND improves significantly on the MMD.This provides evidence that the proposed method can also deal with high-dimensional datasets.In case the ND does not perform the best, it ties with the MMD or Discrepancy.The ND never performs significantly worse.This ranking of the methods exactly corresponds to the order of the bounds given by Theorem 7 under our optimistic probabilistic assumptions.This supports our hypothesis that we find ourselves more often in a more optimistic average-case scenario.

Decomposition of probabilistic bounds
Since we are in the realizable setting we can compute u = h − f with the true labeling function f and our trained model h.Thus we can compute each term in the sum of G(u, M) . ., λ 9 , etc. Averaged over 100 runs of the random active learner.λ 1 in most cases contributes little and in general all λ i contribute to G(u, M).This supports the optimistic probabilistic assumptions in (7) during the experiments. 6We show the contribution of each eigenvalue to G(u, M).In Fig. 2 we show this decomposition using a stacked bar chart during several active learning experiments of the baseline active learner 'Random'. 7Here EV1 indicates the largest absolute eigenvalue, its contribution is given by ū2 1 |λ 1 | (see also ( 7)).EV 2 -9 to indicate the summed contribution: 9 i=2 ū2 i |λ i |, etc.The mean contributions over 100 runs are shown.Observe that the contribution of |λ 1 | to G(u, M) is often small, it is shown by the small white bar at the bottom of the barchart.Therefore the Discrepancy active learner chooses suboptimal samples: its strategy is optimal for a worst-case scenario G(u, M) = 4Λ 2 |λ 1 | that is very rare.We observe that typically all λ i contribute to G(u, M) supporting our probabilistic assumption.

Agnostic setting
For completeness, we briefly mention the agnostic setting, for all details see Appendix F. In the agnostic setting the rankings of methods can change and performance differences become less significant.The ND still improves more upon the MMD than the reverse, however, the trend is less significant.Because our assumption η = 0 is violated our theoretical analysis is less applicable.
For the MNIST experiments we however find that the results for some subproblems almost coincides with the realizeable setting: apparently, for the MNIST dataset the model misspefication is very small.This may be because the dataset is of relatively high dimensionalion.

Influence of subsampling
We briefly mention an additional experiment that we have performed on the splice dataset to see how subsampling affects performance.To this end we measure the performance while we vary the pool size P by changing the amount of subsampling.This to investigate how the proposed methods would perform for problems with a larger scale.For all details please see Appendix G, here we will summarize our findings.
For small pool sizes all active learners experience a drop in performance.We find the larger the pool, the better the performance, up until some point at which the performance levels off.The experiment provides evidence that if finer subsampling is used or larger datasets are used, methods typically improve in performance up to a point where performance levels off.

Discussion
In the experiments we have observed that in the realizable setting the order of the bounds under our more optimistic probabilistic assumptions give the best indication of active learning performance.The empirical decomposition of G(u, M) during experiments also supports our hypothesis that we generally find ourselves in a more optimistic scenario instead of a worst case scenario.
Still it is meaningful to look at worst-case guarantees, though the worst-case should be expected to occur.The worst-case assumed by the Discrepancy can never occur in the realizable setting, and we believe it is also highly unlikely in the agnostic setting.The strength of our probabilistic approach is that it considers all scenarios equally and does not focus too much on specific scenarios, making the strategy more robust.
Our work illustrates that the order of bounds can change under varying conditions and thus tightness of bounds is not the whole story.The conditions under which the bounds hold are equally important, and should reflect the mathematical setting as much as possible.For example, in a different setting where an adversary would pick u, the Discrepancy active learner would be most appropriate.This insight illustrates that not only by obtaining tighter bounds active learning performance can be improved, but by finding more appropriate assumptions (bound-based) active learners can be improved as well.
Our work supports the idea of Germain et al. (2013) who introduce a probabilistic version of the Discrepancy bound for the zero-one loss (Ben-David et al. 2010).Our conclusions also support that the direct Cortes et al. (2019) takes: by using more accurate assumptions to better characterize the the worst case scenario, performance may be improved.
In our study we have focused on minimizing the mean squared error.It would be interesting to investigate the extension of the Nuclear Discrepancy to other loss functions, in particular the zero-one loss.As far as we can see, however, such an extension is not trivial.The above mentioned probabilistic version of the Discrepancy by Germain et al. (2013) may provide some inspiration to achieve this, but they offer a PAC Bayes approach that cannot be easily adapted to the probabilistic setting we consider.
Where the experiments in the realizable setting provide clear insights, the results concerning the agnostic setting are not fully understood.A more in depth experimental study of the agnostic setting is complicated by unexpected effects of η.Since probabilistic bounds are the most informative in the realizable setting, it is of interest to consider probabilistic bounds for the agnostic setting as well.
In our experiments we have used greedy optimization to compute the batch Qn .It is theoretically possible to optimize a whole batch of queries in one global optimization step.However, for the MMD this problem is known to be NP-hard (Chattopadhyay et al. 2012).Minimizing the Discrepancy is also non-trivial, as illustrated by the involved optimization procedure required by Cortes and Mohri (2014) for domain adaptation.Note that their optimization problem is easier than the optimization problem posed by active learning, where binary constraints are necessary.Since the objective value of the Nuclear Discrepancy is given by an expectation which can be approximated using sampling, we believe it may be possible to speed up the optimization by using approximations.
In this work we have only considered single-shot batch active learning.In regular batchmode active learning label information of previously selected samples can be used to improve query selection.This can be accommodated in our active learner by refining p(u) using label information.Our results have implications for adaptation as well.We suspect our suggested choice of σ L may improve the MMD domain adaptation method (Huang et al. 2007).Furthermore, our results suggest that the ND is a promising objective for adaptation.

Conclusion
To investigate the relation between generalization bounds and active learning performance, we gave several theoretical results concerning the bound of the MMD active learner and the Discrepancy bound.In particular, we showed that the Discrepancy provides the tightest worst-case bound.We introduced a novel quantity; Nuclear Discrepancy, motivated from optimistic probabilistic assumptions derived from the principle of maximum entropy.Under these probabilistic assumptions the ND provides the tightest bound on the expected loss, followed by the MMD, and the Discrepancy provides the loosest bound.
Experimentally, we observed that in the realizable setting the Discrepancy performs the worst, illustrating that tighter worst-case bounds do not guarantee improved active learning performance.Our optimistic probabilistic analysis clearly matches the observed behavior in the realizable setting: the proposed ND active learner improves upon the MMD, and the MMD improves upon the Discrepancy active learner.We find that even on the high-dimensional image dataset MNIST our method is competitive.A similar, weaker, trend is observed in the agnostic case.One of our key conclusions is that not only bound tightness is important for active learning performance, but that appropriate assumptions are equally important.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/),which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

A.1 MMD
The MMD quantity can be computed in practice by rewriting it as follows: In the first step we used that l(x) = l, ψ K L (x) K L due to the reproducing property (Mohri et al. 2012, p. 96).Here ψ K L is the featuremap from X → H L .The second step follows from the linearity of the inner product.In (10) we defined μ P = 1 n P x∈ P ψ K L (x) and similarly for μ Q , note that μ Q , μ P ∈ H L .The last step follows from the fact that the vector in H L maximizing the term in ( 10) is Because of the symmetry of ||μ P − μ Q || K L with respect to P and Q, this derivation also holds if we switch P and Q. Therefore: Therefore for all l ∈ H L the following holds We can compute the MMD quantity in practice by working out the norm with kernel products:

A.2 Discrepancy
In this section we calculate the discrepancy analytically for the squared loss in the linear kernel as in Mansour et al. (2009).We then extend the computation to any arbitrary kernel as in Cortes and Mohri (2014).Finally, we prove the agnostic generalization bound in terms of the Discrepancy (Theorem 3).The theorems and proofs here were first given by Mansour et al. (2009), Cortes and Mohri (2014), and Cortes et al. (2019) but we repeat them here for completeness.
Lemma 2 (Mansour et al. 2009) For h, h ∈ H we have Proof We can show using some algebra, where u = h − h .Rewrite L Q (h, h ) similarly and subtract them to find Since M is a real symmetric matrix, M is a normal matrix and admits an orthonormal eigendecomposition with real eigenvalues Here λ i is the ith eigenvalue and e i is the corresponding orthonormal eigenvector.Since M is normal its eigenvectors form an orthonormal basis for R d .Therefore we can express u in terms of e: Where ūi is the projection of u on e i , ūi = e T i u.Note ū is a rotated version of u and therefore both have the same norm, ||u|| 2 = || ū|| 2 .Now we can rewrite ( 14) as Note that M has r = rank(M) non-zero eigenvalues.Combining ( 14) and ( 15) and taking the absolute value on both sides shows the result.Now we are ready to compute the Discrepancy for the linear kernel.
Theorem 8 (Discrepancy computation (Mansour et al. 2009)) Assume K is the linear kernel, K (x i , x j ) = x T i x j , and l is the squared loss, then where λ i are the eigenvalues of M P, Q = M.
Proof First we use Lemma 2.
Now we solve the left term in the maximization.Observe that this is a weighted sum where each ūi weighs each eigenvalue λ i .To maximize this quantity we put as much weight as possible on the largest postive eigenvalue: u = e i max 2Λ, where i max = arg max i λ i .We find max To solve the second maximization, introduce λi = −λ i .Then we maximize the same quantity as before but now λ replaced by λ.It follows that the maximum is attained for u = e i min 2Λ, where i min = arg min i λ i .We find eliminating the maximum proves the result.Now we will describe how to compute the Discrepancy in case we work with an arbitrary kernel K .In this case we have to work in the RKHS H of the kernel K .Define z(x) = ψ K (x), and let Z P be the datamatrix where each row is given by z(x) : x ∈ P. Define Z Q in the analogously.In this case Theorem 8 still holds, and the Discrepancy is given by the eigenvalues of M Z : However, now we run into problems, since for an arbitrary kernel K the dimensions of H can be very large or infinite, such as the case for the Gaussian kernel.Then we clearly cannot compute the matrix M Z or its eigenvalues.
In the following we show that M Z and M K have the same eigenvalues.Then, to compute the Discrepancy with any kernel K , we can simply use the eigenvalues of M K .First, let us define M K .
where K P P is the n P × n P matrix where entry i, j is given by K (x i , x j ), and where D is a diagonal matrix where Lemma 3 (Cortes and Mohri 2014) The eigenvalues of M Z and M K are the same.
Proof Recalling that Q ∈ P, and using some algebra, it can be shown that M Z can be written as Here we used the fact that K (x i , x j ) = ψ K (x i ), ψ K (x j ) K (kernel trick) to rewrite M K .Since the matrix product AB and B A have the same eigenvalues (Cortes and Mohri 2014), M K and M Z have the same eigenvalues.
Theorem 9 Let K be any arbitrary PSD kernel.Then where λ is the vector of eigenvalues of the matrix M K or M Z , where M K was defined in ( 17) and M Z was defined in (16).
Proof First observe that Theorem 8 still holds, but we have to replace M by M Z in case we use any arbitrary PSD kernel K .Then the result follows from Lemma 3.
Finally, we give the proof of the generalization bound in terms of the Discrepancy for the agnostic setting.Note that this proof was already given by Cortes et al. (2019), here we repeat their proof for completeness.

Proof of Theorem 3 Since l(h(x), f (x))
≤ C, we have that the squared loss l is μ-admissible (Cortes et al. 2019) with μ = 2C, meaning that holds for all h, h ∈ H and any f : X → Y. Let f be any arbitrary element from H .By adding and subtracting terms and applying the triangle inequality, we can show that The first term on the right hand side is by definition bounded by the Discrepancy.For the second term we can show The first inequality follows from applying (19) to each summand.We can bound the third term in the same way, since Q ∈ P. Bounding the first term using the Discrepancy and the last two terms with the bound above we find holds for all f ∈ H .The result follows from minimizing the right hand side with respect to f , bounding ) with its absolute value and reordering terms.

B.1 Proof of agnostic MMD worst case bound (Proposition 1)
Proof of Proposition 1 Let l be any element from H L and define g P = 1 n P x∈ P g(h, f )(x) and define g Q similarly.Define l P = 1 n P x∈ P l(x) and l Q analogously.Using the triangle inequality we can show The first term is bounded by the MMD, see (12).For the second term we have |g P − l P | ≤ max x∈ P |g(h, f )(x) − l(x)|.This bound also holds also for the third term since Q ∈ P. Bounding the second and third term and maximizing over h ∈ H we find that holds for any l ∈ H L and any h ∈ H .The result follows by choosing l to minimize the right hand side, bounding L P (h, f ) − L Q (h, f ) by the bound on its absolute value and reordering terms.

B.2 Adjusting the MMD to the loss and hypothesis set (Theorem 2, Corollary 1)
In the main text we have given a sketch of the proof for the linear kernel.Here, we show a rigorous proof for the linear kernel, and afterward we give the proof for any arbitrary kernel K .The technique of the proof stays the same for any arbitrary kernel K , however, we have to do more bookkeeping.
Theorem 10 (Adjusted MMD linear kernel) Let l be the squared loss and assume f ∈ H (realizable setting), furthermore assume K is the linear kernel, K

Proof of Theorem 10
Let u = h − f .Fix h and f , then we will write g(x) as shorthand for g(h, f )(x) = l(h(x), f (x)).Then g(x) = u(x) 2 .Since K is the linear kernel, we have Furthermore, we have H = X , and thus ψ K (x) = x.Furthermore, ψ K L : H → H L is given by (Shawe-Taylor and Cristianini 2004, chap.9.1)8 : Then we can write u).Now what remains to show is that g ∈ H L , or with other words, that ||g|| K L ≤ 4Λ 2 .We can show that where the last step follows from that h, f ∈ H , and therefore ||u|| K ≤ 2Λ.This shows that g(h, f ) ∈ H L , therefore η MMD = 0.
Before we prove the more general case for any kernel K , let us introduce some additional notation.Also, before we show the proof for the kernel K L , we first do the proof for the kernel K which is slightly simpler, later we extend the result to K L .We define the squared kernel K as: Where f ∈ H and g ∈ H, where H is the RKHS of K .We indicate H as the RKHS of K .We assume K is a PSD kernel.By definition of K the kernel K is a PSD kernel since a squared kernel of a PSD kernel is known to be PSD (Mohri et al. 2012, Theorem 5.3).Now we have two kernels we have two featuremaps: ψ K (x) : X → H and ψ K : H → H .Note that the second featuremap can still be computed with (20).See Table 3 for an overview of the notation used.
Recall that because K is PSD kernel we have that: For x, x ∈ X .Similarly for the kernel K which is also PSD we have that: For f , g ∈ H. Again we define u as: Now using ( 23) we can write this as: in other words, we see that the kernel product of K L can be computed in the RKHS of the kernel K .Thus, the RKHS of K and K L coincide!Thus we have that H = H L .Therefore, Theorem 11 implies that we can generalize all results in terms of K to the kernel K L .Therefore, g is also in the RKHS of K L , and in particular we have that g ∈ H L , and therefore η MMD = 0.
Remark 1 Another way to understand this is to see that the featuremap of K L is given by ) and thus maps to the space H , and from this it follows that H = H L .

Proof of Corollary 1
Theorem 2 tells us to choose where we absorbed the factor of 2 in the exponent in σ L , so σ L = σ √ 2 .

B.3 MMD Computation (Theorem 5 and Corollary 2)
First we prove the Theorem 5 in case K is the lineair kernel for d = 2, afterward we extend the proof to any dimension, and finally we prove Theorem 5 for any PSD kernel.
Theorem 12 (MMD Computation linear kernel d = 2) Let K L (x i , x j ) = K (x i , x j ) 2 and Λ L = 4Λ 2 .Furthermore, assume K is the linear kernel, K (x i , x j ) = x T i x j and d = 2, then Proof If K is the linear kernel, H = X and K L defines a featuremap ψ K L (x) : X → H L which is given by ψ , see (Shawe-Taylor and Cristianini 2004, chap.9.1).From (11) we find 12 ) T and some algebra we can show where the second equality follows since M is a real symmetric matrix and thus its eigenvalues are equal to its singular values (up to a sign change).
Let us first generalize Theorem 12 to any arbitrary dimension d, afterward we extend the result to any kernel K .
Theorem 13 (MMD Computation linear kernel arbitrary dimension) If K L (x i , x j ) = K (x i , x j ) 2 and Λ L = 4Λ 2 .Furthermore, assume K is the linear kernel, K (x i , x j ) = x T i x j and the dimension of the input space X is d, then Proof We can show using (20) that the vector μ P − Q becomes Observe that the entry M ii appears only once in μ P − μ Q , and any entry M i j where i = j, appears as √ 2M i j .Furthermore, note that the diagonal M ii only occurs once in the matrix M.However, any element M i j appears twice in M (since M is symmetric).Therefore is a sum of all entries (squared) of M, and therefore as before.The rest of the proof is identical to the proof of Theorem 12 (d = 2).
Note that due to our careful ordering of the featuremap of the squared kernel, given in (20), we have that this featuremap is still properly defined even if the dimension of d → ∞, such as for a Gaussian kernel.Now we are ready to prove Theorem 5.

Proof of Theorem 5
To show the result holds for any arbitrary kernel K , we have to work in the RKHS of K , thus everywhere x needs to be replaced by z(x) as in Appendix A.2, and then we replace M with M Z .Then Theorem 13 still holds, since the featuremap ψ(x) K is still given by the featuremap of the squared kernel, (20), however in this case the featuremap is with respect z instead of x.This does not influence the proof.We showed that this holds for any dimension d, and the featuremap still exists if d → ∞, thus our results hold for a kernel K with arbitrary dimension of the RKHS H.
For the Gaussian kernel we cannot compute the matrix M. Instead, since M and M K have the same eigenvalues, see Lemma 3, we can compute the MMD instead using the eigenvalues of M K , as in Appendix A.2.
These are not equal, since the matrix M K is not symmetric.Therefore, the eigenvalues of M K are not the same as the singular values of M K (as was the case for M, which is symmetric).

B.4 Probabilistic analysis (Lemma 1, Proposition 2 and Theorem 6)
Proof of Lemma 1 We can show that: where the equality follows from Eq. 13 and the inequality follows from the triangle inequality.
Next, bound L P (h, h ) − L Q (h, h ) using the bound above and reorder terms.The result is obtained after applying the expectation w.r.t.u on both sides and applying the linearity of the expectation.

Proof of Proposition 2
Computing G(2Λe , M) which will be found to be exactly equal to the Discrepancy.After combining this fact with Lemma 1 will result in the desired equality, the inequality follows from Corollary 2.

Proof of Theorem 6
We can show that This equality can be shown by working out the expectation and canceling terms and recognizing the definition of the MMD from ( 6).The inequality follows from ||λ|| 2 ≤ √ r ||λ|| ∞ .The final result follows by combining Eq. 27 with Lemma 1.

B.5 Proof of nuclear discrepancy bound (Theorem 7)
Before we can show the proof of the Nuclear Discrepancy bound, we need the following lemma: Lemma 4 Let p(u) be uniform over all u ∈ U .Then Proof By comparing the volume of a sphere of radius 2Λ and the volume of a sphere of radius w = ||u|| 2 , we can show that for this distribution p(u) we have that Then it is straightforward to show that by integration of p(w).From the symmetry of p(u) it follows that E u ū2 1 = E u ū2 i for all i.From this fact and the linearity of the expectation the result follows.

Proof of Theorem 7
We can show that The first equality follows from switching expectation and sum.The second equality follows from symmetry of p(u).The last equality follows from Lemma 4. The bound can be obtained by combining with Theorem 8.The inequalities follow from the vector norm inequalities

C Remark on probabilistic analysis and choice of U s
The remark in this section will explain why instead of U , we need to take U s (to be defined below).The problem stems from the fact that if we choose p(u) uniform on U , it may seem unclear what it means for u to be randomly sampled from an infinite dimensional sphere uniformly.
We will use the notation of Appendix A.2, since we will work a kernel K with a high-dimensional H, in order to highlight the problem that a lot of eigenvalues may be zero.We are analyzing what happens to for arbitrary distributions, for example, for the uniform distribution p(u), in case the length of a vector z(x) is infinite such as with a Gaussian kernel (then M has infinite eigenvalues, but only r are non-zero).
The RKHS of K , H, is then of infinite dimension.We split H in two parts: H s = span(Z P ), and its orthogonal complement H ⊥ s .Then for any vector in a ∈ H, a = a s + a ⊥ s , where a s ∈ H s and a ⊥ s ∈ H ⊥ s .In particular we have that after training a kernel regularized model, we have h ∈ a s due to the regularization term in the training procedure.Furthermore, for any observed f , we have since Z P f ⊥ s = 0.The same thing holds for L Q (h, f ).Therefore, we may consider it redundant to consider f , h, and we may limit our analysis to f s and h s .In addition we have that since u is projected on eigenvectors of M Z , and only eigenvectors in span(Z P ) have non-zero eigenvalue, and thus u ⊥ s only has components that correspond to eigenvalues that are zero.Thus any u has the same objective as the corresponding u s .Therefore, instead of defining a pdf over U , we define a pdf over U s = {u ∈ H s : ||u|| K ≤ 2Λ}.Then by construction the dimension of u is at most r = rank(M Z ) ≤ n P , which is always finite.Then sampling u is a well defined procedure even in infinite dimensional RKHS.

D Computation of the decomposition of the probabilistic bounds
To compute each term of G(u, M Z ), we can compute the eigendecomposition of M K to compute the eigenvalues, however we also need to know ūi for each i.The computation of ūi , the projection of u onto the eigenvector v i of M Z is non-trivial to compute in case kernels are used.Observe that here v i is the ith eigenvector and not a component.Here we assume v i is not normalized to unit norm (which is why we write it differently from e i ).We give a detailed description in this appendix how to compute ūi .In this case the equation for ūi is: The difficulty in this derivation is finding the vector v i in case kernels are used.Then we need to find v i expressed in terms of the datamatrix Z .Then we can apply the 'kernel trick' to compute (31).By the eigenvalue equation of M Z we have: In case of the linear kernel it is straightforward to compute v i .However, to compute v i for any K , we have to take extra steps.First we show that v i can be expressed in terms of the datamatrix Z , and afterward we find this expression of v i in terms of Z .Note that: Here f is given as a linear combination of Z D , we define as the complete datamatrix.This datamatrix includes the training and test set, since f in our experiments was obtained by training on the whole dataset where the original binary labels of the dataset are used (in the realizeable setting).However note that for any f ∈ H the model f can be written in this way.Similarly, since h is trained on the dataset Q, we can write h as a linear combination of objects Z Q .Combining ( 37) and ( 38) with (31) we find that:

E Experimental settings and dataset characteristics
The active learning methods are evaluated on the datasets shown in Table 4.The datasets marked with * were provided by Cawley and Talbot (2004).Other datasets originate from the UCI Machine Learning repository (Lichman 2013), except the MNIST dataset (LeCun et al. 1998) which is a standalone dataset.
The parameter settings used are displayed in Table 5.To obtain these hyperparameters we repeated the following procedure multiple times.We randomly select 25 examples from the dataset and label these.We train a KRLS model on these samples and evaluate the MSE on all unselected objects.The hyperparameters that result in the best performance after averaging are used in the active learning experiments.

F Results of the agnostic setting
For completeness we discuss the results of the agnostic setting where the original binary labels are used.In this setting η = 0, but η will be small due to our choice of hyperparameters, and therefore we ignore it during active learning (since we also cannot estimate it unless we  The curves are less smooth and have larger standard errors compared to the realizable setting.Therefore the active learning methods are harder to distinguish which is reflected in Table 6 by larger standard deviations and more bold numbers in a single row.Observe that the ranking of the methods can also change, see for example the learning curve on ringnorm: in the realizable setting the ND improved upon the Discrepancy, while in the agnostic setting the reverse is the case for large budgets.In this setting, sometimes the Discrepancy performs the best.From Table 6 we can see that the trends observed in the realizable setting are still observed in the agnostic setting: the ND improves more upon the MMD than the reverse, however, the trend is weaker.This is likely the case because for this setting η MMD and η disc are non-zero, and therefore our theoretical analysis is weakened.Finally, observe that for the MNIST dataset, the learning curves and results as summarized by the AULC for 5vs8 are almost completely identical as in the realizeable setting.Similarly, for 3vs5 differences are also quite small.This indicates that MNIST is very close to realizeable with these found hyperparameter settings.

G Influence of subsampling on performance
We perform an additional experiment on the splice dataset to see how subsampling affects performance.To this end we measure the performance while we vary the pool size by changing the amount of subsampling, The subsampled pool is used as training set (this is the pool from which active learners can select queries, P), all remaining samples are used as testset T .Furthermore we use the same experimental protocol as for the other experiments.

Fig. 1
Fig. 1Learning curves for several datasets for the realizeable setting.Results are averaged over 100 runs.The MSE is measured with respect to random sampling (lower is better)

Fig. 2
Fig.2Decomposition of the sum G(u, M) during active learning for several datasets.EV1 indicates the contribution of λ 1 , EV2-9 indicate the summed contributions of λ 2 , . . ., λ 9 , etc. Averaged over 100 runs of the random active learner.λ 1 in most cases contributes little and in general all λ i contribute to G(u, M).This supports the optimistic probabilistic assumptions

Fig. 3
Fig. 3Learning curves for several datasets for the agnostic setting.Results are averaged over 100 runs.Observe that compared to the realizable setting the variability of the performance increases and therefore performance differences become less significant.Due to unexpected effects of η > 0 the ranking of the methods may change

Fig. 5 Fig. 6 Fig. 7
Fig. 5Results on all benchmark datasets for the realizable setting

Table 1
Visual summary of our work

Table 2
Area Under the mean squared error Learning Curve (AULC) for the strategies in the realizable setting, averaged over 100 runs Bold indicates the best result, or results that are not significantly worse than the best result, according to a paired t-test ( p = 0.05).Parenthesis indicate standard deviation an active learner quickly learns a model with low MSE.If a method in the table is bold, it either means it is the best method (as judged by the mean), or if it is not significantly worse than the best method (as judged by the t-test).

Table 3
This table illustrates the notation used when 2 kernels are involved Transformation

Table 6
Area Under the mean squared error Learning Curve (AULC) for the strategies in the agnostic setting, averaged over 100 runs Bold indicates the best result, or results that are not significantly worse than the best result, according to a paired t-test ( p = 0.05).Parenthesis indicate standard deviation have the labels of P).Several illustrating learning curves are shown in Fig.3, all results are summarized in Table6, all learning curves can be found in Appendix H.3.