Covering Numbers for Convex Functions
Adityanand Guntuboyina, Bodhisattva Sen
I Introduction
Ever since the work of , covering numbers (and their logarithms, known as metric entropy numbers) have been studied extensively in a variety of disciplines. For a subset of a metric space , the -covering number is defined as the smallest number of balls of radius whose union contains . Covering numbers capture the size of the underlying metric space and play a central role in a number of areas in information theory and statistics, including nonparametric function estimation, density estimation, empirical processes and machine learning.
In recent years there has been an upsurge of interest in nonparametric function estimation under convexity based constraints, especially in multi-dimension. In general function estimation, it is well-known (see e.g., ) that the covering numbers of the underlying function space can be used to characterize optimal rates of convergence. They are also useful for studying the rates of convergence of empirical minimization procedures (see e.g., ). Our results have direct implications in this regard in the context of understanding the rates of convergence of the numerous convexity constrained function estimators, e.g., the nonparametric least squares estimator of a convex regression function studied in ; the maximum likelihood estimator of a log-concave density in multi-dimension studied in . Also, similar problems that crucially use convexity/concavity constraints to estimate sets have also received recent attention in the statistical and machine learning literature, see e.g., , and our results can be applied in such settings.
The paper is organized as follows. In Section II, we set up notation and provide motivation for our main results, which are proved in Section III. In Section IV, we draw some connections to previous results on covering numbers for convex functions and prove a related auxiliary result along with some inequalities of possible independent interest.
II Motivation
Bronshtein worked with the class where the functions are uniformly Lipschitz with constant . However, in convexity-based function estimation problems, one usually does not have a known uniform Lipschitz bound on the unknown function class. This leads to difficulties in the analysis of empirical minimization procedures via Bronshtein’s result. To the best of our knowledge, there does not exist any other result on the covering numbers of convex functions that deals with all and does not require the Lipschitz constraint.
In the absence of the uniformly Lipschitz constraint (i.e., if one works with the class instead of ), the covering numbers under the metric are infinite. In other words, the space is not totally bounded under the metric. This can be seen, for example, by noting that the functions
are in , for all , and satisfy
This motivated us to study the covering numbers of the class under a different metric, namely the -metric for . We recall that under the -metric, , the distance between two functions and on is defined as
Our main result in this paper shows that if one works with the -metric as opposed to , then the covering numbers of are finite. Moreover, they are bounded from above and below by constant multiples of for sufficiently small .
where .
Fix . There exist positive constants and , depending only on the dimension and , such that, for every and , we have
for every .
The main ingredient in our proof of the above theorem is an extension of Bronshtein’s theorem to uniformly bounded convex functions having different Lipschitz constraints in different directions. Specifically, for , and for , let denote the set of all real-valued convex functions on the rectangle that are uniformly bounded by and satisfy:
for every ; and for . In other words, the function is Lipschitz on with constant for all .
Clearly, the class that Bronshtein studied is contained in . Also, it is easy to check that every function in is Lipschitz with respect to the Euclidean norm on with Lipschitz constant .
Note that for , the inequality (III-A) is satisfied by every function . As a result, we have the equality . The following result gives an upper bound for the -covering number of and is the main ingredient in the proof of Theorem III.1. Its proof is similar to Bronshtein’s proof [3, Proof of Theorem 6] of his upper bound on and is included in Section IV.
There exist positive constants and , depending only on the dimension , such that for every positive and rectangle , we have
for all .
Note that the right hand side of (III.2) equals unless for all . Thus, Theorem III.2 is only meaningful when for all .
Because is contained in , Theorem III.2 includes Bronshtein’s upper bound on as a special case. Moreover, it gives explicit dependence of the upper bound on the constants and . Bronshtein did not state the dependence on these constants.
We are now ready to prove Theorem III.1 using Theorem III.2. Here is the intuition behind the proof. The class can be thought of as an expansion of the class formed by the removal of the Lipschitz constraints (or equivalently, by setting ). Instead of removing all these Lipschitz constraints at the same time, we remove them sequentially one at a time. This is formally accomplished by induction on the number of indices for which . Each step of the induction argument focuses on the removal of one finite and is thus like solving the one-dimensional problem. We consequently use Dryanov’s ideas from [2, Theorem 3.1] to solve this quasi one-dimensional problem which allows us to complete the induction step.
The scaling identity (1) lets us take and .
We shall prove that there exist positive constants and , depending only on and , such that for every , we have
for . Note that this proves the theorem because we can set for all . Our proof will involve induction on : the number of indices for which .
For , i.e., when for all , (III-A) is a direct consequence of Theorem III.2. In fact, in this case, (III-A) also holds for . Suppose now that (III-A) holds for all for some . We shall then verify it for . Fix such that exactly of them equal infinity. Without loss of generality, we assume that and for . For every sufficiently small , we shall exhibit an -cover of in the -metric whose cardinality has logarithm bounded from above by a constant multiple of . Note that for , the term equals zero. For convenience, let us denote the class by in the rest of this proof.
Fix and choose an integer and such that
For every two functions and on , we can obviously decompose the integral as
For a fixed , consider the problem of covering the functions in on the rectangular strip . Clearly,
where, for ,
By convexity, the restriction of every function in to belongs to the class:
Because , we can use the induction hypothesis to assert the existence of positive constants and , depending only on and , such that for every positive real number , there exists an -cover of in the -metric on of size smaller than
By covering the functions in by the constant function 0 on and up to in the -metric on for , we obtain a cover of the restriction of the functions in to the set in -metric having coverage and cardinality bounded from above by where
for , where is the largest integer such that
Note that if , then which implies . Also, for , we have
where we have used and the fact that has the expression (5). Therefore which can be rewritten as
Using this for and , we deduce that
An exactly similar analysis can be done now to cover the restrictions of the functions in to the set having the same coverage and same cardinality bounded by . For , we note, by convexity, that the restrictions of functions in to the set belong to . By the induction hypothesis, there exist constants and , depending only on and , such that for all , one can get a -cover of in the -metric having cardinality smaller than
Observe that only depends on . By combining the covers of the restrictions of functions in to these three strips , and , we obtain, for , a cover of in the -metric having coverage at most
By relabelling as , we have proved that for ,
This proves (III-A) for all such that exactly of them equal . The proof is complete by induction. ∎
The argument used in the induction step above involved splitting the interval $[0,u],[u,v][v,1][0,u]S_{1}S_{2}S_{1}S_{2}$ is much simpler which shortens the argument considerably.
There exist positive constants and , depending only on the dimension , such that for every , and , we have
for .
As before, by the scaling identity (1), we take , and . For functions defined on , the -metric, , is larger than . We will thus take in the rest of this proof. We prove that for sufficiently small, there exists an -packing subset of , under the -metric, of cardinality larger than a constant multiple of . By a packing subset of , we mean a subset satisfying whenever with .
Fix and let be the positive integer satisfying
Consider the intervals for , such that
,
, for ,
for .
Let denote the set of all -dimensional cubes of the form where . The cardinality of , denoted by , is clearly .
where , for . The functions have the following four key properties:
For every , we have .
For every , we have . This is because whenever , we have for each , which implies .
Let with . For every , we have . To see this, let with . Let and fix . If , then and hence
If and , then
The same above bound holds if . Because , at least one of and will be different. Consequently,
Let denote the collection of all -valued functions on . The cardinality of clearly equals (recall that ).
For each , let
The first two properties of ensure that . The last two properties imply that
We now bound from below the distance between and for . Because the interiors of the cubes in are all disjoint, we can write
Note that from (9) and by symmetry, the value of integral
is the same for all . We have thus shown that
where denotes the Hamming distance.
The quantity can be computed in the following way. Let where . We write
By the change of variable for , we get
Recalling that for all , we get where
Note that is a constant that depends on the dimension alone. Thus, from (10), we deduce
for all . We now use the Varshamov-Gilbert lemma (see e.g., [18, Lemma 4.7]) which asserts the existence of a subset of with cardinality, such that for all with . Thus, from (11) and (8), we get that for every with ,
where . Taking , we have obtained for , an -packing subset of of size where
where depends only on the dimension . This completes the proof. ∎
The explicit packing subset constructed in the above proof consists of functions that can be viewed as perturbations of the quadratic function . Previous lower bounds on the covering numbers of convex functions in [3, Proof of Theorem 6] and [2, Section 2] (for ) are based on perturbations of a function whose graph is a subset of a sphere; a more complicated convex function than . The perturbations of in the above proof can also be used to simplify the lower bound arguments in those papers.
IV Distances between convex functions, and their epigraphs
One of the aims of this section is to provide the proof of Theorem III.2. Our strategy for the proof of Theorem III.2 is similar to Bronshtein’s proof of the upper bound on . The proof involves the following ingredients:
An inequality between the distance between two convex functions and the Hausdorff distance between their epigraphs.
The result of Bronshtein for the covering numbers of convex sets in the Hausdorff metric.
For a convex function on and , let us define the epigraph of by
If , then clearly
for every . Therefore, for every , its epigraph is contained in the -dimensional ball of radius centered at the origin. The following inequality relates the distance between two functions in to the Hausdorff distance between their epigraphs. The Hausdorff distance between two compact, convex sets and in Euclidean space is defined by
where denotes Euclidean distance.
For every pair of functions and in , we have
where the second last inequality follows from the Cauchy-Scwarz (C-S) inequality. Lemma IV.1 now follows because is arbitrary in the above argument. ∎
A more detailed account of Bronshtein’s proof of (12) can be found in Section 8.4 of .
The conclusion of the theorem is clearly only meaningful in the case when for all . We therefore assume this in the rest of this proof.
For every , let us define the function on by
for . Clearly the function belongs to the class and covering to within in the -metric is equivalent to covering . Thus
We thus take, without loss of generality, and for all .
From Lemma IV.1 and the observation that for all , it follows that
Thus from (12), we deduce the existence of two positive constants and , depending only on , such that
if . By the scaling inequality (IV), we obtain
if . By another scaling argument, it follows that
for every and, as a consequence, we get, for every ,
if . Choosing (by differentiation)
if . The proof of the theorem will now be complete by noting that
The terms involving can be absorbed in the constants and . ∎
One might wonder if a version of Lemma IV.2 can be proved for the -metric instead of the -metric, and without any Lipschitz constraints. Such an inequality would, in particular, yield an alternative simpler proof of Theorem III.1. It turns out that one can prove such a bound for the -metric but not for for any . The inequality for is presented next. This inequality could possibly be of independent interest. The reason why such an inequality can not be proved for , is explained in Remark IV.1.
For every pair of functions and in , we have
Note that the Cauchy-Schwarz inequality has been used twice in the above chain of inequalities. We have thus shown that in the case when . One would have a similar inequality in the case when . Combining these two, we obtain (15).
where we have used the inequality .
for sufficiently small, where is the unit vector in the th coordinate direction i.e., if and otherwise. Dividing both sides by and letting , we would get (we use to denote the directional derivative of in the direction ; directional derivatives exist as is convex). Using (16) for , we get . Combining these two inequalities, we get
We now show that for each , both the integrals and are bounded from above by 4. Assume, without loss of generality, that and notice
We fix and focus on the inner integral. Let for . Clearly is a convex function on $v_{r}^{\prime}(x_{1})z=x_{1}\in(0,1)f^{\prime}(x;e_{1})x=(x_{1},\dots,x_{d})\int_{\rho}^{1-\rho}|v_{r}^{\prime}(z)|dzvv_{r}^{\prime}(z)$ is non-decreasing and satisfies
The function clearly satisfies because . This implies that . The identity (IV) therefore gives
Similarly, by working with left derivatives of as opposed to right, we can prove that
Therefore, the integral is at most because it is less than or equal to
This completes the proof of Lemma IV.2. ∎
Lemma IV.2 is not true if is replaced by , for . Indeed, if and for and for all , then it can be easily checked that for ,
As can be arbitrarily close to zero, this clearly rules out any inequality of the form (14) with the -metric replaced by , for .
Lemma IV.2 and Bronshtein’s result (12) can be used to give an alternative proof of Theorem III.1 for the special case . Indeed, the scaling identity (1) lets us take , and . Inequality (14) implies that the covering number is less than or equal to
Thus from (12), we deduce the existence of two positive constants and , depending only on , such that
whenever . Note that, by Remark IV.1, this method of proof does not work in the case of , for .