Sensitivity Analysis on the Sphere and a Spherical ANOVA Decomposition
Pith reviewed 2026-05-16 19:17 UTC · model grok-4.3
The pith
A function on the sphere decomposes into a sum of terms each depending on a coordinate subset and a parity vector.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We establish sensitivity analysis on the sphere by presenting formulas that decompose a function f from the d-sphere to the reals into a sum of terms f_{u,ξ} indexed by subsets u of the d+1 coordinates and a parity vector ξ. Each such term depends only on the variables in u, and the decomposition is enabled by orthogonal basis functions suited to the sphere.
What carries the argument
The spherical ANOVA decomposition with parity, which splits functions using subsets of coordinates and even/odd parity to capture geometry-induced dependencies.
Load-bearing premise
The natural geometry on the sphere leads to specific dependencies between the input variables for which suitable orthogonal basis functions permit the decomposition and approximation.
What would settle it
A calculation showing that the proposed decomposition formula does not sum back to the original function for some test case on the sphere would disprove the claim.
Figures
read the original abstract
We establish sensitivity analysis on the sphere. We present formulas that allow us to decompose a function $f\colon \mathbb S^d\rightarrow \mathbb R$ into a sum of terms $f_{\boldsymbol u,\boldsymbol \xi}$. The index $\boldsymbol u$ is a subset of $\{1,2,\ldots,d+1\}$, where each term $f_{\boldsymbol u,\boldsymbol \xi}$ depends only on the variables with indices in $\boldsymbol u$. In contrast to the classical analysis of variance (ANOVA) decomposition, we additionally use the decomposition of a function into functions with different parity, which adds the additional parameter $\boldsymbol \xi$. The natural geometry on the sphere naturally leads to the dependencies between the input variables. Using certain orthogonal basis functions for the function approximation, we are able to model high-dimensional functions with low-dimensional variable interactions.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript claims to establish sensitivity analysis on the sphere by providing formulas to decompose a function f : S^d → R into a sum of terms f_{u,ξ}, where u is a subset of {1,2,…,d+1} and ξ encodes parity. This extends classical ANOVA by incorporating parity decomposition, using orthogonal basis functions to model high-dimensional functions with low-dimensional variable interactions induced by the sphere's geometry.
Significance. If the central decomposition holds with complete orthogonal projections, the result would offer a novel framework for sensitivity analysis in spherical settings, potentially useful for high-dimensional approximation where variable dependencies arise naturally from the quadratic constraint. The use of parity adds an interesting dimension to the decomposition.
major comments (2)
- [Abstract and §2] The completeness of the decomposition, i.e., that the sum of all f_{u,ξ} recovers f exactly under the spherical measure, is asserted but requires an explicit proof that the chosen orthogonal basis induces projections that are complete in L²(S^d) and that no cross terms arise from the coordinate coupling; without this, the claim that each term depends only on variables in u may not hold uniquely.
- [§4] The specific orthogonal basis functions are referred to but not constructed in detail; an example or definition (such as restrictions of multivariate polynomials or spherical harmonics symmetrized by parity) is needed to verify they permit the subset decomposition while respecting ∑x_i²=1.
minor comments (2)
- Clarify the notation for the parity vector ξ early in the paper.
- [Introduction] Add a reference to classical ANOVA decomposition for comparison in the introduction.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments on our manuscript. We address each major point below and have revised the manuscript to incorporate the requested clarifications and proofs.
read point-by-point responses
-
Referee: [Abstract and §2] The completeness of the decomposition, i.e., that the sum of all f_{u,ξ} recovers f exactly under the spherical measure, is asserted but requires an explicit proof that the chosen orthogonal basis induces projections that are complete in L²(S^d) and that no cross terms arise from the coordinate coupling; without this, the claim that each term depends only on variables in u may not hold uniquely.
Authors: We agree that an explicit proof of completeness strengthens the manuscript. In the revised version, Section 2 now contains a detailed proof that the parity-augmented orthogonal projections are complete in L²(S^d) with respect to the spherical measure. The argument proceeds by verifying that the chosen basis (restrictions of multivariate polynomials orthogonalized on the sphere) spans the full space and that the parity decomposition eliminates cross terms arising from the quadratic constraint ∑x_i²=1, ensuring each f_{u,ξ} depends only on the coordinates in u. revision: yes
-
Referee: [§4] The specific orthogonal basis functions are referred to but not constructed in detail; an example or definition (such as restrictions of multivariate polynomials or spherical harmonics symmetrized by parity) is needed to verify they permit the subset decomposition while respecting ∑x_i²=1.
Authors: We have expanded Section 4 with an explicit construction of the basis functions as restrictions to the sphere of even- and odd-degree multivariate polynomials, orthogonalized via Gram-Schmidt with respect to the surface measure. We also provide a concrete low-dimensional example (d=2) using parity-symmetrized spherical harmonics to demonstrate how the subset decomposition is realized while automatically satisfying the sphere constraint. revision: yes
Circularity Check
No circularity: decomposition grounded in spherical geometry and orthogonal bases
full rationale
The provided abstract and description present formulas for an ANOVA-style decomposition on the sphere that incorporates parity via an additional index ξ, relying on the natural geometry of S^d and unspecified but standard orthogonal basis functions for approximation. No equations or steps are exhibited that reduce the claimed sum of f_{u,ξ} terms to a fitted parameter, self-definition, or self-citation chain; the completeness claim is asserted via the choice of bases without reducing to the input data by construction. This is the expected non-finding for a paper whose central construction is an explicit orthogonal expansion under the spherical measure.
Axiom & Free-Parameter Ledger
axioms (2)
- standard math Orthogonal basis functions exist on the sphere that allow the stated decomposition
- domain assumption Sphere geometry induces natural dependencies among coordinates
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/AlexanderDuality.lean, IndisputableMonolith/Cost/FunctionalEquation.leanalexander_duality_circle_linking, washburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We propose a spherical ANOVA decomposition ... using the decomposition of a function into functions with different parity, which adds the additional parameter ξ ... orthogonal basis functions ... Gegenbauer polynomials ... Jacobi polynomials
-
IndisputableMonolith/Foundation/ArithmeticFromLogic.leanembed_injective, embed_strictMono_of_one_lt unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The terms A_u f and A_v f ... are orthogonal with respect to the surface measure on S^d
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
M. Abramowitz and I. A. Stegun.Handbook of Mathematical Functions, With Formulas, Graphs, and Mathe- matical Tables. National Bureau of Standards, 1964
work page 1964
-
[2]
H.-J. Bungartz and M. Griebel. Sparse grids.Acta Numer., 13:147–269, 2004.doi:10.1017/ S0962492904000182
work page 2004
-
[3]
R. Caflisch, W. Morokoff, and A. Owen. Valuation of mortgage-backed securities using Brownian bridges to reduce effective dimension.J. Comput. Finance, 1(1):27–46, 1997.doi:10.21314/jcf.1997.005
-
[4]
M. Chu, N. Del Buono, L. Lopez, and T. Politi. On the low-rank approximation of data on the unit sphere.SIAM Journal on Matrix Analysis and Applications, 27(1):46–60, Jan. 2005.doi:10.1137/s0895479803433295
-
[5]
S. da Veiga. Kernel-based ANOV A decomposition and shapley effects – application to global sensitivity analysis, 2021.arXiv:2101.05487,doi:10.48550/ARXIV.2101.05487. 24 Sensitivity Analysis on the Sphere and a Spherical ANOV A DecompositionA PREPRINT
- [6]
-
[7]
C. F. Dunkl and Y . Xu.Orthogonal Polynomials of Several Variables. Encyclopedia of Mathematics and its Applications. Cambridge University Press, 2nd edition, 2014.doi:10.1017/CBO9781107786134
-
[8]
The Annals of Statistics , author =
B. Efron and C. Stein. The Jackknife Estimate of Variance.Ann. Stat., 9(3):586–596, 1981.doi:10.1214/ aos/1176345462
-
[9]
F. Filbir and W. Themistoclakis. Polynomial approximation on the sphere using scattered data.Math. Nachr., 281:650–668, 2008.doi:10.1002/mana.200710633
-
[10]
M. Griebel. Sparse grids and related approximation schemes for higher dimensional problems. InFOCM 2005, pages 106–161. Cambridge University Press, 2006
work page 2005
-
[11]
M. Griebel and M. Holtz. Dimension-wise integration of high-dimensional functions with applications to finance. J. Complex., 26(5):455–489, 2010.doi:10.1016/j.jco.2010.06.001
-
[12]
M. Hegland. Adaptive sparse grids. InCTAC2001, volume 44, pages C335–C353, 2003.doi:10.21914/ anziamj.v44i0.685
work page 2003
-
[13]
30 Leland McInnes, John Healy, and Steve Astels
W. Hoeffding. A class of statistics with asymptotically normal distribution.Ann. Math. Statist., 19(3):293–325, 1948.doi:10.1214/aoms/1177730196
-
[14]
M. Holtz.Sparse grid quadrature in high dimensions with applications in finance and insurance, volume 77 ofLecture Notes in Computational Science and Engineering. Springer-Verlag, Berlin, 2011.doi:10.1007/ 978-3-642-16004-2
work page 2011
- [15]
-
[16]
T. Jiang and A. B. Owen. Quasi-regression with shrinkage.Math. Comput. Simul., 62(3):231–241, 2003. 3rd IMACS Seminar on Monte Carlo Methods.doi:10.1016/S0378-4754(02)00253-7
-
[17]
S. Kucherenko, S. Tarantola, and P. Annoni. Estimation of global sensitivity indices for models with dependent variables.Comput. Phys. Commun., 183(4):937–946, 2012.doi:10.1016/j.cpc.2011.12.020
-
[18]
F. Y . Kuo, I. H. Sloan, G. W. Wasilkowski, and H. Wo´zniakowski. On decompositions of multivariate functions. Math. Comp., 79(270):953–966, 2009.doi:10.1090/s0025-5718-09-02319-9
-
[19]
C. Labreuche.Contribution of Subsets of Variables in Global Sensitivity Analysis with Dependent Variables, pages 233–248. Springer Nature Switzerland, Nov. 2024.doi:10.1007/978-3-031-76235-2_18
-
[20]
C. Lemieux and A. B. Owen. Quasi-regression and the relative importance of the ANOV A components of a function. InMCQMC 2000, pages 331–344. Springer Berlin Heidelberg, 2002.doi:10.1007/ 978-3-642-56046-0_22
work page 2000
-
[21]
G. Li and H. Rabitz. General formulation of HDMR component functions with independent and correlated variables.J. Math. Chem., 50:99–130, 2012.doi:10.1007/s10910-011-9898-0
-
[22]
L. Lippert, D. Potts, and T. Ullrich. Fast hyperbolic wavelet regression meets ANOV A.Numer. Math., 154:155– 207, 2023.doi:10.1007/s00211-023-01358-8
-
[23]
R. Liu and A. B. Owen. Estimating mean dimensionality of analysis of variance decompositions.J. Amer. Statist. Assoc., 101(474):712–721, 2006.doi:10.1198/016214505000001410
-
[24]
E. Novak and H. Wo´zniakowski.Tractability of Multivariate Problems Volume I: Linear Information. Eur. Math. Society, EMS Tracts in Mathematics V ol 6, 2008
work page 2008
-
[25]
C. C. Paige and M. A. Saunders. LSQR: An algorithm for sparse linear equations and sparse least squares.ACM Trans. Math. Software, 8:43–71, 1982.doi:10.1145/355984.355989
-
[26]
D. Potts and M. Schmischke. Approximation of high-dimensional periodic functions with Fourier-based meth- ods.SIAM J. Numer. Anal., 59(5):2393–2429, 2021.doi:10.1137/20M1354921
-
[27]
D. Potts and M. Schmischke. Learning multivariate functions with low-dimensional structures using polynomial bases.J. Comput. Appl. Math., 403:113821, 2022.doi:10.1016/j.cam.2021.113821
-
[28]
D. Potts and L. Weidensager. ANOV A-boosting for random Fourier features.Appl. Comput. Harmon. Anal., 79:101789, 2025.doi:10.1016/j.acha.2025.101789
-
[29]
H. Rabitz and ¨O. F. Alis ¸. General foundations of high-dimensional model representations.J. Math. Chem., 25:197–233, 1999.doi:10.1023/A:1019188517934. 25 Sensitivity Analysis on the Sphere and a Spherical ANOV A DecompositionA PREPRINT
-
[30]
S. Rahman. A generalized ANOV A dimensional decomposition for dependent probability measures.SIAM-ASA J. Uncertain. Quantif., 2(1):670–697, 2014.doi:10.1137/120904378
-
[31]
A. Saltelli and I. Sobol. Sensitivity analysis for nonlinear mathematical models: Numerical experience.Mat. Model., 7, 1995
work page 1995
-
[32]
I. M. Sobol. On sensitivity estimation for nonlinear mathematical models.Keldysh Applied Mathematics Institute, 1:112–118, 1990
work page 1990
-
[33]
I. M. Sobol. Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates.Math. Comput. Simul., 55(1-3):271–280, 2001.doi:10.1016/s0378-4754(00)00270-6
-
[34]
Wahba.Spline Models for Observational Data
G. Wahba.Spline Models for Observational Data. Springer, SIAM, 1990.doi:10.1137/1.9781611970128
-
[35]
L. Weidensager.Efficient High-Dimensional Approximation: ANOVA Decomposition Meets Wavelets and Random Fourier Features. PhD thesis, Universit ¨atsverlag Technische Universit ¨at Chemnitz, 2025.doi: 10.51382/978-3-96100-272-6
-
[36]
Y . Xie, R. Shi, H. Schaeffer, and R. Ward. SHRIMP: sparser random feature models via iterative magnitude pruning. InMSML, volume 190, pages 303–318. PMLR, 2022. URL:https://proceedings.mlr.press/ v190/xie22a.html. A Calculation of the spherical ANOV A decompositions in Section 6 In this section we calculate the spherical ANOV A decomposition from Definit...
work page 2022
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.