← All papers
First page of Closed-form fractional radial links for elliptical Mahalanobis discriminant analysis

Closed-form fractional radial links for elliptical Mahalanobis discriminant analysis

Serhii Zabolotnii

math.ST Jul 7, 2026 · v1 stat.ME stat.ML
Lean 4 certifies the structural algebra: the elliptical log-likelihood-ratio radial-link bridge, GAM membership, and the affine-link collapse to QDA.
We study binary classification under shared-generator elliptical class-conditional distributions. The log-likelihood ratio is an additive function of the two squared Mahalanobis radii, with radial link $\varphi=\log g$; QDA is recovered only when this link is affine. We derive the Bayes radial-link family from the within-class radius law and estimate it by a finite fractional-power stochastic-polynomial projection instead of tuning a generic spline. The link is identifiable from the radius law, the plug-in estimator is $\sqrt{n}$-consistent and asymptotically normal under finite-moment regularity conditions, and the induced classifier is asymptotically Bayes-optimal in an iterated sieve limit. The structural bridge, GAM membership, and identity-link/affine-generator dichotomy are verified in Lean 4 without unproven placeholders. Against the global Mahalanobis-GAM of Ghosh et al. (2025), reimplemented with mgcv REML splines at equal input budget, the derived link is never significantly worse on three UCI benchmarks and is decisively better on breast_cancer ($[+0.009,+0.021]$ global, $[+0.109,+0.136]$ global+local). Across six real financial series under temporal-dependence-robust validation, it is never significantly worse than the fitted GAM and is significantly better on three of five heavy-tailed series plus the light-tailed control. Relative to QDA, it improves the heaviest-tailed series (oil $[+0.024,+0.070]$, S&P 500 $[+0.038,+0.126]$, JPY/USD $[+0.009,+0.047]$) and ties elsewhere. A closed-form rate simulation corroborates the $\sqrt{n}$ rate and the predicted excess-risk dichotomy between QDA's approximation-limited floor and the derived link's vanishing excess risk. The contribution is no significant loss relative to a tuned global GAM without spline smoothing-parameter selection, plus improved accuracy over QDA where generator curvature matters.

Binary classification under shared-generator elliptical class-conditional distributions, where QDA is Bayes-optimal only in the Gaussian case. The question is what function of the Mahalanobis radii is Bayes-optimal for non-Gaussian elliptical classes, and whether it should be fitted nonparametrically or derived.

The log-likelihood ratio is written as an additive function of the two squared Mahalanobis radii with radial link phi = log g. The link is derived from the within-class radius law and estimated via a finite fractional-power stochastic-polynomial (PATP) projection fed to ridge logistic regression, rather than a tuned spline GAM. Identifiability, sqrt(n)-consistency, asymptotic normality, and iterated-sieve Bayes-optimality are established. The structural algebra (LLR bridge, GAM membership, affine-link/QDA dichotomy) is verified in Lean 4 with a sorry-free development.

The derived closed-form link is never significantly worse than a REML-tuned Mahalanobis-GAM at equal input budget across UCI and financial benchmarks, and beats QDA on heavy-tailed series with the advantage tracking tail-heaviness. A closed-form rate simulation corroborates the sqrt(n) rate and the predicted excess-risk dichotomy.

Lean theoremStatement
ellLLR_eq_radial_differencethe LLR bridge (any generator, any dimension)
ellLLR_mem_span_radialLambda in span{1, phi(D0^2), phi(D1^2)} (GAM with link phi)
affine_link_mem_identity_spanaffine phi => Gaussian collapse to QDA
sq_not_mem_affine_spanquadratic radial link t->t^2 not in identity span
tGen_not_affineStudent-t generator is a concrete non-affine instance
Lean 4 theorems certifying the structural algebra
dataset (n,d,C)qdaidentitykunchenkoghosh
wine (178,13,3)0.9920.9650.9750.970
breast_cancer (569,30,2)0.9640.9010.9590.960
UCI accuracy: qda / identity / kunchenko (derived) / ghosh (fitted)