A nonlinear population Monte Carlo scheme for the Bayesian estimation of parameters of -stable distributions

DOI10.1016/J.CSDA.2015.09.007MaRDI QIDQ1659482zbMATH OpenOpenAlexFDO

Authors Eugenia Koblents, Joaquín Míguez, M. Rodríguez, Alexandra M. Schmidt

Publication date 15 August 2018

Published in Computational Statistics and Data Analysis (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1510.02702, http://hdl.handle.net/10016/25855

zbMATH Keywords

Bayesian inferenceimportance samplinganimal movement\(\alpha\)-stable distributionsLévy process

Mathematics Subject Classification ID

Infinitely divisible distributions; stable distributions (60E07) Processes with independent increments; Lévy processes (60G51) Computational methods for problems pertaining to statistics (62-08) Bayesian inference (62F15) Monte Carlo methods (65C05) Applications of statistics to environmental and related topics (62P12)

Abstract: The class of $alpha$-stable distributions enjoys multiple practical applications in signal processing, finance, biology and other areas because it allows to describe interesting and complex data patterns, such as asymmetry or heavy tails, in contrast with the simpler and widely used Gaussian distribution. The density associated with a general $alpha$-stable distribution cannot be obtained in closed form, which hinders the process of estimating its parameters. A nonlinear population Monte Carlo (NPMC) scheme is applied in order to approximate the posterior probability distribution of the parameters of an $alpha$-stable random variable given a set of random realizations of the latter. The approximate posterior distribution is computed by way of an iterative algorithm and it consists of a collection of samples in the parameter space with associated nonlinearly-transformed importance weights. A numerical comparison of the main existing methods to estimate the $alpha$-stable parameters is provided, including the traditional frequentist techniques as well as a Markov chain Monte Carlo (MCMC) and a likelihood-free Bayesian approach. It is shown by means of computer simulations that the NPMC method outperforms the existing techniques in terms of parameter estimation error and failure rate for the whole range of values of $alpha$, including the smaller values for which most existing methods fail to work properly. Furthermore, it is shown that accurate parameter estimates can often be computed based on a low number of observations. Additionally, numerical results based on a set of real fish displacement data are provided.

Recommendations

• Parameter estimation of \(\alpha\)-stable distributions using a DRAM algorithm
• Bayesian inference for \(\alpha \)-stable distributions: a random walk MCMC approach
• Likelihood-free Bayesian inference for \(\alpha\)-stable models
• Bayesian Inference for Stable Distributions
• Parameter Estimation of Stable Distributions

Cites work

• scientific article; zbMATH DE number 1639862 (Why is no real title available?)
• scientific article; zbMATH DE number 1301883 (Why is no real title available?)
• scientific article; zbMATH DE number 2117879 (Why is no real title available?)
• A Method for Simulating Stable Random Variables
• A tutorial on approximate Bayesian computation
• Adaptive approximate Bayesian computation
• An adaptive sequential Monte Carlo method for approximate Bayesian computation
• An iterative procedure for the estimation of the parameters of stable laws
• Bayesian Inference for Stable Distributions
• Bayesian inference for \(\alpha \)-stable distributions: a random walk MCMC approach
• Calibrated FFT-based density approximations for \(\alpha\)-stable distributions
• Curse-of-dimensionality revisited: Collapse of the particle filter in very large scale systems
• Density parameter estimation of skewed α-stable distributions
• Handbook of computational statistics. Concepts and methods.
• Likelihood-free Bayesian inference for \(\alpha\)-stable models
• Monte Carlo sampling methods using Markov chains and their applications
• Numerical calculation of stable densities and distribution functions
• Parameter Estimates for Symmetric Stable Distributions
• Particle-kernel estimation of the filter density in state-space models
• Sequential Monte Carlo without likelihoods
• Simple consistent estimators of stable distribution parameters

Cited in

(9)

• On variance stabilisation in population Monte Carlo by double Rao-Blackwellisation
• Likelihood-free Bayesian inference for \(\alpha\)-stable models
• Bayesian inference for \(\alpha \)-stable distributions: a random walk MCMC approach
• Flexible two-point selection approach for characteristic function-based parameter estimation of stable laws
• Fast parallel \(\alpha \)-stable distribution function evaluation and parameter estimation using OpenCL in GPGPUs
• Parameter estimation of \(\alpha\)-stable distributions using a DRAM algorithm
• Goodness-of-fit test for \(\alpha\)-stable distribution based on the quantile conditional variance statistics
• Nonparametric estimation of the kernel function of symmetric stable moving average random functions
• Bayesian computation methods for inference in stochastic kinetic models

Describes a project that uses

Uses Software

• STABLE
• Matlab

This page was built for publication: A nonlinear population Monte Carlo scheme for the Bayesian estimation of parameters of \(\alpha\)-stable distributions

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1659482)