Upper bounds on product and multiplier empirical processes

Shahar Mendelson

Introduction

Empirical processes appear frequently in diverse branches of Mathematics, Statistics and Computer Science.

In its most standard form, an empirical process is indexed by a class of functions defined on a probability space (Ω,μ)(\Omega,\mu). If X1,...,XNX_{1},...,X_{N} are independent and distributed according to μ\mu, the centred empirical process indexed by FF is

One would like to obtain upper and lower bounds on

either in high probability or in expectation. The hope is that the supremum (1.1) may be controlled using some geometric features of the class FF, similar to the ones used in the theory of gaussian processes (for more information on gaussian processes and their connection to the geometry of the underlying class see the books and ).

Let F⊂L2F\subset L_{2} and set {Gf:f∈F}\{G_{f}:f\in F\} to be the centred canonical gaussian process indexed by FF; that is, the gaussian process indexed by FF whose covariance structure is endowed by the inner product in L2L_{2}.

To avoid the (well understood) issue of measurability, set

(1) Let ξ∈Lq\xi\in L_{q} for some q>2q>2 (ξ\xi need not be independent of XX) and set ξ1,...,ξN\xi_{1},...,\xi_{N} to be independent copies of ξ\xi. Consider

which is the supremum of the multiplier process indexed by FF and associated with the multiplier ξ\xi.

Note that unlike the standard notion of a multiplier process, here (ξi)i=1N(\xi_{i})_{i=1}^{N} need not be independent of (Xi)i=1N(X_{i})_{i=1}^{N}.

Besides being a natural object from the theoretical point of view, the significance of multiplier processes may be seen in numerous applications. For example, multiplier processes play a central role in Statistics, when studying prediction and estimation problems (see, e.g. and references therein, and also ), but that is only the tip of the iceberg as far as applications go.

(2) Let FF and HH be classes of functions defined on the probability space (Ω,μ)(\Omega,\mu) and consider the supremum of the product process indexed by FF and HH,

Clearly, (1.3) is a natural object when trying to analyze, for example, empirical correlation, or, what is arguably the most important process as far as applications are concerned, the quadratic empirical process

Although the product process may be viewed as the standard empirical process indexed by the product class F⋅H={f⋅h:f∈F,h∈H}F\cdot H=\{f\cdot h:f\in F,h\in H\}, and the multiplier process as the standard empirical process indexed by the class ξ⋅F={ξf:f∈F}\xi\cdot F=\{\xi f:f\in F\} (at least when (ξi,Xi)i=1N(\xi_{i},X_{i})_{i=1}^{N} are independent), the type of result one is looking for here is rather different. When studying the two processes, one would like to bound the supremum of (1.2) and of (1.3) using some geometric structures of the indexing classes FF and F,HF,H respectively, rather than the structures of ξ⋅F\xi\cdot F and F⋅HF\cdot H, which, in most cases, are hard to handle.

Before we continue exploring the two problems, let us explain what is meant by “some geometric structures of the indexing classes”.

Motivated by chaining methods that have had tremendous impact on the theory of gaussian processes, the geometric parameters we shall focus on are ‘relatives’ of Talagrand’s γ\gamma functionals.

Let (T,d)(T,d) be a metric space. An admissible sequence of TT is a collection of subsets, Ts⊂TT_{s}\subset T, whose cardinality satisfies ∣Ts∣≤22s|T_{s}|\leq 2^{2^{s}} for s≥1s\geq 1, and ∣T0∣=1|T_{0}|=1. For α≥1\alpha\geq 1 and s0≥0s_{0}\geq 0 set

where the infimum is taken with respect to all admissible sequences of TT. When s0=0s_{0}=0 we shall write γα(T,d)\gamma_{\alpha}(T,d) instead of γ0,α(T,d)\gamma_{0,\alpha}(T,d).

For more information on chaining methods we refer the reader to M. Talagrand’s book , which contains an extensive and illuminating survey on the topic.

There exist absolute constants cc and CC for which the following holds. Let F⊂L2F\subset L_{2}, ∣F∣>1|F|>1, and consider {Gf:f∈F}\{G_{f}:f\in F\}, the centred, canonical gaussian process indexed by FF.

The upper bound in Theorem 1.3 is due to Fernique while the lower one is Talagrand’s Majorizing Measures Theorem . The proof of both parts may be found in . It should be noted that the original proofs of both results are based on the majorizing measures mechanism which preceded the modern generic chaining scheme.

Chaining arises as a way of relating the supremum of a random process {Zf:f∈F}\{Z_{f}:f\in F\} to the structure of the indexing set FF. An upper estimate is obtained by combining individual tail bounds that ensure that if ff and hh are close in some sense, the probability that ZfZ_{f} is very different from ZhZ_{h} is small. For example, if {Gf:f∈F}\{G_{f}:f\in F\} is the centred, canonical gaussian process indexed by F⊂L2F\subset L_{2}, and if f,h∈Ff,h\in F, then

which is precisely the sort of tail estimate one would like to have (but unfortunately, analogous versions of (1.4) are not as simple for many interesting processes).

The increment condition (1.4) hints towards a fundamental fact: the supremum of a gaussian process {Gf:f∈F}\{G_{f}:f\in F\} is determined by a single metric, endowed by the L2L_{2} norm In fact, Theorem 1.3 implies a two-sided control using the L2L_{2} metric, but since our focus is on obtaining upper bounds, we will focus only on that direction.. However, when it comes to empirical processes, the situation is rather different. To explain this substantial difference between empirical and gaussian processes it is convenient to use the notion of an Orlicz norm.

Let ff be defined on the probability space (Ω,μ)(\Omega,\mu). For α≥1\alpha\geq 1 set

and let LψαL_{\psi_{\alpha}} be the space of functions for which ∥f∥ψα<∞\|f\|_{\psi_{\alpha}}<\infty.

It is well known that ∥f∥ψα\|f\|_{\psi_{\alpha}} is equivalent to the smallest constant cc for which Pr(∣f∣>t)≤2exp⁡(−tα/cα)Pr(|f|>t)\leq 2\exp(-t^{\alpha}/c^{\alpha}) for every t>0t>0, and also to sup⁡p≥1∥f∥Lp/p1/α\sup_{p\geq 1}\|f\|_{L_{p}}/p^{1/\alpha}.

A class FF is LL-subgaussian if for every f,h∈F∪{0}f,h\in F\cup\{0\},

Using the moment characterization of the ψ2\psi_{2} norm, it is evident that if FF is LL-subgaussian then

for every f,h∈F∪{0}f,h\in F\cup\{0\} and for a suitable absolute constant cc.

Observe that if F⊂L2F\subset L_{2}, the centred, canonical gaussian process {Gf:f∈F}\{G_{f}:f\in F\} satisfies that ∥Gf−Gh∥L2=∥f−h∥L2\|G_{f}-G_{h}\|_{L_{2}}=\|f-h\|_{L_{2}} and

but there is no reason why ∥f−h∥ψ2\|f-h\|_{\psi_{2}} should be equivalent to ∥f−h∥L2\|f-h\|_{L_{2}} unless the class is subgaussian. And though one may show that there is an absolute constant c4c_{4} for which

the underlying metric in (1.5) is the ψ2\psi_{2} metric, which is simply too large to be of any use in most applications.

It turns out that this rather unsatisfactory upper bound may be improved by examining the chaining process more closely.

Let {Zf:f∈F}\{Z_{f}:f\in F\} be a random process, set s0≥0s_{0}\geq 0 and consider an admissible sequence (Fs)s≥0(F_{s})_{s\geq 0} and a collection of functions πs:F→Fs\pi_{s}:F\to F_{s}. Under mild assumptions (for example, that for every f∈Ff\in F, πsf→f\pi_{s}f\to f in an appropriate sense), each ZfZ_{f} may be written as a ‘chain’, which is simply a telescopic sum of the ‘links’ Zπs+1f−ZπsfZ_{\pi_{s+1}f}-Z_{\pi_{s}f}:

Obtaining uniform control over all chains requires balancing the tail estimates for each link that appears at the ss-stage with the total number of links that are involved at that stage. And, since there are at most

links at the ss-stage, it suffices to find levels t(f,s)t(f,s) for which

and obtain a similar estimate for the ‘starting points’ of each chain, Zπs0fZ_{\pi_{s_{0}}f}, to ensure the required uniform control over all possible chains. Such a uniform control results in a high probability upper estimate on sup⁡f∈FZf\sup_{f\in F}Z_{f}.

A trivial yet crucial observation is that by Chebychev’s inequality, a possible choice of t(f,s)t(f,s) is

for p∼2sp\sim 2^{s}, where here and throughout the article we write a∼ba\sim b if there are absolute constants c1c_{1} and c2c_{2} for which c1a≤b≤c2ac_{1}a\leq b\leq c_{2}a.

Thus, an obvious alternative to the γ\gamma-functionals is

where the infimum is taken with respect to all admissible sequences of FF and for πs:F→Fs\pi_{s}:F\to F_{s} which is the nearest point map with respect to the norm ∥ ∥Lu22s\|\ \|_{L_{u^{2}2^{s}}}.

It immediately follows from the decomposition to chains in (1.6) that for u≥4u\geq 4, with probability at least 1−2exp⁡(−c0u22s0)1-2\exp(-c_{0}u^{2}2^{s_{0}}),

where c0c_{0} and c1c_{1} are absolute constants.

The idea of using a complexity parameter that takes into account all the LpL_{p} structures endowed by the process has been introduced in and independently by R. Latała (see, for example, and Exercise 2.2.15 in ).

Let us study two examples in which the way the complexity parameter in (1.7) depends on FF takes a rather simple form.

∙\bullet When considering the centred, canonical gaussian process indexed by FF, the parameter in (1.7) is not new. Indeed, recall that

hence, for an almost optimal admissible sequence (1.7) becomes

and the relations between (1.8) and the upper estimate in Theorem 1.3 via the γ2\gamma_{2} functional are clear.

∙\bullet Next, one may consider (1.7) for the standard empirical process,

Applying Latała’s sharp bound on the moments of sums of independent random variables (see also Theorem 3.5, below), one may show that for every f,h∈Ff,h\in F and every p≥2p\geq 2,

The last example leads to the introduction of the following norms and to a ‘graded version’ of the γ\gamma functionals.

For a random variable ZZ and p≥1p\geq 1, set

Thus, if ZfZ_{f} is the empirical process from (1.10), it follows that

Given a class of functions FF, u≥1u\geq 1 and s0≥0s_{0}\geq 0, put

where the infimum is taken with respect to all admissible sequences (Fs)s≥0(F_{s})_{s\geq 0}, and πsf\pi_{s}f is the nearest point in FsF_{s} to ff with respect to the (u22s)(u^{2}2^{s}) norm.

Again, a straightforward chaining argument shows that for every s0≥0s_{0}\geq 0 and u≥4u\geq 4, with probability at least 1−2exp⁡(−c0u22s0)1-2\exp(-c_{0}u^{2}2^{s_{0}}),

Observe that if FF happens to be LL-subgaussian and u≥2u\geq 2 then Λs0,u(F)\Lambda_{s_{0},u}(F) is equivalent to γs0,2(F,L2)\gamma_{s_{0},2}(F,L_{2}). Indeed, by the moment characterization of the ψ2\psi_{2} norm, there is an absolute constant cc for which, for every p≥2p\geq 2, ∥f∥Lp≤cp∥f∥ψ2\|f\|_{L_{p}}\leq c\sqrt{p}\|f\|_{\psi_{2}}. Therefore,

This leads to the next corollary, which may also be established directly, using a standard chaining argument.

Let FF be an LL-subgaussian class. Then for every u≥4u\geq 4 and s0≥0s_{0}\geq 0, with probability at least 1−2exp⁡(−c02s0u2)1-2\exp(-c_{0}2^{s_{0}}u^{2}),

As an example, let d2(F)=sup⁡f∈F∥f∥L2d_{2}(F)=\sup_{f\in F}\|f\|_{L_{2}} and set s0≥0s_{0}\geq 0 be the largest integer for which

However, and unlike subgaussian examples, there are natural examples in which Λs0,u(F){\Lambda}_{s_{0},u}(F) may be significantly smaller than its ψ2\psi_{2} counterpart. This should not come as a surprise, as both ∥ ∥L2s/2s/2\|\ \|_{L_{2^{s}}}/2^{s/2} and ∥ ∥(2s)\|\ \|_{(2^{s})} are ‘local’ versions of the ψ2\psi_{2} metric: they measure the subgaussian behaviour of the functions involved, but only up to a fixed level, rather than at every level.

Clearly, if TT contains any one of the coordinate directions {e1,...,en}\{e_{1},...,e_{n}\}, then FTF_{T} is not a subset of Lψ2L_{\psi_{2}} and γ2(F,ψ2)\gamma_{2}(F,\psi_{2}) is not even well defined.

Talagrand showed in (see also ) that (1.1) has a geometric interpretation:

2 The main results

Up to now, we have only examined the standard empirical process. Unfortunately, the simple chaining argument used in (1.12) to control that process is rather useless when it comes to dealing with multiplier processes or with product processes. We will show in what follows how the suprema of these two types of processes may be bounded from above in terms of the Λ\Lambda-functionals.

Let us begin by formulating the estimate for multiplier processes.

For q>2q>2, there are constants c0c_{0}, c1,c2c_{1},c_{2} and c3c_{3} that depend only on qq, for which the following holds. Let ξ∈Lq\xi\in L_{q} and set ξ1,...,ξN\xi_{1},...,\xi_{N} to be independent copies of ξ\xi. Fix an integer s0≥0s_{0}\geq 0 and w,u>c0w,u>c_{0}. Then, with probability at least

One simple outcome of Theorem 1.9 is when the class FF happens to be LL-subgaussian.

Recall that d2(F)=sup⁡f∈F∥f∥L2d_{2}(F)=\sup_{f\in F}\|f\|_{L_{2}} and set s0≥0s_{0}\geq 0 to satisfy that

As noted previously, since FF is LL-subgaussian,

Thus, it follows from Theorem 1.9 that with probability at least

Turning to the question of product processes, recall the following fact from (see also Theorem 9.3.1 in ).

For every q>4q>4 there exists a constant c(q)c(q) that depends only on qq for which the following holds. Let F⊂LqF\subset L_{q} be a class of functions on (Ω,μ)(\Omega,\mu). Assume that for every t>0t>0 and every f,h∈F∪{0}f,h\in F\cup\{0\},

where d1d_{1} and d2d_{2} are metrics on F∪{0}F\cup\{0\}. If dq(F)=sup⁡f∈F∥f∥Lqd_{q}(F)=\sup_{f\in F}\|f\|_{L_{q}} and γ=γ2(F,d2)+γ1(F,d1)\gamma=\gamma_{2}(F,d_{2})+\gamma_{1}(F,d_{1}), then

Theorem 1.11 is demonstrated by integrating a high-probability bound. However, the probability estimate established in is far from optimal. Recently, Dirksen obtained the optimal probability estimate under the assumption that F⊂Lψ2F\subset L_{\psi_{2}}, improving earlier results from , and a few months later, Bednorz gave a different proof of the same fact. The following formulation is from .

There exist absolute constants c1c_{1} and c2c_{2} for which the following holds. Let FF be a class of functions on (Ω,μ)(\Omega,\mu) and set γ=γ2(F,ψ2)\gamma=\gamma_{2}(F,\psi_{2}) and dψ2(F)=sup⁡f∈F∥f∥ψ2d_{\psi_{2}}(F)=\sup_{f\in F}\|f\|_{\psi_{2}}. For every u>0u>0, with probability at least 1−2exp⁡(−c1min⁡{Nu,u2})1-2\exp(-c_{1}\min\{\sqrt{N}u,u^{2}\}),

The probability estimate in Theorem 1.12 is indeed optimal, as may be seen by setting t=udψ22(F)/Nt=ud_{\psi_{2}}^{2}(F)/\sqrt{N} in Bernstein’s inequality applied to f2f^{2}.

In comparison, below is our estimate on the supremum of a product process, and in particular, on the supremum of the quadratic process.

There exists an absolute constant c0c_{0} and for every q>4q>4 there exists a constant c1(q)c_{1}(q) that depends only on qq for which the following holds. Let FF and HH be classes of functions on (Ω,μ)(\Omega,\mu), set u≥max⁡{8,q}u\geq\max\{8,\sqrt{q}\} and consider an integer s0≥0s_{0}\geq 0. Then, with probability at least 1−2exp⁡(−c0u22s0)1-2\exp(-c_{0}u^{2}2^{s_{0}}), for every f∈Ff\in F and h∈Hh\in H,

In particular, if F=HF=H, then with probability at least 1−2exp⁡(−c0u22s0)1-2\exp(-c_{0}u^{2}2^{s_{0}}),

Let us present some of the outcomes of Theorem 1.13 for the quadratic process.

∙\bullet Since ∥ ∥(p)≤c∥ ∥ψ2\|\ \|_{(p)}\leq c\|\ \|_{\psi_{2}}, it is evident that Λs0,u(F)≤cγ2(F,ψ2)\Lambda_{s_{0},u}(F)\leq c\gamma_{2}(F,\psi_{2}). Thus,

if ∣F∣>1|F|>1 and one sets 2s0∼(γ2(F,ψ2)/dψ2(F))22^{s_{0}}\sim(\gamma_{2}(F,\psi_{2})/d_{\psi_{2}}(F))^{2}. Also, since ∥f∥ψ2∼sup⁡q≥2∥f∥Lq/q\|f\|_{\psi_{2}}\sim\sup_{q\geq 2}\|f\|_{L_{q}}/\sqrt{q}, it is evident that dq(F)≤c2qdψ2(F)d_{q}(F)\leq c_{2}\sqrt{q}d_{\psi_{2}}(F), and (1.17) recovers Theorem 1.12.

The proofs of Theorem 1.9 and Theorem 1.13 are based on symmetrization, which has been one of the most influential tools in empirical processes theory. The most well-known symmetrization inequalities for empirical processes are the celebrated Giné-Zinn inequalities , but we will use an earlier, “in-probability” version of those inequalities (see, e.g., ).

Thanks to Theorem 1.14, one may prove Theorem 1.9 via a high-probability upper bound on the supremum of the Bernoulli process

where the set VV is a typical coordinate projection of the class FF; that is, for σ=(X1,...,XN)\sigma=(X_{1},...,X_{N}),

In a similar fashion, Theorem 1.13 follows from a high-probability upper bound on

for typical coordinate projections V=PσFV=P_{\sigma}F and W=PσHW=P_{\sigma}H.

This observation dictates the structure of the article. We will first study

Finally, a word about notation. Throughout, c0,c1,...c_{0},c_{1},... denote absolute constants. Their value may change from line to line. c(q)c(q) or cqc_{q} are constants that depend only on the parameter qq, and a≲qba\lesssim_{q}b means that a≤c(q)ba\leq c(q)b.

Chaining and Bernoulli processes

As we noted earlier, our method of analysis consists of two main components. First, gathering accurate information on the structure of a typical coordinate projection of FF (and in the case of the product process, of a typical coordinate projection of HH as well); and second, for a typical σ=(X1,...,XN)\sigma=(X_{1},...,X_{N}) and (ξi)i=1N(\xi_{i})_{i=1}^{N}, analyzing the suprema of the conditioned Bernoulli processes

Moreover, by , this estimate is optimal when k∼t2k\sim t^{2}.

Therefore, the effect xx has on the Bernoulli process depends on

When dealing with products, as we have to, xi=wivix_{i}=w_{i}v_{i} and the decomposition requires additional care. Let II be the union of the sets of the kk largest coordinates of (∣wi∣)i=1N(|w_{i}|)_{i=1}^{N} and the kk largest coordinates of (∣vi∣)i=1N(|v_{i}|)_{i=1}^{N}; thus ∣I∣≤2k|I|\leq 2k. Fix r>1r>1, set r′r^{\prime} to be its conjugate index (i.e., 1/r+1/r′=11/r+1/r^{\prime}=1) and observe that with probability at least 1−2exp⁡(−t2/2)1-2\exp(-t^{2}/2),

As (2.1) is meant to play a part in a chaining argument, the bound should hold uniformly for all the links that appear at the ss-stage; hence, a likely choice in (2.1) at the ss-stage is t∼2s/2t\sim 2^{s/2}.

Let us present a chaining process aimed at bounding

Let (Vs)s≥0(V_{s})_{s\geq 0} be an admissible sequence of VV, and note that under very mild assumptions (for example, that for every v∈Vv\in V, πsv→v\pi_{s}v\to v in an appropriate sense as ss tends to infinity), for every s0≥0s_{0}\geq 0,

Set Δsv=πs+1v−πsv\Delta_{s}v=\pi_{s+1}v-\pi_{s}v; thus

Let r,r′>1r,r^{\prime}>1 be conjugate indices, for every s≥s0s\geq s_{0} set an integer js≥1j_{s}\geq 1 to be determined later and assume that (js)s≥s0(j_{s})_{s\geq s_{0}} is non-decreasing in ss. Finally, let I=Iv,sI=I_{v,s} be the union of the js−1j_{s}-1 largest coordinates of (∣Δsv∣i)i=1N(|\Delta_{s}v|_{i})_{i=1}^{N} and the js−1j_{s}-1 largest coordinates of (∣zi∣)i=1N(|z_{i}|)_{i=1}^{N}. Applying (2.1) for t≥4t\geq 4, it is evident that with probability at least 1−2exp⁡(−t22s/2)1-2\exp(-t^{2}2^{s}/2)

Repeating this argument for the vectors πs0v\pi_{s_{0}}v and summing the probabilities, it follows that with probability at least 1−2exp⁡(−ct22s0)1-2\exp(-ct^{2}2^{s_{0}}), for every v∈Vv\in V,

Motivated by this chaining argument, consider the following structural assumption on VV:

∙\bullet Let (js)s≥s0(j_{s})_{s\geq s_{0}} be a non-decreasing sequence of integers, and for every s≥s0s\geq s_{0}, 1≤js≤N+11\leq j_{s}\leq N+1.

Assume that for every s≥s0s\geq s_{0} and every v∈Vv\in V,

∙\bullet (∑i<js((Δsv)i∗)2)1/2≤∥Δsv∥[s]\left(\sum_{i<j_{s}}\left((\Delta_{s}v)_{i}^{*}\right)^{2}\right)^{1/2}\leq\|\Delta_{s}v\|_{[s]}, (∑i≥js((Δsv)i∗)2p)1/2p≤∥Δsv∥N1/2p\left(\sum_{i\geq j_{s}}\left((\Delta_{s}v)_{i}^{*}\right)^{2p}\right)^{1/2p}\leq\|\Delta_{s}v\|N^{1/2p},

∙\bullet (∑i<js((πsv)i∗)2)1/2≤∥πsv∥[s]\left(\sum_{i<j_{s}}\left((\pi_{s}v)_{i}^{*}\right)^{2}\right)^{1/2}\leq\|\pi_{s}v\|_{[s]}, (∑i≥js((πsv)i∗)2p)1/2p≤∥πsv∥N1/2p\left(\sum_{i\geq j_{s}}\left((\pi_{s}v)_{i}^{*}\right)^{2p}\right)^{1/2p}\leq\|\pi_{s}v\|N^{1/2p}.

Observe that if VV satisfies Assumption 2.1 for p=r′p=r^{\prime}, then by (2.2), with probability at least 1−2exp⁡(−ct22s0)1-2\exp(-ct^{2}2^{s_{0}}), for every v∈Vv\in V,

Using the same notation as above, with probability at least 1−2exp⁡(−ct22s0)1-2\exp(-ct^{2}2^{s_{0}}), for every v∈Vv\in V

Seemingly, there is plenty of freedom in the choices of jsj_{s}, rr and the admissible sequence of VV. However, our main interest is when zz and VV are typical realizations of (ξi)i=1N(\xi_{i})_{i=1}^{N} and V=PσFV=P_{\sigma}F respectively, and the natural choice of an admissible sequence of VV should be endowed by an admissible sequence of the underlying class FF. Thus, one must have adequate control on all the ‘monotone sums’

where Δsf=πs+1f−πsf\Delta_{s}f=\pi_{s+1}f-\pi_{s}f, IsI_{s} is the set of the js−1j_{s}-1 largest coordinates of (∣Δsf∣(Xi))i=1N\left(|\Delta_{s}f|(X_{i})\right)_{i=1}^{N}. In a similar fashion, one must be able to control

Since {Δsf:f∈F}\{\Delta_{s}f:f\in F\} contains at most 22s+22^{2^{s+2}} points that must be controlled uniformly, the individual probability estimate that is required in (2.4) and in (2.5) is exp⁡(−u22s)\exp(-u^{2}2^{s}) for a large enough uu (a choice of u≥4u\geq 4 will do). In what follows, we will show that this almost forces the choice of jsj_{s}.

In addition, obtaining sufficient control on (∑i=jsN(ξi∗)2r)1/2r\left(\sum_{i=j_{s}}^{N}(\xi_{i}^{*})^{2r}\right)^{1/2r} for every jsj_{s} restricts the choice of rr; it will depend on the LqL_{q} space to which ξ\xi belongs.

for NN independent copies of a random variable ZZ will be derived in Section 3.

2 Bernoulli product processes

Chaining for a Bernoulli product process is more involved than the one outlined above. However, the two share a common feature: structural assumptions on the indexing sets.

∙\bullet Let (js)s≥s0(j_{s})_{s\geq s_{0}} be an non-decreasing sequence of integers and for every s≥s0s\geq s_{0}, 1≤js≤N+11\leq j_{s}\leq N+1.

∙\bullet Let (νs)s≥s0(\nu_{s})_{s\geq s_{0}} be a sequence of nonnegative numbers.

Assume that for every s≥s0s\geq s_{0} and every v∈Vv\in V,

∙\bullet (∑i<js((Δsv)i∗)2)1/2≤∥Δsv∥[s]\left(\sum_{i<j_{s}}\left((\Delta_{s}v)_{i}^{*}\right)^{2}\right)^{1/2}\leq\|\Delta_{s}v\|_{[s]}, (∑i≥js−1((Δsv)i∗)2p)1/2p≤∥Δsv∥N1/2p\left(\sum_{i\geq j_{s-1}}\left((\Delta_{s}v)_{i}^{*}\right)^{2p}\right)^{1/2p}\leq\|\Delta_{s}v\|N^{1/2p},

∙\bullet (∑i<js((πsv)i∗)2)1/2≤∥πsv∥[s]\left(\sum_{i<j_{s}}\left((\pi_{s}v)_{i}^{*}\right)^{2}\right)^{1/2}\leq\|\pi_{s}v\|_{[s]}, (∑i≥js−1((πsv)i∗)2p)1/2p≤∥πsv∥N1/2p\left(\sum_{i\geq j_{s-1}}\left((\pi_{s}v)_{i}^{*}\right)^{2p}\right)^{1/2p}\leq\|\pi_{s}v\|N^{1/2p}.

∙\bullet (∑i=js−1js−1((πsv)i∗)2)1/2≤d(V)νs\left(\sum_{i=j_{s-1}}^{j_{s}-1}\left((\pi_{s}v)_{i}^{*}\right)^{2}\right)^{1/2}\leq d(V)\nu_{s}, where we set js0−1=js0j_{s_{0}-1}=j_{s_{0}} and recall that d(V)=sup⁡v∈V∥v∥d(V)=\sup_{v\in V}\|v\|.

There is a slight overlap between the sets of large and small coordinates: one set contains the js−1j_{s}-1 largest coordinates, while the other contains the N−js−1+1N-j_{s-1}+1 smallest coordinates – rather than the more natural set, consisting of the N−js+1N-j_{s}+1 smallest coordinates. The reason for this overlap is a minor technicality that will be used in the proof of Lemma 2.3.

In particular, if ∑s=s0s1−1νs≤cN\sum_{s=s_{0}}^{s_{1}-1}\nu_{s}\leq c\sqrt{N} then

Proof. We will only consider the case s1>s0+1s_{1}>s_{0}+1, as the proof of the case s1=s0+1s_{1}=s_{0}+1 is simpler, and is actually contained in the proof of the former.

Fix πs1v\pi_{s_{1}}v and let Is1−1I_{s_{1}-1} be the set of the js1−1−1j_{s_{1}-1}-1 largest coordinates of (∣(πs1v)i∣)i=1N(|(\pi_{s_{1}}v)_{i}|)_{i=1}^{N}. Therefore,

Repeating this argument for s0<s<s1s_{0}<s<s_{1}, it follows that for every v∈Vv\in V, and for the choice of IsI_{s} as the set of the js−1j_{s}-1 largest coordinates of (∣πs+1v∣i)i=1N(|\pi_{s+1}v|_{i})_{i=1}^{N},

Next, recall that v=∑s≥s1Δsv+πs1vv=\sum_{s\geq s_{1}}\Delta_{s}v+\pi_{s_{1}}v and that js=N+1j_{s}=N+1 for s≥s1s\geq s_{1}. Therefore, by Assumption 2.2 for p=1p=1,

There exist absolute constants c1c_{1} and c2c_{2} for which the following holds. If VV, WW, jsj_{s} and s0s_{0} are as above, then for every t>4t>4, with probability at least 1−2exp⁡(−c1t22s0)1-2\exp(-c_{1}t^{2}2^{s_{0}}),

Note that by Assumption 2.2, Lemma 2.3 and since js−1≤jsj_{s-1}\leq j_{s},

Let II be the union of the js−1j_{s}-1 largest coordinates of (∣Δsv∣i)i=1N(|\Delta_{s}v|_{i})_{i=1}^{N} and the js−1j_{s}-1 largest coordinates of (∣πs+1w∣i)i=1N(|\pi_{s+1}w|_{i})_{i=1}^{N}. Since ∣{Δsv:v∈V}∣≤22s+2|\{\Delta_{s}v:v\in V\}|\leq 2^{2^{s+2}} and ∣{πs+1w:w∈W}∣≤22s+1|\{\pi_{s+1}w:w\in W\}|\leq 2^{2^{s+1}}, it follows that for t≥4t\geq 4, with probability at least 1−2exp⁡(−ct22s)1-2\exp(-ct^{2}2^{s}), for every v∈Vv\in V and w∈Ww\in W,

Repeating the argument for ((πsv)i⋅(Δsw)i)i=1N((\pi_{s}v)_{i}\cdot(\Delta_{s}w)_{i})_{i=1}^{N} and summing over s≥s0s\geq s_{0}, one has that with probability at least 1−2exp⁡(−c1t22s0)1-2\exp(-c_{1}t^{2}2^{s_{0}}), for every v∈Vv\in V and w∈Ww\in W

it is evident that with probability at least 1−2exp⁡(−c2t22s)1-2\exp(-c_{2}t^{2}2^{s}),

Corollary 2.1 and Theorem 2.4 are the first component in the proofs of Theorem 1.9 and Theorem 1.13, respectively. For the other component, one has to show that typical coordinate projections of the indexing classes are well-behaved in the sense of Assumption 2.1 or of Assumption 2.2. To that end, one has to identify the norms ∥ ∥[s]\|\ \|_{[s]} and ∥ ∥\|\ \|, the sequences (js)s≥s0(j_{s})_{s\geq s_{0}} and (νs)s≥s0(\nu_{s})_{s\geq s_{0}} and estimate the resulting complexity terms, Λ(PσF)\Lambda(P_{\sigma}F), Θ(PσF)\Theta(P_{\sigma}F) and d(PσF)d(P_{\sigma}F). The main step towards that goal is presented in the next section.

Structural results - preliminary estimates

Let ZZ be a random variable. Our primary goal is to study the monotone nonincreasing rearrangement of NN independent copies of ZZ. We will show that the norms

There exist absolute constants c0,c1,c2c_{0},c_{1},c_{2} for which the following holds. Let 1≤r<q1\leq r<q and set 0<β<(q/r)−10<\beta<(q/r)-1. Consider Z∈LqZ\in L_{q}, and let Z1,...,ZNZ_{1},...,Z_{N} be independent copies of ZZ. Put 1≤p≤N1\leq p\leq N and set

∙\bullet For every t>2t>2 with probability at least 1−t−2pexp⁡(−p)1-t^{-2p}\exp(-p),

∙\bullet For every t>2t>2, with probability at least 1−t−j0qexp⁡(−p)1-t^{-j_{0}q}\exp(-p),

∙\bullet For t>2t>2 with probability at least 1−c2t−qN−β1-c_{2}t^{-q}N^{-\beta},

We will show that the vector UU consists of the j0−1j_{0}-1 largest coordinates of (∣Zi∣)i=1N(|Z_{i}|)_{i=1}^{N} and that V=(Zi)i=1N−UV=(Z_{i})_{i=1}^{N}-U. The key is to determine correctly the ‘cut-off’ point in that decomposition – which happens to be j0j_{0}.

As the formulation of Theorem 3.1 indicates, the treatment depends on whether or not the decomposition is trivial (i.e., if U=0U=0), and on the required probability estimate.

We begin with the smaller coordinates of (∣Zi∣)i=1N(|Z_{i}|)_{i=1}^{N}, which will be used to define the vector VV in the decomposition of (Zi)i=1N(Z_{i})_{i=1}^{N}.

There exist absolute constants c0,c1c_{0},c_{1} and c2c_{2} for which the following holds. Let 1≤r<q1\leq r<q, set Z∈LqZ\in L_{q} and put Z1,...,ZNZ_{1},...,Z_{N} to be independent copies of ZZ. Fix 1≤p≤N1\leq p\leq N, let

If j0>1j_{0}>1, then with probability at least 1−2t−j0qexp⁡(−p)1-2t^{-j_{0}q}\exp(-p),

And, if j0=1j_{0}=1 and 0<β<(q/r)−10<\beta<(q/r)-1 then with probability at least 1−c2t−qN−β1-c_{2}t^{-q}N^{-\beta},

Proof. Fix ρ≥1\rho\geq 1 to be named later and observe that by a binomial estimate and Markov’s inequality, for every u>0u>0,

Therefore, if u=t∥Z∥Lq(eN/j)ρ/qu=t\|Z\|_{L_{q}}(eN/j)^{\rho/q} then for any j≤Nj\leq N,

Note that (eN/j)−j(ρ−1)≤exp⁡(−p)\left({eN}/{j}\right)^{-j(\rho-1)}\leq\exp(-p) when

Set α=ρr/q\alpha=\rho r/q and observe that if α<1\alpha<1, that is, if ρ<q/r\rho<q/r, then

for an absolute constant c3c_{3}. The first part of the claim now follows by setting ρ=1+(q/r−1)/2\rho=1+(q/r-1)/2, implying that 1−α=(q−r)/2q1-\alpha=(q-r)/2q.

The second part follows an identical path to the first one, by setting 0<β<(q/r)−10<\beta<(q/r)-1, ρ=β+1\rho=\beta+1 and α=(β+1)r/q<1\alpha=(\beta+1)r/q<1.

Let q>4q>4 and r=2r=2. Set β=1/2<(q/r)−1\beta=1/2<(q/r)-1 and thus α<3/4\alpha<3/4. As noted in the proof of Lemma 3.2, with probability at least 1−exp⁡(−p)1-\exp(-p), for every i≥j0i\geq j_{0}

This simple fact will play a role in verifying the third part of Assumption 2.2 for a typical coordinate projection.

Next, let us show how the norms ∥ ∥(p)\|\ \|_{(p)} may be used to upper bound the larger coordinates in a monotone rearrangement of Z1,...,ZNZ_{1},...,Z_{N}.

Let Z1,...,ZNZ_{1},...,Z_{N} be independent copies of a random variable ZZ, set p≥log⁡Np\geq\log N and put 1≤m≤N/2e1\leq m\leq N/2e for which (Nm)≤exp⁡(p)\binom{N}{m}\leq\exp(p). Then, for every t>1t>1, with probability at least 1−t−2pexp⁡(−p)1-t^{-2p}\exp(-p), one has

The proof of Lemma 3.4 is based on the following fact, due to Latała .

Let WW be a nonnegative random variable. If W1,...,WmW_{1},...,W_{m} are independent copies of WW, then

Proof of Lemma 3.4. Since (Nm)≥(Nm−1)m≥exp⁡(m)\binom{N}{m}\geq(\frac{N}{m}-1)^{m}\geq\exp(m), it follows that p≥mp\geq m; thus, by Theorem 3.5 for W=Z2W=Z^{2} and r=pr=p,

and ∥∑i=1mZi2∥Lp≲2p∥Z∥(2p)2\|\sum_{i=1}^{m}Z_{i}^{2}\|_{L_{p}}\lesssim 2p\|Z\|_{(2p)}^{2}. Therefore,

for the choice of u∼t2p∥Z∥(2p)u\sim t\sqrt{2p}\|Z\|_{(2p)} and since (Nm)≤exp⁡(p)\binom{N}{m}\leq\exp(p).

The proof of Theorem 3.1 is simply the combination of Lemma 3.2 and Lemma 3.4.

Let us turn to the main application of Theorem 3.1.

From here on, fix u≥4u\geq 4 and for every s≥0s\geq 0 and r<qr<q set

for a suitable absolute constant c0c_{0} as in Theorem 3.1. Consider a finite class of functions HH, whose cardinality is at most 22s+32^{2^{s+3}}.

Writing jsj_{s} instead of js(r,q)j_{s}(r,q), let us examine three different cases: js=N+1j_{s}=N+1, 2≤js≤N2\leq j_{s}\leq N and js=1j_{s}=1. Motivated by the requirements of the chaining arguments outlined earlier, in all three cases one would like to obtain uniform control over all the functions in HH; thus, the probability estimate with which one must control the decomposition from Theorem 3.1 for each individual function should be at least 1−exp⁡(−2s+3)1-\exp(-2^{s+3}).

∙\bullet When js=N+1j_{s}=N+1, the decomposition is trivial, in the sense that for each function hh, U=(h(Xi))i=1NU=(h(X_{i}))_{i=1}^{N} and V=0V=0. Hence, setting 2p=u22s2p=u^{2}2^{s}, it follows that with probability at least 1−exp⁡(−u22s/2)1-\exp(-u^{2}2^{s}/2), for every h∈Hh\in H, (h(Xi))i=1N=U(h(X_{i}))_{i=1}^{N}=U, and

∙\bullet When 1<js≤N1<j_{s}\leq N, and setting 2p=u22s2p=u^{2}2^{s} once again, it follows that with probability at least 1−2exp⁡(−u22s/2)1-2\exp(-u^{2}2^{s}/2), for every h∈Hh\in H, (h(Xi))i=1N=U+V(h(X_{i}))_{i=1}^{N}=U+V, where

∙\bullet When js=1j_{s}=1, (h(Xi))i=1N=V(h(X_{i}))_{i=1}^{N}=V and U=0U=0. Moreover, because js=1j_{s}=1,

for constants c1c_{1} and c2c_{2} that depend only on rr and qq.

Let 0<β≤(q/r)−10<\beta\leq(q/r)-1 and note that by Theorem 3.1, with probability at least 1−c3N−β1-c_{3}N^{-\beta},

Therefore, (3.4) holds with probability of 1−2exp⁡(−θu22s)1-2\exp(-\theta u^{2}2^{s}) for every h∈Hh\in H if

which is the case when u≥c4(q,r)/βu\geq c_{4}(q,r)/\sqrt{\beta} and θ≤c5(q,r)β\theta\leq c_{5}(q,r)\beta.

Combining these observations yields the following outcome:

There exist absolute constants c0,c1c_{0},c_{1} and for every 1≤r<q1\leq r<q there exist constants c2c_{2} and c3c_{3} that depend only on qq and rr for which the following holds. Set

for 0<β≤(q/r)−10<\beta\leq(q/r)-1 and u≥c2/βu\geq c_{2}/\sqrt{\beta}. If H⊂LqH\subset L_{q} is of cardinality at most 22s+22^{2^{s+2}}, then with probability at least 1−2exp⁡(−c3βu22s)1-2\exp(-c_{3}\beta u^{2}2^{s}), for every h∈Hh\in H, (h(Xi))i=1N=Uh+Vh(h(X_{i}))_{i=1}^{N}=U_{h}+V_{h}; the support of each UhU_{h} is the set of the largest js−1j_{s}-1 coordinates of (∣h(Xi)∣)i=1N(|h(X_{i})|)_{i=1}^{N} while VhV_{h} is supported on its complement;

Proofs of the main results

Let us turn to the implications of Corollary 3.6 in the contexts of Assumption 2.1 and Assumption 2.2.

Let FF be a class of functions and set (Fs)s≥0(F_{s})_{s\geq 0} to be an admissible sequence of FF. Fix s0s_{0} to be named later and set s≥s0s\geq s_{0}. Clearly,

and thus Corollary 3.6 holds for both FsF_{s} and {Δsf:f∈F}\{\Delta_{s}f:f\in F\}.

For every s≥s0s\geq s_{0}, denote by [(Δsf)(Xi)]∗\left[(\Delta_{s}f)(X_{i})\right]^{*} the ii-th largest coordinate of the vector (∣Δsf∣(Xi))i=1N\left(|\Delta_{s}f|(X_{i})\right)_{i=1}^{N}, and put [(πsf)(Xi)]∗\left[(\pi_{s}f)(X_{i})\right]^{*} to be the ii-th largest coordinate of the vector (∣πsf∣(Xi))i=1N\left(|\pi_{s}f|(X_{i})\right)_{i=1}^{N}. Finally, recall that

Let us see how Corollary 3.6 may be used to prove Theorem 1.9.

Let q>2q>2, r=min⁡{1/2+q/4,2}r=\min\{1/2+q/4,2\} and r1=2r′r_{1}=2r^{\prime}, where r′r^{\prime} is the conjugate index of rr. Put q1=2r1q_{1}=2r_{1} and set β=1/2<(q1/r1)−1\beta=1/2<(q_{1}/r_{1})-1. Also, let

Below is the summary of the outcome of Corollary 3.6 when applied to the classes {Δsf:f∈F}\{\Delta_{s}f:f\in F\} and FsF_{s} for s≥s0s\geq s_{0}, q1q_{1} and r1r_{1} as above, and for u>max⁡{q1,4}u>\max\{\sqrt{q_{1}},4\}.

There are constants c1c_{1} and c2c_{2} that depend only on qq and an event of probability at least 1−2exp⁡(−c1u22s0)1-2\exp(-c_{1}u^{2}2^{s_{0}}) on which the following holds. For every f∈Ff\in F and s≥s0s\geq s_{0},

Observe that if q≥4q\geq 4 then r=2r=2, and thus 2r=2r′=42r=2r^{\prime}=4 and q1=8q_{1}=8. Therefore, all the constants in Corollary 4.1 are absolute constants and one may take any u≥8u\geq 8.

∙\bullet ∥ ∥[s]=c2u2s/2∥ ∥(u22s)\|\ \|_{[s]}=c_{2}u2^{s/2}\|\ \|_{(u^{2}2^{s})},

∙\bullet ∥ ∥=c2∥ ∥Lq1\|\ \|=c_{2}\|\ \|_{L_{q_{1}}},

Hence, for an almost optimal admissible sequence in FF,

Recall that u2≥q1u^{2}\geq q_{1}, and thus, for every s≥0s\geq 0,

Combining these observations with Corollary 2.1, it follows that for every (Xi)i=1N(X_{i})_{i=1}^{N} for which Corollary 4.1 holds (i.e., with probability at least 1−2exp⁡(−c1u22s0)1-2\exp(-c_{1}u^{2}2^{s_{0}}) with respect to μN\mu^{N}), and for every (ξi)i=1N(\xi_{i})_{i=1}^{N}, one has that with (εi)i=1N(\varepsilon_{i})_{i=1}^{N} probability at least 1−2exp⁡(−c1t22s0)1-2\exp(-c_{1}t^{2}2^{s_{0}}), for every f∈Ff\in F

Thus, to conclude the proof of Theorem 1.9, one has identify AA and BB for which, with high probability,

and then apply the symmetrization argument of Theorem 1.14.

By Lemma 3.2 and since 2r<1+q/2<q2r<1+q/2<q, one has that with probability at least 1−2exp⁡(−c0u22s0)1-2\exp(-c_{0}u^{2}2^{s_{0}}),

Let q>2q>2 and assume that ξ∈Lq\xi\in L_{q}. If ξ1,...,ξN\xi_{1},...,\xi_{N} are independent copies of ξ\xi and z=(ξi)i=1Nz=(\xi_{i})_{i=1}^{N}, then for every w>1w>1, with probability at least 1−c0w−qN−((q/2)−1)log⁡qN1-c_{0}w^{-q}N^{-((q/2)-1)}\log^{q}N,

where c0c_{0} and c1c_{1} depend only on qq.

Proof. Let η=(q/2)−1\eta=(q/2)-1, fix 1≤k≤N/21\leq k\leq N/2 and set v>0v>0 to be named later. A binomial estimate implies that

For every 1≤k≤N1\leq k\leq N, set uk=w/log⁡(eN/k)u_{k}=w/\log(eN/k) and observe that

The claim follows by summing the probability estimates.

Lemma 4.3 implies that one may select A∼qw∥ξ∥LqA\sim_{q}w\|\xi\|_{L_{q}} and with probability at least 1−c0(q)w−qN−((q/2)−1)log⁡qN−2exp⁡(−c1u22s0)1-c_{0}(q)w^{-q}N^{-((q/2)-1)}\log^{q}N-2\exp(-c_{1}u^{2}2^{s_{0}}),

Let us turn to a version of Theorem 1.9 when ξ∈Lψ2\xi\in L_{\psi_{2}}.

There exist absolute constants c1c_{1} and c2c_{2} for which the following holds. If ξ∈Lψ2\xi\in L_{\psi_{2}} then for every u,w≥8u,w\geq 8, with probability at least 1−2exp⁡(−c1u22s0)−2exp⁡(−c1Nw2)1-2\exp(-c_{1}u^{2}2^{s_{0}})-2\exp(-c_{1}Nw^{2}),

The proof follows a similar path to the proof of Theorem 1.9 with a minor modification in the last step – the bounds on (ξi)i=1N(\xi_{i})_{i=1}^{N}.

By Bernstein’s inequality, with probability at least 1−2exp⁡(−c0Nmin⁡{w2,w4})1-2\exp(-c_{0}N\min\{w^{2},w^{4}\}),

Therefore, if u≥8u\geq 8 and w≥1w\geq 1, then with probability at least

for absolute constants c0,c1c_{0},c_{1} and c2c_{2}.

The rest of the proof is unchanged, for the choices of r=r′=2r=r^{\prime}=2 and q1=4r′=8q_{1}=4r^{\prime}=8, as noted in Remark 4.2.

2 The quadratic process

Following the same path as in the previous section, and thanks to Theorem 2.4, one has to show that typical coordinate projections PσFP_{\sigma}F and PσHP_{\sigma}H satisfy Assumption 2.2 for p=1p=1 and p=2p=2.

Fix q>4q>4 and let js=js(2,q)j_{s}=j_{s}(2,q) be as in (3.3). Set

for α<3/4\alpha<3/4 as in Remark 3.3. It is straightforward to verify that if s1s_{1} is the smallest integer for which js=N+1j_{s}=N+1 then

To handle the first and second parts of Assumption 2.2, one may apply Corollary 3.6 to the classes {Δsf:f∈F}\{\Delta_{s}f:f\in F\}, {Δsh:h∈H}\{\Delta_{s}h:h\in H\}, FsF_{s} and HsH_{s} for s≥s0s\geq s_{0}.

There exists an absolute constant c1c_{1}, a constant c2c_{2} that depends only on qq and an event of probability at least 1−2exp⁡(−c1u22s0)1-2\exp(-c_{1}u^{2}2^{s_{0}}) on which the following holds. Consider u≥qu\geq\sqrt{q} and set js=js(2p,q)j_{s}=j_{s}(2p,q) as in (3.3), for p=1p=1 and p=2p=2. For every f∈Ff\in F and every s≥s0s\geq s_{0},

Therefore, with probability at least 1−2exp⁡(−c1u22s0)1-2\exp(-c_{1}u^{2}2^{s_{0}}), the coordinate projection PσFP_{\sigma}F satisfies Assumption 2.2 for p=1p=1 and p=2p=2, with the choices of

∙\bullet ∥ ∥[s]=c2u2s/2∥ ∥(u22s)\|\ \|_{[s]}=c_{2}u2^{s/2}\|\ \|_{(u^{2}2^{s})};

∙\bullet ∥ ∥=c2∥ ∥Lq\|\ \|=c_{2}\|\ \|_{L_{q}} (and, in particular, d(PσF)≤c2sup⁡f∈F∥f∥Lqd(P_{\sigma}F)\leq c_{2}\sup_{f\in F}\|f\|_{L_{q}});

∙\bullet νs=(∑i=jsjs+1−1(eNi)α)1/2\nu_{s}=\left(\sum_{i=j_{s}}^{j_{s+1}-1}\left(\frac{eN}{i}\right)^{\alpha}\right)^{1/2}, implying that ∑s=s0s1−1νs≤cN\sum_{s=s_{0}}^{s_{1}-1}\nu_{s}\leq c\sqrt{N}.

Moreover, if (Fs)s≥s0(F_{s})_{s\geq s_{0}} is an almost optimal admissible sequence,

This observation, together with Theorem 2.4 and the symmetrization argument of Theorem 1.14 completes the proof of Theorem 1.13.

3 Unconditional log-concave ensembles

Therefore, with probability at least 1−2exp⁡(−c3u22s0)1-2\exp(-c_{3}u^{2}2^{s_{0}}),

which improves the probability estimate from .

Acknowledgements

I am indebted to Vladimir Koltchinskii, Joe Neeman and Dong Xia for their careful reading of this manuscript and the many valuable comments and suggestions they have made.

References