Alphabeta Math
Pipeline-generated
How statement and proof provenance work

The first chip identifies the source of the statement or construction; the second identifies the source of its local proof or verification.

  • Literature-sourced: the exact statement appears in a cited source; only wording and notation differ.
  • AI-adapted: a semantically identical restatement of literature-sourced material, modulo indexing, notation, and boundary cases adopted by the library.
  • AI-generated: a genuinely novel statement formulated by AI, with no source for the claim itself.

These labels describe origin, not correctness: citations and verification chips remain separate evidence.

✓ 13 results · all verified · 10 also independently AI-judged
Every result on this page is machine-checked by a proof checker and read in full and owner-audited; the judge is an additional, independent cross-model AI review of the proofs. The 3 not AI-judged were verified by owner audit (typically over a confirmed judge false positive), not failures.

Studentized permutation tests under weak parameter nulls

1 · Prerequisites

2 · Summary

This page develops two-sample order-(r,r) U-statistics, their first projections, and second-moment control of the higher-order remainder. A finite permutation variance estimate strengthens the exact triangular-mixture coupling and proves the conditional permutation CLT under the weak parameter null, block antisymmetry, mixed second moments and positive pooled projection variance. No assignment-wise centering or sample-fraction convergence rate is required. The all-assignment centered theorems remain useful special cases.

The studentized pivot divides by the estimated standard error r{σ^m2/m+σ^n2/n}1/2; that ratio itself has the standard-normal limit. General joint variance estimates require direct permutation consistency. For estimates depending on their own group alone, IID triangular-mixture consistency is sufficient by the composition-count transfer lemma. The finite-support rank example verifies permutation consistency of its joint estimates directly.

When P=Q, finite-group invariance gives finite-sample exactness for every measurable studentized statistic, with the stated tie conventions. That exactness does not establish validity of unstudentized Wilcoxon statistics under a weak mean or median null.

The explicit finite-enumeration test recomputes both variance estimates on every permuted split and supplies inclusive or tie-randomized p-values. The heterogeneous-law Wilcoxon specialization uses tie-corrected pairwise comparisons, verifies its empirical projection variance estimates under both the true samples and permutation, and tests the probabilistic-index null. It does not replace that null with equality of means or medians.

3 · Logical flowchart

4 · Definitions, theorems and proofs

DefinitionDefinition: AI-adaptedProof: Not applicablejudge pass (gpt-6-luna)Open item page →

Two-sample U-statistic parameter

Definition

Let (X,A) be a measurable space, let P,Q be probability measures on it, and let r≥1 be a fixed integer. Give X2r the product sigma-algebra and the product probability measure P⊗r⊗Q⊗r. A two-sample order-(r,r) kernel is a measurable function h:X2r→R that is symmetric in its first r arguments and separately symmetric in its last r arguments. If h is integrable under this product law, its two-sample U-statistic parameter is

θ(P,Q;h):=∫X2rh d(P⊗r⊗Q⊗r).

The value is finite and real. This definition imposes no null hypothesis and no antisymmetry condition; those are hypotheses of later results.

LemmaStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6-luna)Open item page →

Conditional CLT for a random subset of an IID finite population

Statement

Assume the Axiom of Choice. For each N, let WN,1,…,WN,N be IID real variables with mean 0, variance σN2→σ2>0, and uniformly integrable squares:

lim⁡K→∞sup⁡NE ⁣[WN,121{∣WN,1∣>K}]=0.

Let integer sample sizes mN,nN satisfy mN+nN=N and mN/N→p∈(0,1), and put q=1−p. It is enough to consider sufficiently large N, for which N>1 and 1≤mN,nN≤N−1; for each such N, choose a uniformly random mN-element subset AN⊆{1,…,N} independently of the row. Set

WˉN=N−1∑i=1NWN,i,SN2=(N−1)−1∑i=1N(WN,i−WˉN)2.

Define

TN=∑i∈AN(WN,i−WˉN)mNnNSN2/N,

and set TN=0 if SN2=0. Then the conditional distribution of TN given the finite population converges in probability to N(0,1), uniformly in its real argument.

Facts & Assumptions

Given: IID rows with mean zero, variances σN2→σ2>0, uniformly integrable squares, and a uniform sample without replacement of size mN with mN/N→p∈(0,1).

[F1]

Martingale central limit theorem: a square-integrable martingale-difference array whose predictable variance clock tends to 1 in probability and satisfies Lindeberg has a standard-normal terminal limit, assuming AC.

[F2]

Conditional Hoeffding bound for bounded martingale differences: if a conditionally centered increment has conditional range of width at most c, its conditional exponential moment is at most eλ2c2/8.

[F3]

Square-integrable martingale-difference array and variance clock: for martingale differences Dk, the predictable clock is ∑kE[Dk2∣Fk−1].

[F4]

Simple random sampling without replacement: every subset of the given size is sampled with equal probability.

[F5]

SRSWOR inclusion covariance and sample-mean variance: for fixed population values with finite-population variance SN2, a size-mN sample sum has variance mNnNSN2/N.

[F6]

Independent random elements: the random subset is independent of the IID row, so conditional on the row its design law is uniform SRSWOR.

[F7]

Tower property of conditional expectation: iterated conditioning preserves conditional-expectation identities.

[A1]

The Axiom of Choice is assumed exactly for the martingale CLT and the conditional-expectation interfaces in [F1]–[F3]. The Axiom of Choice

Proof

Given: The hypotheses in the Statement. Prove conditional normal convergence for the displayed random-subset statistic.

1.1givenalgebra

Since EWN,1=0 and σN2 is bounded, Chebyshev gives WˉN→0 in probability. To show N−1∑iWN,i2−σN2→0, first truncate W2 at a fixed K: the truncated variables are bounded, so Chebyshev gives convergence of their row average to its expectation. The expected tail average is at most sup⁡NE[WN,121{∣WN,1∣>K}], which tends to zero as K→∞ by uniform integrability. Hence the untruncated average converges to σ2, and SN2→σ2.

1.2F3F4F6givenconstruct

Conditional on the row, list its centered values in a uniformly random order aN,π(1),…,aN,π(N) and write Sj=∑i≤jaN,π(i). For j<N, E[aN,π(j+1)∣Fj]=−Sj/(N−j), since the remaining values sum to −Sj. Therefore Mj:=Sj/(N−j) is a martingale for j≤mN, and its increments are bounded by 2max⁡i∣aN,i∣/(N−j−1). For j≤mN−1, the denominator is at least nN, which is bounded below by a positive multiple of N.

2.1givenstep 1.1algebra

For every ε>0, a union bound and W2≥ε2N on {∣W∣>εN} give P ⁣(max⁡i≤N∣WN,i∣>εN)≤ε−2E ⁣[WN,121{∣WN,1∣>εN}]⟶0. Together with WˉN=OP(1), this implies max⁡i∣aN,i∣/N→0 for aN,i=WN,i−WˉN.

2.2F2F7givenstep 1.1step 1.2algebra

First consider any fixed finite population c1,…,cN with ∑ici=0 and ∣ci∣≤L. For its random order, the martingale Sj/(N−j) from step 1.2 has increments of absolute value at most CpL/N through j≤mN, for all sufficiently large N. Iterating the conditional Hoeffding bound [F2] with [F7] gives, for each fixed u>0, Pπ(∣Sj∣/N>u)≤2exp⁡(−cpu2N/L2) uniformly in j≤mN; here ∣Sj∣/N>u implies ∣Sj/(N−j)∣>u because N−j≤N. A union bound over at most N prefixes makes the maximum divided by N tend to zero. Apply this argument first to the clipped values WN,i(K)−W‾N(K), and then to the centered clipped squares (WN,i(K)−W‾N(K))2−(WN(K)−W‾N(K))2‾. For fixed K, both are bounded arrays. The average absolute first-power truncation error is at most twice the average of ∣WN,i∣1{∣WN,i∣>K}; its expectation is at most 2K−1sup⁡NE[WN,121{∣WN,1∣>K}]. The average absolute difference of the centered squares tends to zero in probability as K→∞: Cauchy--Schwarz bounds it by the square root of the average squared truncation error times the square root of the average (∣aN,i∣+∣WN,i(K)−W‾N(K)∣)2, whose first factor is small by uniform integrability and whose second is bounded in probability by step 1.1. Thus, taking N→∞ and then K→∞, we obtain both max⁡j≤mN∣Sj∣N→0,max⁡j≤mN∣∑i≤j(aN,π(i)2−aN2‾)∣N→0 in probability.

2.3F3F4step 1.2algebra

Let RN,k be the mean of the values remaining before draw k, let sN,k−12 be their finite-population variance, and let FN,k reveal the first k draws. The Doob martingale for the final sample sum has increments DN,k=nNN−k(aN,π(k)−RN,k),1≤k≤mN. Indeed, after draw k the conditional expected contribution of the remaining mN−k draws is (mN−k) times the new remaining mean; subtracting the previous conditional expectation gives the displayed coefficient. The increments are centered and square-integrable.

3.1step 1.1step 2.2algebra

Before draw j+1, the remaining set has size N−j≥nN≍N. Its mean is −Sj/(N−j), and its raw second moment is aN2‾−1N−j∑i≤j(aN,π(i)2−aN2‾). Both differ uniformly from 0 and aN2‾, respectively, by oP(1) by step 2.2. The remaining finite-population variance therefore differs uniformly from SN2=NaN2‾/(N−1) by oP(1); the factors (N−j)/(N−j−1) tend uniformly to 1 because N−j≥nN≍N. Hence max⁡j≤mN∣sN,j2−SN2∣→0.

4.1F3F5step 1.1step 3.1algebra

A uniform draw from the N−k+1 remaining values has variance ((N−k)/(N−k+1))sN,k−12. On {SN2>0}, the predictable clock of ∑k≤mNDN,k/mNnNSN2/N is 1mNnNSN2/N∑k=1mNnN2sN,k−12(N−k)(N−k+1). Define the normalized increments and their clock to be zero on {SN2=0}. By step 1.1, SN2→Pσ2>0, so this zero-variance event has probability tending to zero. On its complement, step 3.1 permits replacing each sN,k−12 uniformly by SN2 with relative oP(1) error. The remaining deterministic sum telescopes: nN2∑k=1mN1(N−k)(N−k+1)=nN2(1nN−1N)=mNnNN. Thus the defined normalized clock tends to 1 in probability, in agreement with the finite-population variance in [F5]. The normalized terminal sum equals TN also on {SN2=0} by the statement's convention.

5.1givenstep 1.1step 2.1step 4.1algebra

On {SN2>0}, since nN/(N−k)≤1 and ∣RN,k∣≤max⁡i∣aN,i∣, step 2.1 gives max⁡k≤mN∣DN,k∣mNnNSN2/N≤2max⁡i∣aN,i∣mNnNSN2/N⟶0 in probability. Under the zero-variance convention of step 4.1, the normalized increments are zero on {SN2=0} as well; since SN2→Pσ2>0, the maximum normalized increment tends to zero in probability. Therefore, for each fixed η>0, every term in the normalized martingale Lindeberg sum is zero with probability tending to one.

6.1A1F1step 1.1step 2.1step 4.1step 5.1algebra

Suppose conditional convergence failed along a subsequence. From steps 1.1 and 2.1, choose a further subsequence on which SN2→σ2>0 and max⁡i∣aN,i∣/N→0 almost surely. For each positive rational η, the probability that the normalized clock differs from 1 by more than η tends to zero by step 4.1. Its conditional probability given the row therefore tends to zero in probability by Markov's inequality. A diagonal further subsequence makes all those conditional probabilities tend to zero almost surely. For almost every resulting sequence of rows, the conditional clock converges in probability to 1; the normalized increments have a deterministic uniform bound tending to zero by step 5.1, so for every fixed η>0 the Lindeberg sum on that row's permutation space is exactly zero for all sufficiently large rows. The conditional permutation spaces may vary with N, whereas the array definition [F3] is stated on one fixed probability space. To justify using the martingale CLT in this setting, repeat the proof of [F1] with the row-specific conditional laws PNπ and expectations ENπ. Every conditional expectation and variance estimate is formed within one row, and no joint event involving different rows is used. In particular, the only bounded-convergence step in [F1]'s proof can be read rowwise: if ΓN→1 in probability under laws PN, then for bounded continuous f and any δ>0, ∣ENf(ΓN)−f(1)∣≤sup⁡∣x−1∣≤δ∣f(x)−f(1)∣+2∥f∥∞PN(∣ΓN−1∣>δ), which tends to zero after N→∞ and then δ↓0. All other steps of the proof of [F1], including the Lindeberg estimate and the predictable-clock truncation, are already rowwise; hence that characteristic-function proof yields the martingale CLT for these finite permutation rows. Pad each row with zero increments. The conditional CDF converges to Φ at every t, contradicting the assumed failure. Thus it converges in probability at each t.

7.1step 1.1step 6.1algebra∎

Finally, choose T so that Φ(−T) and 1−Φ(T) are small, and partition [−T,T] into finitely many points where adjacent Φ values differ by at most ε. Monotonicity of each conditional CDF bounds its deviation between grid points by the endpoint deviations plus ε. Pointwise convergence at the finite grid and the tail bounds yield uniform convergence in probability. The convention when SN2=0 affects an event whose probability tends to zero by step 1.1 and σ2>0.

LemmaStatement: Literature-sourcedProof: AI-adaptedprecheck passOpen item page →

Hypergeometric and binomial composition counts are contiguous

Statement

Assume AC. For each N, let MN of N objects be marked, with pN=MN/N→p∈(0,1), and let kN/N→s∈(0,1). Let BN∼Bin⁡(kN,pN) and let HN count marked objects in a uniform size-kN sample without replacement from the N objects. Then the law of HN is contiguous to the law of BN: for every sequence of sets AN, if P(BN∈AN)→0, then P(HN∈AN)→0.

More precisely, the likelihood ratio LN(j)=P(HN=j)/P(BN=j) satisfies

LN(BN)⇒(1−s)−1/2exp⁡ ⁣(−s2(1−s)Z2),Z∼N(0,1),

and this nonnegative limit has mean 1.

Facts & Assumptions

Given: AC, marked-object proportions pN→p∈(0,1), sampling proportions kN/N→s∈(0,1), a binomial count, and a uniform sample without replacement.

[F1]

Lindeberg-Feller central limit theorem: sufficiency gives a standard-normal limit for a centered independent triangular array under the Lindeberg condition, assuming AC.

[F2]

Stirling's formula for factorials states j!∼2πj(j/e)j as positive integers j→∞.

[F3]

Identical distribution and IID families defines IID as mutual independence together with a common law.

[F4]

Moments, variance, and covariance on a probability space defines variance as the centered second moment.

[F5]

For sigma-finite factors, the product measure exists, has the rectangle formula, is sigma-finite, and is unique supplies the finite product probability measure used for that Bernoulli row.

[F6]

Expectations factor over finite products of independent random variables factors expectations of products of integrable functions of independent variables.

[A1]

AC is used exactly through [F1]. The Axiom of Choice

Proof

Given: The hypotheses in the Statement. Prove the likelihood-ratio limit and the stated direction of contiguity.

1.1givenalgebra

The exact mass ratio, with qN=1−pN, is LN(j)=(MNj)(N−MNkN−j)(NkN)(kNj)pNjqNkN−j on the hypergeometric support, and is zero off that support.

1.2F3F5givenconstructalgebra

By [F5], form the kN-fold product of the law assigning masses pN and qN to 1 and 0. The rectangle formula makes its coordinate indicators IN,i independent with that common Bernoulli law, hence IID in the sense of [F3]. Each binary vector with j ones has mass pNjqNkN−j, and there are (kNj) such vectors. Their sum therefore has the specified binomial law, so we may realize BN=∑i=1kNIN,i on this space.

2.1F4F6step 1.2algebra

Each indicator has mean pN and variance pNqN by [F4]. Factorization [F6] gives E[IN,iIN,j]=pN2 for i≠j, so the row sum has mean kNpN and variance kNpNqN.

3.1A1F1step 2.1givenalgebra

Each centered summand has magnitude at most 1, while the row standard deviation is kNpNqN≍N. Thus for each fixed ε>0 all Lindeberg tail terms are zero for large N; under AC [A1], [F1] gives ZN:=BN−kNpNkNpNqN⇒Z∼N(0,1).

4.1F2step 1.1step 3.1algebra

Put sN=kN/N and δN=BN−kNpN. By [step 3.1], δN=OP(N), so the event ∣δN∣≤N2/3 has probability tending to one. On this event, all factorial arguments in [step 1.1] are bounded below by a positive multiple of N, since pN→p and sN→s with p,s∈(0,1). Let H(a,x)=alog⁡a−xlog⁡x−(a−x)log⁡(a−x) and H(t)=H(1,t). Stirling's formula [F2], uniformly for these arguments, expands each log-binomial coefficient as H(a,x)−12log⁡(2πx(a−x)/a)+o(1). The order-N terms at δ=0 cancel because H(pNN,sNpNN)=pNNH(sN), H(qNN,sNqNN)=qNNH(sN), and H(sNN,sNpNN)=sNNH(pN). For FN(δ)=H(pNN,sNpNN+δ)+H(qNN,sNqNN−δ)−H(N,sNN)−H(sNN,sNpNN+δ)−(sNpNN+δ)log⁡pN−(sNqNN−δ)log⁡qN, use Ha′(x)=log⁡((a−x)/x) and Ha′′(x)=−1/x−1/(a−x); direct substitution gives FN(0)=0, FN′(0)=0 and FN′′(0)=−1pNNsN(1−sN)−1qNNsN(1−sN)+1sNNpNqN=−1pNqNN(1−sN). The third derivative is bounded by C/N2 throughout this event because every argument remains at least cN. The square-root terms at δ=0 sum to −12log⁡(2πpNNsN(1−sN))(2πqNNsN(1−sN))(2πNsN(1−sN))(2πsNNpNqN)=−12log⁡(1−sN), and their change is O(∣δN∣/N)=oP(1). Taylor's formula therefore gives log⁡LN(BN)=−12log⁡(1−sN)−δN22pNqNN(1−sN)+OP(∣δN∣3/N2)+oP(1)=−12log⁡(1−sN)−sN2(1−sN)ZN2+oP(1). The Taylor remainder is oP(1) because δN=OP(N), and the omitted event has probability tending to zero by [step 3.1].

5.1step 3.1step 4.1algebra

Since sN→s, [step 3.1] and [step 4.1] imply the displayed likelihood-ratio limit. Its mean is (1−s)−1/2Eexp⁡ ⁣(−s2(1−s)Z2)=(1−s)−1/21−s=1, by direct integration against the standard-normal density.

6.1step 5.1algebra

Each LN(BN)≥0 and E[LN(BN)]=1 by the likelihood-ratio definition. Nonnegative convergence in distribution to a variable with the same mean implies uniform integrability: for fixed K, E[LN∧K]→E[L∧K], so lim sup⁡NE[(LN−K)+]=1−E[L∧K]→0 as K→∞. This also controls E[LN1{LN>2K}].

7.1step 6.1algebra∎

If P(BN∈AN)→0, then for every K>0, P(HN∈AN)=E[LN(BN)1{BN∈AN}]≤KP(BN∈AN)+E[LN(BN)1{LN(BN)>K}]. The first term tends to zero for fixed K and the second is uniformly small for large K by step 6.1. Therefore P(HN∈AN)→0. This proves the required direction of contiguity; the reverse direction is not needed here.

LemmaStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6-luna)Open item page →

Finite-group rank validity on a measurable sample space

Statement

Let a finite group G act on a measurable space (S,S) by measurable maps, and let P be an invariant probability law: P(g−1A)=P(A) for every g∈G and A∈S. For a measurable statistic T:S→R, define

p0(x)=#{g∈G:T(g⋅x)≥T(x)}∣G∣,pU(x,u)=#{g:T(g⋅x)>T(x)}+u#{g:T(g⋅x)=T(x)}∣G∣.

If X∼P and U is independent uniform on [0,1], then p0(X) is superuniform and pU(X,U) is uniform on [0,1]. The sample space need not be finite; only the acting group is finite.

Facts & Assumptions

Given: A finite group action by measurable maps, an invariant probability law, a measurable real statistic, and an independent uniform tie randomizer.

[F1]

Group and abelian group and Left group actions, transitive actions, and faithful actions: G is a group acting with e⋅x=x and (gh)⋅x=g⋅(h⋅x); hence right multiplication by a fixed element is a bijection of G.

[F2]

A measurable function between measurable spaces: the action maps and T are measurable. The product sigma-algebra and its finite iterates makes the coordinate projections on S×[0,1] measurable, and Arithmetic and lattice operations preserve measurability whenever they are defined preserves measurability under the finite sums and products in the rank formulas.

[F3]

Tied upper ranks on a finite multiset: for a finite list of real numbers, an independent uniform index gives a superuniform inclusive upper rank and a uniform tie-randomized upper rank.

[F4]

Independent random elements: the uniform group element and the tie randomizer are independent of the data and of each other.

Proof

Given: The hypotheses in the Statement. Prove both rank conclusions.

1.1F1F3construct

Fix x∈S and form the finite list (T(g⋅x))g∈G. For any h∈G, the list (T(g⋅(h⋅x)))g∈G has exactly the same multiset of values, because right multiplication g↦gh is a bijection of G by [F1]. Therefore, when H is uniform on G, p0(H⋅x) and pU(H⋅x,U) are exactly the inclusive and tie-randomized ranks of a uniform index in this list.

1.2F2givenalgebra

For each g, x↦T(g⋅x) is measurable by [F2]. Its difference from T(x) is measurable, so the indicators of the strict and equal comparisons are measurable; finite sums and products with the second coordinate u remain measurable by [F2]. Thus p0 and the joint map (x,u)↦pU(x,u) are measurable.

1.3F4givenalgebra

Let H be uniform on G, independent of (X,U). For every measurable A, P(H⋅X∈A)=∣G∣−1∑h∈GP(X∈h−1A)=P(A) by invariance. Since H is independent of U, the pair (H⋅X,U) has the same law as (X,U).

2.1F3step 1.1step 1.2step 1.3algebra∎

For each fixed x, steps 1.1–1.2 and [F3] give PH(p0(H⋅x)≤t)≤t for 0≤t≤1, and pU(H⋅x,U) is uniform. Integrating these statements over X, then using equality in law from step 1.3, proves the two conclusions for p0(X) and pU(X,U). No choice principle is used: all group averages and ranks are finite.

LemmaStatement: AI-adaptedProof: AI-adaptedjudge pass (gpt-6-luna)Open item page →

A successive-revelation variance bound for a uniform permutation

Statement

Let π be uniform on the permutations of {1,…,N}, where N≥1, and let f be any real function on this finite set. Write π(k,l) for the permutation obtained by exchanging its entries in positions k,l, with (k,k) the identity. Then

Var⁡(f(π))≤12∑k=1N1N−k+1∑l=kNE[(f(π)−f(π(k,l)))2].

In particular, if every transposition difference has second moment at most b, then Var⁡(f(π))≤Nb/2. These are finite-sum statements and require no choice axiom.

Facts & Assumptions

Given: A positive integer N, a uniform permutation, and a real-valued function on its finite state space.

[F1]

The uniform probability space on a nonempty finite set assigns equal probability to every permutation.

[F2]

Moments, variance, and covariance on a probability space defines variance and covariance by centered products.

[F3]

The Cauchy-Schwarz inequality for finite sums implies that the square of an average of finitely many real numbers is at most the average of their squares.

Proof

Given: The hypotheses in the Statement. Prove the variance estimate by revealing successive entries.

1.1F1F2givenalgebra

Let Fk reveal the first k entries of π, and let Mk be the average of f over permutations agreeing with those entries. Thus M0=Ef, MN=f, and the average of Mk conditional on Fk−1 is Mk−1. For j<k, averaging (Mj−Mj−1)(Mk−Mk−1) first over completions of the first k−1 entries gives zero. Expanding the square of f−Ef=∑k(Mk−Mk−1) therefore gives Var⁡(f)=∑kE[(Mk−Mk−1)2]. All these conditional averages are finite sums.

2.1F1F2F3step 1.1algebra

Fix a prefix of length k−1 and put q=N−k+1. For each remaining image a, let ua be the average of f over completions whose kth entry is a. The conditional variance of Mk is 12q2∑a,b(ua−ub)2. Choose a uniform completion π of the prefix and independently choose l uniformly from {k,…,N}. Given π(k)=a, the image b=π(l) is uniform among the q remaining images. For fixed a,b, swapping position k with the location of b is a bijection from completions with kth image a onto completions with kth image b. Hence E[f(π)−f(π(k,l))∣a,b,Fk−1]=ua−ub. By [F3], the conditional variance just displayed is at most 12E[(f(π)−f(π(k,l)))2∣Fk−1].

3.1step 1.1step 2.1algebra∎

Average step 2.1 over prefixes and sum over k using step 1.1. The independent uniform choice of l gives the displayed estimate. If every squared difference has expectation at most b, each inner average is at most b, so the sum is at most Nb/2. For N=1 both sides of the displayed estimate are zero.

DefinitionDefinition: Literature-sourcedProof: Not applicablejudge pass (gpt-6-luna)Open item page →

Two-sample order-r U-statistic and first projections

Definition

Fix m,n≥r and independent samples X1,…,Xm∼P and Y1,…,Yn∼Q, each IID within its group. For a kernel from Two-sample U-statistic parameter, define

Um,n:=1(mr)(nr)∑I∈({1,…,m}r)∑J∈({1,…,n}r)h(XI;YJ).

Here I,J are unordered index subsets, and h(XI;YJ) means that the selected observations are listed in increasing index order in their respective blocks. Block symmetry makes the value independent of that ordering. Under E∣h∣2<∞ for the P⊗r⊗Q⊗r law, write

gX(X1):=E[h(X1,X2,…,Xr;Y1,…,Yr)∣X1]−θ(P,Q;h), gY(Y1):=E[h(X1,…,Xr;Y1,Y2,…,Yr)∣Y1]−θ(P,Q;h).

These are the centered first projections under P and Q; they are square-integrable equivalence classes, so changes on their respective null sets do not change any later claim.

Remarks

The displayed conditional-expectation notation is shorthand for the completed section-integral versions under the product laws. The corresponding functions are measurable and square-integrable by First projections are measurable and square-integrable ↗.

The finite sum defining Um,n is a measurable real statistic: the finite sample-vector coordinate maps are measurable by Coordinate maps and permutations are measurable on finite products, composition with h is measurable, and its finite sum and constant scaling remain measurable by Arithmetic and lattice operations preserve measurability whenever they are defined. This is a statistic in the sense of Statistics and estimators.

LemmaStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6-luna)Open item page →

Consistency transfers from IID mixture samples to permuted groups

Statement

Assume AC. Let P,Q be probability laws and let mN+nN=N, with pN=mN/N→p∈(0,1). Write PˉN=pNP+(1−pN)Q. For each kN∈{mN,nN}, let TN be a measurable statistic of kN observations and let cN be a deterministic target. If

TN(V1,…,VkN)−cN→0in probability when Vi are IID PˉN,

then TN applied to either group of a uniform random permutation of the pooled independent sample X1,…,XmN∼P, Y1,…,YnN∼Q also converges in probability to cN. The statement does not require pN−p=O(N−1/2) because the IID comparison law uses the exact weight pN.

Facts & Assumptions

Given: Independent IID samples from P and Q, a uniform permutation of their pooled N observations, and an estimator/statistic consistent for cN on IID PˉN samples of the same group size.

[F1]

Hypergeometric and binomial composition counts are contiguous: the hypergeometric label count in a size-kN permuted group is contiguous to the binomial label count in a size-kN IID mixture sample.

[F2]

Identical distribution and IID families: each source group is IID with its stated common law.

[F3]

Independent random elements: the source groups and the random subset selection are independent as specified.

[F4]

Statistics and estimators: a statistic is a measurable function of the ordered observations; an estimator is such a statistic with a specified target.

[A1]

AC is inherited from the likelihood-ratio contiguity result [F1]. The Axiom of Choice

Proof

Given: The assumptions in the Statement. Prove consistency for both permuted group sizes.

1.1F1F2F3

Fix kN equal to either mN or nN. In an IID sample from PˉN, the number BN of observations whose mixture component is P has law Bin⁡(kN,pN). In a uniformly chosen size-kN group from the pooled sample with exactly mN P-observations among N, the corresponding count HN is hypergeometric.

1.2F2F3F4

Conditional on either count being j, the ordered group has the same distribution: its P-positions are uniform among the j-subsets of the kN positions, the values at those positions are IID P, and the remaining values are IID Q. Therefore there is a function fN(j)∈[0,1] giving the conditional probability that ∣TN−cN∣>ε under either sampling scheme.

2.1A1F1step 1.2algebra

The IID mixture consistency assumption gives E[fN(BN)]→0. Markov's inequality implies fN(BN)→0 in probability under the binomial law. By contiguity [F1], fN(HN)→0 in probability under the hypergeometric law. Since 0≤fN≤1, this also gives E[fN(HN)]→0. AC is used through [F1] only.

3.1step 1.1step 1.2step 2.1∎

By step 1.2, E[fN(HN)] is exactly the probability that the statistic differs from cN by more than ε for the randomly permuted size-kN group. This tends to zero. The argument applies separately to kN=mN and kN=nN, proving the assertion.

LemmaStatement: AI-adaptedProof: AI-adaptedprecheck passOpen item page →

Mixture coupling for permuted two-sample U-statistics

Statement

Use the setup of Two-sample order-r U-statistic and first projections with fixed r, m+n=N, and m/N→p∈(0,1). Assume h is symmetric within each block and antisymmetric under exchange of its two blocks. For every assignment of independent P- and Q-distributed arguments to its 2r positions, assume Eh2<∞.

Let PˉN=(m/N)P+(n/N)Q. Couple an IID row from PˉN to the original labelled sample by drawing mixture labels and matching observations until one labelled group is exhausted, then drawing fresh observations for the remaining labels. Reorder the coupled row to agree with the original PmQn ordering except at D unmatched positions. For a uniform permutation π independent of this coupling,

m (Um,n(Zπ)−Um,n(Zˉππ0))⟶0in L2.

Neither the weak null nor assignment-wise kernel centering is needed for this coupling bound. No rate for m/N−p is assumed.

Facts & Assumptions

Given: Fixed order, proportional sample sizes, a block-symmetric and block-antisymmetric kernel, all mixed second moments finite, and the exact triangular mixture.

[F1]

Two-sample order-r U-statistic and first projections defines the kernel average Um,n over M=(mr)(nr) pairs of subsets.

[F2]

Identical distribution and IID families defines the IID convention used for each source group and the constructed mixture row.

[F4]

Cauchy-Schwarz inequality for L2 bounds expected absolute deviations and yields the triangle inequality for L2 norms.

[F5]

Independent random elements defines independence of random elements. The source groups are independent by hypothesis; independence of the auxiliary draws and permutation is imposed by the product construction [F6].

[F7]

The uniform probability space on a nonempty finite set gives uniform permutations and finite conditional averaging.

[F8]

A successive-revelation variance bound for a uniform permutation bounds a permutation function's variance by half the sum of its successive-revelation averaged squared transposition differences.

[A1]

B=max⁡η∈{P,Q}2rEηh2<∞ by the finite family of moment assumptions.

Proof

Given: The hypotheses in the Statement. Prove the L2 coupling bound.

1.1F2F5F6F7givenconstruct

Enlarge the source sample space by the finite product [F6] with N Bernoulli(pN) label factors, N fresh P-draw factors, N fresh Q-draw factors, and one uniform permutation factor [F7], where pN=m/N. Its rectangle formula makes all auxiliary factors mutually independent and independent of the source samples, in the sense of [F5]. Match each P label to the next unused original X and each Q label to the next unused original Y, using distinct fresh draws of the requested law after exhaustion. Conditional on the labels the constructed coordinates are independent with their respective label laws, so their unconditioned row Zˉ is IID PˉN. Align every shared observation with its original position, putting the fresh observations in the remaining positions by a fixed index order; call the resulting array Z~=Zˉπ0. With BN the number of P labels, its mismatch set J has size D=∣BN−m∣. Conditional on the labels and alignment, the pairs (Zi,Z~i) are independent across i; at a shared position both entries coincide, and at a mismatch the original and fresh entries are independent.

1.2F1F7givenalgebra

For any deterministic array a of length N, the permutation average of Um,n(aπ) is zero: the arguments of a fixed selected kernel term form a uniform ordered distinct 2r-tuple, and exchanging its two r-blocks preserves that distribution and negates the kernel. Therefore f(π):=Um,n(Zπ)−Um,n(Z~π) has permutation mean zero for every realized pair of arrays. This is conditional centering of the full statistic, not centering of each mixed-assignment kernel.

2.1F1F3F4step 1.1algebra

Since BN is binomial with parameters N,pN, EBN=m and ED2=NpN(1−pN)≤N/4. Thus ED≤ED2=O(N) by [F4]. Also m,n≍N and M≍N2r. Constants below depend only on r and bounds keeping pN away from 0,1, not on the labels or values.

3.1A1F1F4step 1.1step 2.1algebra

Fix two positions k,l and a permutation π, conditioning on labels and alignment. A selected term in U(aπ)−U(aπ(k,l)) vanishes unless its subset pair meets {k,l}; there are at most CN2r−1 such pairs. Each term in f(π)−f(π(k,l)) is a signed sum of four kernel evaluations, each on distinct independent population coordinates and hence with conditional second moment at most B by [A1]. The L2 triangle inequality [F4], followed by division by M, bounds the norm of this difference by CB/N. If neither π(k) nor π(l) is in J, a contribution also vanishes unless its selected positions meet π−1J outside {k,l}. There are at most CDN2r−2 such pairs, giving the sharper norm bound CBD/N2. Overlap among evaluations causes no independence requirement for this triangle inequality.

4.1F7F8step 1.2step 3.1algebra

Averaging the squared bounds in step 3.1 over uniform π, the probability that either swapped entry is a mismatch is at most 2D/N. Thus, for every fixed position transposition, Eπ,values[(f(π)−f(π(k,l)))2∣labels]≤CB(D/N3+D2/N4). First apply [F8] for fixed realized arrays, where its expectations and variance are only over π, and then average over values conditional on labels. Its sum of inner-average coefficients is N/2. By step 1.2 the conditional permutation mean is identically zero, so Eπ,values[f(π)2∣labels]≤CB(D/N2+D2/N3).

5.1step 2.1step 4.1algebra∎

Averaging step 4.1 over labels and using m≍N and step 2.1 yields E[mf(π)2]≤CB(ED/N+ED2/N2)=O(N−1/2)→0. Since Z~π=Zˉππ0 with the alignment convention in the Statement, this is the required L2 convergence.

Remarks

The 2011 report's Lemma A.1 uses assignment-wise zero means and a fixed-mixture fraction rate. The local proof instead centers the complete statistic by antisymmetry, controls its permutation variance by transpositions, and uses the exact triangular mixture. It does not assert those extra hypotheses follow from the weak null.

LemmaStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6-luna)Open item page →

First projections are measurable and square-integrable

Statement

For the setup of Two-sample order-r U-statistic and first projections, suppose h∈L2(P⊗r⊗Q⊗r). Then the section integrals defining gX and gY are finite and measurable almost everywhere, can be assigned measurable values on exceptional null sets, and define square-integrable centered projections. In particular, EP[gX(X1)]=EQ[gY(Y1)]=0, EP[gX(X1)2]≤EP⊗r⊗Q⊗r[h2],EQ[gY(Y1)2]≤EP⊗r⊗Q⊗r[h2].

Facts & Assumptions

Given: Probability measures P,Q, a product-measurable kernel h with finite second moment under P⊗r⊗Q⊗r, and the projection definitions.

[F1]

Two-sample order-r U-statistic and first projections defines gX and gY as the corresponding conditional section integrals minus θ(P,Q;h).

[F2]

Fubini's theorem for L^1 functions on a sigma-finite product gives integrable sections almost everywhere, measurable section-integral versions, and equality of iterated and product integrals.

[F3]

Cauchy-Schwarz inequality for L2 gives ∣∫ab∣≤∥a∥2∥b∥2.

[F4]

The product sigma-algebra and its finite iterates supplies the finite product sigma-algebras on which the section maps and product integrals are defined.

Proof

Given: h∈L2(P⊗r⊗Q⊗r) and the definitions in [F1]. Prove both projections are measurable, centered, and square-integrable with the displayed bounds.

1.1F1F3

The product measure is a probability measure, so applying [F3] to h and the constant function 1 gives E∣h∣≤(Eh2)1/2<∞. Thus h is integrable and θ(P,Q;h)=Eh is finite.

1.2F1F2F4

By [F4], the product decomposition uses the stated finite product sigma-algebra. Apply [F2] with first coordinate P-distributed and the remaining 2r−1 coordinates distributed according to P or Q. For P-almost every x, the section integral fX(x)=∫h(x,x2,…,xr;y1,…,yr) dPr−1dQr is finite; [F2] supplies a measurable version after setting it to zero on the exceptional null set. Reordering the coordinates and using the Y1 section gives a finite measurable fY for Q-almost every y.

2.1F2F3step 1.2

For each good x, Cauchy–Schwarz [F3] on the remaining probability product measure gives ∣fX(x)∣2≤∫h(x,x2,…,xr;y1,…,yr)2 dPr−1dQr. Integrating in P and using [F2] yields EP[fX(X1)2]≤E[h2]. The same argument for the Y1 section gives EQ[fY(Y1)2]≤E[h2].

3.1F1F2step 1.1step 2.1∎

By [F2], EPfX=Eh=θ(P,Q;h) and EQfY=Eh=θ(P,Q;h). Subtracting this finite mean gives the centered projections in [F1]; step 2.1 and finiteness of θ show both remain square-integrable.

LemmaStatement: Literature-sourcedProof: AI-adaptedprecheck passOpen item page →

Pooled-mixture projection and variance

Statement

Use the setup of Two-sample order-r U-statistic and first projections and suppose h is antisymmetric under exchange of its two blocks and has finite second moment under every independent P/Q assignment to its 2r positions. For t∈(0,1) put Pˉt=tP+(1−t)Q and define

gt(z)=∫h(z,z2,…,zr;w1,…,wr) dPˉt⊗(2r−1)−θ(Pˉt,Pˉt;h).

Then θ(Pˉt,Pˉt;h)=0, the first projection in the second block is −gt, and σt2:=EPˉt[gt2] is finite. The map t↦σt2 is continuous on (0,1). In particular, with N=m+n and pN=m/N→p∈(0,1),

r2(1+mn)σpN2⟶r21−pσp2.

Facts & Assumptions

Given: A block-symmetric kernel, its antisymmetry under exchange of blocks, and finite second moments under all mixed P/Q assignments.

[F1]

Two-sample order-r U-statistic and first projections defines the block kernel and its first projections.

[F2]

First projections are measurable and square-integrable gives the section-integral versions and their L2 bounds.

[F3]

Moments, variance, and covariance on a probability space defines the second-moment variance E[gt2] when Egt=0.

[F4]

Cauchy-Schwarz inequality for L2 makes cross products of the finitely many square-integrable mixed-assignment section functions integrable.

[F5]

Fubini's theorem for L^1 functions on a sigma-finite product supplies measurable integrable section functions and equality of iterated integrals on the mixed product probability spaces.

Proof

Given: The hypotheses in the Statement. Prove the pooled mean, sign relation, finiteness, and continuity claims.

1.1F1givenalgebra

Under Pˉt⊗2r, swapping the first and second blocks preserves the product law and changes h to −h. Therefore θ(Pˉt,Pˉt;h)=−θ(Pˉt,Pˉt;h), so it is zero.

1.2F2F3F4F5givenalgebra

Expand each factor Pˉt=tP+(1−t)Q in the section integral. This writes gt as a finite sum ∑ηcη(t)fη(z), where η ranges over assignments of P and Q to the remaining 2r−1 coordinates, cη(t) is a polynomial, and fη is the corresponding section integral. For either first-coordinate law R∈{P,Q}, write ER,η for product expectation with independent coordinates having law R first and the laws specified by η in the remaining positions. The mixed-moment assumption gives both EP,η[h2]<∞ and EQ,η[h2]<∞. Cauchy--Schwarz [F4] also makes h integrable with first-coordinate law (P+Q)/2. Fubini [F5] on that product space supplies one measurable version of fη, valid (P+Q)-almost everywhere after setting it to zero on the common exceptional set. Cauchy--Schwarz on the remaining probability product then gives ∫∣fη∣2 dR≤ER,η[h2]<∞ for each R∈{P,Q}, as in the section-integral argument of [F2]. Consequently every cross integral ∫fηfη′ dR is finite by [F4]. Expanding EPˉtgt2=tEP[gt2]+(1−t)EQ[gt2] now gives a finite sum of polynomial coefficients times these finite cross integrals, so it is finite and continuous in t.

2.1F1F2step 1.1algebra

Holding the first coordinate of h fixed and integrating the other 2r−1 coordinates with law Pˉt defines gt, since the kernel mean in [F1] is zero by step 1.1. For the projection in the second block, exchange the blocks inside the integral. Antisymmetry changes the sign and the common product law is unchanged, so the second-block projection is −gt.

3.1step 1.2algebra∎

Since m/n=pN/(1−pN) and pN→p, step 1.2 gives σpN2→σp2, while 1+m/n=1/(1−pN)→1/(1−p). Multiplication by r2 proves the displayed limit.

LemmaStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6-luna)Open item page →

Two-sample Hoeffding decomposition

Statement

In the setup of Two-sample order-r U-statistic and first projections, assume h∈L2(P⊗r⊗Q⊗r). For each X-coordinate u let Pu integrate that coordinate under P; for each Y-coordinate let Pu integrate under Q, and put Qu=I−Pu. These coordinate operators commute on h by Fubini. For 0≤a,b≤r, choose any a of the X-coordinates and any b of the Y-coordinates, apply Qu in those coordinates and Pu in all others, and call the resulting (reindexed) kernel ha,b. Separate block symmetry makes it independent of the chosen coordinate subsets. It is canonical in each of its a+b arguments: integrating it in any such argument gives zero. Let Ua,b(ha,b) be its average over all a-subsets of the X sample and all b-subsets of the Y sample, with the zero-order convention U0,0=θ(P,Q;h). Then

Um,n−θ(P,Q;h)=∑0≤a,b≤ra+b≥1(ra)(rb) Ua,b(ha,b).

In particular,

Um,n−θ(P,Q;h)=rm∑i=1mgX(Xi)+rn∑j=1ngY(Yj)+Rm,n,

where Rm,n is the sum of the terms with a+b≥2.

Facts & Assumptions

Given: Independent IID samples from P and Q, a symmetric order-(r,r) kernel with finite second moment, and the projection definitions.

[F1]

Two-sample order-r U-statistic and first projections defines Um,n, gX, and gY.

[F2]

First projections are measurable and square-integrable supplies measurable square-integrable coordinate integrals and the Fubini equalities needed for each projection operator.

[F3]

Fubini's theorem for L^1 functions on a sigma-finite product permits the coordinate-integration operators on different coordinates to commute.

[F4]

Cauchy-Schwarz inequality for L2 gives the L2 contraction bound for integration over a probability coordinate.

Proof

Given: The hypotheses in the Statement. Prove the finite decomposition and identify its order-one terms.

1.1F2F3F4algebra

For each coordinate u, Cauchy–Schwarz [F4] on its probability marginal and Fubini [F3] give ∥Puh∥2≤∥h∥2. The pointwise inequality ∣h−Puh∣2≤2∣h∣2+2∣Puh∣2 then gives ∥Quh∥2≤2∥h∥2. Thus every finite product of the coordinate operators is a well-defined L2 class. The operator Pu is idempotent because integrating an already coordinate-independent function changes nothing. For distinct coordinates u,v, Fubini [F3] gives PuPv=PvPu on each integrable section, so Pu,Qu commute across distinct coordinates.

2.1step 1.1algebra

Since I=Pu+Qu for every one of the 2r coordinates, multiplying these finitely many identities and using commutativity expands h as the sum over all subsets of coordinates of the corresponding products of Pu and Qu. This is a finite algebraic expansion; its empty selected subset is P1⋯P2rh=θ(P,Q;h).

3.1step 1.1step 2.1algebra

If coordinate u is selected, applying Pu to that component gives zero because PuQu=Pu−Pu2=0; the other commuting factors do not alter this identity. Thus each component is canonical in each selected coordinate. Block symmetry identifies the (ra)(rb) components with a selected X-coordinates and b selected Y-coordinates as reindexed copies of ha,b.

4.1F1step 2.1step 3.1algebra

Apply the expansion in step 2.1 to each summand h(XI;YJ) in Um,n, summing over sample subsets and selected kernel coordinates together. Each fixed a-subset of the X sample is contained in exactly (m−ar−a) sample r-subsets; in each, exactly one choice of a kernel coordinates selects it. Thus its coefficient after averaging is (m−ar−a)/(mr)=(ra)/(ma). Similarly, each fixed b-subset of the Y sample has coefficient (rb)/(nb). These identities also hold when a=0 or b=0. Block symmetry identifies the resulting functions with ha,b, so the combined contribution is (ra)(rb)Ua,b(ha,b). Removing the constant component gives the first displayed decomposition.

5.1F1step 4.1algebra∎

For (a,b)=(1,0), the integrated kernel is exactly gX from [F1], and its subset average is m−1∑igX(Xi); its coefficient is (r1)=r. For (a,b)=(0,1) the average is n−1∑jgY(Yj) with the same coefficient r. All remaining pairs have a+b≥2 and are, by definition, the remainder Rm,n.

LemmaStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6-luna)Open item page →

Second-moment order of the degenerate remainder

Statement

In the setup of Two-sample Hoeffding decomposition, fix r and suppose m/N→p∈(0,1), n/N→1−p, with N=m+n. At row N, let the two groups be independent IID samples from probability measures PN,QN and use the same block-symmetric kernel h. If B:=sup⁡Nmax⁡η∈{PN,QN}2rEη[h2]<∞, where Eη integrates independent arguments with the assigned marginal laws. Write EN for expectation under the row-N sample law. Then the remainder in the first-projection decomposition satisfies

EN[Rm,n2]=O(N−2),NRm,n→0in L2.

Facts & Assumptions

Given: Fixed order r, proportionally growing m,n, independent IID groups with row laws PN,QN, and the canonical Hoeffding components from the preceding lemma. Fixed P,Q is the special case PN=P, QN=Q.

[F1]

Two-sample Hoeffding decomposition: Rm,n is the finite sum of the terms (ra)(rb)Ua,b(ha,b) with a+b≥2; every selected-coordinate kernel is canonical.

[F2]

Moments, variance, and covariance on a probability space: variance is the second centered moment and covariance is the expectation of the product of centered variables.

[F3]

Cauchy-Schwarz inequality for L2: ∣E[AB]∣≤(EA2)1/2(EB2)1/2.

Proof

Given: The hypotheses in the Statement. Prove the remainder is OL2(N−1).

1.1F1F2algebra

Fix (a,b) with a+b≥1. Expand Ua,b(ha,b) as the average over its (ma)(nb) pairs of index subsets. Two summands with different X-index sets or different Y-index sets have covariance zero: choose an observation appearing in one selected-coordinate set but not the other, condition on all other observations, and integrate its canonical coordinate to zero by [F1]. The identical-index terms each have second moment at most EN[ha,b2], where EN is expectation under the row-N sample law. Therefore EN[Ua,b(ha,b)2]≤EN[ha,b2](ma)(nb), with the convention (m0)=(n0)=1.

1.2F1F3algebra

The coordinate operators in the decomposition have L2 norm at most 1 for Pu by Cauchy–Schwarz on a probability section. The pointwise inequality ∣f−Puf∣2≤2∣f∣2+2∣Puf∣2 then bounds the L2 norm of I−Pu by 2. Hence each EN[ha,b2]≤42rEN[h2]≤42rB, uniformly over N and the finitely many pairs (a,b). Since m,n≍N and a,b are fixed, (ma)(nb)≍Na+b.

2.1F1step 1.1step 1.2algebra∎

Distinct Hoeffding components in [F1] are orthogonal by the same canonical-coordinate argument used in step 1.1. Summing their variances with the fixed coefficients (ra)2(rb)2 gives EN[Rm,n2]=∑a+b≥2(ra)2(rb)2EN[Ua,b(ha,b)2]≤Cr,pBN−2, because there are finitely many pairs and every denominator has order at least N2. Multiplying by N proves EN[(NRm,n)2]→0.

DefinitionDefinition: Literature-sourcedProof: Not applicableOpen item page →

Studentized two-sample U-statistic and permutation distribution

Definition

For the two-sample statistic in Two-sample order-r U-statistic and first projections, let σ^m2 and σ^n2 be nonnegative measurable statistics of the two samples intended to estimate the variances of gX(X1) and gY(Y1). Each may depend on both groups; separate-group dependence is an additional property, not part of this definition. Define the estimated standard error and studentized statistic by

Vm,nSE:=rσ^m2m+σ^n2n,Sm,n:={Um,n/Vm,nSE,Vm,nSE>0,0,Vm,nSE=0.

For a uniform random permutation π of the pooled N=m+n observations, recompute both group statistics and both variance estimates on the permuted groups; the permutation distribution of Sm,n is the conditional law of Sm,nπ given the pooled observations.

Remarks

The report's scale is Vm,nCR=rσ^m2+(m/n)σ^n2 and its statistic is Um,n/Vm,nCR, with a m multiplier. Since Vm,nSE=Vm,nCR/m, the usual standard-error studentization is Sm,n=m Um,n/Vm,nCR. The standard-normal limits for studentized statistics use Sm,n itself as the pivot, under the null, nondegeneracy, and variance-consistency hypotheses of the corresponding limit theorem. An additional m multiplying this Sm,n is not part of the studentized pivot.

Consistency and permutation consistency are hypotheses of asymptotic validity results, not consequences of this definition. For an estimate depending on both groups, consistency on an IID pooled-mixture row alone does not imply permutation consistency; that condition must be proved for the actual estimate.

TheoremStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6-luna)Open item page →

Projection CLT for a two-sample U-statistic

Statement

Assume AC and the setup of Two-sample order-r U-statistic and first projections, r is fixed, m/N→p∈(0,1), θ(P,Q;h)=0, and 0<EP[gX2]=σX2<∞ and 0<EQ[gY2]=σY2<∞. If the fixed-order kernel has finite second moments under all mixed assignments, then

m Um,n⇒N ⁣(0, r2(σX2+p1−pσY2)).

Facts & Assumptions

Given: Independent IID samples from P and Q, fixed r, m/N→p∈(0,1), a zero parameter, and positive finite first-projection variances.

[F1]

Two-sample order-r U-statistic and first projections defines Um,n and its centered projections gX,gY.

[F2]

Second-moment order of the degenerate remainder gives NRm,n→0 in L2 under the mixed second-moment premises.

[F3]

Lindeberg-Levy iid central limit theorem: under AC, each IID centered real sample mean with positive finite variance has its normalized normal limit.

[F5]

Slutsky's theorem for real random variables transfers a weak limit after adding a term converging to zero in probability.

[A1]

AC is assumed exactly for the CLT and characteristic-function suppliers in [F3]–[F4]. The Axiom of Choice

Proof

Given: The assumptions in the Statement. Prove the centered normal limit at scale m.

1.1F1F2givenalgebra

By [F1] and [F2], m Um,n=r1m∑i=1mgX(Xi)+rmn(1n∑j=1ngY(Yj))+mRm,n, and the final term tends to zero in probability because m/N is bounded and E[(NRm,n)2]→0.

1.2F1F3A1given

The summands gX(Xi) are IID, centered, and have variance σX2; [F3] gives m−1/2∑igX(Xi)⇒N(0,σX2). The gY(Yj) are IID, centered, and have variance σY2, so the same result gives n−1/2∑jgY(Yj)⇒N(0,σY2). The two sums are independent because the two samples are independent.

2.1F1F4step 1.1step 1.2algebra

Since m/n→p/(1−p), the characteristic function of the sum of the first two terms in step 1.1 is the product of the two independent characteristic functions. By [F4] this product converges to the characteristic function of a centered normal with variance r2(σX2+(p/(1−p))σY2); the positivity assumptions make that variance nonzero.

3.1F5step 1.1step 2.1∎

The first two terms in step 1.1 therefore have the stated weak limit. The last term converges to zero in probability by step 1.1, so [F5] gives the claimed limit for m Um,n. All divisions are by positive finite variances, and no claim is made when either first-projection variance is zero.

CorollaryStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6-luna)Open item page →

Finite-sample exactness under equal distributions

Statement

For the studentized statistic in Studentized two-sample U-statistic and permutation distribution, assume P=Q, m,n≥r, and use the convention there when a permuted denominator is zero. For every such finite m,n, the inclusive permutation p-value is superuniform under this null. If an independent uniform [0,1] variable randomizes within the tied boundary rank, the resulting permutation p-value is exactly uniform.

Facts & Assumptions

Given: Independent IID samples of sizes m,n from a common probability law, the measurable coordinate-permutation action on their pooled vector, and the studentized statistic with its zero-denominator convention.

[F1]

Identical distribution and IID families: the pooled observations are IID when P=Q.

[F3]

Coordinate maps and permutations are measurable on finite products: each coordinate permutation is measurable.

[F4]

Studentized two-sample U-statistic and permutation distribution: the statistic is a measurable real function, including on zero-denominator inputs.

[F5]

Exact finite-group permutation p-values and tests: under an invariant law, every measurable statistic on a finite measurable group action has a superuniform inclusive orbit-rank p-value and an exactly uniform tie-randomized p-value.

Proof

Given: P=Q, independent IID samples with m,n≥r, all coordinate permutations of the pooled sample, and the rank conventions in [F5].

1.1F1F2givenalgebra

By [F1], the pooled vector has product law P⊗N. On every measurable rectangle A1×⋯×AN, its probability is ∏i=1NP(Ai); permuting coordinates merely reorders these factors. By uniqueness in [F2], the product law is invariant under every coordinate permutation.

2.1F3F4F5step 1.1

The symmetric group on N coordinates is finite and acts by measurable maps by [F3]. The statistic is measurable on the entire pooled sample space by [F4], including any zero-denominator permutation. Thus the hypotheses of [F5] hold for this action and invariant law.

3.1F5step 2.1∎

Apply [F5] to conclude that the inclusive p-value is superuniform and the independent-uniform tie-randomized p-value is uniform on [0,1]. The result is finite-sample and does not require a nondegenerate first projection.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Permutation CLT for a two-sample U-statistic under a weak null

Statement

Assume AC. Let X1,…,Xm be IID P and independently let Y1,…,Yn be IID Q, with N=m+n, m/N→p∈(0,1), and θ(P,Q;h)=0. Fix r and a measurable kernel h that is symmetric in each block and antisymmetric under block exchange. Assume for every assignment of independent P- and Q-distributed arguments to the 2r positions that 0<E[h2]<∞. Let σp2=EPˉ[gPˉ2]>0 for Pˉ=pP+(1−p)Q, where gPˉ is its first projection. Then the conditional permutation CDF of m Um,n satisfies

sup⁡t∈R∣R^m,nmU(t)−Φ ⁣(tτp)∣→P0,τp2=r21−pσp2.

Facts & Assumptions

Given: The sample and kernel hypotheses in the Statement, including finite positive second moments under every mixed argument assignment.

[F1]

Two-sample order-r U-statistic and first projections defines Um,n and its two first projections.

[F2]

Two-sample Hoeffding decomposition gives the first-order projection plus the finite sum of canonical higher-order terms.

[F3]

Second-moment order of the degenerate remainder makes the scaled higher-order remainder negligible.

[F4]

Pooled-mixture projection and variance identifies the pooled projection gPˉ and its variance scale τp2.

[F5]

Conditional CLT for a random subset of an IID finite population gives conditional normal convergence for a uniform sample without replacement from an IID pooled row.

[F6]

Mixture coupling for permuted two-sample U-statistics couples the original permuted U-statistic to the IID-mixture permuted U-statistic under antisymmetry and all mixed second moments.

[A1]

AC is used through [F5]. The Axiom of Choice

Proof

Given: The weak-null and mixed-moment kernel hypotheses and m/N→p∈(0,1). Prove the conditional permutation CDF limit for mUm,n.

1.1F1F2F3F4given

Let pN=m/N and PˉN=pNP+(1−pN)Q. Under an IID PˉN row, antisymmetry gives pooled parameter zero and first projections gN and −gN by [F4]. Expanding the 2r mixture factors shows that EPˉN⊗2rh2 is a finite convex combination of the mixed-assignment second moments in the hypotheses, hence is uniformly bounded in N. The same bound applies to every mixed assignment when both row laws are PˉN, so the triangular version of [F3] applies to the two groups formed by a uniform permutation and makes m times the degenerate remainder tend to zero in joint probability.

2.1F1F4F5givenalgebra

The linear projection for a uniform size-m subset A of the pooled row is mLNπ=rNnm(∑i∈AgN(Zˉi)−mN∑i=1NgN(Zˉi)). Expanding each factor of PˉN shows that gN is a finite linear combination of the section-integral functions for the finitely many mixed P/Q assignments. The assignment-wise second-moment bounds make these functions square-integrable under P and Q; their coefficients are bounded as pN→p. If H is the sum of their squared absolute values, then H is integrable under both P and Q and gN2≤CH for one constant C. Hence, for K>0, EPˉN[gN21{gN2>K}]≤C(EP[H1{H>K/C}]+EQ[H1{H>K/C}])⟶0, so the row gN(Zˉi) has uniformly integrable squares. Its variance is σN2→σp2. Let Sg,N2=(N−1)−1∑i=1N(gN(Zˉi)−gˉN)2 be its empirical finite-population variance. Apply step 1.1 of [F5] to this IID row: the uniformly integrable squares give N−1∑igN(Zˉi)2−σN2→P0, while EgˉN=0 and Var⁡(gˉN)=σN2/N→0. Therefore Sg,N2−σN2→P0 and Sg,N2→Pσp2.

3.1A1F5step 1.1step 2.1algebra

Put TNπ=∑i∈A(gN(Zˉi)−gˉN)mnSg,N2/N when Sg,N2>0, and set it to zero otherwise. Apply [F5] to this random subset. The exact identity on {Sg,N2>0} is mLNπ=cNTNπ, where cN=rN/n Sg,N. Since Sg,N2→Pσp2>0, cN→Prσp/1−p=τp and the event Sg,N2=0 has probability tending to zero. The conditional CDF of TNπ converges uniformly in probability to Φ by [F5]; thus, for each fixed t, the conditional CDF of cNTNπ converges in probability to Φ(t/τp). Indeed, on cN>0 its conditional CDF is FN,Zˉ(t/cN), whose difference from Φ(t/cN) is bounded by sup⁡u∣FN,Zˉ(u)−Φ(u)∣, and Φ(t/cN)→PΦ(t/τp). The zero-variance event is negligible. The conditional probability of the remainder exceeding any fixed threshold tends to zero in probability, since its expectation tends to zero by step 1.1 and Markov's inequality. Another threshold squeeze and continuity of the normal CDF give the same pointwise conditional limit for mUm,n(Zˉπ).

4.1F6step 3.1algebra

Let π,π′ be conditionally independent uniform permutations. By step 3.1, the two mixture-row statistics have asymptotically independent N(0,τp2) limits: conditionally on Zˉ, their joint CDF is the product of the two conditional CDFs, and each factor converges in probability to a constant. Apply the coupling [F6] to both permutations; the union bound and its L2 conclusion preserve this joint limit for the original sample.

5.1step 3.1step 4.1algebra∎

Write RN(t) for the original conditional permutation CDF. The one-permutation convergence gives ERN(t)→Φ(t/τp); the two independent permutations in step 4.1 give E[RN(t)2]→Φ(t/τp)2. Therefore Var⁡(RN(t))→0, so RN(t)→PΦ(t/τp) at every t. Since the limit CDF is continuous, choose finite tail cutoffs and a finite grid with arbitrarily small Φ increments; monotonicity bounds the supremum between grid values. This proves the stated uniform convergence.

Remarks

The proof uses the exact triangular mixture and the locally strengthened coupling estimate. It requires neither zero kernel means for other assignments nor a convergence rate for m/N−p. The weak null is needed for the accompanying true-sampling limit; antisymmetry alone centers the permutation distribution.

TheoremStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6-luna)Open item page →

Permutation CLT for a two-sample U-statistic under mixed-assignment centering

Statement

Assume AC. Let X1,…,Xm be IID P and independently let Y1,…,Yn be IID Q, with N=m+n, m/N→p∈(0,1), and θ(P,Q;h)=0. Fix r and a measurable kernel h that is symmetric in each block and antisymmetric under block exchange. Assume for every assignment of independent P- and Q-distributed arguments to the 2r positions that E[h]=0 and 0<E[h2]<∞. Let σp2=EPˉ[gPˉ2]>0 for Pˉ=pP+(1−p)Q, where gPˉ is its first projection. Then the conditional permutation CDF of m Um,n satisfies

sup⁡t∈R∣R^m,nmU(t)−Φ ⁣(tτp)∣→P0,τp2=r21−pσp2.

Facts & Assumptions

Given: The sample, all-assignment centering and moment assumptions, and positive pooled projection variance in the Statement.

[F1]

Permutation CLT for a two-sample U-statistic under a weak null gives the conditional permutation limit under the weak null, block antisymmetry, finite positive mixed second moments and positive pooled projection variance.

[A1]

AC is inherited through [F1]. The Axiom of Choice

Proof

Given: The hypotheses in the Statement. Apply the general weak-null result.

1.1F1given

The stated all-assignment assumptions include the weak null and every finite positive mixed second moment required by [F1]. Block symmetry, antisymmetry, proportional growth and positive pooled projection variance are the same hypotheses as in [F1].

2.1A1F1step 1.1algebra∎

Under AC [A1], apply [F1] to obtain the stated uniform conditional CDF limit and the identical scale τp2=r2σp2/(1−p). The additional assignment-wise centering assumptions cause no change in that scale.

Remarks

This is the all-assignment-centered specialization of the general weak-null theorem. Its extra centering assumptions are sufficient but are not needed by the locally proved permutation coupling.

TheoremStatement: AI-adaptedProof: AI-adaptedjudge pass (gpt-6-luna)Open item page →

Studentized permutation validity under a weak two-sample parameter null

Statement

Assume AC. Let independent IID samples X1,…,Xm∼P and Y1,…,Yn∼Q satisfy N=m+n, m/N→p∈(0,1), and the weak null θ(P,Q;h)=0. Fix r and a measurable kernel symmetric within each block and antisymmetric under block exchange. Assume 0<Eηh2<∞ for every independent P/Q assignment η to the 2r positions. Suppose the direct first-projection variances σX2,σY2 and limiting pooled-mixture first-projection variance σp2 are positive and finite. Put pN=m/N, PˉN=pNP+(1−pN)Q, and σN2=EPˉNgN2→σp2.

Use the nonnegative measurable variance estimates and zero-denominator convention in Studentized two-sample U-statistic and permutation distribution. Assume σ^m2→PσX2, σ^n2→PσY2 under the original samples. For general estimates depending on both groups, assume also permutation consistency: for each ε>0 and k∈{m,n},

Pπ(∣σ^k,π2−σN2∣>ε∣Z)→P0.

A sufficient condition for this last hypothesis is that each estimate is a statistic of its own group alone, and is consistent for σN2 on an IID PˉN sample of that group's size. Joint estimates may instead establish permutation consistency directly.

Then the conditional permutation CDF of Sm,n=Um,n/Vm,nSE converges uniformly in probability to Φ, and Sm,n⇒N(0,1) under the true sampling law, where

Vm,nSE=rσ^m2/m+σ^n2/n.

Facts & Assumptions

Given: The weak-null, moment, positive-projection and estimator-consistency hypotheses in the Statement.

[F1]

Permutation CLT for a two-sample U-statistic under a weak null gives the conditional permutation limit of mU with variance τp2=r2σp2/(1−p). Its proof also shows σN2→σp2 from the finite family of mixed second moments.

[F2]

Consistency transfers from IID mixture samples to permuted groups transfers consistency of a statistic of one IID triangular-mixture group to that permuted group.

[F3]

Projection CLT for a two-sample U-statistic gives the true-sampling CLT with variance a2=r2(σX2+(p/(1−p))σY2)>0.

[F4]

Studentized two-sample U-statistic and permutation distribution defines VSE and sets the ratio to zero if that standard error is zero.

[F5]

Slutsky's theorem for real random variables divides a weak limit by a positive probability limit.

[A1]

AC is inherited through [F1]–[F3]. The Axiom of Choice

Proof

Given: The hypotheses in the Statement. Prove both standard-normal limits and the sufficient consistency criterion.

1.1A1F2givenalgebra

If an estimate depends on its own group alone and has the stated IID triangular-mixture consistency, [F2] gives its unconditional convergence under a random permutation. The expectation of its conditional failure probability given Z is that unconditional failure probability, which tends to zero. Markov's inequality makes the conditional failure probability tend to zero in probability, proving the sufficient condition. For general joint estimates, this conditional convergence is assumed directly. A union bound treats the two estimates simultaneously.

2.1F4step 1.1givenalgebra

Thus mVSE,π=rσ^m,π2+(m/n)σ^n,π2 converges in conditional probability, in probability over Z, to τp=rσp2/(1−p)>0. In particular, the conditional probability of a zero denominator tends to zero.

3.1A1F1step 2.1algebra

For fixed t and 0<δ<τp/2, restrict to ∣mVSE,π−τp∣<δ. There the event Sπ≤t is bounded between the numerator events with thresholds t(τp−δ) and t(τp+δ), reversing their order for t<0. The excluded event has conditional probability tending to zero by step 2.1. Apply [F1] at those thresholds and then let δ↓0 to obtain conditional CDF convergence to Φ(t). A finite grid on a bounded interval with small normal tail probabilities and small adjacent normal CDF increments, together with monotonicity, makes this convergence uniform in t.

4.1A1F3F4F5givenalgebra∎

Under the original sampling law, estimator consistency and [F4] give mVSE→Pa>0, while [F3] gives mU⇒N(0,a2). By [F5], their ratio converges to N(0,1). The zero-denominator convention changes an event of probability tending to zero.

Remarks

The standard-error ratio S is the pivot. The report's normalization is VCR=mVSE, so its m(U/VCR) equals S. An extra factor m on the present S would be erroneous.

Assignment-wise centering is not required by the locally proved coupling and permutation CLT. The report's stronger coupling assumptions are not inferred from the weak null, and no unavailable 2013 supplementary proof is used. IID-mixture consistency alone transfers by [F2] only for separate-group statistics; general joint estimates require the explicit permutation condition above.

TheoremStatement: AI-adaptedProof: AI-adaptedprecheck passjudge pass (gpt-6-luna)Open item page →

Studentized permutation validity under mixed-assignment centering

Statement

Assume AC and the sampling and kernel conditions of Permutation CLT for a two-sample U-statistic under mixed-assignment centering. Let the direct first-projection variances σX2,σY2 and pooled variance σp2 be positive and finite, and put σN2=EPˉNgN2→σp2. Use the studentized statistic in Studentized two-sample U-statistic and permutation distribution. Suppose both variance estimates are consistent for their direct projection variances under the original samples and are permutation-consistent for σN2 in the conditional sense of Studentized permutation validity under a weak two-sample parameter null. A sufficient condition for the latter is that each estimate depends on its own group alone and is consistent for σN2 on an IID PˉN sample of that size.

Then the conditional permutation CDF of Sm,n=Um,n/Vm,nSE converges uniformly in probability to Φ, and under the true sampling law Sm,n⇒N(0,1).

Facts & Assumptions

Given: All-assignment centering and moments, positive direct and pooled projection variances, direct consistency, and either direct permutation consistency or the sufficient separate-group triangular-mixture condition.

[F1]

Studentized permutation validity under a weak two-sample parameter null proves both studentized limits and the sufficient separate-group consistency condition under the weak null and mixed moments.

[F2]

Permutation CLT for a two-sample U-statistic under mixed-assignment centering states the sampling and all-assignment centered kernel conditions inherited here.

[A1]

AC is inherited through [F1]. The Axiom of Choice

Proof

Given: The hypotheses in the Statement. Apply the general studentized theorem.

1.1F1F2given

The all-assignment centered conditions in [F2] include the weak null, block antisymmetry, fixed order, proportional sample growth and all mixed moments needed by [F1]. The remaining projection and estimator conditions here are exactly the corresponding conditions of [F1].

2.1A1F1step 1.1∎

Under AC [A1], apply [F1] to obtain both standard-normal limits and, when used, its separate-group sufficient consistency criterion.

Remarks

This specialization retains the original all-assignment centered branch. Joint two-group variance estimates satisfy the theorem when their permutation consistency is established directly; the separate-group transfer lemma does not by itself prove consistency of arbitrary joint estimates.

AlgorithmAlgorithm: Literature-sourcedProof: Not applicablejudge pass (gpt-6-luna)Open item page →

Studentized two-sample U-statistic permutation test

Definition

Inputs: Two samples of sizes m,n≥r, a measurable order-(r,r) kernel symmetric within each group, two supplied nonnegative measurable variance-estimation routines, and a significance level α∈(0,1). Use the U-statistic in Two-sample order-r U-statistic and first projections and the standard-error pivot S in Studentized two-sample U-statistic and permutation distribution, including its value zero when the estimated standard error is zero. Choose the upper alternative statistic T=S, lower alternative statistic T=−S, or two-sided statistic T=∣S∣. Optionally supply an independent uniform random variable U∈[0,1] for ties.

  1. Pool the observations into z=(X1,…,Xm,Y1,…,Yn) and compute T(z).
  2. Enumerate all N! coordinate permutations π of N=m+n. For each permutation, treat its first m observations as the first group and its last n as the second. Recompute the U-statistic and both variance estimates from those groups, apply the same zero-denominator rule, and compute T(zπ).
  3. Return the inclusive p-value pinc=N!−1#{π:T(zπ)≥T(z)} and reject when pinc≤α. If tie randomization was requested, instead return pU=N!−1{#{π:T(zπ)>T(z)}+U#{π:T(zπ)=T(z)}} and reject when pU≤α.

This is complete finite enumeration, with permutations counted with their multiplicities even when observations repeat. It does not define a Monte Carlo approximation using a subset of permutations. No generic weak-null validity is asserted by the procedure alone; that requires the applicable studentization and estimator-consistency theorem.

TheoremStatement: AI-adaptedProof: AI-adaptedjudge pass (gpt-6-luna)Open item page →

Studentized Wilcoxon validity under heterogeneous laws

Statement

Assume AC. Let independent IID real samples X1,…,Xm∼P and Y1,…,Yn∼Q satisfy m,n≥2 and m/(m+n)→p∈(0,1). Define the tie-corrected Wilcoxon kernel and weak null by h(x,y)=1{y>x}−1{y<x},θ=Eh(X,Y)=0. Equivalently, P(X<Y)+12P(X=Y)=1/2 for independent X∼P,Y∼Q. Put a(x)=EQh(x,Y), b(y)=EPh(X,y) and gp(z)=EpP+(1−p)Qh(z,V). Assume EPa2>0, EQb2>0 and EpP+(1−p)Qgp2>0.

Set U=(mn)−1∑i,jh(Xi,Yj), Ri=n−1∑jh(Xi,Yj) and Cj=m−1∑ih(Xi,Yj). Use the measurable nonnegative projection variance estimates σ^m2=1m−1∑i(Ri−U)2,σ^n2=1n−1∑j(Cj−U)2. Recompute them in every permutation. With S=U/σ^m2/m+σ^n2/n and the zero-denominator convention of Studentized two-sample U-statistic and permutation distribution, both S⇒N(0,1) and its conditional permutation CDF converges uniformly in probability to Φ. Consequently, each of the upper, lower and two-sided tests in Studentized two-sample U-statistic permutation test has rejection probability tending to α under this weak null. If P=Q, inclusive p-values are finite-sample superuniform and tie-randomized p-values are exactly uniform, without a nondegeneracy requirement for this finite-sample assertion.

The laws may be different and may have atoms. This null is a probabilistic-index null, not an assertion that equality of means or medians suffices for Wilcoxon validity.

Facts & Assumptions

Given: AC and the laws, positive projection variances, samples and estimators in the Statement.

[F1]

For a fixed antisymmetric two-sample U-kernel with positive direct and pooled first-projection variances, finite positive mixed second moments, sampling variance consistency and direct conditional permutation consistency of the estimates imply both standard-normal limits. Studentized permutation validity under a weak two-sample parameter null

[F2]

The standard-error pivot uses r=1 here and has value zero when its standard error vanishes. Studentized two-sample U-statistic and permutation distribution

[F3]

The test recomputes both variances for every permutation and uses inclusive or randomized orbit ranks, for S, −S or ∣S∣. Studentized two-sample U-statistic permutation test

[F4]

When P=Q, the pooled product law is invariant under coordinate permutations, as proved by the exactness corollary. For any measurable score on that orbit, inclusive ranks are superuniform and independent-uniform tie ranks are uniform, including its zero-denominator rule. Finite-sample exactness under equal distributions Finite-group rank validity on a measurable sample space

[F5]

For a square-integrable variable Z, P(∣Z−EZ∣>t)≤Var⁡(Z)/t2. Chebyshev's inequality for random variables

Proof

1.1F1givenalgebra

The bounded measurable kernel is antisymmetric and h(x,x)=0. All mixed second moments are finite. Positivity of EPa2 excludes a point-mass P because a is constant then; likewise positivity of EQb2 excludes a point-mass Q. A non-point-mass law has positive probability that two independent draws differ: otherwise every event would have probability zero or one, and on the real line its distribution function would identify a point mass. Thus the P/P and Q/Q second moments of h are positive. The P/Q second moment is positive because its vanishing would give a=0 P-almost surely; the Q/P case follows by antisymmetry. The first projections under the null are exactly a,b. For μN=(m/N)P+(n/N)Q, its pooled projection gN is centered by antisymmetry; its second moment σN2 converges to the positive displayed pooled variance because it is a polynomial in m/N with bounded coefficients.

2.1F5step 1.1algebra

Conditional on Xi, Ri averages n independent variables bounded by one, so E(Ri−a(Xi))2≤1/n. Consequently m−1∑i(Ri−a(Xi))2→0 in probability by the elementary nonnegative bound P(Z>t)≤EZ/t. The corresponding bound for Cj−b(Yj) is 1/m. The averages of a(Xi) and a(Xi)2 converge to their expectations by [F5], since their variances are bounded by 1/m; similarly for b. For bounded lists in [−1,1], replacing a list by one with mean squared error e changes its mean by at most e and its mean square by at most 2e, by the finite-sum Cauchy–Schwarz inequality. Since the two row/column means equal U, the stated sample variances are therefore consistent for EPa2 and EQb2.

3.1F5step 2.1algebra

Fix a deterministic pooled list z1,…,zN. Put di=N−1∑kh(zi,zk) and AN=N−1∑idi2. Antisymmetry gives ∑idi=0. For a uniform split into groups A and B of sizes m,n, condition on i∈A. Then B is a uniform n-subset of the other N−1 indices, so its average of h(zi,zj) has mean Ndi/(N−1) and variance at most 4/n. Indeed, sampling k of M bounded numbers without replacement gives variance (M−k)/(k(M−1)) times their population variance when M>1, by expanding the indicators with means k/M and pair probabilities k(k−1)/(M(M−1)); their population variance is at most one. The loose bound four covers every permitted size. Thus the expected within-A average of (Ri−di)2 is at most 4/n+(N−1)−2. For B's column averages the same argument approximates −dj with error at most 4/m+(N−1)−2. Finally, the sample means of di and di2 in either split group converge to the full-list means in conditional probability, with variance at most 4/min⁡(m,n). The finite-list bound from step 2.1 shows both recomputed variance estimates converge to AN in conditional permutation probability, uniformly over every pooled list. The factors m/(m−1) and n/(n−1) tend to one.

4.1F5step 3.1algebra

Under the original independent pooled observations, AN=N−3∑i,j,kh(Zi,Zj)h(Zi,Zk) is a bounded ordered-triple average. Terms with repeated indices comprise at most 3N2 terms. For distinct indices, expectations factor according to their fixed P/Q labels; allowing all index triples in this deterministic label average is exactly integration of h(z,v)h(z,w) under μN⊗3. Removing the repeated indices changes either expectation by at most 6/N, hence ∣EAN−σN2∣≤6/N. Two triple terms with disjoint index sets are independent. At most 9N5 ordered pairs of triples share an index, and each covariance has absolute value at most two, giving Var⁡(AN)≤18/N. By [F5], AN−σN2→0 in probability. Combining with the uniform conditional bound of step 3.1 proves the exact direct joint-estimator permutation consistency required by [F1], without asserting that joint-estimator consistency follows automatically from IID-mixture consistency.

5.1F1F2step 1.1step 2.1step 3.1step 4.1algebra

Steps 1.1, 2.1, 3.1, 4.1 verify every premise of [F1] for r=1 and the actual empirical estimates. Its two limits apply to [F2], including its vanishing-probability zero-denominator event. This proves the two asserted standard-normal limits for heterogeneous laws, with or without ties.

6.1F3F4step 5.1algebra∎

For T=S, −S or ∣S∣, the two limits in step 5.1 imply convergence of the sampling CDF and the conditional permutation CDF to the same continuous limit H, respectively the standard normal or absolute-standard-normal law. Uniform convergence for the absolute value follows by subtracting the two normal CDF tails; the limiting CDF is continuous. The maximal conditional atom is at most twice the sup-norm CDF error from H, and therefore vanishes in probability. Both inclusive and randomized ranks in [F3] consequently differ from 1−H(T) by a quantity tending to zero in probability. Since T converges to the law with continuous CDF H, 1−H(T) converges to uniform: at each u in (0,1) this follows by the continuous strictly increasing quantile of H on its support. Thus rejection probability tends to α. For P=Q, [F4] instead gives the stated exact finite-sample conclusions directly, whether or not the positive-variance asymptotic hypotheses hold.

5 · Examples, counterexamples and false statements

None yet.

Sources