Testing hypotheses simultaneously is a frequent issue in statistical practice, e.g. in genomic research. A widely used criterion for deciding which of these hypotheses should be rejected is the so-called ”false discovery rate” (FDR) promoted by Benjamini and Hochberg . The FDR is the expectation of the ”false discovery proportion” (FDP), the ratio of the number of false rejections and all rejections . Let a level be given. Under the so-called basic independence (BI) assumption we have FDR for the classical Benjamini and Hochberg linear step-up test, briefly denoted by BH test. Here, is the unknown amount of true null hypotheses. To achieve a higher power it is of great interest to exhaust the FDR as good as possible. Especially, if is not close to there is space for improvement of the BH test. That is why since the beginning of this century the interest of adaptive tests grows. The idea is to estimate by an appropriate estimator in a first step and to apply the BH test for the (data depended) level
in the second step. Heuristically, we obtain for a good estimatorthat FDR. Benjamini and Hochberg  suggested an early estimator for leading to an FDR controlling test. Before them Schweder and Spjøtvoll  already discussed estimators for using plots of the empirical distribution function of the -values. The number of estimators suggested in the literature is huge, here only a short selection: Benjamini et al. , Blanchard and Roquain [7, 8] and Zeisel et al. . We want to emphasize the papers of the Storey  and Storey et al.  and, in particular, the Storey estimator based on a tuning parameter . We refer to Storey and Tibshirani  for a discussion of the adjustment of the tuning parameter . Generalized Storey estimators with data dependent weights, which were already discussed by Heesen and Janssen , will be our prime example for our general results. A nice property of them is that we have finite sample FDR control, see . Sufficient conditions for finite sample FDR control on general estimators can be found in Sarkar  and Heesen and Janssen [20, 21].
Beside the FDR control there are also other control criteria, for example the family-wise error rate FWER. Also for the control of FWER adaptive tests, i.e. tests using an (plug-in) estimator for , are used and discussed in the literature, see e.g. Finner and Gontscharuk  and Sarkar et al. .
Stochastic process methods were applied to study the asymptotic behavior of the FDP, among others to calculate asymptotic confidence intervals, and the familywise error rate (FWER) in detail, seeGenovese and Wassermann , Meinshausen and Bühlmann , Meinshausen and Rice  and Neuvial . Dealing with a huge amount of -values the fluctuation of the FDP becomes, of course, relevant. Ferreira and Zwinderman  presented formulas for higher moments of FDP for the BH test and Roquain and Villers  did so for step-up and step-down tests with general (but data independent) critical values. We generalize these formulas to adaptive step-up tests using general estimators for . In particular, we derive an exact finite sample formula for the variability of FDP. As an application of this we discuss the consistency of FDP and present sufficient and necessary conditions for it. We also discuss the more challenging case of sparsity in the sense that as . This situation can be compared to the one of Abramovich et al. 
, who derived an estimator of the (sparse) mean of a multivariate normal distribution using FDR procedures.
Outline of the results. In Section 2 we introduce the model as well as the adaptive step-up tests, and in particular the generalized Storey estimators which serve as prime examples. Section 3 provides exact finite sample variance formulas for the FDP under the BI model. Extensions to higher moments can be found in the appendix, see Section 9. These results apply to the variability and the consistency of FDP, see Section Section 4. Roughly speaking we have consistency if we have stable estimators and the number of rejections tends to infinity. Section 5 is devoted to concrete adaptive step-up tests mainly based on the convex combinations of generalized Storey estimators with data dependent weights. We will see that consistency cannot be achieved in general. Under mild assumptions the adaptive tests based on the estimators mentioned above are superior compared to the BH test: The FDR is more exhausted but remains finitely controlled by the level . Furthermore, they are consistent at least when the BH test is consistent. In Section 6 we discuss least favorable configurations which serve as useful technical tools. For the reader’s convenience we add a discussion and summary of the paper in Section 7. All proofs are collected in Section 8.
2.1 The model and general step-up tests
Let us first describe the model and the procedures. A multiple testing problem consists on null hypotheses with associated -value on a common probability space . We will always use the basic independence (BI) assumption given by
The set of hypotheses can be divided in the disjoint union of unknown portions of true null and false null , respectively. Denote by the cardinality of for .
The vectors of-values and are independent, where each dependence structure is allowed for the -values for the false hypotheses.
of the true null are independent and uniformly distributed on, i.e. for all .
Throughout the paper let be nonrandom. As in Heesen and Janssen  the results can be extended to more general models with random by conditioning under . By using this modification the results easily carry over to familiar mixture models discussed, for instance, by Abramovich et al.  and Genovese and Wassermann . We study adaptive multiple step-up tests with estimated critical values extending the famous Benjamini and Hochberg  step-up test, briefly denoted by BH test. In the following we recall the definition of this kind of tests. Let
denote possibly data dependent critical values. As an example for the critical values we recall the ones for the BH test, which do not depend on the data:
If denote the ordered -value then the number of rejections is given by
and the multiple procedure acts as follows:
be the number of falsely rejected null hypothesis. Then the false discovery rateand the false discovery proportion are given by
Good multiple tests like the BH test or the frequently applied adaptive test of Storey et al.  control the FDR at a pre-specified acceptance error bound at least under the BI assumption. Besides the control, two further aspects are of importance and discussed below:
To make the test sensitive for signal detection the FDR should exhaust the level as best as possible.
On the other hand the variability of the FDP, see (2.3), is of interest in order to judge the stability of the test.
For a large class of adaptive tests exact FDR formulas were established in Heesen and Janssen . Here the reader can find new -controlling tests with high FDR. These formulas are now completed by formulas for exact higher FDP moments and, in particular, for the variance. These results open the door for a discussion about the consistency of multiple tests, i.e.
As already stated we can not expect consistency in general. In the following we discuss the BH test for two extreme case.
Let be fixed. Then is minimal for the so-called Dirac uniform configuration DU, where all entries of are equal to zero. Under this configuration in distribution with
. The limit variable belongs to the class of linear Poisson distributions, seeFinner et al. , Jain  and Consul and Famoye . Hence, the BH-test is never consistent under BI for fixed since (2.5) is violated.
More information about DU and least favorable configurations can be found in Section 6. The requirement for consistency will be somehow in between these two extreme cases where the assumption will be always needed.
2.2 Our step-up tests and underlying assumptions
In the following we introduce the adaptive step-up tests we consider in this paper. Let be a fixed level and let , , be a tuning parameter and we agree that no null with should be rejected. The latter is not restrictive for practice since it is very unusual to reject a null if the corresponding -value exceeds, for instance, . We divide the range of the -values in a decision region , where all with have a chance to be rejected, and an estimation region , where are used to estimate , see Figure 1.
To be more specific we consider estimators of the form
for estimating , which are measurable functions depending only on . As usual we denote by the empirical distribution function of the -values . As motivated in the introduction we now plug-in these estimators in the BH test. Doing this we obtain the data driven critical values
where we promote to use the upper bound as Heesen and Janssen  already did. The following two quantities will be rapidly used: Through
Throughout this paper, we investigate different mild assumptions. For our main results we fix the following two:
Suppose that is always positive and
If only is valid then our results apply to appropriate subsequences. The most interesting case is since otherwise (if ) the FDR can be controlled, i.e. , by rejecting everything.
A prominent example for an adaptive test controlling the FDR by is given by the Storey estimator (2.11):
A refinement was established by Heesen and Janssen . They introduced a couple of inspection points , where is estimated on each interval . As motivation for this idea observe that the Storey estimator can be rewritten as the following linear combination
with weights , where . The ingredients
are also estimators for , which were used by Liang and Nettleton  in another context. Under BI the following theorem was proved by Heesen and Janssen . A discussion of their consistency is given in Section 5.
Theorem 2.3 (cf. Thm 10 in )
Let be random weights for with . The adaptive step-up tests using the estimator
controls the FDR, i.e. .
Finally, we want to present a necessary condition of asymptotic FDR control. It was proven by Heesen and Janssen  for a greater class than the BI models, namely reverse martingals. The same condition was already used by Finner and Gontscharuk  for asymptotic FWER control.
This section provides exact second moment formulas of for our adaptive step-up tests for a fixed regime . Our method of proof relies on conditioning with respect to the -algebra
Conditional under the (non-observable) -algebra the quantities , and are fixed values. But only and are given by the data and observable. The FDR formula (2.9) is now completed by an exact variance formula. The proof offers also a rapid approach to the known variance formula of Ferreira and Zwinderman  for the Benjamini and Hochberg test (with and ). Without loss of generality we can assume that is ordered by
Now, we introduce a new -value vector . If then set equal to but replace one -value with by for one , for convenience take the smallest integer with this property. If then set . Moreover, let be the number of rejections of the adaptive test for substituted vector of -value. Note that remains unchanged when considering instead of .
Suppose that our assumptions (A2) are fulfilled:
The second moment of is given by
The variance of FDP fulfils
Exact higher moment formulas are established in the appendix, see Section 9.
4 The variability of and the Consistency of adaptive multiple tests
The exact variance formula applies to the stability of the FDP and its consistency if tends to infinity. If not stated otherwise, all limits are meant as . In the following we need a further mild assumption:
There is some such that for all .
Clearly, (A3) is fulfilled for the trivial estimator and for all generalized weighted estimators of the form (2.13) with . Note that (A1) and (A3) imply and, hence, (2.5) is a necessary condition for consistency in this case. In the following we give boundaries for the variance of depending on the leading term in the variance formula of Theorem 3.5:
Since under (A1) we have consistency iff . In the following we present sufficient and necessary conditions for this.
Under mild additional assumptions we can improve the convergence in expectation from Remark 4.8(b). Recall that and depend, of course, on the pre-specified level . In comparison to the rest of this paper we consider in the following theorem we consider more than one level. That is why we prefer (only) for this theorem the notation and
The next example points out that consistency may depend on the level and adaptive tests may be consistent while the BH test is not so. A proof of the statements is given in Section 8.
Let be i.i.d. uniformly distributed on . Consider , and -values from the false null given by with , . The BH test BH() with level is not consistent while BH is consistent. But the adaptive test Stor using the classical Storey estimator (2.10) is consistent.
5 Consistent and inconsistent regimes
Below we will exclude the ugly estimator which could lead to rejecting all hypotheses with . To avoid this effect let us introduce:
There exists a constant with
Consistency for the case was already discussed by Genovese and Wassermann , who used a stochastic process approach. Also Ferreira and Zwinderman  used their formulas for the moments of to discuss the consistency for the BH test. By their Proposition 2.2 or our Theorem 4.7 it is sufficient to show for that in -probability. For this purpose Ferreira and Zwinderman  found conditions such that in -probability. The sparse signal case is more delicate since always tends to even for adaptive tests. Recall for the following lemma that is the largest intersection point of and the random Simes line , observe .
Besides the result of Theorem 5.11 we already know that a further necessary condition for consistency is (4.3) which is assumed to be fulfilled in the following. Turning to convergent subsequence we can assume without loss of generality under (A1), (A3) and (A4) that . In this case (4.3) is equivalent to
In the following we will rather work with (5.1) instead of (4.3). Due to Lemma 5.12 the question about consistency can be reduced in the sparse signal case to the comparison of the random Simes line defined above and close to .
In the case we need a sequence tending to .
For the DU-configuration the assumption (5.2) is fulfilled for as long as the necessary condition holds.
As already stated, consistency only holds under certain additional assumptions. In the following we compare consistency of the classical BH test and adaptive tests with an appropriate estimator.
5.1 Combination of generalized Storey estimators
In the following we become more concrete by discussing the combined Storey estimators (2.13) introduced in Section 1. For this purpose we need the following assumption to ensure that (A4) is fulfilled.
Suppose that .
Suppose that . Then (5.1) holds with and
Suppose that and we have with probability one that
for every . Moreover, assume that
and all . If there is some and such that
where , then we have an improvement of the asymptotically compared to the Benjamini-Hochberg procedure, i.e.
(Consistency) Suppose that the weights are asymptotically constant, i.e. a.s. for all , and fulfil (5.6). Assume that
and for all . Moreover, suppose that
for some , where . Additionally, assume if . Then (5.1) holds for some and consistency of the BH test implies consistency of the adaptive test. Moreover, if (5.2) holds for and a sequence with then we have always consistency of the adaptive test.
It is easy to see that the assumptions of (c) imply the ones of (b). Typically, the -values , , from the false null are stochastically smaller than the uniform distribution, i.e. for all (with strict inequality for some ). This may lead to (5.7) or (5.10).
5.2 Asymptotically optimal rejection curve
Our results can be transferred to general deterministic critical values (2.1), which are not of the form (2.7) and do not use a plug-in estimator for . To label this case we use . Analogously to Section 4 and Section 9 we define for by setting -values from the true null to . By the same arguments as in the proof of Theorem 3.5 we obtain
The first formula can also be found in Benditkis et al. , see the proof of Theorem 2 therein. The proof of the second one is left to the reader. By these formulas we can now treat an important class of critical values given by
A necessary condition for the valid step-up tests is . This condition holds for the critical values (5.12) if
These critical values are closely related to
of the asymptotically optimal rejection curve introduced by Finner et al. . Note that the case is excluded on purpose because it would lead to . The remaining coefficient has to be defined separately such that , see Finner et al.  and Gontscharuk  for a detailed discussion. It is well-known that neither for (5.12) with and nor for (5.14) we have control of the FDR by over all BI models simultaneously. This follows from Lemma 4.1 of Heesen and Janssen  since . However, Heesen and Janssen  proved that for all fixed , and there exists a unique parameter such that
where the supremum is taken over all BI models at sample size . The value may be found under the least favorable configuration DU using numerical methods.
By transferring our techniques to this type of critical values we get the following sufficient and necessary conditions for consistency.
If , and then (5.15) is sufficient for consistency and, moreover, in this case.
6 Least favorable configurations and consistency
Below least favorable configurations (LFC) are derived for the -value of the false portion. When deterministic critical values are increasing then the FDR is decreasing in each argument , , for fixed , see Benjamini and Yekutieli  or Benditkis et al.  for a short proof. Here and subsequently, we use ”increasing” and ”decreasing” in their weak form, i.e. equality is allowed, whereas other authors use ”nondecreasing” and ”nonincreasing” for this purpose. In that case the Dirac uniform configuration DU, see Example 2.1, has maximum FDR, i.e. it is LFC. LFC are sometimes useful tools for all kind of proofs.
Our exact moment formulas provide various LFC-results which are collected below. To formulate these we introduce a new assumption
Let be increasing for each coordinate .
Below we are going to condition on . By (BI2) we may write , where represents the distribution of under for , and .
Theorem 6.20 (LFC for adaptive tests)
Suppose that (A2) is fulfilled. Define the vector .
The conditional FDR conditioned on
only depends on the portion