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 . If are independent and distributed according to , the centred empirical process indexed by 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 , 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 and set to be the centred canonical gaussian process indexed by ; that is, the gaussian process indexed by whose covariance structure is endowed by the inner product in .
To avoid the (well understood) issue of measurability, set
(1) Let for some ( need not be independent of ) and set to be independent copies of . Consider
which is the supremum of the multiplier process indexed by and associated with the multiplier .
Note that unlike the standard notion of a multiplier process, here need not be independent of .
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 and be classes of functions defined on the probability space and consider the supremum of the product process indexed by and ,
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 , and the multiplier process as the standard empirical process indexed by the class (at least when 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 and respectively, rather than the structures of and , 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 functionals.
Let be a metric space. An admissible sequence of is a collection of subsets, , whose cardinality satisfies for , and . For and set
where the infimum is taken with respect to all admissible sequences of . When we shall write instead of .
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 and for which the following holds. Let , , and consider , the centred, canonical gaussian process indexed by .
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 to the structure of the indexing set . An upper estimate is obtained by combining individual tail bounds that ensure that if and are close in some sense, the probability that is very different from is small. For example, if is the centred, canonical gaussian process indexed by , and if , 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 is determined by a single metric, endowed by the norm In fact, Theorem 1.3 implies a two-sided control using the 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 be defined on the probability space . For set
and let be the space of functions for which .
It is well known that is equivalent to the smallest constant for which for every , and also to .
A class is -subgaussian if for every ,
Using the moment characterization of the norm, it is evident that if is -subgaussian then
for every and for a suitable absolute constant .
Observe that if , the centred, canonical gaussian process satisfies that and
but there is no reason why should be equivalent to unless the class is subgaussian. And though one may show that there is an absolute constant for which
the underlying metric in (1.5) is the 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 be a random process, set and consider an admissible sequence and a collection of functions . Under mild assumptions (for example, that for every , in an appropriate sense), each may be written as a ‘chain’, which is simply a telescopic sum of the ‘links’ :
Obtaining uniform control over all chains requires balancing the tail estimates for each link that appears at the -stage with the total number of links that are involved at that stage. And, since there are at most
links at the -stage, it suffices to find levels for which
and obtain a similar estimate for the ‘starting points’ of each chain, , to ensure the required uniform control over all possible chains. Such a uniform control results in a high probability upper estimate on .
A trivial yet crucial observation is that by Chebychev’s inequality, a possible choice of is
for , where here and throughout the article we write if there are absolute constants and for which .
Thus, an obvious alternative to the -functionals is
where the infimum is taken with respect to all admissible sequences of and for which is the nearest point map with respect to the norm .
It immediately follows from the decomposition to chains in (1.6) that for , with probability at least ,
where and are absolute constants.
The idea of using a complexity parameter that takes into account all the 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 takes a rather simple form.
When considering the centred, canonical gaussian process indexed by , 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 functional are clear.
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 and every ,
The last example leads to the introduction of the following norms and to a ‘graded version’ of the functionals.
For a random variable and , set
Thus, if is the empirical process from (1.10), it follows that
Given a class of functions , and , put
where the infimum is taken with respect to all admissible sequences , and is the nearest point in to with respect to the norm.
Again, a straightforward chaining argument shows that for every and , with probability at least ,
Observe that if happens to be -subgaussian and then is equivalent to . Indeed, by the moment characterization of the norm, there is an absolute constant for which, for every , . Therefore,
This leads to the next corollary, which may also be established directly, using a standard chaining argument.
Let be an -subgaussian class. Then for every and , with probability at least ,
As an example, let and set be the largest integer for which
However, and unlike subgaussian examples, there are natural examples in which may be significantly smaller than its counterpart. This should not come as a surprise, as both and are ‘local’ versions of the metric: they measure the subgaussian behaviour of the functions involved, but only up to a fixed level, rather than at every level.
Clearly, if contains any one of the coordinate directions , then is not a subset of and 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 -functionals.
Let us begin by formulating the estimate for multiplier processes.
For , there are constants , and that depend only on , for which the following holds. Let and set to be independent copies of . Fix an integer and . Then, with probability at least
One simple outcome of Theorem 1.9 is when the class happens to be -subgaussian.
Recall that and set to satisfy that
As noted previously, since is -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 there exists a constant that depends only on for which the following holds. Let be a class of functions on . Assume that for every and every ,
where and are metrics on . If and , 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 , 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 and for which the following holds. Let be a class of functions on and set and . For every , with probability at least ,
The probability estimate in Theorem 1.12 is indeed optimal, as may be seen by setting in Bernstein’s inequality applied to .
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 and for every there exists a constant that depends only on for which the following holds. Let and be classes of functions on , set and consider an integer . Then, with probability at least , for every and ,
In particular, if , then with probability at least ,
Let us present some of the outcomes of Theorem 1.13 for the quadratic process.
Since , it is evident that . Thus,
if and one sets . Also, since , it is evident that , 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 is a typical coordinate projection of the class ; that is, for ,
In a similar fashion, Theorem 1.13 follows from a high-probability upper bound on
for typical coordinate projections and .
This observation dictates the structure of the article. We will first study
Finally, a word about notation. Throughout, denote absolute constants. Their value may change from line to line. or are constants that depend only on the parameter , and means that .
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 (and in the case of the product process, of a typical coordinate projection of as well); and second, for a typical and , analyzing the suprema of the conditioned Bernoulli processes
Moreover, by , this estimate is optimal when .
Therefore, the effect has on the Bernoulli process depends on
When dealing with products, as we have to, and the decomposition requires additional care. Let be the union of the sets of the largest coordinates of and the largest coordinates of ; thus . Fix , set to be its conjugate index (i.e., ) and observe that with probability at least ,
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 -stage; hence, a likely choice in (2.1) at the -stage is .
Let us present a chaining process aimed at bounding
Let be an admissible sequence of , and note that under very mild assumptions (for example, that for every , in an appropriate sense as tends to infinity), for every ,
Set ; thus
Let be conjugate indices, for every set an integer to be determined later and assume that is non-decreasing in . Finally, let be the union of the largest coordinates of and the largest coordinates of . Applying (2.1) for , it is evident that with probability at least
Repeating this argument for the vectors and summing the probabilities, it follows that with probability at least , for every ,
Motivated by this chaining argument, consider the following structural assumption on :
Let be a non-decreasing sequence of integers, and for every , .
Assume that for every and every ,
, ,
, .
Observe that if satisfies Assumption 2.1 for , then by (2.2), with probability at least , for every ,
Using the same notation as above, with probability at least , for every
Seemingly, there is plenty of freedom in the choices of , and the admissible sequence of . However, our main interest is when and are typical realizations of and respectively, and the natural choice of an admissible sequence of should be endowed by an admissible sequence of the underlying class . Thus, one must have adequate control on all the ‘monotone sums’
where , is the set of the largest coordinates of . In a similar fashion, one must be able to control
Since contains at most points that must be controlled uniformly, the individual probability estimate that is required in (2.4) and in (2.5) is for a large enough (a choice of will do). In what follows, we will show that this almost forces the choice of .
In addition, obtaining sufficient control on for every restricts the choice of ; it will depend on the space to which belongs.
for independent copies of a random variable 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.
Let be an non-decreasing sequence of integers and for every , .
Let be a sequence of nonnegative numbers.
Assume that for every and every ,
, ,
, .
, where we set and recall that .
There is a slight overlap between the sets of large and small coordinates: one set contains the largest coordinates, while the other contains the smallest coordinates – rather than the more natural set, consisting of the 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 then
Proof. We will only consider the case , as the proof of the case is simpler, and is actually contained in the proof of the former.
Fix and let be the set of the largest coordinates of . Therefore,
Repeating this argument for , it follows that for every , and for the choice of as the set of the largest coordinates of ,
Next, recall that and that for . Therefore, by Assumption 2.2 for ,
There exist absolute constants and for which the following holds. If , , and are as above, then for every , with probability at least ,
Note that by Assumption 2.2, Lemma 2.3 and since ,
Let be the union of the largest coordinates of and the largest coordinates of . Since and , it follows that for , with probability at least , for every and ,
Repeating the argument for and summing over , one has that with probability at least , for every and
it is evident that with probability at least ,
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 and , the sequences and and estimate the resulting complexity terms, , and . The main step towards that goal is presented in the next section.
Structural results - preliminary estimates
Let be a random variable. Our primary goal is to study the monotone nonincreasing rearrangement of independent copies of . We will show that the norms
There exist absolute constants for which the following holds. Let and set . Consider , and let be independent copies of . Put and set
For every with probability at least ,
For every , with probability at least ,
For with probability at least ,
We will show that the vector consists of the largest coordinates of and that . The key is to determine correctly the ‘cut-off’ point in that decomposition – which happens to be .
As the formulation of Theorem 3.1 indicates, the treatment depends on whether or not the decomposition is trivial (i.e., if ), and on the required probability estimate.
We begin with the smaller coordinates of , which will be used to define the vector in the decomposition of .
There exist absolute constants and for which the following holds. Let , set and put to be independent copies of . Fix , let
If , then with probability at least ,
And, if and then with probability at least ,
Proof. Fix to be named later and observe that by a binomial estimate and Markov’s inequality, for every ,
Therefore, if then for any ,
Note that when
Set and observe that if , that is, if , then
for an absolute constant . The first part of the claim now follows by setting , implying that .
The second part follows an identical path to the first one, by setting , and .
Let and . Set and thus . As noted in the proof of Lemma 3.2, with probability at least , for every
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 may be used to upper bound the larger coordinates in a monotone rearrangement of .
Let be independent copies of a random variable , set and put for which . Then, for every , with probability at least , one has
The proof of Lemma 3.4 is based on the following fact, due to Latała .
Let be a nonnegative random variable. If are independent copies of , then
Proof of Lemma 3.4. Since , it follows that ; thus, by Theorem 3.5 for and ,
and . Therefore,
for the choice of and since .
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 and for every and set
for a suitable absolute constant as in Theorem 3.1. Consider a finite class of functions , whose cardinality is at most .
Writing instead of , let us examine three different cases: , and . 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 ; thus, the probability estimate with which one must control the decomposition from Theorem 3.1 for each individual function should be at least .
When , the decomposition is trivial, in the sense that for each function , and . Hence, setting , it follows that with probability at least , for every , , and
When , and setting once again, it follows that with probability at least , for every , , where
When , and . Moreover, because ,
for constants and that depend only on and .
Let and note that by Theorem 3.1, with probability at least ,
Therefore, (3.4) holds with probability of for every if
which is the case when and .
Combining these observations yields the following outcome:
There exist absolute constants and for every there exist constants and that depend only on and for which the following holds. Set
for and . If is of cardinality at most , then with probability at least , for every , ; the support of each is the set of the largest coordinates of while 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 be a class of functions and set to be an admissible sequence of . Fix to be named later and set . Clearly,
and thus Corollary 3.6 holds for both and .
For every , denote by the -th largest coordinate of the vector , and put to be the -th largest coordinate of the vector . Finally, recall that
Let us see how Corollary 3.6 may be used to prove Theorem 1.9.
Let , and , where is the conjugate index of . Put and set . Also, let
Below is the summary of the outcome of Corollary 3.6 when applied to the classes and for , and as above, and for .
There are constants and that depend only on and an event of probability at least on which the following holds. For every and ,
Observe that if then , and thus and . Therefore, all the constants in Corollary 4.1 are absolute constants and one may take any .
,
,
Hence, for an almost optimal admissible sequence in ,
Recall that , and thus, for every ,
Combining these observations with Corollary 2.1, it follows that for every for which Corollary 4.1 holds (i.e., with probability at least with respect to ), and for every , one has that with probability at least , for every
Thus, to conclude the proof of Theorem 1.9, one has identify and for which, with high probability,
and then apply the symmetrization argument of Theorem 1.14.
By Lemma 3.2 and since , one has that with probability at least ,
Let and assume that . If are independent copies of and , then for every , with probability at least ,
where and depend only on .
Proof. Let , fix and set to be named later. A binomial estimate implies that
For every , set and observe that
The claim follows by summing the probability estimates.
Lemma 4.3 implies that one may select and with probability at least ,
Let us turn to a version of Theorem 1.9 when .
There exist absolute constants and for which the following holds. If then for every , with probability at least ,
The proof follows a similar path to the proof of Theorem 1.9 with a minor modification in the last step – the bounds on .
By Bernstein’s inequality, with probability at least ,
Therefore, if and , then with probability at least
for absolute constants and .
The rest of the proof is unchanged, for the choices of and , 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 and satisfy Assumption 2.2 for and .
Fix and let be as in (3.3). Set
for as in Remark 3.3. It is straightforward to verify that if is the smallest integer for which then
To handle the first and second parts of Assumption 2.2, one may apply Corollary 3.6 to the classes , , and for .
There exists an absolute constant , a constant that depends only on and an event of probability at least on which the following holds. Consider and set as in (3.3), for and . For every and every ,
Therefore, with probability at least , the coordinate projection satisfies Assumption 2.2 for and , with the choices of
;
(and, in particular, );
, implying that .
Moreover, if 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 ,
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.