Smoothed empirical likelihood for GARCH models with heavy-tailed errors

ABSTRACT

This paper proposes an empirical likelihood (EL) method for estimating the GARCH(p, q) models with heavy-tailed errors. Using the kernel smoothing method, we derive a smoothed EL ratio statistic, which yields a smoothed EL estimator. Moreover, we derive a profile EL for the partial parameters in the presence of nuisance parameters. Simulations and empirical results are conducted to illustrate our proposed method.     

Adaptive LASSO estimation for ARDL models with GARCH innovations

ABSTRACT

In this paper, we show the validity of the adaptive least absolute shrinkage and selection operator (LASSO) procedure in estimating stationary autoregressive distributed lag(p,q) models with innovations in a broad class of conditionally heteroskedastic models. We show that the adaptive LASSO selects the relevant variables with probability converging to one and that the estimator is oracle efficient, meaning that its distribution converges to the same distribution of the oracle-assisted least squares, i.e., the least square estimator calculated as if we knew the set of relevant variables beforehand. Finally, we show that the LASSO estimator can be used to construct the initial weights. The performance of the method in finite samples is illustrated using Monte Carlo simulation.     

Robust estimates for GARCH models

In this paper we present two robust estimates for GARCH models. The first is defined by the minimization of a conveniently modified likelihood and the second is similarly defined, but includes an additional mechanism for restricting the propagation of the effect of one outlier on the next estimated conditional variances. We study the asymptotic properties of our estimates proving consistency and asymptotic normality. A Monte Carlo study shows that the proposed estimates compare favorably with respect to other robust estimates. Moreover, we consider some real examples with financial data that illustrate the behavior of these estimates.     

Robust and efficient estimation of GARCH models based on Hellinger distance

It is well known that financial data frequently contain outlying observations. Almost all methods and techniques used to estimate GARCH models are likelihood-based and thus generally non-robust against outliers. Minimum distance method, as an important tool for statistical inferences and a competitive alternative for achieving robustness, has surprisingly not been well explored for GARCH models. In this paper, we proposed a minimum Hellinger distance estimator (MHDE) and a minimum profile Hellinger distance estimator (MPHDE), depending on whether the innovation distribution is specified or not, for estimating the parameters in GARCH models. The construction and investigation of the two estimators are quite involved due to the non-i.i.d. nature of data. We proved that the MHDE is a consistent estimator and derived its bias in explicit expression. For both of the proposed estimators, we demonstrated their finite-sample performance through simulation studies and compared with the well-established methods including MLE, Gaussian Quasi-MLE, Non-Gaussian Quasi-MLE and Least Absolute Deviation estimator. Our numerical results showed that MHDE and MPHDE have much better performance than MLE-based methods when data are contaminated while simultaneously they are very competitive when data is clean, which testified to the robustness and efficiency of the two proposed MHD-type estimations.     

Spectral analysis of Markov switching GARCH models with statistical inference

We derive matrix expressions in closed form for the autocovariance function and the spectral density of Markov switching GARCH models and their powers. For this, we apply the Riesz–Fischer theorem which defines the spectral representation as the Fourier transform of the autocovariance function. Under suitable assumptions, we prove that the sample estimator of the spectral density is consistent and asymptotically normally distributed. Further statistical implications in terms of order identification and parameter estimation are discussed. A simulation study confirms the validity of the asymptotic properties. These methods are also well suited for financial market applications, and in particular for the analysis of time series in the frequency domain, as shown in some proposed real-world examples.     

The use of GARCH models in VaR estimation

We evaluate the performance of an extensive family of ARCH models in modeling the daily Value-at-Risk (VaR) of perfectly diversified portfolios in five stock indices, using a number of distributional assumptions and sample sizes. We find, first, that leptokurtic distributions are able to produce better one-step-ahead VaR forecasts; second, the choice of sample size is important for the accuracy of the forecast, whereas the specification of the conditional mean is indifferent. Finally, the ARCH structure producing the most accurate forecasts is different for every portfolio and specific to each equity index.     

Empirical characteristic function tests for GARCH innovation distribution using multipliers

Goodness-of-fit tests for the innovation distribution in GARCH models based on measuring deviations between the empirical characteristic function of the residuals and the characteristic function under the null hypothesis have been proposed in the literature. The asymptotic distributions of these test statistics depend on unknown quantities, so their null distributions are usually estimated through parametric bootstrap (PB). Although easy to implement, the PB can become very computationally expensive for large sample sizes, which is typically the case in applications of these models. This work proposes to approximate the null distribution through a weighted bootstrap. The procedure is studied both theoretically and numerically. Its asymptotic properties are similar to those of the PB, but, from a computational point of view, it is more efficient.     

Self-weighted recursive estimation of GARCH models

The generalized autoregressive conditional heteroscedasticity (GARCH) processes are frequently used to investigate and model financial returns. They are routinely estimated by computationally complex off-line estimation methods, for example, by the conditional maximum likelihood procedure. However, in many empirical applications (especially in the context of high-frequency financial data), it seems necessary to apply numerically more effective techniques to calibrate and monitor such models. The aims of this contribution are: (i) to review the previously introduced recursive estimation algorithms and to derive self-weighted alternatives applying general recursive identification instruments, and (ii) to examine these methods by means of simulations and an empirical application.     

Time series models with asymmetric Laplace innovations

We propose autoregressive moving average (ARMA) and generalized autoregressive conditional heteroscedastic (GARCH) models driven by asymmetric Laplace (AL) noise. The AL distribution plays, in the geometric-stable class, the analogous role played by the normal in the alpha-stable class, and has shown promise in the modelling of certain types of financial and engineering data. In the case of an ARMA model we derive the marginal distribution of the process, as well as its bivariate distribution when separated by a finite number of lags. The calculation of exact confidence bands for minimum mean-squared error linear predictors is shown to be straightforward. Conditional maximum likelihood-based inference is advocated, and corresponding asymptotic results are discussed. The models are particularly suited for processes that are skewed, peaked, and leptokurtic, but which appear to have some higher order moments. A case study of a fund of real estate returns reveals that AL noise models tend to deliver a superior fit with substantially less parameters than normal noise counterparts, and provide both a competitive fit and a greater degree of numerical stability with respect to other skewed distributions.     

The finite-time ruin probability of a discrete-time risk model with GARCH discounted factors and dependent risks

The finite-time ruin probability of a discrete-time risk model with dependent stochastic discount factors and dependent insurance and financial risks is investigated in this paper. Assume that the stochastic discount factors follow a GARCH process and the one-period insurance and financial risks form a sequence of independent and identically distributed random pairs, which are the copies of a random pair with a bivariate Sarmanov dependent distribution. When the common distribution of claim-sizes is heavy-tailed, we establish an asymptotic estimate for the finite-time ruin probability. Applying the result to a special case, we also get conservative asymptotic bounds. A numerical simulation is given at the end of the paper.     

Zero-inflated compound Poisson distributions in integer-valued GARCH models

In this paper we introduce a wide class of integer-valued stochastic processes that allows to take into consideration, simultaneously, relevant characteristics observed in count data namely zero inflation, overdispersion and conditional heteroscedasticity. This class includes, in particular, the compound Poisson, the zero-inflated Poisson and the zero-inflated negative binomial INGARCH models, recently proposed in literature. The main probabilistic analysis of this class of processes is here developed. Precisely, first- and second-order stationarity conditions are derived, the autocorrelation function is deduced and the strict stationarity is established in a large subclass. We also analyse in a particular model the existence of higher-order moments and deduce the explicit form for the first four cumulants, as well as its skewness and kurtosis.     

Analysis of autoregressive models with symmetric stable innovations

This paper develops algorithms for fitting autoregressive models with symmetric stable innovations using auto-covariation function. A recursive algorithm is proposed for generalized Yule-Walker estimation of autoregressive coefficients and partial auto-covariation function. It also introduces a new information criterion, useful for consistent order selection. Applications of the proposed methods are illustrated using observations simulated from autoregressive models with symmetric stable innovations as well as by analysing a set of real data.     

Data cloning estimation of GARCH and COGARCH models

GARCH models include most of the stylized facts of financial time series and they have been largely used to analyse discrete financial time series. In the last years, continuous-time models based on discrete GARCH models have been also proposed to deal with non-equally spaced observations, as COGARCH model based on Lévy processes. In this paper, we propose to use the data cloning methodology in order to obtain estimators of GARCH and COGARCH model parameters. Data cloning methodology uses a Bayesian approach to obtain approximate maximum likelihood estimators avoiding numerically maximization of the pseudo-likelihood function. After a simulation study for both GARCH and COGARCH models using data cloning, we apply this technique to model the behaviour of some NASDAQ time series.     

The Minimum Density Power Divergence Estimation for the Lognormal Density

In this article, we implement the minimum density power divergence estimation for estimating the parameters of the lognormal density. We compare the minimum density power divergence estimator (MDPDE) and the maximum likelihood estimator (MLE) in terms of robustness and asymptotic distribution. The simulations and an example indicate that the MDPDE is less biased than MLE and is as good as MLE in terms of the mean square error under various distributional situations.     

Monitoring distributional changes of squared residuals in GARCH models

Change point monitoring for distributional changes in time-series models is an important issue. In this article, we propose two monitoring procedures to detect distributional changes of squared residuals in GARCH models. The asymptotic properties of our monitoring statistics are derived under both the null of no change in distribution and the alternative of a change in distribution. The finite sample properties are investigated by a simulation.     

On a stochastic process with a heavy-tailed distributed component describing inventory model type of (s,S)

In the present study, the stochastic process X(t) describing inventory model type of (s, S) with a heavy-tailed distributed demands is considered. The asymptotic expansions at sufficiently large values of parameter β = S ? s for the ergodic distribution and nth-order moment of the process X(t) based on the main results of the studies Teugels (1968 Teugels, J.L. (1968). Renewal theorems when the first or the second moment is infinite. Ann. Math. Stat. 39(4):1210–1219.[Crossref] , [Google Scholar]) and Geluk and Frenk (2011 Geluk, J.L., Frenk, J.B.G. (2011). Renewal theory for random variables with a heavy tailed distribution and finite variance. Stat. Probab. Lett. 81:77–82.[Crossref], [Web of Science ®] , [Google Scholar]) are obtained.     

Sensitivity analysis on ruin probabilities with heavy-tailed claims

In this note, we consider the classical insurance risk model with heavy-tailed claim distributions. By using the Pollaczek–Khinchin Formula, we provide some sensitivity analysis on the ruin probability.     

Zero-inflated Poisson and negative binomial integer-valued GARCH models

Zero inflation means that the proportion of 0's of a model is greater than the proportion of 0's of the corresponding Poisson model, which is a common phenomenon in count data. To model the zero-inflated characteristic of time series of counts, we propose zero-inflated Poisson and negative binomial INGARCH models, which are useful and flexible generalizations of the Poisson and negative binomial INGARCH models, respectively. The stationarity conditions and the autocorrelation function are given. Based on the EM algorithm, the estimating procedure is simple and easy to be implemented. A simulation study shows that the estimation method is accurate and reliable as long as the sample size is reasonably large. A real data example leads to superior performance of the proposed models compared with other competitive models in the literature.     

Estimation of the parameters of symmetric stable ARMA and ARMA–GARCH models

In this article, we first propose the modified Hannan–Rissanen Method for estimating the parameters of autoregressive moving average (ARMA) process with symmetric stable noise and symmetric stable generalized autoregressive conditional heteroskedastic (GARCH) noise. Next, we propose the modified empirical characteristic function method for the estimation of GARCH parameters with symmetric stable noise. Further, we show the efficiency, accuracy and simplicity of our methods with Monte-Carlo simulation. Finally, we apply our proposed methods to model the financial data.