Abstract:The results of a series of theoretical studies are reported, examining the convergence rate for different approximate representations of $alpha $ -stable distributions. Although they play a key role in modelling random processes with jumps and discontinuities, the use of $alpha $ -stable distributions in inference often leads to analytically intractable problems. The LePage series, which is a probabilistic representation employed in this work, is used to transform an intractable, infinite-dimensional inference problem into a finite-dimensional (conditionally Gaussian) parametric problem. A major component of our approach is the approximation of the tail of this series by a Gaussian random variable. Standard statistical techniques, such as Expectation-Maximization (EM), Markov chain Monte Carlo, and Particle Filtering, can then be readily applied. In addition to the asymptotic normality of the tail of this series, we establish explicit, nonasymptotic bounds on the approximation error. Their proofs follow classical Fourier-analytic arguments, using Esséen's smoothing lemma. Specifically, we consider the distance between the distributions of: $(i)$ the tail of the series and an appropriate Gaussian; $(ii)$ the full series and the truncated series; and $(iii)$ the full series and the truncated series with an added Gaussian term. In all three cases, sharp bounds are established, and the theoretical results are compared with the actual distances (computed numerically) in specific examples of symmetric $alpha-$ -stable distributions. This analysis facilitates the selection of appropriate truncations in practice and offers theoretical guarantees for the accuracy of resulting estimates. One of the main conclusions obtained is that, for the purposes of inference, the use of a truncated series together with an approximately Gaussian error term has superior statistical properties and is likely a preferable choice in practice.

Non-asymptotic Error Bounds for Sequential MCMC and Stability of Feynman-Kac Propagators

Sequential Monte Carlo samplers: error bounds and insensitivity to initial conditions

On the Behaviour of the Backward Interpretation of Feynman-Kac Formulae under Verifiable Conditions

Linear Variance Bounds for Particle Approximations of Time-Homogeneous Feynman-Kac Formulae

Non-asymptotic entropic bounds for non-linear kinetic Langevin sampler with second-order splitting scheme

Convergence Speed and Approximation Accuracy of Numerical MCMC

Asymptotic Stability of a Class of Nonlinear Stochastic Systems Undergoing Markovian Jumps

Convergence of sequential Markov Chain Monte Carlo methods: I. Nonlinear flow of probability measures

On the Particle Approximation of Lagged Feynman-Kac Formulae

Spectral Independence via Stability and Applications to Holant-Type Problems

Information metrics for long-time errors in splitting schemes for stochastic dynamics and parallel KMC

Explicit error bounds for lazy reversible Markov Chain Monte Carlo

Asymptotic bias of inexact Markov Chain Monte Carlo methods in high dimension

Bounding Wasserstein distance with couplings

Exponential Stability of Slowly Decaying Solutions to the Kinetic-Fokker-Planck Equation

A non-asymptotic error analysis for parallel Monte Carlo estimation from many short Markov chains

Convergence Bounds for Sequential Monte Carlo on Multimodal Distributions using Soft Decomposition

Short-time Fokker-Planck propagator beyond the Gaussian approximation

Statistical guarantees for stochastic Metropolis-Hastings

Quasi-ergodic theorems for Feynman-Kac semigroups and large deviation for additive functionals

Nonasymptotic Gaussian Approximation for Inference With Stable Noise