REVIEW 2 major objections 4 minor 40 references
A new class of tests for convex-ordered families based on expected order statistics
T0 review · 2 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read A general class of nonparametric tests uses expected order statistics to decide whether an unknown distribution lies in a convex-ordered family relative to a known reference, with unbiasedness, monotone power, and consistency for heavy…
desk verdict Solid new tests for convex-ordered families with an honest heavy-tail extension; the finite-sample proof of monotone power has a real but repairable gap that needs a referee's attention. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The engine is the identity $\pi^G_{j:m} = G(E(G^{-1}(B_{j:m})))$, where $B_{j:m}$ is a $\beta$$(j, m-j+1)$ random variable, together with the Jensen bound of Proposition 1 relating $P(X \leq \mu_{j:m})$ to this quantity. The expected order statistics $\mu_{j:m}(F)$ are estimated by L-estimators, weighted linear combinations of sample order statistics, and Theorem 1 gives an almost-sure strong law for these estimators when $E|X|$ is infinite, under maximum-domain-of-attraction and regular-variation conditions. The finite-sample theory rests on Lemma 1, which uses the linearly interpolated empirical quantile composition to show that $F \le_c H$ implies $\widetilde F_n(\mu_{j:m}(F_n)) \le_{\mathrm{st}} \widetilde H_n(\mu_{j:m}(H_n))$; this stochastic monotonicity is what carries the unbiasedness and monotone-power results for every sample size.
What would settle it
Simulate from a pair $F \le_c H$ (for example exponential versus Weibull with shape 1.5) and check every sample path for whether $\widetilde H_n^{-1}\circ\widetilde F_n$ is convex between consecutive order statistics; a single sample with a downward kink disproves the Lemma 1 premise, and Monte Carlo can then test whether rejection probabilities under $F$ and $H$ violate the claimed monotonicity.
Extended reading notes
Core claim
The central claim is that, for any absolutely continuous reference distribution $G$, the test statistics $T^{G+}_{m,p}$ and $T^{G-}_{m,p}$ --- the positive and negative parts of the gap between the bounds $\pi^G_{j:m}$ and the linearly interpolated empirical CDF evaluated at estimated expected order statistics --- solve the problem of testing $H_0: F$ belongs to the location-scale family of $G$ against the alternatives that $F$ is strictly convex- or concave-ordered with respect to $G$. Proposition 1 provides the load-bearing bound: $F \in \mathcal{F}^{cx}_G$ implies $P(X \leq \mu_{j:m}) \leq \pi^G_{j:m}$, and $F \in \mathcal{F}^{cv}_G$ implies the reverse inequality. The paper proves that these tests are unbiased and have monotone power at every fixed sample size, and that they are consistent both in the finite-mean case and, through a constrained version that discards non-converging ranks, for heavy-tailed distributions with infinite or undefined mean. It also reports simulations indicating that the tests work for the IHR/DHR, IOR/DOR, and DRHR/IRHR families, including detection of non-monotone hazard rates where the normalized-spacings test can be misleading.
Load-bearing premise
The finite-sample results depend on the claim that linearly interpolating the empirical distribution preserves the convexity of the transformed quantile function in Lemma 1; if a sample can make the smoothed quantile composition non-convex where the true composition is convex, unbiasedness and monotone power no longer follow.
Editorial extensions
If this is right
- The same testing recipe applies to every absolutely continuous reference $G$, including uniform, exponential, negative exponential, log-logistic, Frechet, and Cauchy, with no support restrictions.
- Families such as IOR, DOR, and DRHR, for which the paper says no comparable tests currently exist, become testable with guaranteed size control and consistency.
- The constrained version of the test remains consistent when the mean is infinite, so a Cauchy or other very heavy-tailed $F$ can still be tested using ranks near the median, such as $\hat\mu_{2:3}$.
- For the exponential reference, running the IHR and DHR versions together can detect a non-monotone hazard rate: the proposed tests reject in favour of both alternatives, while normalized-spacings tests may incorrectly report a single monotone direction.
- The tests complement goodness-of-fit procedures: rejecting $H_0$ while a goodness-of-fit test accepts the shape class provides evidence that $F$ lies in the convex-ordered family outside the location-scale family.
Reading between the lines
- A natural extension, not developed in the paper, is to automate the choice of safe ranks by estimating the tail index from the sample (for example with a Hill-type estimator) and then selecting $m$ and $\ell$ data-adaptively rather than from a prior lower bound on $\alpha$.
- The same bounds could support confidence statements about shape rather than only tests: the distance between $\pi^G_{j:m}$ and $F(\mu_{j:m})$ is a population measure of how far $F$ is from the location-scale family, and its L-estimator could be used to construct an effect-size summary.
- Because the tests are location- and scale-invariant and do not restrict support, they could be applied iteratively to narrow down the shape of a distribution, as the paper does for river-flow data; this suggests a model-selection workflow where families are accepted or rejected one by one.
- The strictness of Jensen's inequality in the consistency proof means the tests distinguish the location-scale family from the rest of the convex-ordered family; they do not by themselves certify membership in the shape class, so a combined use with existing goodness-of-fit tests is the intended reading.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper proposes a class of nonparametric tests for convex-ordered families F^cx_G and F^cv_G, with null hypothesis H^G_0: F belongs to the location-scale family of a known absolutely continuous reference distribution G. The test statistics are L_p-norms of the positive or negative parts of the vector pi^G_{j:m} - \tilde F_n(\hat\mu_{j:m}(F_n)), where pi^G_{j:m} = G(E(G^{-1}(B_{j:m}))) and \hat\mu_{j:m} is an L-estimator of the expected order statistic. The paper proves an almost-sure limit theorem for these L-estimators under extreme-value domain-of-attraction conditions, including cases where the mean is infinite, and claims finite-sample unbiasedness and monotone power through Lemma 1 and Theorem 2, plus consistency through Propositions 5 and 6. Simulations compare the method with the Proschan-Pyke test and apply it to river flow data.
Significance. If the finite-sample claims are properly established, this is a useful contribution: it provides a unified testing method for several shape-constrained families, including IOR, DOR, DRHR and IRHR, for which direct tests are largely unavailable, and it does so without support restrictions and with robustness to infinite means. Theorem 1, the strong law for L-estimators under heavy tails, is a valuable standalone result, and the paper is transparent about code and simulations. The central inequality Proposition 1 is taken from the authors' prior published work, but it is a Jensen-type bound rather than a fitted quantity, so I do not see a circularity problem. The main caveat is that the proof of the finite-sample stochastic-order result is incomplete as written, so the advertised unbiasedness and monotone power properties are not yet logically guaranteed by the text.
major comments (2)
- [Section 4.1, proof of Theorem 2] The proof passes from componentwise stochastic inequalities \tilde U_j ≤_st \tilde Z_j, obtained from Lemma 1, to stochastic ordering of the increasing function ψ_p(z) = ||z||_p by invoking Theorem 1.A.3 of Shaked and Shanthikumar. That theorem requires the coordinates within each vector to be independent. Here the coordinates are all functions of the same empirical CDF H_n (or F_n), so they are dependent, and componentwise stochastic dominance does not imply multivariate stochastic dominance for dependent coordinates. This is a load-bearing gap because Corollary 2 and the advertised finite-sample unbiasedness and monotone power rest on Theorem 2. The gap appears repairable: the proof of Lemma 1 constructs a coupling F_n = H_n ∘ H^{-1} ∘ F for which the inequalities hold pointwise for all j simultaneously, and applying ψ_p to both sides of that pointwise vector inequality would establish the desired stochastic order directly. The repair should be written into the proof.
- [Section 4.1, Theorem 2(2)] The second displayed inequality in Theorem 2 has the wrong direction. Under F ≤_c H, Lemma 1 gives \tilde F_n(μ_{j:m}(F_n)) ≤ \tilde H_n(μ_{j:m}(H_n)). Since x ↦ (π^G_{j:m} - x)_- is decreasing, the componentwise inequality reverses for the statistic T^{G-}_{m,p}. The paragraph after Lemma 1 also states that under F ≤_c H the rejection probability for H^{G}_{1-} is smaller under F than under H, and Corollary 2(2) requires P(T^{G-}_{m,p}(F_n) ≥ c^-_{α,n}) ≤ P(T^{G-}_{m,p}(H_n) ≥ c^-_{α,n}) when F ≤_c H. Therefore part 2 of Theorem 2 should read ≤, not ≥, or should be reformulated with the order reversed. As written, the theorem contradicts Corollary 2(2) when H = G and F ≤_c G.
minor comments (4)
- [Section 4.1, proof of Lemma 1] The statement that \tilde H_n^{-1} ∘ \tilde F_n 'is a convex function, interpolating H_n^{-1} ∘ F_n' needs one explicit justification: this composition is the piecewise-linear interpolant of the convex function H^{-1} ∘ F at the increasing knots x_i, and a piecewise-linear interpolant of a convex function with increasing knots is convex because its successive slopes are nondecreasing. Adding this sentence would make the proof easier to follow.
- [Section 4.2, proof of Theorem 1] Several displays in the proof of Theorem 1 have lost their superscript formatting. In particular, the expression X_{n-k+1:n} n^{-1-(m-j)} k^{m-j} and the conditions E(h,j,m) and Q(h,j,m) are hard to parse; the exponents should be typeset explicitly so the reader can verify the convergence conditions.
- [Section 5.4 and Section 6] The references to 'Table 2b-(a)' and 'Table 2b-(b)' are confusing; the table should be given distinct sublabels or separate tables. There are also small typos, such as 'well-kown' in Section 2.1 and 'Kologorov-Smirnov' in Section 6, which should be corrected.
- [Section 4.2] The constrained test statistic T^{G+}_{m,ℓ,p} is introduced with notation that is not fully defined before first use; please state explicitly that ℓ is the number of selected indices j_k and that the indices are chosen so that the corresponding L-estimators converge under Theorem 1.
Circularity Check
No circular derivation: the central inequality is a parameter-free Jensen bound derived in the text, and the test construction does not fit or reuse its own output. A minor self-citation is not load-bearing; a separate proof gap in Theorem 2 is a correctness concern, not circularity.
full rationale
The paper's central claim rests on Proposition 1: if G^{-1}∘F is convex then P(X ≤ μ_{j:m}) ≤ G(E(G^{-1}(B_{j:m}))). The paper itself derives this from Jensen's inequality in Section 2.2 ('Jensen's inequality implies that EX_{j:m} ≤ F^{-1}∘G(E(G^{-1}∘F(X_{j:m}))) = F^{-1}∘G(E(G^{-1}(B_{j:m})))'), then cites Arab et al. (2025) only as a summary of the property. This is a parameter-free bound whose assumptions do not include the test's conclusions, so the self-citation is not load-bearing. The test statistics T^{G+}_{m,p} and T^{G-}_{m,p} compare empirical non-exceedance probabilities with these fixed bounds; no parameter is fitted to data and then renamed as a prediction, and the consistency proofs (Propositions 5 and 6) use standard Glivenko-Cantelli and L-estimator convergence arguments rather than assuming the claimed result. Lemma 1 and Theorem 2 attempt a stochastic-monotonicity proof of unbiasedness and monotone power; here the text contains a genuine correctness gap: the proof of Theorem 2 invokes Theorem 1.A.3 of Shaked and Shantikumar (2007) while stating only componentwise stochastic order, omitting that theorem's independence condition, which the dependent coordinates eU_j and eZ_j (all functions of the same empirical CDF) do not satisfy. This is repairable via the coupling in Lemma 1, but as written it is a missing-support issue, not a circular equation. Similarly, the finite-sample properties of the constrained statistics T^{G+}_{m,ℓ,p} are asserted without proof ('It is easy to check...'), another omitted-support issue. No step in the derivation reduces a claimed prediction to an input by construction, so the circularity score is 0.
Assumptions & free parameters
free parameters (3)
- m (number of order statistics in test statistic) =
user-selected; simulations suggest 10-20% of n for IHR/IOR, larger for heavy-tailed DOR
- p (norm order) =
p=1 chosen after simulations (Table 1)
- ℓ (number of order statistics used in constrained heavy-tailed test) =
user-selected, e.g., ℓ=m-2 for IOR; for DOR m=25,ℓ=5 etc.
assumptions (5)
- domain assumption G is absolutely continuous with known density g
- standard math F is absolutely continuous or has a density a.e. under the alternative
- standard math Jensen's inequality applied to convex/concave transforms to obtain Proposition 1 bounds
- standard math Mason (1982) strong laws for linear functions of order statistics
- domain assumption F ∈ D+(Φα) or D-(Φβ) domain of attraction conditions for heavy-tailed consistency
Cite this review
Pith. "Pith review of A new class of tests for convex-ordered families based on expected order statistics." pith.science (2026). https://pith.science/paper/MLFYK2MA
@misc{pith2026250114075,
author = {Pith},
title = {Pith review of: A new class of tests for convex-ordered families based on expected order statistics},
year = {2026},
howpublished = {\url{https://pith.science/paper/MLFYK2MA}},
note = {Machine review of arXiv:2501.14075}
}
abstract
Consider a pair of cumulative distribution functions $F$ and $G$, where $F$ is unknown and $G$ is a known reference distribution. Given a sample from $F$, we propose tests to detect the convexity or the concavity of $G^{-1}\circ F$ versus equality in distribution (up to location and scale transformations). This framework encompasses well-known cases, including increasing hazard rate distributions, as well as some other relevant families that have garnered attention more recently, for which no tests are currently available. We introduce test statistics based on the estimated probability that the random variable of interest does not exceed a given expected order statistic, which, in turn, is estimated via L-estimation. The tests are unbiased, consistent, and exhibit monotone power with respect to the convex transform order. To ensure consistency, we show that our L-estimators satisfy a strong law of large numbers, even when the mean is not finite, thereby making the tests suitable for heavy-tailed distributions. Unlike other approaches, these tests are broadly applicable, regardless of the choice of $G$ and without support restrictions. The performance of the method under various conditions is demonstrated via simulations, and its applicability is illustrated through a concrete example.
Figures
Figures from the paper (2 more)
Reference graph
Works this paper leans on
-
[1]
author Anis, M. , year 2013 . title A family of tests for exponentiality against ifr alternatives . journal Journal of Statistical Planning and Inference volume 143 , pages 1409--1415
work page 2013
-
[2]
author Anis, M. , author Basu, K. , year 2014 . title Tests for exponentiality against NBUE alternatives: a M onte C arlo comparison . journal Journal of Statistical Computation and Simulation volume 84 , pages 231--247
work page 2014
-
[3]
author Arab, I. , author Lando, T. , author Oliveira, P. , year 2025 . title Inequalities and bounds for expected order statistics from convex ordered families . journal Journal of Applied Probability , pages 1--19
work page 2025
-
[4]
author Barlow, R.E. , author Marshall, A.W. , author Proschan, F. , year 1963 . title Properties of probability distributions with monotone hazard rate . journal The Annals of Mathematical Statistics , pages 375--389
work page 1963
-
[5]
author Barlow, R.E. , author Van Zwet, W.R. , year 1969 . title Asymptotic properties of isotonic estimators for the generalized failure rate function. P art 1: strong consistency . type Technical Report . California Univ Berkeley O perations R esearch center
work page 1969
-
[6]
author Beare, B.K. , year 2021 . title Least favourability of the uniform distribution for tests of the concavity of a distribution function . journal Stat volume 10 , pages e376
work page 2021
-
[7]
author Bickel, P.J. , author Doksum, K.A. , year 1969 . title Tests for monotone failure rate based on normalized spacings . journal The Annals of Mathematical Statistics , pages 1216--1235
work page 1969
-
[8]
author Billingsley, P. , year 2013 . title Convergence of probability measures . publisher John Wiley & Sons
work page 2013
Show all 40 references
-
[9]
, author Goldie, C.M
author Bingham, N.H. , author Goldie, C.M. , author Teugels, J.L. , year 1989 . title Regular variation . volume volume 27 . publisher Cambridge university press
1989
-
[10]
, author Savits, T.H
author Block, H.W. , author Savits, T.H. , author Singh, H. , year 1998 . title The reversed hazard rate function . journal Probability in the Engineering and informational Sciences volume 12 , pages 69--90
1998
-
[11]
, year 1974
author Bryson, M.C. , year 1974 . title Heavy-tailed distributions: properties and tests . journal Technometrics volume 16 , pages 61--68
1974
-
[12]
, year 2002
author Carolan, C.A. , year 2002 . title The least concave majorant of the empirical distribution function . journal Canadian Journal of Statistics volume 30 , pages 317--328
2002
-
[13]
, author Embrechts, P
author Chen, Y. , author Embrechts, P. , author Wang, R. , year 2024 . title An unexpected stochastic dominance: Pareto distributions, dependence, and diversification . journal Operations Research
2024
-
[14]
, author Shneer, S
author Chen, Y. , author Shneer, S. , year 2024 . title Stochastic dominance for super heavy-tailed random variables . journal arXiv preprint arXiv:2408.15033
2024 arXiv
-
[15]
, author Heckman, N
author Gijbels, I. , author Heckman, N. , year 2004 . title Nonparametric testing for a monotone hazard function via normalized spacings . journal Journal of Nonparametric Statistics volume 16 , pages 463--477
2004
-
[16]
, year 1956
author Grenander, U. , year 1956 . title On the theory of mortality measurement: part II . journal Scandinavian Actuarial Journal volume 1956 , pages 125--153
1956
-
[17]
, author Jongbloed, G
author Groeneboom, P. , author Jongbloed, G. , year 2012 . title Isotonic L _2 -projection test for local monotonicity of a hazard . journal Journal of Statistical Planning and Inference volume 142 , pages 1644--1658
2012
-
[18]
, author Jongbloed, G
author Groeneboom, P. , author Jongbloed, G. , year 2014 . title Nonparametric estimation under shape constraints . number 38 , publisher Cambridge University Press
2014
-
[19]
, author Ferreira, A
author de Haan, L. , author Ferreira, A. , year 2006 . title Extreme value theory: an introduction . volume volume 3 . publisher Springer
2006
-
[20]
, author Van Keilegom, I
author Hall, P. , author Van Keilegom, I. , year 2005 . title Testing for monotone increasing hazard rate . journal Annals of Statistics volume 33 , pages 1109--1137
2005
-
[21]
, year 2023
author Lando, T. , year 2023 . title Testing departures from the increasing hazard rate property . journal Statistics & Probability Letters volume 193 , pages 109736
2023
-
[22]
, author Arab, I
author Lando, T. , author Arab, I. , author Eduardo Oliveira, P. , year 2024 . title Nonparametric inference about increasing odds rate distributions . journal Journal of Nonparametric Statistics volume 36 , pages 435--454
2024
-
[23]
, author Arab, I
author Lando, T. , author Arab, I. , author Oliveira, P.E. , year 2022 . title Properties of increasing odds rate distributions with a statistical application . journal Journal of Statistical Planning and Inference volume 221 , pages 313--325
2022
-
[24]
, author Arab, I
author Lando, T. , author Arab, I. , author Oliveira, P.E. , year 2023 . title Transform orders and stochastic monotonicity of statistical functionals . journal Scandinavian Journal of Statistics volume 50 , pages 1183--1200
2023
-
[25]
, author Olkin, I
author Marshall, A.W. , author Olkin, I. , year 2007 . title Life Distributions . publisher Springer, New York
2007
-
[26]
, year 1982
author Mason, D.M. , year 1982 . title Some characterizations of strong laws for linear functions of order statistics . journal The Annals of Probability , pages 1051--1057
1982
-
[27]
, author Anis, M
author Mitra, M. , author Anis, M. , year 2008 . title An L -statistic approach to a test of exponentiality against IFR alternatives . journal Journal of Statistical Planning and Inference volume 138 , pages 3144--3148
2008
-
[28]
, year 2024
author M \"u ller, A. , year 2024 . title Some remarks on the effect of risk sharing and diversification for infinite mean risks . journal arXiv:2411.10139
2024
-
[29]
, year 1981
author Oja, H. , year 1981 . title On location, scale, skewness and kurtosis of univariate distributions . journal Scandinavian Journal of statistics , pages 154--168
1981
-
[30]
, author Pyke, R
author Proschan, F. , author Pyke, R. , year 1967 . title Tests for monotone failure rate , in: booktitle Fifth Berkley Symposium , pp. pages 293--313
1967
-
[31]
, year 2008
author Resnick, S.I. , year 2008 . title Extreme values, regular variation, and point processes . volume volume 4 . publisher Springer Science & Business Media
2008
-
[32]
, author Wright, F.T
author Robertson, T. , author Wright, F.T. , author Dykstra, R.L. , year 1988 . title Order restricted statistical inference . publisher Wiley
1988
-
[33]
, year 2009
author Serfling, R.J. , year 2009 . title Approximation theorems of mathematical statistics . publisher John Wiley & Sons
2009
-
[34]
, author Shantikumar, J.G
author Shaked, M. , author Shantikumar, J.G. , year 2007 . title Stochastic Orders . publisher Springer, New York
2007
-
[35]
, year 2017 a
author Staudte, R.G. , year 2017 a. title Inference for quantile measures of kurtosis, peakedness, and tail weight . journal Communications in Statistics-Theory and Methods volume 46 , pages 3148--3163
2017
-
[36]
, year 2017 b
author Staudte, R.G. , year 2017 b. title The shapes of things to come: Probability density quantiles . journal Statistics volume 51 , pages 782--800
2017
-
[37]
, year 1964
author Van Zwet, W.R. , year 1964 . title Convex transformations of random variables . journal MC Tracts
1964
-
[38]
, year 1980
author Van Zwet, W.R. , year 1980 . title A strong law for linear functions of order statistics . journal The Annals of Probability volume 5 , pages 986--990
1980
-
[39]
, year 1977
author Wellner, J.A. , year 1977 . title A G livenko- C antelli theorem and strong laws of large numbers for functions of order statistics . journal The Annals of Statistics volume 5 , pages 473--480
1977
-
[40]
, author Wang, Y
author Zimmer, W.J. , author Wang, Y. , author Pathak, P.K. , year 1998 . title Log-odds rate and monotone log-odds rate distributions . journal Journal of quality technology volume 30 , pages 376--385
1998
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.