REVIEW 3 major objections 4 minor 24 references
Gaussian Integral based Bayesian Smoother
T0 review · 3 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read This paper proposes a Gaussian-integral RTS smoother that evaluates polynomial moments exactly and reports lower RMSE than three standard smoothers on a Van der Pol oscillator.
desk verdict A sensible smoothing extension with a misprinted core formula: as printed, the exact-integration theorem gives the wrong moments, though the fix looks trivial. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The machinery is the Gaussian-integral moment identity for multivariate polynomials (Theorem 1, Eq. (12)), together with its extension to $x f(x)^\top$-type integrals (Theorem 2, Eq. (15)). The identity diagonalizes the covariance $P = S\,\mathrm{diag}(c_1,\ldots,c_n)\,S^{-1}$, expands each monomial about $\hat{x}$ with multinomial coefficients, and uses the fact that a one-dimensional Gaussian integral of $x^m$ vanishes for odd $m$ and equals a Gamma-function expression for even $m$. This converts every integral in the forward filter and in the backward smoother cross-covariance (8) into a finite combinatorial sum, so the smoother's only Gaussian assumption is the moment-matched posterior; no quadrature points, $\sigma$ points, or Taylor linearization are introduced. The RTS update equations (10)-(11) themselves remain unchanged.
What would settle it
For a scalar Gaussian $x \sim N(\mu,\sigma^2)$, evaluate Eq. (12) for $f(x)=x^2$ and $f(x)=x^3+2x$ and compare with the closed forms $\mu^2+\sigma^2$ and $\mu^3+3\mu\sigma^2$; any mismatch invalidates the formula, and the same check can be run on the multivariate example used in the Van der Pol simulation.
Extended reading notes
Core claim
The central claim is that every integral the RTS smoother needs, both in the forward filtering pass and in the backward cross-covariance calculation, can be written as a finite sum when the process and measurement functions are polynomials. For a polynomial $f(x)$ and a Gaussian pdf $N(x;\hat{x}, P)$, the integral $\int f(x)N(x;\hat{x},P)\,dx$ is evaluated exactly by diagonalizing $P$ into an orthogonal matrix $S$ and eigenvalues $c_i$, expanding the shifted monomials multinomially, and applying the one-dimensional odd-even moment identity, so each term is a product of Gamma functions, powers of $\hat{x}$, and entries of $S$. Theorem 2 extends the same identity to cross-integrals $\int x_g f_{g'}(x)N(x;\hat{x},P)\,dx$, which supply the cross-covariance $P_{x_k,x_{k+1}}$ in the backward step. On the Van der Pol oscillator, whose dynamics contain $x_1^2x_2$ and $x_3(1-x_1^2)x_2$ terms, the resulting GIRTSS reports lower RMSE than ERTSS, URTSS, and CRTSS at essentially every backward time step, with the same estimate as the Gaussian-integral filter at the final time. The paper's own caveat is that for non-polynomial nonlinearities the exactness is lost once the function is replaced by a truncated Taylor expansion.
Load-bearing premise
The load-bearing premise is that the combinatorial Gaussian-integral formula of Theorem 1, whose proof is cited to earlier work rather than given here and whose transcription into Theorem 2 contains notation that is hard to verify as printed, correctly computes polynomial moments under a Gaussian density; the smoother's exactness claim collapses if that identity is wrong.
Editorial extensions
If this is right
- For any state-space model with polynomial process and measurement functions, GIRTSS evaluates the Gaussian moments appearing in the RTS recursion exactly, so the remaining error comes only from the Gaussian approximation of the posterior, not from numerical integration.
- At the final time step the GIRTSS estimate coincides with the Gaussian-integral filter's, and moving backward the smoother error decreases, so the backward pass strictly adds information rather than re-approximating it.
- In the Van der Pol simulation, GIRTSS achieves lower average RMSE than ERTSS, URTSS, and CRTSS on all three states; the margin over ERTSS is the largest, and the method's relative execution time is about 6.3 times that of the EKF-based smoother.
- For non-polynomial nonlinearities the paper's method still applies only after a Taylor expansion, and then the exactness property no longer holds; the accuracy advantage will depend on how well the truncated polynomial represents the function.
Reading between the lines
- If the combinatorial identity is verified and numerically stable, the same exact-moment substitution could be dropped into other Gaussian-assumed estimators, including fixed-lag smoothers and Gaussian filters, wherever the integrand is a polynomial times a Gaussian density.
- A natural stress test beyond the paper's example is a higher-degree polynomial or stiff Van der Pol dynamics with a smaller sampling interval; the index-combination sum in Eq. (15) grows quickly with degree and state dimension, so the computational cost may erode the accuracy gain in that regime.
- The method's reliance on an eigendecomposition of the covariance every step suggests near-degenerate covariances, where eigenvalues $c_i$ are very small or close together, could be the practical limit on the formula's numerical reliability.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes a Gaussian-integral RTS smoother (GIRTSS) for discrete-time nonlinear state-space models with Gaussian noise. The key idea is to evaluate the Gaussian integrals appearing in the forward-filtering and backward-smoothing equations exactly when the nonlinear functions are polynomials, using a multivariate Gaussian integral formula from the authors' earlier work [19]. The resulting smoother is then inserted into the standard RTS smoothing recursions. The algorithm is tested on a Van der Pol oscillator model and compared with ERTSS, URTSS, and CRTSS in terms of RMSE and relative execution time.
Significance. If the central formulas were correct, the contribution would be a clean and potentially useful specialization: for polynomial nonlinearities, the filtering and smoothing integrals are computed exactly rather than by linearization or cubature, which should improve accuracy when the Gaussian assumption is adequate. The manuscript is honest about the restriction to polynomial nonlinearities and notes that non-polynomial functions require Taylor approximation. However, the paper does not prove its central Theorem 1 (only cites [19]), and the theorem as printed is false; moreover, the pseudocode contains a sign error in the Kalman update. The numerical validation is a single Monte Carlo study without statistical error bars or code release. With a corrected theorem, a self-contained derivation, and a corrected algorithm, the idea is defensible but incremental.
major comments (3)
- [Section III, Theorem 1, Eq. (12)] The Gamma-function argument in Eq. (12) is incorrect. For the univariate case n=1, xhat=0, c=1, f1(x)=x^2, the formula as printed returns (1/sqrt(pi))*2*Gamma(2.5)=1.5, whereas E[x^2]=1; for x^4 it returns 26.25 instead of 3. The standard zero-mean Gaussian moment requires Gamma((m+1)/2), so the printed argument should be (Σ_j a_{j,i+1}+1)/2, not Σ_j a_{j,i+1}+1/2. Because Theorem 2 (Eq. (15)) and Algorithm 1 (lines 4-8) inherit this expression, the central exact-integration claim is not supported by the manuscript as written.
- [Section IV, Algorithm 1, line 10] Algorithm 1 line 10 writes xhat_{k|k} = xhat_{k|k-1} - K_k(y_k - yhat_{k|k-1}); the Kalman update requires a plus sign before the Kalman gain term. As printed, the filtering recursion is incorrect. If the simulations used the correct plus sign, the pseudocode must be corrected; if the simulations used the printed minus sign, the reported RMSE values are not produced by the algorithm described.
- [Section III, Theorem 1 and Theorem 2] The proof of Theorem 1 is only cited to [19], and Theorem 2 follows directly from Theorem 1. Because Eq. (12) is the load-bearing exact-integration identity, the manuscript should give a self-contained derivation or at least a complete, correct statement with all indices carefully defined. The current notation is ambiguous: the role of the matrix S in the factor ∏_{j=1}^n S_{i,j}^{a_{j,i+1}} and the summation convention in the Gamma argument are not clear enough for a reader to implement Algorithm 1 without the external source.
minor comments (4)
- [Section II, Eq. (6)] The joint density in Eq. (6) is written as p(x_k, x_{k+1}|y_{1:k}) = p(x_{k+1}|x_{1:k}) p(x_k|y_{1:k}); the first factor should be p(x_{k+1}|x_k), and the conditioning on y_{1:k} should be made explicit.
- [Throughout] There are several typographical issues: Eq. (11) contains 'Pxk.xk+1' with a dot instead of a proper subscript, Note 1 writes 'I = [I1 I2 ... I3]' where the last entry should be I_n, and Table I reports the URTSS RET as 1.962 in one row and 1.96 in the neighboring block.
- [Section IV, Tables I and Figures 1-2] The numerical comparison is based on a single Van der Pol scenario with 1000 Monte Carlo runs but no standard errors, confidence intervals, or statistical significance tests; the average RMSE differences are small for some states, so the superiority claim would be more convincing with additional scenarios and error bars.
- [Reproducibility] No code or pseudocode-verification details are provided; given the discrepancy in Algorithm 1 line 10, the authors should state explicitly whether the simulations were run with the plus or minus sign and ideally release the implementation.
Circularity Check
No significant circularity: the GI-based smoother applies a parameter-free moment formula (cited from prior work) to standard RTS smoothing equations; no fitted input is relabeled as a prediction.
full rationale
The derivation chain is: (i) assume Gaussian filtering densities; (ii) evaluate the integrals in the forward filter and in the cross-covariance (8) using the multivariate Gaussian moment formula of Theorem 1; (iii) feed those moments into the standard RTS update equations (10)-(11). Theorem 1 is quoted from [19] and Theorem 2 is a direct row-wise application of Theorem 1 to x_g f_{g'}(x). Neither theorem is defined in terms of the smoother output, and no parameter is fitted to the Van der Pol data and then renamed a prediction. The exactness claim is a mathematical property of polynomial moments under a Gaussian, and the simulations are independent comparisons against ERTSS, URTSS, and CRTSS. The only self-citation is [19], whose overlapping authors (Kumar and Bhaumik) supply the cited moment formula; because that cited result is a parameter-free, externally checkable mathematical identity whose assumptions do not include the proposed smoother, this citation does not make the argument circular. Concerns in the manuscript that are correctness issues rather than circularity: Eq. (12)'s Gamma-function argument appears misprinted (a moment of order r should scale as (r+1)/2 rather than r+1/2 as printed), Algorithm 1 line 10 has a minus sign in the Kalman mean update, and the proof of Theorem 1 is deferred to [19] rather than proved in the paper. These omissions and typos affect reliability, but they do not amount to the derivation being equivalent to its inputs by construction.
Assumptions & free parameters
assumptions (3)
- standard math Theorem 1 from [19] provides the exact Gaussian integral of polynomial functions over a Gaussian pdf.
- domain assumption State distributions remain Gaussian under the Bayesian recursion.
- domain assumption Process and measurement noises are Gaussian, uncorrelated, and with known covariances Q and R.
Cite this review
Pith. "Pith review of Gaussian Integral based Bayesian Smoother." pith.science (2026). https://pith.science/paper/FOVTN6FS
@misc{pith2026250106881,
author = {Pith},
title = {Pith review of: Gaussian Integral based Bayesian Smoother},
year = {2026},
howpublished = {\url{https://pith.science/paper/FOVTN6FS}},
note = {Machine review of arXiv:2501.06881}
}
read the original abstract
This work introduces the Gaussian integration to address a smoothing problem of a nonlinear stochastic state space model. The probability densities of states at each time instant are assumed to be Gaussian, and their means and covariances are evaluated by utilizing the odd-even properties of Gaussian integral, which are further utilized to realize Rauch-Tung-Striebel (RTS) smoothing expressions. Given that the Gaussian integration provides an exact solution for the integral of a polynomial function over a Gaussian probability density function, it is anticipated to provide more accurate results than other existing Gaussian approximation-based smoothers such as extended Kalman, cubature Kalman, and unscented Kalman smoothers, especially when polynomial types of nonlinearity are present in the state space models. The developed smoothing algorithm is applied to the Van der Pol oscillator, where the nonlinearity associated with their dynamics is represented using polynomial functions. Simulation results are provided to demonstrate the superiority of the proposed algorithm.
Reference graph
Works this paper leans on
-
[19]
A new method for nonlin ear state estimation problem,
K. Kumar, S. Das, and S. Bhaumik, “A new method for nonlin ear state estimation problem,” Digital Signal Processing , vol. 132, p. 103788, 2023
work page 2023
-
[1]
Multiple object trajectory estima tion using backward simulation,
Y . Xia, L. Svensson, ´A. F. Garc´ ıa Fern´ andez, J. L. Williams, D. Svens- son, and K. Granstr¨ om, “Multiple object trajectory estima tion using backward simulation,” IEEE Transactions on Signal Processing , vol. 70, pp. 3249–3263, 2022
work page 2022
-
[2]
GVINS: Tightly coupled GNSS–v isual– inertial fusion for smooth and consistent state estimation ,
S. Cao, X. Lu, and S. Shen, “GVINS: Tightly coupled GNSS–v isual– inertial fusion for smooth and consistent state estimation ,” IEEE Trans- actions on Robotics , vol. 38, no. 4, pp. 2004–2021, 2022
work page 2004
-
[3]
R ao- Blackwellized posterior linearization backward SLAM,
´A. F. Garc´ ıa Fern´ andez, R. Hostettler, and S. S¨ arkk¨ a, “R ao- Blackwellized posterior linearization backward SLAM,” IEEE Trans- actions on V ehicular Technology, vol. 68, no. 5, pp. 4734–4747, 2019
work page 2019
-
[4]
C. Ritter, T. Wollmann, J. Y . Lee, A. Imle, B. M¨ uller, O. T . Fackler, R. Bartenschlager, and K. Rohr, “Data fusion and smoothing f or proba- bilistic tracking of viral structures in fluorescence micro scopy images,” Medical Image Analysis , vol. 73, p. 102168, 2021
work page 2021
-
[5]
Efficient particle smoothing for Bayesian inference in dynamic survival models,
P . Munezero, “Efficient particle smoothing for Bayesian inference in dynamic survival models,” Computational Statistics , vol. 37, no. 2, pp. 975–994, 2022
work page 2022
-
[6]
J. Zhang, Q. Zhang, Y . Li, and Q. Wang, “Sequential Bayesi an inference for agent-based models with application to the Chinese busi ness cycle,” Economic Modelling , vol. 126, p. 106381, 2023
work page 2023
-
[7]
S. S¨ arkk¨ a and L. Svensson, Bayesian filtering and smoothing . Cam- bridge University Press, 2023, vol. 17
work page 2023
Show all 24 references
-
[8]
Bar Shalom, X
Y . Bar Shalom, X. R. Li, and T. Kirubarajan, Estimation with applica- tions to tracking and navigation: theory algorithms and sof tware. John Wiley & Sons, 2004
2004
-
[9]
Maximum likelih ood estimates of linear dynamic systems,
H. E. Rauch, F. Tung, and C. T. Striebel, “Maximum likelih ood estimates of linear dynamic systems,” AIAA journal, vol. 3, no. 8, pp. 1445–1450, 1965
1965
-
[10]
Gaussian filters for nonlinear filte ring problems,
K. Ito and K. Xiong, “Gaussian filters for nonlinear filte ring problems,” IEEE Transactions on Automatic Control , vol. 45, no. 5, pp. 910–927, 2000
2000
-
[11]
The unscented Kalman filte r for nonlinear estimation,
E. A. Wan and R. V an Der Merwe, “The unscented Kalman filte r for nonlinear estimation,” in Proceedings of the IEEE 2000 Adaptive Sys- tems for Signal Processing, Communications, and Control Sy mposium. IEEE, 2000, pp. 153–158
2000
-
[12]
Unscented Rauch–Tung–Striebel smoothe r,
S. S¨ arkk¨ a, “Unscented Rauch–Tung–Striebel smoothe r,” IEEE transac- tions on Automatic Control , vol. 53, no. 3, pp. 845–849, 2008
2008
-
[13]
Cubature Kalman filters,
I. Arasaratnam and S. Haykin, “Cubature Kalman filters, ” IEEE Trans- actions on Automatic Control , vol. 54, no. 6, pp. 1254–1269, 2009
2009
-
[14]
Monte Carlo filters for non-linear state estimation,
E. Bølviken, P . J. Acklam, N. Christophersen, and J. M. S tørdal, “Monte Carlo filters for non-linear state estimation,” Automatica, vol. 37, no. 2, pp. 177–183, 2001
2001
-
[15]
Nonlinear s moothing theory,
C. T. Leondes, J. B. Peller, and E. B. Stear, “Nonlinear s moothing theory,” IEEE Transactions on Systems Science and Cybernetics , vol. 6, no. 1, pp. 63–71, 1970
1970
-
[16]
Cubature Kalman smoothe rs,
I. Arasaratnam and S. Haykin, “Cubature Kalman smoothe rs,” Automat- ica, vol. 47, no. 10, pp. 2245–2250, 2011
2011
-
[17]
On Gaussian optimal smoo thing of non- linear state space models,
S. Sarkka and J. Hartikainen, “On Gaussian optimal smoo thing of non- linear state space models,” IEEE Transactions on Automatic Control , vol. 55, no. 8, pp. 1938–1941, 2010
1938
-
[18]
Fourier-Hermite Rauch -Tung-Striebel smoother,
J. Sarmavuori and S. S¨ arkk¨ a, “Fourier-Hermite Rauch -Tung-Striebel smoother,” in 2012 Proceedings of the 20th European Signal Processing Conference (EUSIPCO) . IEEE, 2012, pp. 2109–2113
2012
-
[20]
Impulses and physiological states in the oretical models of nerve membrane,
R. FitzHugh, “Impulses and physiological states in the oretical models of nerve membrane,” Biophysical Journal , vol. 1, no. 6, pp. 445–466, 1961
1961
-
[21]
A novel and computationally efficient joint un scented Kalman filtering scheme for parameter estimation of a class of nonli near sys- tems,
A. Onat, “A novel and computationally efficient joint un scented Kalman filtering scheme for parameter estimation of a class of nonli near sys- tems,” IEEE Access , vol. 7, pp. 31 634–31 655, 2019
2019
-
[22]
Oscillatory Kalman filtering for Duffing, Coulomb, and Van d er Pol oscillators,
V . G. Y amalakonda, G. Kumar, R. B. Pachori, and A. K. Sing h, “Oscillatory Kalman filtering for Duffing, Coulomb, and Van d er Pol oscillators,” Signal Processing, vol. 211, p. 109091, 2023
2023
-
[23]
Stochastic integration filter,
J. Dun´ ık, O. Straka, and M. ˇSimandl, “Stochastic integration filter,” IEEE Transactions on Automatic Control, vol. 58, no. 6, pp. 1561–1566, 2013
2013
-
[24]
A brief look at Gaussian integrals,
W. O. Straub, “A brief look at Gaussian integrals,” Article: Pasadena California, 2009
2009
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.