Track new results across pure and applied mathematics.
Search papers, labs, and topics across Lattice.
The authors derive the exact closed-form fourth cumulant of the limiting null distribution for the Baringhaus–Henze–Epps–Pulley (BHEP) multivariate normality test across arbitrary dimensions and smoothing parameters. This resolves a long-standing analytical bottleneck that previously required costly numerical eigenvalue decompositions or inaccurate lower-moment approximations to compute critical values. Leveraging these first four moments, the resulting Johnson and Pearson system approximations deliver essentially instantaneous, highly accurate analytic $p$-values that significantly outperform traditional lognormal baselines across dimensions and tail probabilities.
Testing high-dimensional multivariate normality no longer requires expensive numerical eigenvalue calculations—a newly derived fourth cumulant unlocks virtually instantaneous, high-precision analytic $p$-values.
The Baringhaus--Henze--Epps--Pulley (BHEP) tests form a widely applicable class of affine invariant and consistent tests for multivariate normality. Under the null hypothesis, the test statistic converges to a weighted sum of independent chi-squared random variables. Although closed-form expressions for the first three cumulants of this limiting distribution are known for arbitrary dimension $d$ and smoothing parameter $β$, an explicit expression for the fourth cumulant has so far been available only in a special univariate case. We derive the fourth cumulant for arbitrary $d\geq1$ and $β>0$ and use the resulting first four moments to construct Johnson and Pearson approximations to the limiting null distribution. These yield essentially instantaneous analytic approximations to BHEP $p$-values, without numerical eigenvalue calculations. Using the recently obtained complete spectrum as a benchmark, we show that both four-moment approximations are highly accurate over a broad range of dimensions, smoothing parameters and upper-tail probabilities, and substantially improve on two- and three-parameter lognormal approximations. Monte Carlo results indicate good finite-sample calibration for many parameter combinations, although convergence to the limiting distribution may be slow in higher dimensions for extreme values of the smoothing parameter.