A Closed-Form Formula for Consistent Lipschitz Regression on Metric Spaces with Sparse Neural Network Realizations
Abstract
Several classical machine-learning methods, such as KRRs and SVRs, are both computationally and analytically tractable since their estimators either admit closed-form expressions or are obtained by minimizing convex training objectives; neither feature is generally available for deep neural networks. We address this by introducing a simple closed-form ``two-stage'' compositional formula fΜ for reconstructing an unknown Lipschitz function f:π³β β on a metric space (π³,Ο) from N i.i.d. noisy observations. Our main result is a high-probability uniform (Lβ) recovery guarantee that jointly controls approximation and statistical errors while enjoying an optimization error of zero; in particular, we do not assume oracle access to an approximate ERM. Our secondary main results establish the optimality of our formula in three complementary senses. 1) Function space: On Ahlfors-regular metric spaces, the hypothesis class parameterized by our formula attains the optimal fat-shattering dimension. 2) Parameter space: Its dependence on the parameters is maximally numerically stable, in the sense that a smaller approximation error cannot be achieved with a smaller Lipschitz dependence on the model parameters. 3) Forward pass: Its dependence on the input is maximally regular, matching the Lipschitz constant of the target function f. When π³=[0,1]d is equipped with the ββ norm, fΜ admits algorithmic ReLU-MLP and exact ReLU-multi-head transformer realizations of depth πͺ(log(N)) with πͺ(N) nonzero parameters.
Related
- Risk Equivalence between RKHS Regression and Sequence Models for Lipschitz Spectral Algorithms
- PINN-LPCDD: Physics-Informed Neural Network with Learnable and Physically Consistent Domain Decomposition for Compressible Multi-Medium Flows
- Substochastic Path Operators for Refinement-Consistent Network Risks
- Linear regression with probabilistic networks for network-linked data
- Popularity Regression-Based Latent Space Model for Large-Scale Sparse Network Analysis
BibTeX
@misc{hong2026closedform,
title = {{A Closed-Form Formula for Consistent Lipschitz Regression on Metric Spaces with Sparse Neural Network Realizations}},
author = {Ruiyang Hong and Hrad Ghoukasian and Anastasis Kratsios},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.03129},
eprint = {2609.03129},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.03129},
}
RIS
TY - UNPB
TI - A Closed-Form Formula for Consistent Lipschitz Regression on Metric Spaces with Sparse Neural Network Realizations
AU - Ruiyang Hong
AU - Hrad Ghoukasian
AU - Anastasis Kratsios
PY - 2026
DA - 2026/09/02
PB - arXiv
DO - 10.48550/arxiv.2609.03129
UR - https://arxiv.org/abs/2609.03129
AB - Several classical machine-learning methods, such as KRRs and SVRs, are both computationally and analytically tractable since their estimators either admit closed-form expressions or are obtained by minimizing convex training objectives; neither feature is generally available for deep neural networks. We address this by introducing a simple closed-form ``two-stage'' compositional formula fΜ for reconstructing an unknown Lipschitz function f:π³β β on a metric space (π³,Ο) from N i.i.d. noisy observations. Our main result is a high-probability uniform (L^β) recovery guarantee that jointly controls approximation and statistical errors while enjoying an optimization error of zero; in particular, we do not assume oracle access to an approximate ERM. Our secondary main results establish the optimality of our formula in three complementary senses. 1) Function space: On Ahlfors-regular metric spaces, the hypothesis class parameterized by our formula attains the optimal fat-shattering dimension. 2) Parameter space: Its dependence on the parameters is maximally numerically stable, in the sense that a smaller approximation error cannot be achieved with a smaller Lipschitz dependence on the model parameters. 3) Forward pass: Its dependence on the input is maximally regular, matching the Lipschitz constant of the target function f. When π³=[0,1]^d is equipped with the β^β norm, fΜ admits algorithmic ReLU-MLP and exact ReLU-multi-head transformer realizations of depth πͺ(log(N)) with πͺ(N) nonzero parameters.
ER -