Testing and estimation of the index of stability of univariate and bivariate symmetric α-stable distributions via modified Greenwood statistic
Pith reviewed 2026-05-10 09:04 UTC · model grok-4.3
The pith
The modified Greenwood statistic enables testing for bivariate symmetric alpha-stable distributions and estimation of the stability index in univariate cases.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We extend the modified Greenwood statistic to a bivariate setting and examine its probabilistic properties within the class of alpha-stable distributions, with a focus on the sub-Gaussian case. We introduce a novel testing approach that considers two variations of the modified Greenwood statistic as test statistics for the bivariate case. In the univariate setting, we adapt the proposed testing methodology for estimating the stability index. The simulation studies presented demonstrate that our proposed methodology outperforms classical approaches previously used in this context and serves as an effective tool for distinguishing between Gaussian and alpha-stable distributions with a index of
What carries the argument
The modified Greenwood statistic for symmetric distributions, extended to bivariate samples to serve as a test statistic and basis for stability index estimation.
If this is right
- Enables construction of tests for whether bivariate observations follow a joint symmetric alpha-stable law.
- Provides a direct estimator for the stability index from univariate samples.
- Yields improved ability to separate alpha-stable laws from Gaussian ones when the index is close to 2.
- Supplies explicit test statistics based on two variants of the modified Greenwood statistic for the bivariate setting.
Where Pith is reading between the lines
- The bivariate extension opens the possibility of joint modeling of heavy-tailed pairs of observations in applications such as portfolio risk or paired sensor data.
- Because the method focuses on the sub-Gaussian regime, it may serve as a template for similar statistics on other symmetric heavy-tailed families.
- Practical deployment would require checking robustness when the data only approximately follows the symmetric alpha-stable assumption.
Load-bearing premise
Data samples are drawn from symmetric alpha-stable distributions and the modified Greenwood statistic has the examined probabilistic properties within this class, particularly in the sub-Gaussian case.
What would settle it
A simulation experiment in which the proposed bivariate tests or univariate estimator fails to outperform classical methods when distinguishing Gaussian data from alpha-stable data with stability index near 2 would falsify the performance claim.
Figures
read the original abstract
We propose a testing and estimation methodology for univariate and bivariate symmatric $\alpha$-stable distributions using a modified version of the Greenwood statistic. Originally designed for positive-valued random variables, the Greenwood statistic, and its modified version tailored for symmetric distributions, have been predominantly applied to univariate random samples. In this paper, we extend the modified Greenwood statistic to a bivariate setting and examine its probabilistic properties within the class of $\alpha$-stable distributions, with a focus on the sub-Gaussian case. Additionally, we introduce a novel testing approach that considers two variations of the modified Greenwood statistic as test statistics for the bivariate case. In the univariate setting, we adapt the proposed testing methodology for estimating the stability index. The simulation studies presented demonstrate that our proposed methodology outperforms classical approaches previously used in this context and serves as an effective tool for distinguishing between Gaussian and $\alpha$-stable distributions with a stability index close to 2. The theoretical and simulation results are further illustrated with practical data examples.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes extending the modified Greenwood statistic to bivariate symmetric α-stable distributions and examines its probabilistic properties (with focus on the sub-Gaussian case). It develops testing procedures using two variants of the statistic for the bivariate setting, adapts the approach for estimating the stability index α in the univariate case, and reports simulation results claiming outperformance over classical methods, particularly for distinguishing Gaussian (α=2) from α-stable laws with α close to 2. Theoretical results and simulations are illustrated on practical data examples.
Significance. If the probabilistic properties are shown to hold uniformly for α ∈ (1,2] and the simulation comparisons are reproducible with full controls, the work could provide a useful addition to the limited set of tools for inference on near-Gaussian stable laws, which arise in financial modeling and signal processing. The bivariate extension and data examples add practical value, though the reliance on simulations for the key distinction claim limits the strength of the contribution.
major comments (3)
- [Section 3] Section 3 (probabilistic properties): The derivations and limiting results for the modified Greenwood statistic are presented with explicit focus on the sub-Gaussian case. However, the central claim requires reliable behavior for α close to but strictly less than 2, where the tails are heavier and quadratic functionals may lack the same integrability. No continuity argument, uniform integrability result, or separate derivation for α ∈ (1,2) is provided to justify extrapolation of the test statistic's distribution or moments to the near-Gaussian alternatives used in the simulations.
- [Section 5] Section 5 (simulation studies): The reported outperformance in distinguishing α=2 from α near 2 and in estimation accuracy lacks accompanying standard errors, confidence intervals, or details on the number of Monte Carlo replications, exact grid of α values tested, and data-generation parameters. Without these, it is impossible to assess whether the claimed superiority is statistically significant or sensitive to simulation design choices.
- [Section 4] Section 4 (bivariate testing): The two proposed test statistics for the bivariate case are introduced without a clear power analysis or comparison to existing bivariate stable tests under local alternatives near α=2. This leaves the practical advantage of the bivariate extension unsubstantiated beyond the univariate estimation results.
minor comments (3)
- [Abstract] The abstract contains a typographical error ('symmatric' instead of 'symmetric').
- [Section 2] Notation for the modified Greenwood statistic in the bivariate setting should be introduced with an explicit definition before its use in the testing procedures.
- [Discussion] The paper would benefit from a brief discussion of how the proposed estimator behaves when the underlying distribution deviates from exact symmetry, even if the main focus is symmetric α-stable laws.
Simulated Author's Rebuttal
We thank the referee for the constructive comments, which help clarify the scope and strengthen the presentation of our results. We respond point by point to the major comments below.
read point-by-point responses
-
Referee: [Section 3] Section 3 (probabilistic properties): The derivations and limiting results for the modified Greenwood statistic are presented with explicit focus on the sub-Gaussian case. However, the central claim requires reliable behavior for α close to but strictly less than 2, where the tails are heavier and quadratic functionals may lack the same integrability. No continuity argument, uniform integrability result, or separate derivation for α ∈ (1,2) is provided to justify extrapolation of the test statistic's distribution or moments to the near-Gaussian alternatives used in the simulations.
Authors: The derivations in Section 3 are deliberately restricted to the sub-Gaussian case because this setting yields explicit limiting distributions via the properties of multivariate Gaussians and the associated quadratic forms. For the univariate estimation procedure and the simulation study of behavior near α = 2, we rely on the empirical performance shown in Section 5. We agree that a formal bridge between the sub-Gaussian limit and α ∈ (1,2) is not supplied. In the revision we will insert a short remark in Section 3 invoking the continuous mapping theorem together with the fact that symmetric α-stable laws converge weakly to the Gaussian law as α → 2^−; under the moment conditions already verified for the Greenwood statistic this implies convergence of the relevant functionals. This addition will be limited to a remark rather than a full re-derivation. revision: partial
-
Referee: [Section 5] Section 5 (simulation studies): The reported outperformance in distinguishing α=2 from α near 2 and in estimation accuracy lacks accompanying standard errors, confidence intervals, or details on the number of Monte Carlo replications, exact grid of α values tested, and data-generation parameters. Without these, it is impossible to assess whether the claimed superiority is statistically significant or sensitive to simulation design choices.
Authors: We accept that the simulation section must be made fully reproducible. The revised manuscript will report: 10 000 Monte Carlo replications for each configuration, the precise grid α = 1.01, 1.05, …, 1.95, 2.00, Monte Carlo standard errors and 95 % confidence intervals for every rejection rate and MSE value, and the exact data-generation protocol (Chambers–Mallows–Stuck algorithm with fixed random seeds). These additions will allow readers to judge both statistical significance and sensitivity to design choices. revision: yes
-
Referee: [Section 4] Section 4 (bivariate testing): The two proposed test statistics for the bivariate case are introduced without a clear power analysis or comparison to existing bivariate stable tests under local alternatives near α=2. This leaves the practical advantage of the bivariate extension unsubstantiated beyond the univariate estimation results.
Authors: Section 4 presents the bivariate extension together with the null limiting distributions under the sub-Gaussian assumption; the data examples illustrate its use. We acknowledge that a dedicated power comparison against existing bivariate procedures (e.g., characteristic-function or likelihood-based tests) under local alternatives is absent. In the revision we will add a concise simulation subsection (or appendix) that evaluates empirical power of the two proposed statistics for α = 1.90, 1.95, 1.99 against a benchmark empirical-characteristic-function test, using the same Monte Carlo settings as Section 5. This will directly address the practical advantage of the bivariate extension. revision: yes
Circularity Check
No significant circularity; derivation and validation are independent
full rationale
The paper extends the modified Greenwood statistic to the bivariate case for symmetric alpha-stable laws, examines its probabilistic properties (with explicit focus on the sub-Gaussian regime), introduces corresponding test statistics and an estimation procedure for the stability index, and reports simulation results showing outperformance versus classical methods. No load-bearing step reduces a claimed prediction or uniqueness result to a fitted parameter, self-definition, or self-citation chain by the paper's own equations. The simulation-based validation constitutes an external check rather than a tautology, rendering the central claims self-contained against the stated assumptions and benchmarks.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Symmetric alpha-stable distributions possess well-defined probabilistic properties that allow extension of the modified Greenwood statistic, especially in the sub-Gaussian case.
Reference graph
Works this paper leans on
-
[1]
K. Skowronek, M. Arendarczyk, R. Zimroz, A. Wyłomańska, Modified Greenwood statistic and its ap- plication for statistical testing, Journal of Computation al and Applied Mathematics 452 (2024) 116122
work page 2024
-
[2]
M. Greenwood, The statistical study of infectious disea ses, Journal of the Royal Statistical Society 109 (2) (1946) 85–110
work page 1946
-
[3]
P. A. P. Moran, The random division of an interval, Supple ment to the Journal of the Royal Statistical Society 9 (1) (1947) 92
work page 1947
-
[4]
R. B. D’Agostino, M. A. Stephens, Goodness-of-fit techni ques, Marcel Dekker, Inc., USA, 1986
work page 1986
-
[5]
H. Albrecher, J. Teugels, Asymptotic analysis of a measu re of variation, Theory of Probability and Mathematical Statistics 74 (2007) 1–10
work page 2007
-
[6]
H. Albrecher, J. L. Teugels, K. Scheicher, A combinatori al identity for a problem in asymptotic statistics, Applicable Analysis and Discrete Mathematics 3 (1) (2009) 6 4–68
work page 2009
-
[7]
H. Albrecher, S. A. Ladoucette, J. L. Teugels, Asymptoti cs of the sample coefficient of variation and the sample dispersion, Journal of Statistical Planning and Inference 140 (2) (2010) 358–368
work page 2010
-
[8]
M. Arendarczyk, T. J. Kozubowski, A. K. Panorska, The Gre enwood statistic, stochastic dominance, clustering and heavy tails, Scandinavian Journal of Statis tics 49 (1) (2022) 331–352
work page 2022
-
[9]
M. Arendarczyk, T. J. Kozubowski, A. K. Panorska, A Compu tational Approach to Confidence Intervals and Testing for Generalized Pareto Index Using the Greenwoo d Statistic, REVSTAT-Statistical Journal 21 (3) (2023) 367–388
work page 2023
- [10]
- [11]
-
[12]
H. Albrecher, B. García Flores, Asymptotic analysis of generalized Greenwood statistics for very heavy tails, Statistics & Probability Letters 185 (2022) 109429
work page 2022
-
[13]
J. S. Rao, M. Kuo, Asymptotic Results on the Greenwood St atistic and Some of its Generalizations, Journal of the Royal Statistical Society: Series B (Methodo logical) 46 (2) (1984) 228–237
work page 1984
-
[14]
H. Scott Hurd, J. B. Kaneene, The application of simulat ion models and systems analysis in epidemi- ology: a review, Preventive Veterinary Medicine 15 (2) (199 3) 81–99
-
[15]
M. C. Riley, A. Clare, R. D. King, Locational distributi on of gene functional classes in arabidopsis thaliana, BMC Bioinformatics 8 (1) (2007) 112
work page 2007
-
[16]
B. D. Peterson-Burch, D. Nettleton, D. F. Voytas, Genom ic neighborhoods for arabidopsisretrotrans- posons: a role for targeted integration in the distribution of the metaviridae, Genome Biology 5 (10) (2004) R78
work page 2004
-
[17]
S. L. DeRuiter, I. L. Boyd, D. E. Claridge, C. W. Clark, C. Gagnon, B. L. Southall, P. L. Tyack, Delphinid whistle production and call matching during play back of simulated military sonar, Marine Mammal Science 29 (2) (2013) E46–E59
work page 2013
-
[18]
F. Moscone, E. Tosetti, Testing for error cross section independence with an application to us health expenditure, Regional Science and Urban Economics 40 (5) (2 010) 283–291, advances In Spatial Econo- metrics. 13
-
[19]
N. Benlagha, W. Hemrit, Does investment in insurance st ocks reap diversification benefits? Static and time varying copula modeling, Communications in Statis tics - Simulation and Computation 52 (4) (2023) 1384–1402
work page 2023
-
[20]
J. del Castillo, I. Serra, Likelihood inference for gen eralized pareto distribution, Computational Statis- tics & Data Analysis 83 (2015) 116–128
work page 2015
-
[21]
D. Vermetten, B. van Stein, F. Caraffini, L. L. Minku, A. V. Kononova, Bias: A toolbox for bench- marking structural bias in the continuous domain, IEEE Tran sactions on Evolutionary Computation 26 (6) (2022) 1380–1393
work page 2022
- [22]
-
[23]
R. Pakyari, N. Balakrishnan, Goodness-of-fit tests for progressively type-ii censored data from loca- tion–scale distributions, Journal of Statistical Computa tion and Simulation 83 (1) (2013) 167–178
work page 2013
-
[24]
O. Shchur, A. C. Turkmen, T. Januschowski, J. Gasthaus, S. Günnemann, Detecting anomalous event sequences with temporal point processes, in: M. Ranzato, A. Beygelzimer, Y. Dauphin, P. Liang, J. W. Vaughan (Eds.), Advances in Neural Information Processing Systems, Vol. 34, Curran Associates, Inc., 2021, pp. 13419–13431
work page 2021
- [25]
-
[26]
L. Henriques-Rodrigues, M. I. Gomes, D. Pestana, Stati stics of extremes in athletics, REVSTAT- Statistical Journal 9 (2) (2011) 127–153
work page 2011
-
[27]
G. Samoradnitsky, M. S. Taqqu, Stable Non-Gaussian Ran dom Processes: Stochastic Models with Infinite Variance, Routledge, 1994
work page 1994
- [28]
-
[29]
A. Janicki, A. Weron, Simulation and Chaotic Behavior o f Alpha-stable Stochastic Processes, Marcel Dekker, Inc., New York, 1994
work page 1994
-
[30]
E. J. Wegman, S. C. Schwartz, J. B. Thomas (Eds.), Topics in Non-Gaussian Signal Processing, New York: Springer, 1989
work page 1989
-
[31]
J. P. Nolan, Univariate Stable Distributions. Models f or Heavy Tailed Data, Springer, 2020
work page 2020
-
[32]
B. Mandelbrot, The Pareto-Lévy Law and the distributio n of income, International Economic Review 1 (2) (1960) 79–106
work page 1960
-
[33]
S. T. Rachev, S. Mittnik, Stable Paretian Models in Fina nce, Wiley, New York, 2000
work page 2000
-
[34]
P. V. Bidarkota, B. V. Dupoyet, J. H. McCulloch, Asset pr icing with incomplete information and fat tails, Journal of Economic Dynamics and Control 33(6) (2009 ) 1314–1331
work page 2009
- [35]
- [36]
-
[37]
P. Barthelemy, J. Bertolotti, D. S. Wiersma, A Lévy fligh t for light, Nature 453 (2008) 495–498. 14
work page 2008
-
[38]
I. M. Sokolov, J. Mai, A. Blumen, Paradoxal diffusion in c hemical space for nearest-neighbor walks over polymer chains, Physical Review Letters 79 (1997) 857
work page 1997
-
[39]
M. A. Lomholt, T. Ambjornsson, R. Metzler, Optimal targ et search on a fast-folding polymer chain with volume exchange, Physical Review Letter 95 (2005) 2606 03
work page 2005
-
[40]
R. Durrett, J. Foo, K. Leder, J. Mayberry, F. Michor, Int ratumor heterogeneity in evolutionary models of tumor progression, Genetics 188 (2011) 1–17
work page 2011
-
[41]
B. L. Lan, M. Toda, Fluctuations of healthy and unhealth y heartbeat intervals, Europhysics Letters 102(1) (2013) 18002–p1–p6
work page 2013
-
[42]
C.-K. Peng, J. Mietus, J. M. Hausdorff, S. Havlin, H. E. St anley, A. L. Goldberger, Long-range anti- correlations and non-Gaussian behavior of the heartbeat, P hysical Review Letters 70 (1993) 1343
work page 1993
-
[43]
P. D. Ditlevsen, Observation of alpha-stable noise ind uced millennial climate changes from an ice-core record, Geophysical Research Letters 26 (1999) 1441–1444
work page 1999
-
[44]
D. Middleton, Non-Gaussian Noise Models in Signal Proc essing for Telecommunications: New Methods and Results for Class A and Class B Noise Models, IEEE Transac tions on Information Theory 45 (1999) 1129–1149
work page 1999
-
[45]
R. Li, Z. Zhao, Y. Zhong, C. Qi, H. Zhang, The stochastic g eometry analyses of cellular networks with alpha-stable self-similarity, IEEE Transactions on Commu nications 67(3) (2019) 2487–2503
work page 2019
-
[46]
J. Janczura, T. Barszcz, R. Zimroz, A. Wyłomańska, Mach ine condition change detection based on data segmentation using a three-regime, α -stable hidden markov model, Measurement in press (2023)
work page 2023
-
[47]
G. Yu, C. Li, J. Zhang, A new statistical modeling and det ection method for rolling element bearing faults based on alpha-stable distribution, Mechanical Sys tems and Signal Processing 41 (2013) 155–175
work page 2013
-
[48]
W. Żuławiński, K. Maraj-Zygmąt, H. Shiri, A. Wyłomańsk a, R. Zimroz, Framework for stochastic modelling of long-term non-homogenous data with non-Gauss ian characteristics for machine condition prognosis, Mechanical Systems and Signal Processing 184 (1 84) (2023) 109677
work page 2023
-
[49]
A. Grzesiek, A. Michalak, A. Wyłomańska, How to describ e the linear dependence for heavy-tailed distributed data, in: F. Chaari, J. Leskow, A. Wylomanska, R . Zimroz, A. Napolitano (Eds.), Nonsta- tionary Systems: Theory and Applications, Springer Intern ational Publishing, Cham, 2022, pp. 69–92
work page 2022
-
[50]
M. Teimouri, S. Rezakhah, A. Mohammadpour, Parameter e stimation using the em algorithm for symmetric stable random variables and sub-Gaussian random vectors, Journal of Statistical Theory and Applications 17 (3) (2018) 439–461
work page 2018
-
[51]
M. Jabłońska-Sabuka, M. Teuerle, A. Wyłomańska, Bivar iate sub-Gaussian model for stock index re- turns, Physica A: Statistical Mechanics and its Applicatio ns 486 (2017) 628–637
work page 2017
- [52]
-
[53]
S. O. Lozza, T. Lando, F. Petronio, T. Tichý, Asymptotic multivariate dominance: A financial appli- cation, Methodology and Computing in Applied Probability 1 8 (4) (2016) 1097–1115
work page 2016
-
[54]
S. Ortobelli, T. Tichý, On the impact of semidefinite pos itive correlation measures in portfolio theory, Annals of Operations Research 235 (1) (2015) 625–652
work page 2015
-
[55]
A. Mahmood, M. Chitre, Generating random variates for s table sub-Gaussian processes with memory, Signal Processing 131 (2017) 271–279. 15
work page 2017
-
[56]
P. Hao, O. Karakuş, A. Achim, Robust Kalman filters based on the sub-Gaussian α -stable distribution, Signal Processing 224 (2024) 109574
work page 2024
-
[57]
A. Guadagnini, M. Riva, S. P. Neuman, Recent advances in scalable non-Gaussian geostatistics: The generalized sub-Gaussian model, Journal of Hydrology 562 ( 2018) 685–691
work page 2018
-
[58]
A. Kabašinskas, L. Sakalauskas, I. Vaičiulyt˙ e, An ana lytical em algorithm for sub-Gaussian vectors, Mathematics 9 (9) (2021)
work page 2021
-
[59]
K. V. Mardia, Measures of multivariate skewness and kur tosis with applications, Biometrika 57 (3) (1970) 519–530
work page 1970
-
[60]
C. M. Jarque, A. K. Bera, A test for normality of observat ions and regression residuals, International Statistical Review / Revue Internationale de Statistique 5 5 (2) (1987) 163–172
work page 1987
- [61]
-
[62]
J. P. Royston, Some Techniques for Assessing Multivara te Normality Based on the Shapiro-Wilk W, Journal of the Royal Statistical Society Series C: Applied S tatistics 32 (2) (2018) 121–133
work page 2018
-
[63]
A. M. Mathai, S. B. Provost, Quadratic forms in random va riables: theory and applications, Marcel Dekker, 1992
work page 1992
-
[64]
Łukasz Bielak, A. Grzesiek, J. Janczura, A. Wyłomańska , Market risk factors analysis for an interna- tional mining company. multi-dimensional, heavy-tailed- based modelling, Resources Policy 74 (2021) 102308
work page 2021
-
[65]
L. Rabiner, A tutorial on hidden markov models and selec ted applications in speech recognition, Pro- ceedings of the IEEE 77 (2) (1989) 257–286
work page 1989
-
[66]
A. Grzesiek, M. Mrozińska, P. Giri, S. Sundar, A. Wyłoma ńska, The covariation-based yule–walker method for multidimensional autoregressive time series wi th α -stable distributed noise, International Journal of Advances in Engineering Sciences and Applied Mat hematics 13 (4) (2021) 394–414
work page 2021
-
[67]
K. Maraj-Zygmat, G. Sikora, M. Pitera, A. Wyłomańska, G oodness-of-fit test for stochastic processes using even empirical moments statistic, Chaos: An Interdis ciplinary Journal of Nonlinear Science 33 (1) (2023) 013128
work page 2023
-
[68]
K. Skowronek, R. Zimroz, A. Wyłomańska, Testing for fini te variance with applications to vi- bration signals from rotating machines, Journal of Mathema tics in Industry 14 (19) (2024) https://doi.org/10.1186/s13362–024–00157–6
-
[69]
K. Skowronek, T. Barszcz, J. Antoni, R. Zimroz, A. Wyłom ańska, Assessment of background noise properties in time and time–frequency domains in the contex t of vibration-based local damage detection in real environment, Mechanical Systems and Signal Process ing 199 (2023) 110465. 16
work page 2023
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.