arxiv: 2604.05974 · v1 · submitted 2026-04-07 · 📊 stat.ME

Recognition: 2 theorem links

· Lean Theorem

Nonparametric Statistical Inference for Multivariate Niche Overlap

Jonas Beck , Solomon Harrar

Authors on Pith no claims yet

Pith reviewed 2026-05-10 18:33 UTC · model grok-4.3

classification 📊 stat.ME

keywords niche overlapnonparametric estimationmultivariate distributionsbootstrap inferenceecological statisticsspecies differentiationstatistical testing

0 comments

The pith

A nonparametric overlap index enables estimation and bootstrap inference for multivariate niche overlap.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper introduces a nonparametric index to quantify overlap between multivariate distributions describing species niches. It proposes estimators for the index and establishes their asymptotic properties under suitable conditions. Bootstrap-based procedures are developed to perform hypothesis tests and build simultaneous confidence intervals that remain reliable even with small samples. Numerical studies show the methods control type I error and maintain power across scenarios, while an application to fish stable-isotope data illustrates how the approach can distinguish species niches without parametric assumptions.

Core claim

We introduce a nonparametric overlap index for multivariate data and develop estimators whose asymptotic behavior is characterized. Bootstrap-based methods are proposed to enable statistical testing and simultaneous confidence intervals in small-sample settings, offering a robust alternative to parametric approaches in ecological studies.

What carries the argument

The nonparametric overlap index, a distribution-free measure of shared probability mass between multivariate distributions, equipped with consistent estimators and valid bootstrap inference procedures.

Load-bearing premise

The underlying multivariate distributions permit a well-defined overlap index and satisfy regularity conditions sufficient for the asymptotic results and bootstrap validity to hold.

What would settle it

A Monte Carlo experiment with known true overlap values in which the bootstrap confidence intervals fail to achieve nominal coverage or the tests exceed the nominal type I error rate.

Figures

Figures reproduced from arXiv: 2604.05974 by Jonas Beck, Solomon Harrar.

**Figure 2.** Figure 2: Empirical power for the tests in Section 2.4 for [PITH_FULL_IMAGE:figures/full_fig_p015_2.png] view at source ↗

**Figure 3.** Figure 3: Empirical power for the tests in Section 2.4 with a lognormal [PITH_FULL_IMAGE:figures/full_fig_p016_3.png] view at source ↗

**Figure 4.** Figure 4: 95% Confidence Intervals: Percentile-based, Normal-based, and [PITH_FULL_IMAGE:figures/full_fig_p018_4.png] view at source ↗

**Figure 5.** Figure 5: Empirical power for the tests from Section S1.4, based on [PITH_FULL_IMAGE:figures/full_fig_p041_5.png] view at source ↗

**Figure 6.** Figure 6: Empirical power for the tests from Section S1.4, for increasing [PITH_FULL_IMAGE:figures/full_fig_p042_6.png] view at source ↗

**Figure 7.** Figure 7: Pairwise Scatterplot Matrix of Fish Data [PITH_FULL_IMAGE:figures/full_fig_p043_7.png] view at source ↗

**Figure 8.** Figure 8: Chi-square Q-Q Plots of Mahalanobis Distances by Species [PITH_FULL_IMAGE:figures/full_fig_p044_8.png] view at source ↗

read the original abstract

In ecological studies niche overlap is often used to quantify species interaction and dynamics. This paper develops a robust, nonparametric statistical framework for quantifying and analyzing multivariate niche overlap. Parametric methods are often constrained by restrictive assumptions and tend to underperform in complex multivariate settings. We introduce a nonparametric overlap index and propose estimators for it. Further, we investigate asymptotic properties of the estimators. We also propose bootstrap-based inference procedures that enable statistical testing and simultaneous confidence intervals in small sample settings. Extensive numerical examples demonstrate that our proposed methods maintain correct size and exhibit robust power across various scenarios. We illustrate the practical utility of our methodology using stable isotope measurements from multiple fish species and provide distinct ecological insights regarding species niche differentiation.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper introduces a nonparametric multivariate niche overlap index with bootstrap inference that simulations and a fish isotope example support as practical.

read the letter

The core advance here is a nonparametric index for overlap in multivariate niche data, plus estimators whose asymptotics are worked out and bootstrap procedures for tests and simultaneous intervals that hold in small samples. That combination targets a real gap where parametric assumptions often fail in ecological settings with multiple dimensions like isotope ratios or abundances. The simulations check size and power across scenarios, and the stable isotope fish data shows it can surface species differentiation that matters for the application. Those pieces give the work concrete grounding rather than pure theory. The construction avoids obvious circularity or self-referential fitting, and the regularity conditions invoked are the usual ones for consistency and bootstrap validity. No load-bearing gaps jump out from the claims and checks provided. One softer spot is that direct head-to-head comparisons with existing multivariate or nonparametric overlap tools could be expanded to quantify the practical gain more sharply. The assumptions on the underlying distributions are standard but could be illustrated with a few edge-case simulations to reassure users. This is for ecologists and applied statisticians who need overlap measures on multivariate data without strong parametric restrictions. Readers working with small-sample niche studies or isotope data will find usable methods and code-level checks. The evidence from numerics and the real example is sufficient to warrant a serious referee rather than a desk reject, even if revisions will likely tighten the comparisons and assumption discussion.

Referee Report

0 major / 2 minor

Summary. The manuscript develops a nonparametric statistical framework for multivariate niche overlap in ecology. It introduces a nonparametric overlap index, proposes corresponding estimators, derives their asymptotic properties, and develops bootstrap-based procedures for hypothesis testing and simultaneous confidence intervals suitable for small samples. Simulations are used to demonstrate that the methods achieve correct size and robust power across scenarios, and the approach is illustrated with stable isotope measurements from multiple fish species to yield ecological insights on niche differentiation.

Significance. If the index definition, estimators, and bootstrap validity hold under the stated regularity conditions, the work supplies a flexible nonparametric alternative to parametric niche-overlap methods that frequently underperform in multivariate settings. The combination of asymptotic theory, small-sample bootstrap inference, simulation validation, and a real-data ecological application constitutes a practical contribution to statistical ecology.

minor comments (2)

Abstract: the claim of 'correct size and robust power' is stated without reference to the specific simulation designs or the form of the overlap index; adding one sentence on the index definition would improve informativeness while remaining within abstract length limits.
The manuscript would benefit from an explicit statement, early in the methods section, of the precise regularity conditions (e.g., smoothness, support positivity) required for the asymptotic normality and bootstrap consistency results.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive and accurate summary of our manuscript, which correctly identifies the nonparametric overlap index, asymptotic theory, bootstrap procedures, simulation studies, and ecological application. We appreciate the recommendation for minor revision and the recognition that the work provides a flexible alternative to parametric methods in multivariate settings.

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper defines a new nonparametric multivariate overlap index, constructs estimators from it, derives asymptotic consistency and normality under standard regularity conditions on the distributions, and develops bootstrap procedures for inference. These steps rely on external nonparametric statistical theory rather than reducing by construction to the paper's own fitted quantities, self-definitions, or self-citation chains. No load-bearing step renames a known result, smuggles an ansatz, or treats a fitted input as a prediction; the numerical simulations and real-data examples function as independent checks of the claimed properties.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The framework rests on standard nonparametric regularity conditions for asymptotics and bootstrap consistency plus the domain assumption that a meaningful multivariate overlap measure exists for the ecological distributions; no free parameters or invented entities are mentioned.

axioms (1)

domain assumption The multivariate distributions admit a well-defined overlap index and satisfy regularity conditions for asymptotic normality and bootstrap validity.
Required for the proposed estimators and inference procedures to have the claimed properties.

pith-pipeline@v0.9.0 · 5402 in / 1165 out tokens · 99116 ms · 2026-05-10T18:33:53.479128+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We introduce a nonparametric overlap index ... I(H(s),F(s)i)=P(X(s)i,low<Z(s)<X(s)i,up) ... rewritten as integrals of H∘F−1 ... Hadamard-differentiable ... functional delta method
IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

bootstrap-based inference ... asymptotic normality with covariance Σ

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

13 extracted references · 13 canonical work pages

[1]

(2003).An Introduction to Multivariate Statistical Analysis

Anderson, T. (2003).An Introduction to Multivariate Statistical Analysis. Wiley Series in Probability and Statistics. Wiley. Audia, P. G., Freeman, J., and Reynolds, P. D. (2006). Organizational found- ings in community context: Instruments manufacturers and their inter- relationship with other organizations.Administrative Science Quarterly, 51:381 –

work page 2003
[2]

B., and Bathke, A

Beck, J., Langthaler, P. B., and Bathke, A. C. (2025). Combining stochastic tendency and distribution overlap towards improved nonparametric effect measures and inference.Scandinavian Journal of Statistics, 52(3):1138–

work page 2025
[3]

B., Maitner, B., Harris, D

Blonder, B., Morrow, C. B., Maitner, B., Harris, D. J., Lamanna, C., Violle, C., Enquist, B. J., and Kerkhoff, A. J. (2018). New approaches for de- lineating n-dimensional hypervolumes.Methods in Ecology and Evolution, 9(2):305–319. Box, G. E. P. (1954). Some Theorems on Quadratic Forms Applied in the Study of Analysis of Variance Problems, I. Effect of I...

work page 2018
[4]

and Hettmansperger, T

Brown, B. and Hettmansperger, T. (2002). Kruskal–Wallis, multiple com- parisons and efron dice.Australian & New Zealand Journal of Statistics, 44(4):427–438. 20 Brunner, E., Bathke, A. C., and Konietschke, F. (2018).Rank and pseudo- rank procedures for independent observations in factorial designs. Springer, Cham, Switzerland. Brunner, E., Dette, H., and ...

work page arXiv 2002
[5]

Hutchinson, G. E. (1957). Population studies: Animal ecology and demography—concluding remarks. InCold Spring Harbor Symposia on Quantitative Biology, volume 22, pages 415–427. Cold Spring Harbor Lab- oratory Press. Junker, R. R., Kuppler, J., Bathke, A. C., Schreyer, M. L., and Trutschnig, W. (2016). Dynamic range boxes – a robust nonparametric approach ...

work page 1957
[6]

B., Gladow, K.-P., Kr¨ uger, O., and Beck, J

Langthaler, P. B., Gladow, K.-P., Kr¨ uger, O., and Beck, J. (2024). A novel method for nonparametric statistical inference for niche overlap in multiple species.Biometrical Journal, 66(7):e202400013. 22 Marcus, R., Eric, P., and Gabriel, K. R. (1976). On closed testing proce- dures with special reference to ordered analysis of variance.Biometrika, 63(3):...

work page 2024
[7]

S1.1 Statistical Model Let the independent and identically distributed random vectors Xk = (X (1) k ,

in the multivariate case for two-samples. S1.1 Statistical Model Let the independent and identically distributed random vectors Xk = (X (1) k , . . . , X(d) k )⊤ ∼F, fork= 1, . . . , nrepresent the data for the first sample, whereX (s) k denotes the observation on thes-th endpoint of thek-th subject in the first sample. Similarly, let Yk = (Y (1) k , . . ...

work page 2018
[8]

The first term equation shows the close relation to the ROC-curve

as I(F (s), G(s)) = Z 1 0 F (s) ◦G (s)−1 (1−α/2)dα− Z 1 0 F (s) ◦G (s)−1 (α/2)dα = 2 Z ∞ G(s)−1(1/2) F (s)dG(s) − Z G(s)−1 (1/2) −∞ F (s)dG(s) . The first term equation shows the close relation to the ROC-curve. Note thatF (s) =G (s) implies the corresponding overlap indexI(F (s), G(s)) = 1/2. However, the converse is not true. In addition, the overlap in...

work page 2018
[9]

, dandσ (s)∗ is the bootstrapped standard deviation of thes-th component

σ(s)∗ fors= 1, . . . , dandσ (s)∗ is the bootstrapped standard deviation of thes-th component. By Theorem S1.1, under the null hypothesis,Thas asymptotic multivariate normal distribution with mean vector0and covariance matrix P ∗, whereP ∗ is the correlation matrix of the empirical bootstrap sample. According to the so-called Max-T test (Konietschke et al...

work page 2012
[10]

,|Td|)≥z(1−α, P ∗), wherez(1−α, P ∗) the equicoordinate quantiles

the null hypothesis (S1.3) may be rejected at levelα, if T0 = max(|T1|, . . . ,|Td|)≥z(1−α, P ∗), wherez(1−α, P ∗) the equicoordinate quantiles. The quantiles can computed with theR-packagemvtnorm(Genz et al., 2021; Genz and Bretz, 2009). This test, however, tends to lack power as shown later in the Simulation section. Notwithstanding this limitation, an ...

work page 2021
[11]

, d, and this holds uniformly over any subin- terval [a, b] of [0,1] (Hsieh and Turnbull, 1996, Theorem 2.2)

almost surely for alls= 1, . . . , d, and this holds uniformly over any subin- terval [a, b] of [0,1] (Hsieh and Turnbull, 1996, Theorem 2.2). Sinceaandb are selected without restriction, the relationship remains valid over all such subintervals. LetD=D[0,1] be the Skorokhod space, i.e., the space of all C` adl` ag functions on [0,1], andE⊂R. Consider the...

work page 1996
[12]

and bIn,m(F,G) = Φ( ˆF (1) ◦( ˆG(1))−1,

Results are shown for dimensionsd= 2 andd= 5, as described in Example S2.5. and bIn,m(F,G) = Φ( ˆF (1) ◦( ˆG(1))−1, . . . , ˆF (d) ◦( ˆG(d))−1) Sinceϕis continuously Hadamard-differentiable (van der Vaart, 1998, Theorem 20.9), so is Φ by chain rule. Then, the desired result follows by the Cr` amer-Wold device and the delta method for empirical processes (...

work page 1998
[13]

for distribution functionsFandGdefined onR

Results are shown for d= 2 andd= 5, as described in Example S2.6. for distribution functionsFandGdefined onR. This transformation is Hadamard-differentiable tangentially toD[a, b]×C[a, b] (van der Vaart and Wellner, 2023, Comment 4 in Section 3.10). Then the compositionϕ◦ψ, whereϕis defined as in (S3.2), is Hadamard-differentiable (van der Vaart, 1998, Th...

work page 2023