A central limit theorem with application to inference in α-stable regression models

It is well known that the α-stable distribution, while having no closed form density function in the general case, admits a Poisson series representation (PSR) in which the terms of the series are a function of the arrival times of a unit rate Poisson process. In our previous work we have shown how to carry out inference for regression models using this series representation, which leads to a very convenient conditionally Gaussian framework, amenable to straightforward Gaussian inference procedures. The PSR has to be truncated to a finite number of terms for practical purposes. The residual terms have been approximated in our previous work by a Gaussian distribution with fully characterised moments. In this paper we present a new Central Limit Theorem (CLT) for the residual terms which serves to justify our previous approximation of the residual as Gaussian. Furthermore, we provide an analysis of the asymptotic convergence rate expressed in the CLT.

[1]  Simon J. Godsill,et al.  Enhanced Poisson sum representation for alpha-stable processes , 2011, 2011 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP).

[2]  Simon J. Godsill,et al.  Fully Bayesian inference for α-stable distributions using a Poisson series representation , 2015, Digit. Signal Process..

[3]  Simon J. Godsill,et al.  Linear gaussian computations for near-exact Bayesian Monte Carlo inference in skewed alpha-stable time series models , 2012, 2012 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP).

[4]  P. Lee,et al.  14. Simulation and Chaotic Behaviour of α‐Stable Stochastic Processes , 1995 .

[5]  R. Katz,et al.  Extreme events in a changing climate: Variability is more important than averages , 1992 .

[6]  E. Fama The Behavior of Stock-Market Prices , 1965 .

[7]  Josiane Zerubia,et al.  SAR image filtering based on the heavy-tailed Rayleigh model , 2006, IEEE Transactions on Image Processing.

[8]  Raoul LePage,et al.  Multidimensional infinitely divisidle variables and processes Part II , 1981 .

[9]  W. Linde STABLE NON‐GAUSSIAN RANDOM PROCESSES: STOCHASTIC MODELS WITH INFINITE VARIANCE , 1996 .

[10]  Friedrich Götze,et al.  Bounds for the accuracy of Poissonian approximations of stable laws , 1996 .

[11]  B. Mandelbrot New Methods in Statistical Economics , 1963, Journal of Political Economy.

[12]  K. Chung,et al.  Limit Distributions for Sums of Independent Random Variables. , 1955 .

[13]  C. Mallows,et al.  A Method for Simulating Stable Random Variables , 1976 .

[14]  Marco J. Lombardi,et al.  On-line Bayesian Estimation of Signals in Symmetric α-Stable Noise , 2004 .

[15]  M. Taqqu,et al.  Stable Non-Gaussian Random Processes : Stochastic Models with Infinite Variance , 1995 .

[16]  Alin Achim,et al.  Novel Bayesian multiscale method for speckle removal in medical ultrasound images , 2001, IEEE Transactions on Medical Imaging.

[17]  Raoul LePage,et al.  Appendix Multidimensional infinitely divisible variables and processes. Part I: Stable case , 1980 .

[18]  Joel Zinn,et al.  Convergence to a Stable Distribution Via Order Statistics , 1981 .

[19]  C. L. Nikias,et al.  Signal processing with alpha-stable distributions and applications , 1995 .

[20]  Simon J. Godsill,et al.  A poisson series approach to Bayesian Monte Carlo inference for skewed alpha-stable distributions , 2014, 2014 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP).

[21]  J. Lindeberg Eine neue Herleitung des Exponentialgesetzes in der Wahrscheinlichkeitsrechnung , 1922 .

[22]  M. Parlange,et al.  Statistics of extremes in hydrology , 2002 .

[23]  P. Levy,et al.  Calcul des Probabilites , 1926, The Mathematical Gazette.