pith. sign in

arxiv: 2511.07059 · v2 · submitted 2025-11-10 · 📊 stat.ME

Applying the Polynomial Maximization Method to Estimate ARIMA Models with Asymmetric Non-Gaussian Innovations

Pith reviewed 2026-05-17 23:55 UTC · model grok-4.3

classification 📊 stat.ME
keywords ARIMA modelspolynomial maximizationnon-Gaussian innovationssemiparametric estimationasymmetric distributionsMonte Carlo simulationtime series estimationhigher-order moments
0
0 comments X

The pith

PMM2 delivers more efficient ARIMA estimates for asymmetric non-Gaussian innovations

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper introduces the second-order polynomial maximization method (PMM2) as a way to estimate ARIMA models when the error terms follow asymmetric distributions rather than the usual Gaussian assumption. Classical estimators lose efficiency in such cases, which are common in financial and economic time series. PMM2 works by using higher-order moments and cumulants in a semiparametric way, without needing the complete probability distribution of the innovations. Large-scale simulations across different sample sizes and distributions show clear gains in precision for skewed innovations, with relative efficiencies up to 1.9, while performing equally well as standard methods when the Gaussian assumption holds. The approach therefore offers a direct improvement for modeling real data that deviates from normality.

Core claim

The second-order polynomial maximization method provides a semiparametric estimator for the parameters of ARIMA(p, d, q) processes that utilizes cumulants beyond the second order to achieve higher efficiency than maximum likelihood, conditional sum of squares, or ordinary least squares when the innovations are asymmetric and non-Gaussian, as validated by Monte Carlo experiments involving 128,000 simulations.

What carries the argument

The second-order polynomial maximization method (PMM2), a semiparametric estimation technique that forms polynomial estimating equations from the observed series to incorporate skewness information from the innovations.

If this is right

  • For ARIMA(1,1,0) with sample size 500, PMM2 achieves 37 to 47 percent variance reduction relative to classical methods for Gamma, lognormal, and chi-squared innovations.
  • PMM2 attains the same statistical efficiency as ordinary least squares when the innovations are Gaussian.
  • The performance advantage appears for sample sizes of at least 200 and skewness levels of 0.5 or higher.
  • Computational requirements remain similar to those of maximum likelihood estimation.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Applying PMM2 to actual financial return series could yield more accurate volatility estimates and better risk assessments.
  • Integrating PMM2 with GARCH-type models might address both innovation asymmetry and time-varying variance in one framework.
  • Developing automatic order selection procedures compatible with PMM2 would simplify its use in practice.

Load-bearing premise

The series is correctly identified as an ARIMA process and the innovations have finite moments of order four or higher that allow the polynomial method to capture asymmetry.

What would settle it

If PMM2 is applied to data generated from an ARIMA model with innovations having infinite fourth moments, the efficiency advantage should disappear or reverse compared to classical estimators.

Figures

Figures reproduced from arXiv: 2511.07059 by Serhii Zabolotnii.

Figure 1
Figure 1. Figure 1: Relative efficiency of PMM2 versus CSS as a function of residual skewness [PITH_FULL_IMAGE:figures/full_fig_p018_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Normalised quality metrics (AIC, BIC, RMSE, MAE) for CSS-ML and PMM2 across [PITH_FULL_IMAGE:figures/full_fig_p020_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Differences between PMM2 and CSS-ML (negative values favour PMM2) for informa [PITH_FULL_IMAGE:figures/full_fig_p021_3.png] view at source ↗
read the original abstract

Classical estimators for ARIMA parameters (MLE, CSS, OLS) assume Gaussian innovations, an assumption frequently violated in financial and economic data exhibiting asymmetric distributions with heavy tails. We develop and validate the second-order polynomial maximization method (PMM2) for estimating ARIMA$(p,d,q)$ models with non-Gaussian innovations. PMM2 is a semiparametric technique that exploits higher-order moments and cumulants without requiring full distributional specification. Monte Carlo experiments (128,000 simulations) across sample sizes $N \in \{100, 200, 500, 1000\}$ and four innovation distributions demonstrate that PMM2 substantially outperforms classical methods for asymmetric innovations. For ARIMA(1,1,0) with $N=500$, relative efficiency reaches 1.58--1.90 for Gamma, lognormal, and $\chi^2(3)$ innovations (37--47\% variance reduction). Under Gaussian innovations PMM2 matches OLS efficiency, avoiding the precision loss typical of robust estimators. The method delivers major gains for moderate asymmetry ($|\gamma_3| \geq 0.5$) and $N \geq 200$, with computational costs comparable to MLE. PMM2 provides an effective alternative for time series with asymmetric innovations typical of financial markets, macroeconomic indicators, and industrial measurements. Future extensions include seasonal SARIMA models, GARCH integration, and automatic order selection.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript develops and validates the second-order polynomial maximization method (PMM2), a semiparametric estimator for ARIMA(p,d,q) parameters that exploits higher-order moments and cumulants without requiring a full distributional specification for the innovations. Monte Carlo experiments comprising 128,000 simulations across N in {100,200,500,1000} and four innovation distributions (Gaussian, Gamma, lognormal, χ²(3)) are used to demonstrate that PMM2 yields relative efficiencies of 1.58–1.90 (37–47% variance reduction) versus classical OLS/MLE/CSS for asymmetric innovations in the ARIMA(1,1,0) case at N=500, while matching OLS efficiency under Gaussian innovations.

Significance. If the efficiency gains hold without offsetting bias, PMM2 would constitute a practical semiparametric alternative for ARIMA estimation in financial, macroeconomic, and industrial time series that routinely exhibit asymmetry and heavy tails. The fact that the method incurs no precision loss under Gaussianity distinguishes it from many robust alternatives and, together with computational cost comparable to MLE, supports its potential adoption for moderate asymmetry (|γ₃| ≥ 0.5) and N ≥ 200.

major comments (2)
  1. [Monte Carlo experiments] Monte Carlo experiments section: only relative-efficiency (variance-ratio) results are presented for the asymmetric cases; no bias, MSE, or coverage diagnostics are reported for Gamma, lognormal, or χ²(3) innovations. Because PMM2 exploits higher cumulants that interact with the ARIMA lag structure after differencing, any finite-sample bias could offset the reported 37–47% variance reduction, rendering the net improvement in estimation accuracy unclear.
  2. [Abstract and results] Abstract and results: the statement that PMM2 “matches OLS efficiency” under Gaussian innovations is given without an accompanying table or explicit variance comparison that would confirm the absence of any small-sample penalty; this comparison is load-bearing for the claim that PMM2 avoids the typical precision loss of robust estimators.
minor comments (2)
  1. [Abstract] The abstract does not define the precise formula used for relative efficiency or state the number of Monte Carlo replications per (N, distribution, model) cell.
  2. [Methods] Notation for the innovation cumulants (γ₃, etc.) should be introduced explicitly in the methods section before being used in the efficiency discussion.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments on our manuscript. We address each major comment below and outline the revisions we will make to improve the clarity and completeness of the Monte Carlo validation.

read point-by-point responses
  1. Referee: [Monte Carlo experiments] Monte Carlo experiments section: only relative-efficiency (variance-ratio) results are presented for the asymmetric cases; no bias, MSE, or coverage diagnostics are reported for Gamma, lognormal, or χ²(3) innovations. Because PMM2 exploits higher cumulants that interact with the ARIMA lag structure after differencing, any finite-sample bias could offset the reported 37–47% variance reduction, rendering the net improvement in estimation accuracy unclear.

    Authors: We agree that a fuller set of diagnostics would strengthen the results. Although the experiments emphasize relative efficiency to quantify the variance reductions achieved by PMM2, we recognize that bias or MSE could offset those gains in finite samples. In the revised manuscript we will add tables reporting bias, MSE, and coverage probabilities for all four innovation distributions (including the three asymmetric cases) across the full range of sample sizes. These additions will allow direct assessment of whether the reported efficiency gains translate into net improvements in estimation accuracy. revision: yes

  2. Referee: [Abstract and results] Abstract and results: the statement that PMM2 “matches OLS efficiency” under Gaussian innovations is given without an accompanying table or explicit variance comparison that would confirm the absence of any small-sample penalty; this comparison is load-bearing for the claim that PMM2 avoids the typical precision loss of robust estimators.

    Authors: We accept that an explicit side-by-side comparison is needed to substantiate the claim. The Monte Carlo design includes the Gaussian case, but the presentation focuses on the asymmetric results. We will insert a concise table in the results section that reports the variance ratios (or relative efficiencies) of PMM2 versus OLS for the Gaussian innovations at each sample size. This will confirm the absence of any small-sample penalty and directly support the statement that PMM2 matches OLS efficiency under Gaussianity. revision: yes

Circularity Check

0 steps flagged

No significant circularity; performance claims rest on independent Monte Carlo validation

full rationale

The paper introduces PMM2 as a semiparametric estimator that exploits higher-order moments and cumulants for ARIMA models with asymmetric innovations, then validates its relative efficiency through 128,000 independent Monte Carlo simulations across multiple distributions, sample sizes, and ARIMA specifications. These simulation-based efficiency ratios (e.g., 1.58–1.90 for N=500) are generated externally to the estimator definition and do not reduce to fitted parameters or self-referential predictions by construction. No load-bearing step equates the claimed outperformance to the method's own inputs, self-citations, or ansatzes; the derivation chain remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The claim rests on standard ARIMA modeling assumptions and the existence of higher moments for semiparametric estimation; no free parameters or new entities are introduced in the abstract.

axioms (2)
  • domain assumption The observed time series follows a correctly specified ARIMA(p,d,q) process that is stationary and invertible after differencing.
    This is the foundational assumption required for any ARIMA parameter estimation procedure.
  • domain assumption Innovations have finite moments of order at least four to permit computation and use of cumulants in the polynomial maximization.
    Necessary for the semiparametric PMM2 to operate without a full distributional specification.

pith-pipeline@v0.9.0 · 5557 in / 1585 out tokens · 40155 ms · 2026-05-17T23:55:06.086674+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Polynomial Maximization Method with Fractional Polynomial Basis: A Frequentist Bridge to Bayesian Fractional Polynomials

    stat.ME 2026-05 unverdicted novelty 6.0 partial

    PMM-FP extends polynomial maximization to fractional polynomial bases and derives a closed-form variance-reduction coefficient g2 for asymmetric non-Gaussian errors, formalized in Lean 4 and checked via Monte Carlo.

  2. Variance-Reduced Manifold Sampling via Polynomial-Maximization Density Estimation

    stat.ME 2026-05 unverdicted novelty 4.0

    PMM-MASEM introduces a gated PMM2/PMM3 density estimator on kNN shell spacings for MASEM, reducing MSE by 22-36% on asymmetric regimes while falling back to MLE on flat Exp(1) spacings and showing mixed results overall.

Reference graph

Works this paper leans on

46 extracted references · 46 canonical work pages · cited by 2 Pith papers

  1. [1]

    George E. P. Box, Gwilym M. Jenkins, Gregory C. Reinsel, and Greta M. Ljung.Time Series Analysis: Forecasting and Control. John Wiley & Sons, Hoboken, NJ, 5th edition, 2015

  2. [2]

    OTexts, 3rd edition, 2021

    RobJ.HyndmanandGeorgeAthanasopoulos.Forecasting: Principles and Practice. OTexts, 3rd edition, 2021. Accessed: 2025-01-15

  3. [3]

    Kantelhardt, Stephan A

    Jan W. Kantelhardt, Stephan A. Zschiegner, Eva Koscielny-Bunde, Shlomo Havlin, Armin Bunde, and H. Eugene Stanley. Multifractal detrended fluctuation analysis of nonstationary time series.Physica A: Statistical Mechanics and Its Applications, 316(1-4):87–114, 2002

  4. [4]

    Young Shin Kim and Frank J. Fabozzi. Approximation of skewed and leptokurtic return distributions.Applied Financial Economics, 22(16):1299–1316, 2012

  5. [5]

    Fat tails in financial return distributions revisited: Evidence from the korean stock market.Physica A: Statistical Mechanics and its Applications, 526:121114, 2019

    Cheoljun Eom, Taisei Kaizoji, and Enrico Scalas. Fat tails in financial return distributions revisited: Evidence from the korean stock market.Physica A: Statistical Mechanics and its Applications, 526:121114, 2019

  6. [6]

    Commodity prices and international inflation, 1851– 1913.Journal of International Money and Finance, 144:103084, 2024

    Stefan Gerlach and Rebecca Stuart. Commodity prices and international inflation, 1851– 1913.Journal of International Money and Finance, 144:103084, 2024. Also available as CEPR Discussion Paper DP16526

  7. [7]

    On the optimal prediction of extreme events in heavy-tailed time series with applications to solar flare forecasting.Journal of Time Series Analysis, 2025

    Victor Verma, Stilian Stoev, and Yang Chen. On the optimal prediction of extreme events in heavy-tailed time series with applications to solar flare forecasting.Journal of Time Series Analysis, 2025

  8. [8]

    Igoris Belovas, Leonidas Sakalauskas, Vadimas Starikovičius, and Edward W. Sun. Mixed- stable models: An application to high-frequency financial data.Entropy, 23(6):739, 2021

  9. [9]

    Modeling and simulation of financial returns under non-gaussian distributions.Physica A: Statistical Mechanics and Its Applications, 622:128886, 2023

    Federica De Domenico, Giacomo Livan, Guido Montagna, and Oreste Nicrosini. Modeling and simulation of financial returns under non-gaussian distributions.Physica A: Statistical Mechanics and Its Applications, 622:128886, 2023

  10. [10]

    Pötscher

    Benedikt M. Pötscher. Noninvertibility and pseudo-maximum likelihood estimation of mis- specified arma models.Econometric Theory, 7(4):435–449, 1991

  11. [11]

    Quasi-maximum likelihood estimation of garch models with heavy-tailed likelihoods.Journal of Business & Economic Statistics, 32(2):178– 191, 2014

    Jianqing Fan, Li Qi, and Dacheng Xiu. Quasi-maximum likelihood estimation of garch models with heavy-tailed likelihoods.Journal of Business & Economic Statistics, 32(2):178– 191, 2014. 31

  12. [12]

    Maximum likelihood estimation for nearly non- stationary stable autoregressive processes.Journal of Time Series Analysis, 33(3):424–439, 2012

    Rongmao Zhang and Chor-Yiu Sin. Maximum likelihood estimation for nearly non- stationary stable autoregressive processes.Journal of Time Series Analysis, 33(3):424–439, 2012

  13. [13]

    On the forecasting of high-frequency finan- cial time series based on arima model improved by deep learning.Journal of Forecasting, 39(7):1081–1097, 2020

    Long Li, Siyu Leng, Jun Yang, and Guang Yu. On the forecasting of high-frequency finan- cial time series based on arima model improved by deep learning.Journal of Forecasting, 39(7):1081–1097, 2020

  14. [14]

    Dowe, Shelton Peiris, and Ellie Kim

    David L. Dowe, Shelton Peiris, and Ellie Kim. A novel arfima-ann hybrid model for fore- casting time series—and its role in explainable ai.Journal of Econometrics and Statistics, 5(1):107–127, 2025

  15. [15]

    Inference robustness of arima models under non-normality.Metrika, 26:43–56, 1989

    Johannes Ledolter. Inference robustness of arima models under non-normality.Metrika, 26:43–56, 1989

  16. [16]

    Peter J. Huber. Robust estimation of a location parameter.The Annals of Mathematical Statistics, 35(1):73–101, 1964

  17. [17]

    Nora Muler, Daniel Peña, and Victor J. Yohai. Robust estimation for arma models.The Annals of Statistics, 37(2):816–840, 2009

  18. [18]

    A robust m-estimator for gaussian arma time series based on the whittle approximation.Applied Mathematical Modelling, 134:666–678, 2024

    Valdério Anselmo Reisen, Céline Lévy-Leduc, and Carlo Corrêa Solci. A robust m-estimator for gaussian arma time series based on the whittle approximation.Applied Mathematical Modelling, 134:666–678, 2024

  19. [19]

    Quantile time series regression models revisited

    Christis Katsouris. Quantile time series regression models revisited. arXiv preprint arXiv:2308.06617, 2023

  20. [20]

    Least absolute deviations estimation for arch and garch models

    Liang Peng and Qiwei Yao. Least absolute deviations estimation for arch and garch models. Biometrika, 90(4):967–975, 2003

  21. [21]

    Self-weighted least absolute deviation estimation for infinite variance autore- gressive models.Journal of the Royal Statistical Society: Series B, 67(3):381–393, 2005

    Shiqing Ling. Self-weighted least absolute deviation estimation for infinite variance autore- gressive models.Journal of the Royal Statistical Society: Series B, 67(3):381–393, 2005

  22. [22]

    Model-based pricing for financial derivatives with negative skew- ness and excess kurtosis.Econometric Theory, 31(6):1199–1242, 2015

    Ke Zhu and Shiqing Ling. Model-based pricing for financial derivatives with negative skew- ness and excess kurtosis.Econometric Theory, 31(6):1199–1242, 2015

  23. [23]

    A student-t mixture autoregressive model with applications to heavy-tailed financial data.Biometrika, 96(3):751–760, 2009

    Chun Shan Wong, Wai Sum Chan, and Pui Lam Kam. A student-t mixture autoregressive model with applications to heavy-tailed financial data.Biometrika, 96(3):751–760, 2009. 32

  24. [24]

    Mensah, Eric N

    Samuel Ampadu, Eric T. Mensah, Eric N. Aidoo, Alexander Boateng, and Daniel Maposa. A comparativestudyoferrordistributionsinthegarchmodelthroughamontecarlosimulation approach.Scientific African, 23:e01988, 2024

  25. [25]

    Gramacy, Christian L

    Timothy Graves, Robert B. Gramacy, Christian L. E. Franzke, and Nicholas W. Watkins. Efficient bayesian inference for natural time series using arfima processes.Nonlinear Pro- cesses in Geophysics, 22:679–700, 2015

  26. [26]

    A comprehensive survey of deep learning for time series forecasting: Architectural diversity and open challenges

    Tian Zhou, Ziqing Niu, Xue Wang, Liang Sun, and Rong Jin. A comprehensive survey of deep learning for time series forecasting: Architectural diversity and open challenges. Artificial Intelligence Review, 58(3):1–45, 2025. Comprehensive review of Transformer ar- chitectures and attention mechanisms for time series

  27. [27]

    Arima and attention-based cnn-lstm hybrid neural network for battery life estimation

    Xiaoming Liu, Wei Chen, and Yue Zhang. Arima and attention-based cnn-lstm hybrid neural network for battery life estimation. InProceedings of the 2024 8th International Conference on Machine Learning and Soft Computing, pages 45–51. ACM, 2024. Hybrid model combining ARIMA preprocessing with attention-based neural networks

  28. [28]

    Advances in deep learning for time series forecasting/classification winter 2024.Medium: Deep Data Science, 2024

    Isaac Godfried. Advances in deep learning for time series forecasting/classification winter 2024.Medium: Deep Data Science, 2024. Review of foundation models including TimeGPT- 1, Chronos, and TimesFM for time series forecasting

  29. [29]

    Kunchenko.Polynomial Parameter Estimations of Close to Gaussian Random Variables

    Yuriy P. Kunchenko.Polynomial Parameter Estimations of Close to Gaussian Random Variables. Shaker Verlag, Aachen, Germany, 2002

  30. [30]

    Kunchenko.Stokhastychni Polinomy

    Yuriy P. Kunchenko.Stokhastychni Polinomy. Naukova Dumka, Kyiv, 2006. In Ukrainian

  31. [31]

    Zabolotnii, Zygmunt L

    Serhii W. Zabolotnii, Zygmunt L. Warsza, and Oleksandr Tkachenko. Polynomial estimation of linear regression parameters for the asymmetric pdf of errors. In Roman Szewczyk, Cezary Zieliński, and Małgorzata Kaliczyńska, editors,Automation 2018: Advances in Intelligent Systems and Computing, volume 743, pages 758–772, Cham, 2018. Springer

  32. [32]

    Serhii Zabolotnii, Volodymyr Khotunov, Andrii Chepynoha, and Oleksandr Tkachenko. Es- timating parameters of linear regression with an exponential power distribution of errors by using a polynomial maximization method.Eastern-European Journal of Enterprise Tech- nologies, 1(4-109):64–73, 2021

  33. [33]

    Joint signal parameter estimation in non-gaussian noise by the method of polynomial maximization.Journal of Electrical Engineering, 67(3):217– 221, 2016

    Volodymyr Palahin and Jozef Juhár. Joint signal parameter estimation in non-gaussian noise by the method of polynomial maximization.Journal of Electrical Engineering, 67(3):217– 221, 2016. 33

  34. [34]

    Warsza and Serhii W

    Zygmunt L. Warsza and Serhii W. Zabolotnii. A polynomial estimation of measurand pa- rameters for samples of non-gaussian symmetrically distributed data. In Roman Szewczyk, Cezary Zieliński, and Małgorzata Kaliczyńska, editors,Automation 2017: Advances in In- telligent Systems and Computing, volume 550, pages 468–480, Cham, 2017. Springer

  35. [35]

    Zabolotnii, Zygmunt L

    Serhii W. Zabolotnii, Zygmunt L. Warsza, and Oleksandr Tkachenko. Estimation of lin- ear regression parameters of symmetric non-gaussian errors by polynomial maximization method. In Roman Szewczyk, Jurek Sąsiadek, and Małgorzata Kaliczyńska, editors,Au- tomation 2019: Advances in Intelligent Systems and Computing, volume 920, pages 636–649, Cham, 2020. Springer

  36. [36]

    Large sample properties of generalized method of moments estimators

    Lars Peter Hansen. Large sample properties of generalized method of moments estimators. Econometrica, 50(4):1029–1054, 1982

  37. [37]

    Econometric issues in the analysis of regressions with generated regressors

    Adrian Pagan. Econometric issues in the analysis of regressions with generated regressors. International Economic Review, 25(1):221–247, 1984

  38. [38]

    Whitney K. Newey. A method of moments interpretation of sequential estimators.Eco- nomics Letters, 14(2-3):201–206, 1984

  39. [39]

    Bassett, Gilbert

    Roger Koenker and Jr. Bassett, Gilbert. Regression quantiles.Econometrica, 46(1):33–50, 1978

  40. [40]

    Hampel, Elvezio M

    Frank R. Hampel, Elvezio M. Ronchetti, Peter J. Rousseeuw, and Werner A. Stahel.Robust Statistics: The Approach Based on Influence Functions. John Wiley & Sons, New York, 1986

  41. [41]

    Harvey.Dynamic Models for Volatility and Heavy Tails: With Applications to Financial and Economic Time Series

    Andrew C. Harvey.Dynamic Models for Volatility and Heavy Tails: With Applications to Financial and Economic Time Series. Cambridge University Press, Cambridge, 2013

  42. [42]

    Model-based clustering of multiple time series.Journal of Business & Economic Statistics, 24(1):78–89, 2006

    Sylvia Frühwirth-Schnatter and Sylvia Kaufmann. Model-based clustering of multiple time series.Journal of Business & Economic Statistics, 24(1):78–89, 2006

  43. [43]

    Stochastic volatility model with leverage and asym- metricallyheavy-tailederrorusingghskewstudent’st-distribution.Computational Statistics & Data Analysis, 56(11):3690–3704, 2012

    Jouchi Nakajima and Yasuhiro Omori. Stochastic volatility model with leverage and asym- metricallyheavy-tailederrorusingghskewstudent’st-distribution.Computational Statistics & Data Analysis, 56(11):3690–3704, 2012

  44. [44]

    Cambridge University Press, Cambridge, 2005

    Roger Koenker.Quantile Regression. Cambridge University Press, Cambridge, 2005

  45. [45]

    Forecasting the real price of oil in a changing world: A forecast combination approach.Journal of Business & Economic Statistics, 34(3):338– 351, 2016

    Christiane Baumeister and Lutz Kilian. Forecasting the real price of oil in a changing world: A forecast combination approach.Journal of Business & Economic Statistics, 34(3):338– 351, 2016. 34

  46. [46]

    Rivera-Castro

    Shawkat Hammoudeh, Juan Carlos Reboredo, and Miguel A. Rivera-Castro. Nonlinear co-movement in oil and stock markets.Energy Economics, 49:84–94, 2015. 35