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 F{\cal F} of a metric space (X,ρ)({\mathcal{X}},\rho), the ϵ\epsilon-covering number M(F,ϵ;ρ)M({\cal F},\epsilon;\rho) is defined as the smallest number of balls of radius ϵ\epsilon whose union contains F{\cal F}. 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 C([a,b]d,B,Γ){\mathcal{C}}([a,b]^{d},B,\Gamma) where the functions are uniformly Lipschitz with constant Γ\Gamma. 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 d≥1d\geq 1 and does not require the Lipschitz constraint.

In the absence of the uniformly Lipschitz constraint (i.e., if one works with the class C([a,b]d,B){\mathcal{C}}([a,b]^{d},B) instead of C([a,b]d,B,Γ){\mathcal{C}}([a,b]^{d},B,\Gamma)), the covering numbers under the L∞L_{\infty} metric are infinite. In other words, the space C([a,b]d,B){\mathcal{C}}([a,b]^{d},B) is not totally bounded under the L∞L_{\infty} metric. This can be seen, for example, by noting that the functions

are in C(,1){\mathcal{C}}(,1), for all j≥1j\geq 1, and satisfy

This motivated us to study the covering numbers of the class C([a,b]d,B){\mathcal{C}}([a,b]^{d},B) under a different metric, namely the LpL_{p}-metric for 1≤p<∞1\leq p<\infty. We recall that under the LpL_{p}-metric, 1≤p<∞1\leq p<\infty, the distance between two functions ff and gg on [a,b]d[a,b]^{d} is defined as

Our main result in this paper shows that if one works with the LpL_{p}-metric as opposed to L∞L_{\infty}, then the covering numbers of C([a,b]d,B){\mathcal{C}}([a,b]^{d},B) are finite. Moreover, they are bounded from above and below by constant multiples of ϵ−d/2\epsilon^{-d/2} for sufficiently small ϵ\epsilon.

where ϵ′:=(b−a)−d/pϵ/B\epsilon^{\prime}:=(b-a)^{-d/p}\epsilon/B.

Fix 1≤p<∞1\leq p<\infty. There exist positive constants cc and ϵ0\epsilon_{0}, depending only on the dimension dd and pp, such that, for every B>0B>0 and b>ab>a, we have

for every ϵ≤ϵ0B(b−a)d/p\epsilon\leq\epsilon_{0}B(b-a)^{d/p}.

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 B∈(0,∞)B\in(0,\infty), Γi∈(0,∞]\Gamma_{i}\in(0,\infty] and ai<bia_{i}<b_{i} for i=1,…,di=1,\dots,d, let C(∏i=1d[ai,bi];B;Γ1,…,Γd){\mathcal{C}}\left(\prod_{i=1}^{d}[a_{i},b_{i}];B;\Gamma_{1},\dots,\Gamma_{d}\right) denote the set of all real-valued convex functions ff on the rectangle [a1,b1]×⋯×[ad,bd][a_{1},b_{1}]\times\dots\times[a_{d},b_{d}] that are uniformly bounded by BB and satisfy:

for every i=1,…,di=1,\dots,d; xi,yi∈[ai,bi]x_{i},y_{i}\in[a_{i},b_{i}] and xj∈[aj,bj]x_{j}\in[a_{j},b_{j}] for j≠ij\neq i. In other words, the function x↦f(x1,…,xi−1,x,xi+1,…,xd)x\mapsto f(x_{1},\dots,x_{i-1},x,x_{i+1},\dots,x_{d}) is Lipschitz on [ai,bi][a_{i},b_{i}] with constant Γi\Gamma_{i} for all xj∈[aj,bj],j≠ix_{j}\in[a_{j},b_{j}],j\neq i.

Clearly, the class C([a,b]d,B,Γ){\mathcal{C}}([a,b]^{d},B,\Gamma) that Bronshtein studied is contained in C([a,b]d;B;Γ,…,Γ){\mathcal{C}}([a,b]^{d};B;\Gamma,\dots,\Gamma). Also, it is easy to check that every function ff in C(∏i[ai,bi];B;Γ1,…,Γd){\mathcal{C}}\left(\prod_{i}[a_{i},b_{i}];B;\Gamma_{1},\dots,\Gamma_{d}\right) is Lipschitz with respect to the Euclidean norm on ∏i[ai,bi]\prod_{i}[a_{i},b_{i}] with Lipschitz constant Γ12+⋯+Γd2\sqrt{\Gamma_{1}^{2}+\dots+\Gamma_{d}^{2}}.

Note that for Γi=∞\Gamma_{i}=\infty, the inequality (III-A) is satisfied by every function ff. As a result, we have the equality C([a,b]d,B)=C([a,b]d;B;∞,…,∞){\mathcal{C}}([a,b]^{d},B)={\mathcal{C}}([a,b]^{d};B;\infty,\dots,\infty). The following result gives an upper bound for the ϵ\epsilon-covering number of C(∏i[ai,bi];B;Γ1,…,Γd){\mathcal{C}}(\prod_{i}[a_{i},b_{i}];B;\Gamma_{1},\dots,\Gamma_{d}) 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 C([a,b]d,B,Γ){\mathcal{C}}([a,b]^{d},B,\Gamma) and is included in Section IV.

There exist positive constants cc and ϵ0\epsilon_{0}, depending only on the dimension dd, such that for every positive B,Γ1,…,ΓdB,\Gamma_{1},\dots,\Gamma_{d} and rectangle [a1,b1]×⋯×[ad,bd][a_{1},b_{1}]\times\dots\times[a_{d},b_{d}], we have

for all 0<ϵ≤ϵ0{B+∑i=1dΓi(bi−ai)}0<\epsilon\leq\epsilon_{0}\{B+\sum_{i=1}^{d}\Gamma_{i}(b_{i}-a_{i})\}.

Note that the right hand side of (III.2) equals ∞\infty unless Γi<∞\Gamma_{i}<\infty for all i=1,…,di=1,\dots,d. Thus, Theorem III.2 is only meaningful when Γi<∞\Gamma_{i}<\infty for all i=1,…,di=1,\dots,d.

Because C([a,b]d,B,Γ){\mathcal{C}}([a,b]^{d},B,\Gamma) is contained in C([a,b]d;B;Γ1,…,Γd){\mathcal{C}}([a,b]^{d};B;\Gamma_{1},\dots,\Gamma_{d}), Theorem III.2 includes Bronshtein’s upper bound on C([a,b]d,B,Γ){\mathcal{C}}([a,b]^{d},B,\Gamma) as a special case. Moreover, it gives explicit dependence of the upper bound on the constants a,b,Ba,b,B and Γ\Gamma. 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 C([a,b]d,B){\mathcal{C}}([a,b]^{d},B) can be thought of as an expansion of the class C([a,b]d;B;Γ1,…,Γd){\mathcal{C}}([a,b]^{d};B;\Gamma_{1},\dots,\Gamma_{d}) formed by the removal of the dd Lipschitz constraints Γ1,…,Γd\Gamma_{1},\dots,\Gamma_{d} (or equivalently, by setting Γ1=⋯=Γd=∞\Gamma_{1}=\dots=\Gamma_{d}=\infty). Instead of removing all these dd 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 ii for which Γi=∞\Gamma_{i}=\infty. Each step of the induction argument focuses on the removal of one finite Γi\Gamma_{i} 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 a=0,b=1a=0,b=1 and B=1B=1.

We shall prove that there exist positive constants cc and ϵ0\epsilon_{0}, depending only on dd and pp, such that for every Γi∈(0,∞]\Gamma_{i}\in(0,\infty], we have

for 0<ϵ≤ϵ00<\epsilon\leq\epsilon_{0}. Note that this proves the theorem because we can set Γi=∞\Gamma_{i}=\infty for all i=1,…,di=1,\dots,d. Our proof will involve induction on ll: the number of indices ii for which Γi=∞\Gamma_{i}=\infty.

For l=0l=0, i.e., when Γi<∞\Gamma_{i}<\infty for all i=1,…,di=1,\dots,d, (III-A) is a direct consequence of Theorem III.2. In fact, in this case, (III-A) also holds for p=∞p=\infty. Suppose now that (III-A) holds for all l<kl<k for some k∈{1,…,d}k\in\{1,\dots,d\}. We shall then verify it for l=kl=k. Fix Γi∈(0,∞]\Gamma_{i}\in(0,\infty] such that exactly kk of them equal infinity. Without loss of generality, we assume that Γ1=⋯=Γk=∞\Gamma_{1}=\dots=\Gamma_{k}=\infty and Γi<∞\Gamma_{i}<\infty for i>ki>k. For every sufficiently small ϵ>0\epsilon>0, we shall exhibit an ϵ\epsilon-cover of C(d;1;∞,…,∞,Γk+1,…,Γd){\mathcal{C}}(^{d};1;\infty,\dots,\infty,\Gamma_{k+1},\dots,\Gamma_{d}) in the LpL_{p}-metric whose cardinality has logarithm bounded from above by a constant multiple of (∑i>kΓi+2)d/2ϵ−d/2(\sum_{i>k}\Gamma_{i}+2)^{d/2}\epsilon^{-d/2}. Note that for k=dk=d, the term ∑i>kΓi\sum_{i>k}\Gamma_{i} equals zero. For convenience, let us denote the class C(d;1;∞,…,∞,Γk+1,…,Γd){\mathcal{C}}(^{d};1;\infty,\dots,\infty,\Gamma_{k+1},\dots,\Gamma_{d}) by G{\mathcal{G}} in the rest of this proof.

Fix η>0\eta>0 and choose an integer AA and δ1,…,δA+1\delta_{1},\dots,\delta_{A+1} such that

For every two functions ff and gg on d^{d}, we can obviously decompose the integral ∫∣f−g∣p\int|f-g|^{p} as

For a fixed m=1,…,Am=1,\dots,A, consider the problem of covering the functions in G{\mathcal{G}} on the rectangular strip [δm,δm+1]×d−1[\delta_{m},\delta_{m+1}]\times^{d-1}. Clearly,

where, for x=(x1,…,xd)∈dx=(x_{1},\dots,x_{d})\in^{d},

By convexity, the restriction of every function ff in G{\mathcal{G}} to [δm,δm+1]×d−1[\delta_{m},\delta_{m+1}]\times^{d-1} belongs to the class:

Because 2(δm+1−δm)/δm<∞2(\delta_{m+1}-\delta_{m})/\delta_{m}<\infty, we can use the induction hypothesis to assert the existence of positive constants ϵ0\epsilon_{0} and cc, depending only on dd and pp, such that for every positive real number αm≤ϵ0\alpha_{m}\leq\epsilon_{0}, there exists an αm\alpha_{m}-cover of C(d;1;2(δm+1−δm)/δm,∞,…,∞,Γk+1,…,Γd)){\mathcal{C}}(^{d};1;2(\delta_{m+1}-\delta_{m})/\delta_{m},\infty,\dots,\infty,\Gamma_{k+1},\ldots,\Gamma_{d})) in the LpL_{p}-metric on d^{d} of size smaller than

By covering the functions in G{\mathcal{G}} by the constant function 0 on [0,δ1]×d−1[0,\delta_{1}]\times^{d-1} and up to αm\alpha_{m} in the LpL_{p}-metric on [δm,δm+1]×d−1[\delta_{m},\delta_{m+1}]\times^{d-1} for m=1,…,Am=1,\dots,A, we obtain a cover of the restriction of the functions in G{\mathcal{G}} to the set [0,u]×d−1[0,u]\times^{d-1} in LpL_{p}-metric having coverage S11/pS_{1}^{1/p} and cardinality bounded from above by exp⁡(S2)\exp(S_{2}) where

for m=1,…,A+1m=1,\dots,A+1, where AA is the largest integer such that

Note that if η≤1\eta\leq 1, then log⁡η≤0\log\eta\leq 0 which implies ζm≤1\zeta_{m}\leq 1. Also, for m=2,…,Am=2,\dots,A, we have

where we have used δA<u\delta_{A}<u and the fact that uu has the expression (5). Therefore ζm≥2ζm−1\zeta_{m}\geq 2\zeta_{m-1} which can be rewritten as

Using this for r=2r=2 and r=dr=d, we deduce that

An exactly similar analysis can be done now to cover the restrictions of the functions in G{\mathcal{G}} to the set [v,1]×d−1[v,1]\times^{d-1} having the same coverage S11/pS_{1}^{1/p} and same cardinality bounded by exp⁡(S2)\exp(S_{2}). For [u,v]×d−1[u,v]\times^{d-1}, we note, by convexity, that the restrictions of functions in G{\mathcal{G}} to the set [u,v]×d−1[u,v]\times^{d-1} belong to C([u,v]×d−1;1;2/u,∞,…,∞,Γk+1,…,Γd){\mathcal{C}}([u,v]\times^{d-1};1;2/u,\infty,\dots,\infty,\Gamma_{k+1},\dots,\Gamma_{d}). By the induction hypothesis, there exist constants cc and ϵ0\epsilon_{0}, depending only on dd and pp, such that for all η≤ϵ0\eta\leq\epsilon_{0}, one can get a ϵ\epsilon-cover of C([u,v]×d−1;1;2/u,∞,…,∞,Γk+1,…,Γd){\mathcal{C}}([u,v]\times^{d-1};1;2/u,\infty,\dots,\infty,\Gamma_{k+1},\dots,\Gamma_{d}) in the LpL_{p}-metric having cardinality smaller than

Observe that uu only depends on pp. By combining the covers of the restrictions of functions in G{\mathcal{G}} to these three strips [0,u]×d−1[0,u]\times^{d-1}, [u,v]×d−1[u,v]\times^{d-1} and [v,1]×d−1[v,1]\times^{d-1}, we obtain, for η≤ϵ0\eta\leq\epsilon_{0}, a cover of G{\mathcal{G}} in the LpL_{p}-metric having coverage at most

By relabelling (17/3)1/pη(17/3)^{1/p}\eta as ϵ\epsilon, we have proved that for ϵ≤(3/17)1/pϵ0\epsilon\leq(3/17)^{1/p}\epsilon_{0},

This proves (III-A) for all Γ1,…,Γd\Gamma_{1},\dots,\Gamma_{d} such that exactly kk of them equal ∞\infty. The proof is complete by induction. ∎

The argument used in the induction step above involved splitting the interval $intothethreeintervalsinto the three intervals[0,u],[u,v]andand[v,1],andthensubsequentlysplittingtheinterval, and then subsequently splitting the interval[0,u]intosmallersubintervals.WehaveborrowedthisideafromDryanov[2,ProofofTheorem3.1].WemustmentionhoweverthatDryanovusesamoreelaborateargumenttoboundsumsoftheforminto smaller subintervals. We have borrowed this idea from Dryanov [2, Proof of Theorem 3.1]. We must mention however that Dryanov uses a more elaborate argument to bound sums of the formS_{1}andandS_{2}.Ourwayofcontrolling. Our way of controllingS_{1}andandS_{2}$ is much simpler which shortens the argument considerably.

There exist positive constants cc and ϵ0\epsilon_{0}, depending only on the dimension dd, such that for every p≥1p\geq 1, B>0B>0 and b>ab>a, we have

for ϵ≤ϵ0B(b−a)d/p\epsilon\leq\epsilon_{0}B(b-a)^{d/p}.

As before, by the scaling identity (1), we take a=0a=0, b=1b=1 and B=1B=1. For functions defined on d^{d}, the LpL_{p}-metric, p>1p>1, is larger than L1L_{1}. We will thus take p=1p=1 in the rest of this proof. We prove that for ϵ\epsilon sufficiently small, there exists an ϵ\epsilon-packing subset of C(d,1){\mathcal{C}}(^{d},1), under the L1L_{1}-metric, of cardinality larger than a constant multiple of ϵ−d/2\epsilon^{-d/2}. By a packing subset of C(d,1){\mathcal{C}}(^{d},1), we mean a subset FF satisfying ∣∣f−g∣∣1≥ϵ||f-g||_{1}\geq\epsilon whenever f,g∈Ff,g\in F with f≠gf\neq g.

Fix 0<η≤4(2+d−1)−20<\eta\leq 4(2+\sqrt{d-1})^{-2} and let k:=k(η)k:=k(\eta) be the positive integer satisfying

Consider the intervals I(i)=[u(i),v(i)]I(i)=[u(i),v(i)] for i=1,…,ki=1,\dots,k, such that

0≤u(1)<v(1)≤u(2)<v(2)≤⋯≤u(k)<v(k)≤10\leq u(1)<v(1)\leq u(2)<v(2)\leq\dots\leq u(k)<v(k)\leq 1,

v(i)−u(i)=ηv(i)-u(i)=\sqrt{\eta}, for i=1,…,ki=1,\dots,k,

u(i+1)−v(i)=12η(d−1)u(i+1)-v(i)=\frac{1}{2}\sqrt{\eta(d-1)} for i=1,…,k−1i=1,\dots,k-1.

Let S{\mathcal{S}} denote the set of all dd-dimensional cubes of the form I(i1)×⋯×I(id)I(i_{1})\times\dots\times I(i_{d}) where i1,…,id∈{1,…,k}i_{1},\dots,i_{d}\in\{1,\dots,k\}. The cardinality of S{\mathcal{S}}, denoted by ∣S∣|{\mathcal{S}}|, is clearly kdk^{d}.

where f0(x):=1d(x12+⋯+xd2)f_{0}(x):=\frac{1}{d}\left(x_{1}^{2}+\dots+x_{d}^{2}\right), for x∈dx\in^{d}. The functions hS,S∈Sh_{S},S\in{\mathcal{S}} have the following four key properties:

For every x∈dx\in^{d}, we have hS(x)≤hS(1,…,1)≤1h_{S}(x)\leq h_{S}(1,\dots,1)\leq 1.

For every x∈Sx\in S, we have hS(x)≥f0(x)h_{S}(x)\geq f_{0}(x). This is because whenever x∈Sx\in S, we have u(ij)≤xj≤v(ij)u(i_{j})\leq x_{j}\leq v(i_{j}) for each jj, which implies {xj−u(ij)}{v(ij)−xj}≥0\{x_{j}-u(i_{j})\}\{v(i_{j})-x_{j}\}\geq 0.

Let S,S′∈SS,S^{\prime}\in{\mathcal{S}} with S≠S′S\neq S^{\prime}. For every x∈S′x\in S^{\prime}, we have hS(x)≤f0(x)h_{S}(x)\leq f_{0}(x). To see this, let S′=I(i1′)×⋯×I(id′)S^{\prime}=I(i^{\prime}_{1})\times\dots\times I(i^{\prime}_{d}) with I(ij′)=[u(ij′),v(ij′)]I(i^{\prime}_{j})=[u(i^{\prime}_{j}),v(i^{\prime}_{j})]. Let x∈S′x\in S^{\prime} and fix 1≤j≤d1\leq j\leq d. If I(ij)=I(ij′)I(i_{j})=I(i_{j}^{\prime}), then xj∈I(ij)=[u(ij),v(ij)]x_{j}\in I(i_{j})=[u(i_{j}),v(i_{j})] and hence

If I(ij)≠I(ij′)I(i_{j})\neq I(i_{j}^{\prime}) and u(ij′)<v(ij′)<u(ij)<v(ij)u(i^{\prime}_{j})<v(i^{\prime}_{j})<u(i_{j})<v(i_{j}), then

The same above bound holds if u(ij)<v(ij)<u(ij′)<v(ij′)u(i_{j})<v(i_{j})<u(i^{\prime}_{j})<v(i^{\prime}_{j}). Because S≠S′S\neq S^{\prime}, at least one of iji_{j} and ij′i^{\prime}_{j} will be different. Consequently,

Let {0,1}S\{0,1\}^{{\mathcal{S}}} denote the collection of all {0,1}\{0,1\}-valued functions on S{\mathcal{S}}. The cardinality of {0,1}S\{0,1\}^{{\mathcal{S}}} clearly equals 2∣S∣2^{|{\mathcal{S}}|} (recall that ∣S∣=kd|{\mathcal{S}}|=k^{d}).

For each θ∈{0,1}S\theta\in\{0,1\}^{{\mathcal{S}}}, let

The first two properties of hS,S∈Sh_{S},S\in{\mathcal{S}} ensure that gθ∈C(d,1)g_{\theta}\in{\mathcal{C}}(^{d},1). The last two properties imply that

We now bound from below the L1L_{1} distance between gθg_{\theta} and gθ′g_{\theta^{\prime}} for θ,θ∈{0,1}S\theta,\theta\in\{0,1\}^{{\mathcal{S}}}. Because the interiors of the cubes in S{\mathcal{S}} are all disjoint, we can write

Note that from (9) and by symmetry, the value of integral

is the same for all S∈SS\in{\mathcal{S}}. We have thus shown that

where Υ(θ,θ′):=∑S∈S{θ(S)≠θ′(S)}{\Upsilon}(\theta,\theta^{\prime}):=\sum_{S\in{\mathcal{S}}}\left\{\theta(S)\neq\theta^{\prime}(S)\right\} denotes the Hamming distance.

The quantity ζ\zeta can be computed in the following way. Let S=I(i1)×⋯×I(id)S=I(i_{1})\times\dots\times I(i_{d}) where I(ij)=[u(ij),v(ij)]I(i_{j})=[u(i_{j}),v(i_{j})]. We write

By the change of variable yj={xj−u(ij)}/{v(ij)−u(ij)}y_{j}=\{x_{j}-u(i_{j})\}/\{v(i_{j})-u(i_{j})\} for j=1,…,dj=1,\dots,d, we get

Recalling that v(i)−u(i)=ηv(i)-u(i)=\sqrt{\eta} for all i=1,…,ki=1,\dots,k, we get ζ=ηd/2ηγd\zeta=\eta^{d/2}\eta\gamma_{d} where

Note that γd\gamma_{d} is a constant that depends on the dimension dd alone. Thus, from (10), we deduce

for all θ,θ′∈{0,1}S\theta,\theta^{\prime}\in\{0,1\}^{{\mathcal{S}}}. We now use the Varshamov-Gilbert lemma (see e.g., [18, Lemma 4.7]) which asserts the existence of a subset WW of {0,1}S\{0,1\}^{{\mathcal{S}}} with cardinality, ∣W∣≥exp⁡(∣S∣/8)|W|\geq\exp(|{\mathcal{S}}|/8) such that Υ(τ,τ′)≥∣S∣/4{\Upsilon}(\tau,\tau^{\prime})\geq|{\mathcal{S}}|/4 for all τ,τ′∈W\tau,\tau^{\prime}\in W with τ≠τ′\tau\neq\tau^{\prime}. Thus, from (11) and (8), we get that for every τ,τ′∈W\tau,\tau^{\prime}\in W with τ≠τ′\tau\neq\tau^{\prime},

where c1:=γd4(2+d−1)−dc_{1}:=\frac{\gamma_{d}}{4}(2+\sqrt{d-1})^{-d}. Taking ϵ:=c1η\epsilon:=c_{1}\eta, we have obtained for ϵ≤ϵ0:=4c1(2+d−1)−2\epsilon\leq\epsilon_{0}:=4c_{1}(2+\sqrt{d-1})^{-2}, an ϵ\epsilon-packing subset of C(d,1){\mathcal{C}}(^{d},1) of size M:=∣W∣M:=|W| where

where cc depends only on the dimension dd. 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 f0f_{0}. Previous lower bounds on the covering numbers of convex functions in [3, Proof of Theorem 6] and [2, Section 2] (for d=1d=1) are based on perturbations of a function whose graph is a subset of a sphere; a more complicated convex function than f0f_{0}. The perturbations of f0f_{0} 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 M(C([a,b]d,B,Γ),ϵ;L∞)M({\mathcal{C}}([a,b]^{d},B,\Gamma),\epsilon;L_{\infty}). The proof involves the following ingredients:

An inequality between the L∞L_{\infty} 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 ff on d^{d} and B>0B>0, let us define the epigraph Vf(B)V_{f}(B) of ff by

If f∈C(d,B)f\in{\mathcal{C}}(^{d},B), then clearly

for every (x1,…,xd+1)∈Vf(B)(x_{1},\dots,x_{d+1})\in V_{f}(B). Therefore, for every f∈C(d,B)f\in{\mathcal{C}}(^{d},B), its epigraph Vf(B)V_{f}(B) is contained in the (d+1)(d+1)-dimensional ball of radius d+B2\sqrt{d+B^{2}} centered at the origin. The following inequality relates the L∞L_{\infty} distance between two functions in C(d;B;Γ1,…,Γd){\mathcal{C}}(^{d};B;\Gamma_{1},\dots,\Gamma_{d}) to the Hausdorff distance between their epigraphs. The Hausdorff distance between two compact, convex sets CC and DD in Euclidean space is defined by

where ∣⋅∣|\cdot| denotes Euclidean distance.

For every pair of functions ff and gg in C(d;B;Γ1,…,Γd){\mathcal{C}}(^{d};B;\Gamma_{1},\dots,\Gamma_{d}), we have

where the second last inequality follows from the Cauchy-Scwarz (C-S) inequality. Lemma IV.1 now follows because x∈dx\in^{d} 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 Γi<∞\Gamma_{i}<\infty for all i=1,…,di=1,\dots,d. We therefore assume this in the rest of this proof.

For every f∈C(∏i=1d[ai,bi];B;Γ1,…,Γd)f\in{\mathcal{C}}\left(\prod_{i=1}^{d}[a_{i},b_{i}];B;\Gamma_{1},\dots,\Gamma_{d}\right), let us define the function f^\hat{f} on d^{d} by

for t1,t2,…,td∈t_{1},t_{2},\ldots,t_{d}\in. Clearly the function f^\hat{f} belongs to the class C(d;B;Γ1(b1−a1),…,Γd(bd−ad)){\mathcal{C}}\left(^{d};B;\Gamma_{1}(b_{1}-a_{1}),\dots,\Gamma_{d}(b_{d}-a_{d})\right) and covering f^\hat{f} to within ϵ\epsilon in the L∞L_{\infty}-metric is equivalent to covering ff. Thus

We thus take, without loss of generality, ai=0a_{i}=0 and bi=1b_{i}=1 for all i=1,…,di=1,\dots,d.

From Lemma IV.1 and the observation that Vf(B)∈Kd+1(d+B2)V_{f}(B)\in{\mathcal{K}}^{d+1}(\sqrt{d+B^{2}}) for all f∈C(d,B)f\in{\mathcal{C}}(^{d},B), it follows that

Thus from (12), we deduce the existence of two positive constants cc and ϵ0\epsilon_{0}, depending only on dd, such that

if ϵ≤ϵ0(d+B2)(1+Γ12+⋯+Γd2)\epsilon\leq\epsilon_{0}\sqrt{(d+B^{2})(1+\Gamma_{1}^{2}+\dots+\Gamma_{d}^{2})}. By the scaling inequality (IV), we obtain

if ϵ≤ϵ0(d+B2)(1+∑iΓi2(bi−ai)2)\epsilon\leq\epsilon_{0}\sqrt{(d+B^{2})(1+\sum_{i}\Gamma_{i}^{2}(b_{i}-a_{i})^{2})}. By another scaling argument, it follows that

for every Υ>0\Upsilon>0 and, as a consequence, we get, for every Υ>0\Upsilon>0,

if ϵ≤ϵ0(dΥ2+B2)(1+∑iΓi2(bi−ai)2/Υ2)\epsilon\leq\epsilon_{0}\sqrt{(d\Upsilon^{2}+B^{2})(1+\sum_{i}\Gamma_{i}^{2}(b_{i}-a_{i})^{2}/\Upsilon^{2})}. Choosing (by differentiation)

if ϵ≤ϵ0(B+d∑iΓi2(bi−ai)2)\epsilon\leq\epsilon_{0}\left(B+\sqrt{d\sum_{i}\Gamma_{i}^{2}(b_{i}-a_{i})^{2}}\right). The proof of the theorem will now be complete by noting that

The terms involving dd can be absorbed in the constants cc and ϵ0\epsilon_{0}. ∎

One might wonder if a version of Lemma IV.2 can be proved for the LpL_{p}-metric instead of the L∞L_{\infty}-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 L1L_{1}-metric but not for LpL_{p} for any p>1p>1. The inequality for L1L_{1} is presented next. This inequality could possibly be of independent interest. The reason why such an inequality can not be proved for Lp,p>1L_{p},p>1, is explained in Remark IV.1.

For every pair of functions ff and gg in C(d,1){\mathcal{C}}(^{d},1), we have

Note that the Cauchy-Schwarz inequality has been used twice in the above chain of inequalities. We have thus shown that g(x)−f(x)≤ρ(1+∣mg(x)∣)g(x)-f(x)\leq\rho(1+|m_{g}(x)|) in the case when f(x)<g(x)f(x)<g(x). One would have a similar inequality in the case when f(x)>g(x)f(x)>g(x). Combining these two, we obtain (15).

where we have used the inequality (1−2ρ)d≥1−2dρ(1-2\rho)^{d}\geq 1-2d\rho.

for t>0t>0 sufficiently small, where eie_{i} is the unit vector in the iith coordinate direction i.e., ei(j):=1e_{i}(j):=1 if i=ji=j and otherwise. Dividing both sides by tt and letting t↓0t\downarrow 0, we would get mf(x)(i)≤f′(x;ei)m_{f}(x)(i)\leq f^{\prime}(x;e_{i}) (we use f′(x;v)f^{\prime}(x;v) to denote the directional derivative of ff in the direction vv; directional derivatives exist as ff is convex). Using (16) for t<0t<0, we get mf(x)(i)≥−f′(x;−ei)m_{f}(x)(i)\geq-f^{\prime}(x;-e_{i}). Combining these two inequalities, we get

We now show that for each ii, both the integrals ∫[ρ,1−ρ]d∣f′(x;ei)∣\int_{[\rho,1-\rho]^{d}}|f^{\prime}(x;e_{i})| and ∫[ρ,1−ρ]d∣f′(x;−ei)∣\int_{[\rho,1-\rho]^{d}}|f^{\prime}(x;-e_{i})| are bounded from above by 4. Assume, without loss of generality, that i=1i=1 and notice

We fix u=(x2,…,xd)∈[ρ,1−ρ]d−1u=(x_{2},\dots,x_{d})\in[\rho,1-\rho]^{d-1} and focus on the inner integral. Let v(z):=f(z,x2,…,xd)v(z):=f(z,x_{2},\dots,x_{d}) for z∈z\in. Clearly vv is a convex function on $anditsrightderivative,and its right derivative,v_{r}^{\prime}(x_{1})atthepointat the pointz=x_{1}\in(0,1)equalsequalsf^{\prime}(x;e_{1})wherewherex=(x_{1},\dots,x_{d}).Theinnerintegralthusequals. The inner integral thus equals\int_{\rho}^{1-\rho}|v_{r}^{\prime}(z)|dz.Becauseoftheconvexityof. Because of the convexity ofv,itsrightderivative, its right derivativev_{r}^{\prime}(z)$ is non-decreasing and satisfies

The function v(⋅)v(\cdot) clearly satisfies ∣v(z)∣≤1|v(z)|\leq 1 because f∈C(d,1)f\in{\mathcal{C}}(^{d},1). This implies that ∫ρ1−ρ∣vr′(z)∣dz≤4\int_{\rho}^{1-\rho}|v_{r}^{\prime}(z)|dz\leq 4. The identity (IV) therefore gives

Similarly, by working with left derivatives of vv as opposed to right, we can prove that

Therefore, the integral ∫[ρ,1−ρ]d∣mf∣\int_{[\rho,1-\rho]^{d}}|m_{f}| is at most 8d8d because it is less than or equal to

This completes the proof of Lemma IV.2. ∎

Lemma IV.2 is not true if L1L_{1} is replaced by LpL_{p}, for p>1p>1. Indeed, if d=1d=1 and fα(x):=max⁡(0,1−(x/α))f_{\alpha}(x):=\max(0,1-(x/\alpha)) for 0<α≤10<\alpha\leq 1 and g(x):=0g(x):=0 for all x∈x\in, then it can be easily checked that for 1≤p<∞1\leq p<\infty,

As α\alpha can be arbitrarily close to zero, this clearly rules out any inequality of the form (14) with the L1L_{1}-metric replaced by LpL_{p}, for 1<p≤∞1<p\leq\infty.

Lemma IV.2 and Bronshtein’s result (12) can be used to give an alternative proof of Theorem III.1 for the special case p=1p=1. Indeed, the scaling identity (1) lets us take a=0a=0, b=1b=1 and B=1B=1. Inequality (14) implies that the covering number M(C(d,1),ϵ;L1)M\left({\mathcal{C}}(^{d},1),\epsilon;L_{1}\right) is less than or equal to

Thus from (12), we deduce the existence of two positive constants cc and ϵ0\epsilon_{0}, depending only on dd, such that

whenever ϵ≤ϵ0\epsilon\leq\epsilon_{0}. Note that, by Remark IV.1, this method of proof does not work in the case of LpL_{p}, for 1<p<∞1<p<\infty.

References