Randomized Quadrature with Periodic Kernels:  Applications to Cavalieri Volume Estimation

Francisco Javier Soto Sánchez

doi:10.5566/ias.3810

Authors

Francisco Javier Soto Sánchez University Rey Juan Carlos

DOI:

https://doi.org/10.5566/ias.3810

Keywords:

Cavalieri volume estimation, kernel quadrature, radial basis functions, reproducing kernel Hilbert spaces, stereological methods

Abstract

This paper studies randomized algorithms for unbiased numerical integration of d-dimensional periodic functions using kernel-based quadrature rules, with particular emphasis on rules induced by periodic radial basis function (RBF) kernels. The integration points are either deterministically generated or locally perturbed and then randomly shifted, introducing structured randomness into the scheme. The analysis builds on tools from the theory of reproducing kernel Hilbert spaces (RKHS) and Sobolev interpolation. It is shown that the resulting estimators achieve optimal variance decay rates, effectively capturing the smoothness of the integrand even when the assumed regularity is overestimated. The work is motivated by Cavalieri volume estimation, a classical problem in stereology. The theoretical results generalize this framework to higher dimensions and provide a Fourier-based perspective on smoothness, yielding a flexible and mathematically grounded alternative for randomized quadrature with periodic structure.

References

Aronszajn N (1950). Theory of reproducing kernels. Trans Am Math Soc 68:337–404.

Bach F (2017). On the equivalence between kernel quadrature rules and random feature expansions. J Mach Learn Res 18:1–38.

Baddeley A, Dorph-Petersen KA, Vedel Jensen EB (2006). A note on the stereological implications of irregular spacing of sections. J Microsc 222:177–81.

Baddeley A, Jensen EBV (2004). Stereology for statisticians. Chapman and Hall/CRC.

Briol FX, Oates C, Girolami M, Osborne MA (2015). Frank–Wolfe Bayesian quadrature: Probabilistic integration with theoretical guarantees. Adv Neural Inf Process Syst 28.

Briol FX, Oates CJ, Girolami M, Osborne MA, Sejdinovic D (2019). Probabilistic integration. Stat Sci 34:1–22.

Cobos F, Kühn T, Sickel W (2016). Optimal approximation of multivariate periodic Sobolev functions in the sup-norm. J Funct Anal 270:4196–4212.

Cruz-Orive LM (2024). Stereology: Theory and Applications. Cham: Springer.

Dick J (2007). Explicit constructions of quasi–Monte Carlo rules for the numerical integration of high-dimensional periodic functions. SIAM J Numer Anal 45:2141–76.

Dick J (2008). Walsh spaces containing smooth functions and quasi–Monte Carlo rules of arbitrary high order. SIAM J Numer Anal 46:1519–53.

Dick J, Kuo FY, Sloan IH (2013). High-dimensional integration: the quasi–Monte Carlo way. Acta Numer 22:133–288.

Dick J, Nuyens D, Pillichshammer F (2014). Lattice rules for nonperiodic smooth integrands. Numer Math 126:259–91.

Fuselier E, Hangelbroek T, Narcowich FJ, Ward JD, Wright GB (2014). Kernel based quadrature on spheres and other homogeneous spaces. Numer Math 127:57–92.

Fuselier E, Wright GB (2012). Scattered data interpolation on embedded submanifolds with restricted positive definite kernels: Sobolev error estimates. SIAM J Numer Anal 50:1753–76.

García-Fiñana M (2006). Confidence intervals in Cavalieri sampling. J Microsc 222:146–57.

García-Fiñana M, Cruz-Orive LM (2004). Improved variance prediction for systematic sampling on R. Statistics 38:243–72.

González-Villa J, Cruz M, Cruz-Orive LM (2017). On the precision of the nucleator. Image Anal Stereol 36:121–32.

Gual-Arnau X, Cruz-Orive LM (2002). Variance prediction for pseudosystematic sampling on the sphere. Adv Appl Probab 34:469–83.

Janáček J (2006). Variance of periodic measure of bounded set with random position. Comment Math Univ Carolin 47:443–55.

Janáček J (2008). Asymptotics of variance of the lattice point count. Czechoslovak Math J 58:751–8.

Janacek J, Jirak D (2019). Variance of the isotropic uniform systematic sampling. Image Anal Stereol 38:261–7.

Kanagawa M, Sriperumbudur BK, Fukumizu K (2016). Convergence guarantees for kernel-based quadrature rules in misspecified settings. Adv Neural Inf Process Syst 29.

Kanagawa M, Sriperumbudur BK, Fukumizu K (2020). Convergence analysis of deterministic kernel-based quadrature rules in misspecified settings. Found Comput Math 20:155–94.

Karvonen T, Oates CJ, Särkkä S (2018). A Bayes–Sard cubature method. Adv Neural Inf Process Syst 31.

Kiderlen M, Dorph-Petersen KA (2017). The Cavalieri estimator with unequal section spacing revisited. Image Anal Stereol 36:133–9.

Korobov NM (1963). Number-theoretic methods in approximate analysis.

Krieg D, Pozharska K, Ullrich M, Ullrich T (2023). Sampling recovery in L2 and other norms. arXiv preprint arXiv:2305.07539.

Leobacher G, Pillichshammer F (2014). Introduction to quasi–Monte Carlo integration and applications. Springer.

Lieb EH, Loss M (2001). Analysis, vol. 14. American Mathematical Society.

Matérn B (2013). Spatial variation, vol. 36. Springer.

Narcowich FJ, Schaback R, Ward JD (2002). Approximation in Sobolev spaces by kernel expansions. J Approx Theory 114:70–83.

Narcowich FJ, Ward JD, Wendland H (2006). Sobolev error estimates and a Bernstein inequality for scattered data interpolation via radial basis functions. Constr Approx 24:175–86.

Niederreiter H (1992). Random Number Generation and Quasi–Monte Carlo Methods. SIAM.

Novak E, Woźniakowski H (2012). Tractability of Multivariate Problems. European Mathematical Society.

Oates C, Girolami M (2016). Control functionals for quasi–Monte Carlo integration. In: Artificial Intelligence and Statistics. PMLR.

O’Hagan A (1991). Bayes–Hermite quadrature. J Stat Plan Inference 29:245–60.

Pausinger F, Gomez-Perez D, Gonzalez-Villa J (2019). Estimation of volume using the nucleator and lattice points. Image Anal Stereol 38:141–50.

Sloan IH (1985). Lattice methods for multiple integration. J Comput Appl Math 12:131–43.

Soto FJ (2025). Sharp unified smoothness theory for Cavalieri estimation via Fourier decay. Axioms 14:786.

Stehr M, Kiderlen M (2020a). Asymptotic variance of Newton–Cotes quadratures based on randomized sampling points. Adv Appl Probab 52:1284–307.

Stehr M, Kiderlen M (2020b). Improving the Cavalieri estimator under non-equidistant sampling and dropouts. Image Anal Stereol 39:197–212.

Stehr M, Kiderlen M (2025). Confidence intervals for Newton–Cotes quadratures based on stationary point processes. Stat Methods Appl 1–29.

Stehr M, Kiderlen M, Dorph-Petersen KA (2022). Improving Cavalieri volume estimation based on non-equidistant planar sections: The trapezoidal estimator. J Microsc 288:40–53.

Wendland H (1995). Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree. Adv Comput Math 4:389–96.

Wendland H (2004). Scattered data approximation, vol. 17. Cambridge University Press.

Zhu SX (2012). Compactly supported radial basis functions: How and why? Tech Rep OCCAM Preprint 12/57, Oxford Centre for Collaborative Applied Mathematics, University of Oxford. Eprints Archive, Oxford, UK.

Ziegel J, Baddeley A, Dorph-Petersen KA, Jensen EBV (2010). Systematic sampling with errors in sample locations. Biometrika 97:1–13.

Ziegel J, Vedel Jensen EB, Dorph-Petersen KA (2011). Variance estimation for generalized Cavalieri estimators. Biometrika 98:187–98.