Consistent Estimation of the Value at Risk When the Error Distribution of the Volatility Model is Misspecified

A two‐step approach for conditional value at risk estimation is considered. First, a generalized quasi‐maximum likelihood estimator is employed to estimate the volatility parameter, then the empirical quantile of the residuals serves to estimate the theoretical quantile of the innovations. When the instrumental density h of the generalized quasi‐maximum likelihood estimator is not the Gaussian density, both the estimations of the volatility and of the quantile are generally asymptotically biased. However, the two errors counterbalance and lead to a consistent estimator of the value at risk. We obtain the asymptotic behavior of this estimator and show how to choose optimally h.     

Local Likelihood for non‐parametric ARCH(1) models

Abstract. We propose a non‐parametric local likelihood estimator for the log‐transformed autoregressive conditional heteroscedastic (ARCH) (1) model. Our non‐parametric estimator is constructed within the likelihood framework for non‐Gaussian observations: it is different from standard kernel regression smoothing, where the innovations are assumed to be normally distributed. We derive consistency and asymptotic normality for our estimators and show, by a simulation experiment and some real‐data examples, that the local likelihood estimator has better predictive potential than classical local regression. A possible extension of the estimation procedure to more general multiplicative ARCH(p) models with p > 1 predictor variables is also described.     

On the asymptotic properties of a feasible estimator of the continuous time long memory parameter

This article considers a fractional noise model in continuous time and examines the asymptotic properties of a feasible frequency domain maximum likelihood estimator of the long memory parameter. The feasible estimator is one that maximizes an approximation to the likelihood function (the approximation arises from the fact that the spectral density function involves the finite truncation of an infinite summation). It is of interest therefore to explore the conditions required of this approximation to ensure the consistency and asymptotic normality of this estimator.     

Statistical Inference for Unified Garch–Itô Models with High‐Frequency Financial Data

The existing estimation methods for the model parameters of the unified GARCH–Itô model (Kim and Wang, 2014 ) require long period observations to obtain the consistency. However, in practice, it is hard to believe that the structure of a stock price is stable during such a long period. In this article, we introduce an estimation method for the model parameters based on the high‐frequency financial data with a finite observation period. In particular, we establish a quasi‐likelihood function for daily integrated volatilities, and realized volatility estimators are adopted to estimate the integrated volatilities. The model parameters are estimated by maximizing the quasi‐likelihood function. We establish asymptotic theories for the proposed estimator. A simulation study is conducted to check the finite sample performance of the proposed estimator. We apply the proposed estimation approach to the Bank of America stock price data.     

Minimum α‐divergence estimation for arch models

Abstract. This paper considers a minimum α‐divergence estimation for a class of ARCH(p) models. For these models with unknown volatility parameters, the exact form of the innovation density is supposed to be unknown in detail but is thought to be close to members of some parametric family. To approximate such a density, we first construct an estimator for the unknown volatility parameters using the conditional least squares estimator given by Tjøstheim [Stochastic processes and their applications (1986) Vol. 21, pp. 251–273]. Then, a nonparametric kernel density estimator is constructed for the innovation density based on the estimated residuals. Using techniques of the minimum Hellinger distance estimation for stochastic models and residual empirical process from an ARCH(p) model given by Beran [Annals of Statistics (1977) Vol. 5, pp. 445–463] and Lee and Taniguchi [Statistica Sinica (2005) Vol. 15, pp. 215–234] respectively, it is shown that the proposed estimator is consistent and asymptotically normal. Moreover, a robustness measure for the score of the estimator is introduced. The asymptotic efficiency and robustness of the estimator are illustrated by simulations. The proposed estimator is also applied to daily stock returns of Dell Corporation.     

Maximum likelihood estimation of higher‐order integer‐valued autoregressive processes

Abstract. In this article, we extend the earlier work of Freeland and McCabe [Journal of time Series Analysis (2004) Vol. 25, pp. 701–722] and develop a general framework for maximum likelihood (ML) analysis of higher‐order integer‐valued autoregressive processes. Our exposition includes the case where the innovation sequence has a Poisson distribution and the thinning is binomial. A recursive representation of the transition probability of the model is proposed. Based on this transition probability, we derive expressions for the score function and the Fisher information matrix, which form the basis for ML estimation and inference. Similar to the results in Freeland and McCabe (2004) , we show that the score function and the Fisher information matrix can be neatly represented as conditional expectations. Using the INAR(2) specification with binomial thinning and Poisson innovations, we examine both the asymptotic efficiency and finite sample properties of the ML estimator in relation to the widely used conditional least squares (CLS) and Yule–Walker (YW) estimators. We conclude that, if the Poisson assumption can be justified, there are substantial gains to be had from using ML especially when the thinning parameters are large.     

Least tail‐trimmed squares for infinite variance autoregressions

We develop a robust least squares estimator for autoregressions with possibly heavy tailed errors. Robustness to heavy tails is ensured by negligibly trimming the squared error according to extreme values of the error and regressors. Tail‐trimming ensures asymptotic normality and super‐‐convergence with a rate comparable to the highest achieved amongst M‐estimators for stationary data. Moreover, tail‐trimming ensures robustness to heavy tails in both small and large samples. By comparison, existing robust estimators are not as robust in small samples, have a slower rate of convergence when the variance is infinite, or are not asymptotically normal. We present a consistent estimator of the covariance matrix and treat classic inference without knowledge of the rate of convergence. A simulation study demonstrates the sharpness and approximate normality of the estimator, and we apply the estimator to financial returns data. Finally, tail‐trimming can be easily extended beyond least squares estimation for a linear stationary AR model. We discuss extensions to quasi‐maximum likelihood for GARCH, weighted least squares for a possibly non‐stationary random coefficient autoregression, and empirical likelihood for robust confidence region estimation, in each case for models with possibly heavy tailed errors.     

Robust estimation for copula Parameter in SCOMDY models

In this study, we study the robust estimation for the copula parameter in semiparametric copula‐based multivariate dynamic (SCOMDY) models proposed by Chen and Fan (2006). To this end, instead of the pseudo maximum likelihood estimator in Chen and Fan (2006), we use a minimum density power divergence estimator (MDPDE) proposed by Basu et al. (1998). It is shown that the MDPDE is consistent and asymptotically normal under regularity conditions. We compare the performance between the two estimators when outliers exist through a simulation study.     

Further Results on Pseudo-Maximum Likelihood Estimation and Testing in the Constant Elasticity of Variance Continuous Time Model

Constant elasticity volatility processes have been shown to be useful, for example, to encompass a number of existing models that have closed-form likelihood functions. In this article, we extend the existing literature in two directions: first we find explicit closed form solutions of the pseudo maximum likelihood estimators (MLEs) by discretizing the diffusion function and we provide their asymptotic theory in the context of the constant elasticity of variance (CEV) model characterized by a general CEV parameter ρ ≥ 0. Second we obtain bias expansions for those pseudo MLEs also in terms of ρ ≥ 0. We provide a general framework since only the cases with ρ = 0 and ρ = 0.5 have been considered in the literature so far. When the time series is not positive almost surely, we need to impose the restriction that ρ is a non-negative integer.     

Gaussian Maximum Likelihood Estimation For ARMA Models. I. Time Series

Abstract. We provide a direct proof for consistency and asymptotic normality of Gaussian maximum likelihood estimators for causal and invertible autoregressive moving‐average (ARMA) time series models, which were initially established by Hannan [Journal of Applied Probability (1973) vol. 10, pp. 130–145] via the asymptotic properties of a Whittle's estimator. This also paves the way to establish similar results for spatial processes presented in the follow‐up article by Yao and Brockwell [Bernoulli (2006) in press].     

Small‐b and Fixed‐b Asymptotics for Weighted Covariance Estimation in Fractional Cointegration

In a standard cointegrating framework, Phillips (1991) introduced the weighted covariance (WC) estimator of cointegrating parameters. Later, Marinucci (2000) applied this estimator to fractional circumstances and, like Phillips (1991), analysed the so‐called small‐b asymptotic approximation to its sampling distribution. Recently, an alternative limiting theory (fixed‐b asymptotics) has been successfully employed to approximate sampling distributions. With the purpose of comparing both approaches, we derive here the fixed‐b limit of WC estimators in a fractional setting, filling also some gaps in the traditional (small‐b) theory. We also provide some Monte Carlo evidence that suggests that the fixed‐b limit is more accurate.     

Finite Sample Theory of QMLE in ARCH Models with Dynamics in the Mean Equation

Abstract. We provide simulation and theoretical results concerning the finite‐sample theory of quasi‐maximum‐likelihood estimators in autoregressive conditional heteroskedastic (ARCH) models when we include dynamics in the mean equation. In the setting of the AR(q)–ARCH(p), we find that in some cases bias correction is necessary even for sample sizes of 100, especially when the ARCH order increases. We warn about the existence of important biases and potentially low power of the t‐tests in these cases. We also propose ways to deal with them. We also find simulation evidence that when conditional heteroskedasticity increases, the mean‐squared error of the maximum‐likelihood estimator of the AR(1) parameter in the mean equation of an AR(1)‐ARCH(1) model is reduced. Finally, we generalize the Lumsdaine [J. Bus. Econ. Stat. 13 (1995) pp. 1–10] invariance properties for the biases in these situations.     

ESTIMATION OF THE MULTIVARIATE AUTOREGRESSIVE MOVING AVERAGE HAVING PARAMETER RESTRICTIONS AND AN APPLICATION TO ROTATIONAL SAMPLING

Abstract. The vector autoregressive moving average model with nonlinear parametric restrictions is considered. A simple and easy-to-compute Newton-Raphson estimator is proposed that approximates the restricted maximum likelihood estimator which takes full advantage of the information contained in the restrictions. In the case when there are no parametric restrictions, our Newton-Raphson estimator is equivalent to the estimator proposed by Reinsel et al. (Maximum likelihood estimators in the multivariate autoregressive moving-average model from a generalized least squares view point. J. Time Ser. Anal. 13 (1992), 133–45). The Newton-Raphson estimation procedure also extends to the vector ARMAX model. Application of our Newton-Raphson estimation method in rotational sampling problems is discussed. Simulation results are presented for two different restricted models to illustrate the estimation procedure and compare its performance with that of two alternative procedures that ignore the parametric restrictions.     

Time‐series models with an EGB2 conditional distribution

A time‐series model in which the signal is buried in noise that is non‐Gaussian may throw up observations that, when judged by the Gaussian yardstick, are outliers. We describe an observation‐driven model, based on an exponential generalized beta distribution of the second kind (EGB2), in which the signal is a linear function of past values of the score of the conditional distribution. This specification produces a model that is not only easy to implement but which also facilitates the development of a comprehensive and relatively straightforward theory for the asymptotic distribution of the maximum‐likelihood (ML) estimator. Score‐driven models of this kind can also be based on conditional t distributions, but whereas these models carry out what, in the robustness literature, is called a soft form of trimming, the EGB2 distribution leads to a soft form of Winsorizing. An exponential general autoregressive conditional heteroscedastic (EGARCH) model based on the EGB2 distribution is also developed. This model complements the score‐driven EGARCH model with a conditional t distribution. Finally, dynamic location and scale models are combined and applied to data on the UK rate of inflation.     

Bootstrapping a weighted linear estimator of the ARCH parameters

Abstract.  A standard assumption while deriving the asymptotic distribution of the quasi maximum likelihood estimator in ARCH models is that all ARCH parameters must be strictly positive. This assumption is also crucial in deriving the limit distribution of appropriate linear estimators (LE). We propose a weighted linear estimator (WLE) of the ARCH parameters in the classical ARCH model and show that its limit distribution is multivariate normal even when some of the ARCH coefficients are zero. The asymptotic dispersion matrix involves unknown quantities. We consider appropriate bootstrapped version of this WLE and prove that it is asymptotically valid in the sense that the bootstrapped distribution (given the data) is a consistent estimate (in probability) of the distribution of the WLE. Although we do not show theoretically that the bootstrap outperforms the normal approximation, our simulations demonstrate that it yields better approximations than the limiting normal.     

Robustness of Zero Crossing Estimator

Zero crossing (ZC) statistic is the number of zero crossings observed in a time series. The expected value of the ZC specifies the first‐order autocorrelation of the processes. Hence, we can estimate the autocorrelation by using the ZC estimator. The asymptotic consistency and normality of the ZC estimator for scalar Gaussian processes are already discussed in 1980. In this article, first, we derive the joint asymptotic distribution of the ZC estimator for ellipsoidal processes. Next, we show the variance of the ZC estimator does not attain the Cramer–Rao lower bound (CRLB). However, it is shown that the ZC estimator has robustness when the process is contaminated by an outlier. In contrast with this, we observe that the quasi‐maximum likelihood estimator (QMLE) attains the CRLB. However, we can see that QMLE is sensitive for the outlier.     

NON‐STATIONARITY AND QUASI‐MAXIMUM LIKELIHOOD ESTIMATION ON A DOUBLE AUTOREGRESSIVE MODEL

This article first studies the non‐stationarity of the first‐order double AR model, which is defined by the random recurrence equation , where γ₀ > 0, α₀ ≥ 0, and {η_t}is a sequence of i.i.d. symmetric random variables. It is shown that the double AR(1) model is explosive under the condition . Based on this, it is shown that the quasi‐maximum likelihood estimator of (φ₀,α₀) is consistent and asymptotically normal so that the unit root problem does not exist in the double AR(1) model. Simulation studies are carried out to assess the performance of the quasi‐maximum likelihood estimator in finite samples.     

SOME PROPERTIES OF THE MAXIMUM LIKELIHOOD ESTIMATOR IN THE SIMULTANEOUS SWITCHING AUTOREGRESSIVE MODEL

Abstract. The simultaneous switching autoregressive (SSAR) model proposed by Kunitomo and Sato (A non-linearity in economic time series and disequilibrium econometric models. In Theory and Application of Mathematical Statistics (ed. A. Takemura). Tokyo:University of Tokyo Press (in Japanese), 1994; Asymmetry in economic time series and simultaneous switching autoregressive model. Struct. Change Econ. Dyn. , forthcoming (1994).) is a Markovian non-linear time series model. We investigate the finite sample as well as the asymptotic properties of the least squares estimator and the maximum likelihood (ML) estimator. Due to a specific simultaneity involved in the SSAR model, the least squares estimator is badly biased. However, the ML estimator under the assumption of Gaussian disturbances gives reasonable estimates.     

QMLE for Quadratic ARCH Model with Long Memory

We discuss parametric quasi‐maximum likelihood estimation for quadratic autoregressive conditional heteroskedasticity (ARCH) process with long memory introduced in Doukhan emphet al. (2016) and Grublyt? and ?karnulis (2016) with conditional variance involving the square of inhomogeneous linear combination of observable sequence with square summable weights. The aforementioned model extends the quadratic ARCH model of Sentana ( 1995 ) and the linear ARCH model of Robinson ( 1991 ) to the case of strictly positive conditional variance. We prove consistency and asymptotic normality of the corresponding quasi‐maximum likelihood estimates, including the estimate of long memory parameter 0 < d < 1/2. A simulation study of empirical mean‐squared error is included.     

DIFFERENTIAL GEOMETRY OF ARMA MODELS

Abstract. A general approach for the development of a statistical inference on autoregressive moving-average (ARMA) models is presented based on geometric arguments. ARMA models are characterized as members of the curved exponential family. Geometric properties of ARMA models are computed and used to suggest parameter transformations that satisfy predetermined properties. In particular, the effect on the asymptotic bias of the maximum likelihood estimator of model parameters is illustrated. Hypothesis testing of parameters is discussed through the application of a modified form of the likelihood ratio test statistic.