Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Max-to-Sum Ratio

Statement

For realized nonnegative observations x1,…,xnx_1,\dots,x_n with sn>0s_n>0, define

sn=∑i=1nxi,mn=max⁡1≤i≤nxi,rn=mnsn.s_n=\sum_{i=1}^n x_i,\qquad m_n=\max_{1\le i\le n}x_i,\qquad r_n=\frac{m_n}{s_n}.

The realized max-to-sum ratio rnr_n is the share of the total produced by the largest observed value. Its random-sample counterpart is Rn=Mn/SnR_n=M_n/S_n, where Sn=∑i=1nXiS_n=\sum_{i=1}^n X_i and Mn=max⁡1≤i≤nXiM_n=\max_{1\le i\le n}X_i. When X1,X2,…X_1,X_2,\dots are iid nonnegative variables with finite positive mean μ=E[X]\mu=\mathbb E[X], then

Rn→0almost surely.R_n\to0\qquad\text{almost surely}.

Thus a persistent large realized value rnr_n is a warning sign that the sample sum is being governed by extremes rather than by aggregation. For Pareto Type I samples with α≤1\alpha\le1, the ordinary mean is infinite, so the finite-mean argument below does not apply.

Lowercase symbols describe the observed list; uppercase symbols describe random samples. For a nonnegative list with positive total, 1/n≤rn≤11/n\le r_n\le1. If the total is zero, the ratio is undefined. For the almost-sure limit, Sn/n→μ>0S_n/n\to\mu>0 ensures that Sn>0S_n>0 eventually almost surely; any convention on earlier zero totals gives the same limit.

Extreme-dominance intuition

In a thin-tailed or comfortable finite-mean sample, the largest observation may be memorable, but it should eventually become a negligible fraction of the whole. The total grows like nμn\mu, while the largest draw grows more slowly than nn on the scale needed to dominate the sum.

Fat-tailed samples can look different. One observation can represent a large fraction of all observed mass, especially near or below the Pareto mean boundary α=1\alpha=1. This is the same mechanism behind LLN Failure Under Infinite Mean, but RnR_n makes the dominance visible without plotting the entire running average. It complements the path view in Pre-Asymptotic LLN Behavior.

Finite-mean proof

We prove the finite-mean theorem: if the iid nonnegative observations have a positive finite mean, then the maximum becomes a negligible share of the sum. The proof uses the strong law and Borel--Cantelli. The infinite-mean Pareto regimes are cited or described heuristically, because their exact limits need heavier regular-variation machinery.

Assume Xi≥0X_i\ge0 are iid and 0<μ=E[X]<∞0<\mu=\mathbb E[X]<\infty. The strong law gives Sn/n→μS_n/n\to\mu almost surely. It remains to show that Mn/n→0M_n/n\to0 almost surely.

For any ε>0\varepsilon>0, integrability implies

∑n=1∞P(Xn>εn)<∞,\sum_{n=1}^{\infty}\mathbb P(X_n>\varepsilon n)<\infty,

because the sum is bounded by a constant multiple of ∫0∞P(X>t) dt=E[X]\int_0^\infty\mathbb P(X>t)\,dt=\mathbb E[X]. By Borel--Cantelli, Xn/n→0X_n/n\to0 almost surely. Therefore Mn/n→0M_n/n\to0 almost surely: after a random finite index all new observations are at most εn\varepsilon n, while the finitely many old observations divided by nn vanish. Hence

MnSn=Mn/nSn/n→0\frac{M_n}{S_n}=\frac{M_n/n}{S_n/n}\to0

almost surely.

For a Pareto Type I tail, Pareto Moment Existence shows that the mean exists exactly when α>1\alpha>1. The regimes are qualitatively different:

For iid nonnegative variables with P(X>0)>0P(X>0)>0, the converse to the finite-mean theorem also holds: Rn→0R_n\to0 almost surely if and only if E[X]<∞\mathbb E[X]<\infty. This converse is cited from O'Brien (1980), not proved above. It does not rule out convergence in probability with an infinite mean, as the exact Pareto α=1\alpha=1 case shows.

An exact running-list diagnostic

Consider the illustrative list 1,1,1,1,16,1,1,1,11,1,1,1,16,1,1,1,1. No distribution is fitted to these numbers; they isolate the effect of a record observation.

Prefix length nnTotal sns_nMaximum mnm_nLargest share rnr_n
4411/41/4
520164/54/5
924162/32/3

Adding the record changes the largest share abruptly. Later small observations increase the denominator while leaving the maximum fixed. Reordering the list changes this path but preserves the final ratio.

For an exact Pareto law with α>1\alpha>1, the Frechet maximum scale is xmn1/αx_m n^{1/\alpha}, whereas the sum is asymptotic to nμn\mu. Thus a typical ratio has scale n1/α−1n^{1/\alpha-1}: its exponent is −3/13-3/13 at α=1.3\alpha=1.3 and −2/3-2/3 at α=3\alpha=3. This scale comparison supplies no universal finite-nn quantile or confidence bound. The infinite-mean limits cited above require separate asymptotic results; the simulation below illustrates finite samples.

Paths and repeated samples across the mean boundary

Paired path and sampling-quantile panels show running maximum-to-sum ratios for Pareto exponents 0.8, 1 and 1.3, including large upward record jumps and slower decay at the mean boundary.

Exact iid Pareto draws with xm=1x_m=1 and α=0.8,1,1.3\alpha=0.8,1,1.3. Each path panel retains all 30,00030{,}000 prefixes for three independent paths. Each companion panel shows pointwise 10th and 90th sampling percentiles and the median from 400 independent paths, evaluated at 100 requested log-spaced checkpoints (duplicate integer checkpoints removed). These are sampling quantiles under specified laws, not confidence intervals or simultaneous path bands. Independent PCG64 streams use seeds 20261030–20261032 for the displayed paths and 20261040–20261042 for the ensembles. The broad shares at α=0.8\alpha=0.8 do not imply convergence to one; the slow decay at α=1\alpha=1 is compatible with an infinite mean and convergence of the ratio to zero in probability. Finite plots do not prove these asymptotic statements.

max_to_sum_path returns one ratio per prefix. Zero-total prefixes are undefined and return nan.

from xmath.incerto.estimators import max_to_sum_path

path = max_to_sum_path([1, 1, 1, 1, 16, 1, 1, 1, 1])
for n in (4, 5, 9):
    print(f"n={n:d}: ratio={path[n-1]:.4g}")
print("Zero prefix: [" + ", ".join(f"{v:.4g}" for v in max_to_sum_path([0, 2, 2])) + "]")

Output

n=4: ratio=0.25
n=5: ratio=0.8
n=9: ratio=0.6667
Zero prefix: [nan, 1, 0.5]

Regenerate the figure with uv run python demos/incerto/max_to_sum.py.

Reference distributions for the diagnostic

Use model-based simulation as the primary reference: specify or fit an iid parent distribution on nonnegative values, simulate independent samples of the observed size nn, and compute RnR_n for each sample. Compare the observed rnr_n with that simulated distribution. Report the parent model, fitted parameters (including α\alpha for a Pareto model), sample size, simulation count, and seed. The reference assumes the stated parent law is appropriate; a fitted-model comparison is conditional on the fit and needs a separate assessment of parameter uncertainty before it is used as a calibrated test. A tail-only fit also needs a model for the body when SnS_n includes all observations. A tail exponent alone does not determine a finite-nn reference. If the model permits zero totals, state whether the reference is conditional on Sn>0S_n>0.

As a complement, use mm-out-of-nn subsampling without replacement, with m→∞m\to\infty and m/n→0m/n\to0: compute the ratio on all size-mm subsets, or approximate that distribution using many uniformly sampled subsets. Validity requires a limiting law for the statistic with the normalization appropriate to the parent tail. For example, for a stated exact Pareto exponent, the preceding maximum and sum limits give the following local normalizing factors ckc_k for ckRkc_kR_k:

Parent regimeckc_k
α>1\alpha>1k1−1/αk^{1-1/\alpha}
α=1\alpha=1log⁡k\log k
0<α<10<\alpha<11

The Pareto note proves the boundary normalizations and cites the nondegenerate infinite-mean limit. Compare cnrnc_nr_n with the subsample values cmrmc_mr_m, using the same model parameters. Raw mm-sample quantiles are not nn-sample quantiles when cm≠cnc_m\ne c_n. For a general regularly varying parent, slowly varying factors can change the normalization. If its parameters are estimated, justify that their substitution leaves the normalized limit unchanged; mere consistency of an exponent estimate need not control a power of nn.

For iid data, a nondegenerate weak limit of the normalized statistic and m/n→0m/n\to0 justify convergence of the subsampling distribution at continuity points of the limiting CDF. Estimated centering or normalization must have negligible effect at the subsample scale. We cite this result from Politis & Romano (1994); Bickel, Götze, and van Zwet (1997), Theorem 1, p. 9 gives a proof and an explicit plug-in condition. These are asymptotic conditions; report the chosen mm and sensitivity to it.

Ordinary nn-out-of-nn empirical bootstrap validity is not assumed here. Its maximum cannot exceed the observed maximum, and with a unique sample maximum it equals that maximum with conditional probability 1−(1−1/n)n→1−e−11-(1-1/n)^n\to1-e^{-1}. This can obstruct a continuous extreme-value limit. Heavy-tailed sums can also have non-Gaussian stable limits for which ordinary bootstrap approximations fail; see Bickel, Götze, and van Zwet (1997), §2.3, Example 4, p. 5. The finite-mean theorem for RnR_n proves no resampling validity statement.

Caveats

References

References
  1. Darling, D. A. (1952). The Influence of the Maximum Term in the Addition of Independent Random Variables. Transactions of the American Mathematical Society, 73(1), 95–107. 10.1090/S0002-9947-1952-0048726-0
  2. O’Brien, G. L. (1980). A Limit Theorem for Sample Maxima and Heavy Branches in Galton–Watson Trees. Journal of Applied Probability, 17(2), 539–545. 10.2307/3213043
  3. Politis, D. N., & Romano, J. P. (1994). Large Sample Confidence Regions Based on Subsamples under Minimal Assumptions. The Annals of Statistics, 22(4), 2031–2050. 10.1214/aos/1176325770
  4. Embrechts, P., Klüppelberg, C., & Mikosch, T. (1997). Modelling Extremal Events. Springer. 10.1007/978-3-642-33483-2
  5. Resnick, S. I. (2007). Heavy-Tail Phenomena: Probabilistic and Statistical Modeling. Springer. 10.1007/978-0-387-45024-7