Second Maximum of a Gaussian Random Field and Exact (t-)Spacing test
This provides a novel statistical test for detecting sparse signals in various applications like Gaussian tensors and neural networks, though it is incremental in extending existing maximum-based methods.
The paper tackles the problem of detecting sparse alternatives in Gaussian random fields by introducing the concept of the second maximum and deriving an exact spacing test based on the distribution of the maximum conditioned on this second maximum. The result is a perfectly calibrated test under the null hypothesis with high power for sparse detection, supported by numerical experiments.
In this article, we introduce the novel concept of the second maximum of a Gaussian random field on a Riemannian submanifold. This second maximum serves as a powerful tool for characterizing the distribution of the maximum. By utilizing an ad-hoc Kac Rice formula, we derive the explicit form of the maximum's distribution, conditioned on the second maximum and some regressed component of the Riemannian Hessian. This approach results in an exact test, based on the evaluation of spacing between these maxima, which we refer to as the spacing test. We investigate the applicability of this test in detecting sparse alternatives within Gaussian symmetric tensors, continuous sparse deconvolution, and two-layered neural networks with smooth rectifiers. Our theoretical results are supported by numerical experiments, which illustrate the calibration and power of the proposed tests. More generally, this test can be applied to any Gaussian random field on a Riemannian manifold, and we provide a general framework for the application of the spacing test in continuous sparse kernel regression. Furthermore, when the variance-covariance function of the Gaussian random field is known up to a scaling factor, we derive an exact Studentized version of our test, coined the $t$-spacing test. This test is perfectly calibrated under the null hypothesis and has high power for detecting sparse alternatives.