Robust M-Estimation of Multivariate GARCH Models
|
|
- Hector Dawson
- 5 years ago
- Views:
Transcription
1 Robust M-Estimation of Multivariate GARCH Models Kris Boudt Christophe Croux Faculty of Business and Economics, K.U.Leuven, Belgium Abstract In empirical work on multivariate financial time series, it is common to postulate a Multivariate GARCH model. Due to measurement error or unusual economic events, some observations can be in discordance with the model assumptions. It is desirable that these outliers have little influence on the estimation result. In a Monte Carlo study, we show that the Gaussian quasi-maximum likelihood estimator can be highly affected by outliers. As a more robust alternative, we propose to use M-estimators. By downweighting extreme returns in the loss function of the M-estimator and in the MGARCH equation, we obtain robust parameter estimates and conditional covariance matrix predictions. We prove consistency of a wide class of M-estimators for MGARCH models with elliptical innovations. Simulations and a real data example document the benefits of the robust approach. JEL classification: C13; C32; C51 Keywords: GARCH models, M-estimators, multivariate time series, outliers, robust methods Correspondence to: Kris Boudt, Faculty of Business and Economics, 69 Naamsestraat, B-3000 Leuven, Belgium. Kris.Boudt@econ.kuleuven.be Tel: Fax:
2 1 Introduction It is widely recognized that the volatility and the correlations of daily and weekly financial returns may be time-varying. There is no doubt that incorporating these features into the estimation of the conditional covariance matrix of financial returns can lead to better decisions on portfolio optimization, asset pricing and risk management. One way to do this is to specify the conditional covariance matrix as a deterministic function of the innovations in the mean equation. This approach defines the class of Multivariate GARCH (MGARCH) models (see Bauwens et al., 2006, for a recent survey). We propose a new estimator and a new MGARCH volatility forecasting model that, compared to the classical MGARCH methods, have the advantage of being robust to outliers in the data. It is common to estimate MGARCH models by maximizing a Gaussian likelihood function. Jeantheau (1998) shows that this estimation procedure, called the Gaussian Quasi-Maximum Likelihood (QML) estimation, can produce consistent parameter estimates, even when the true distribution is not Gaussian. This is an important result, since it is common to find that after correcting the returns for the dynamics in the conditional covariance matrix, the marginal distribution of the standardized return series is still heavy tailed. Engle and González-Rivera (1991) show that, for the univariate GARCH model, 2
3 the QML estimator has a degree of inefficiency that increases with the degree of departure from normality. A current practice, therefore, consists in doing ML estimation under the assumption of a multivariate Student t distribution (Fiorentini et al., 2003, and Pesaran and Pesaran, 2007). According to Franses and Ghijsels (1999) the leptokurtosis in the standardized returns can also be interpreted as an indication of outliers in the data set. The effect of outliers on estimating univariate GARCH models has been studied in the literature and robust estimators have been proposed. One approach to estimating GARCH models in the presence of outliers is to prune iteratively the outliers and fit the model to the remaining data until no more outliers are detected (Franses and Ghijsels, 1999). Another approach is to use an estimator that is robust to outliers. Sakata and White (1998), Park (2002) and Peng and Yao (2003) estimate the GARCH model by minimizing a robust measure of scale of the residuals. Mancini et al. (2005) and Muler and Yohai (2002,2007) propose a robust M-estimator that assigns a much lower weight to outliers than the Gaussian ML estimator does. Newey and Steigerwald (1997) prove consistency of the Student t Quasi-Maximum Likelihood (QML) estimator, provided the true innovation density is unimodal and symmetric around zero. All authors find that in the presence of outliers, the Gaussian QML estimator collapses and that robust estimators are more reliable. Van Dijk et 3
4 al. (1999) find that LM-tests for conditional heteroscedasticity perform better when they are based on robust parameter estimates rather than on the Gaussian QML estimates. Muler and Yohai (2002,2007) argue that in the presence of outliers, it may be more prudent to specify the conditional variance as a bounded function of past returns. Doing so, one limits the impact of outliers on subsequent conditional variance predictions which depend in an autoregressive way on past returns. To our best knowledge, we are the first to study the effect of outliers on the Gaussian QML estimator of multivariate GARCH models and to propose robust alternatives. In Section 2 we prove consistency of M-estimators of MGARCH models with elliptical innovations. In Section 3 we argue that the Gaussian QML estimator can be made more robust by replacing its loss function with a loss function whose derivative increases more rapidly, and by bounding the impact of news on the conditional covariance matrix process. Section 4 proposes new diagnostic checks for outlier identification and for correct model specification. Section 5 illustrates the new methodology through a Monte Carlo study and a financial application. We show that our proposal of a robust M-estimator has a reasonably good efficiency at the assumed model and that, in contrast with the Gaussian QML estimator, it remains reliable when the data is outlier contaminated. In the empirical case of the daily returns on the Nasdaq and NYSE indices, we find that our proposal for a robust 4
5 MGARCH volatility prediction model reflects better reality than the classical model, especially regarding the difference in reaction of volatility to the 1987 stock market crash and the burst of the internet bubble in Section 6 summarizes our conclusions and outlines the implications for further research. 2 M-estimation 2.1 Definition We consider a centered, conditionally elliptically distributed N-variate random process {r t } whose conditional covariance matrix process {H t,θ } is an MGARCH process. Under this model, each random variable r t can be expressed as r t = H 1/2 t,θ u t, (1) where u t is a zero mean random variable with an elliptically symmetric distribution. Let I rt 1 be the information incorporated in the past of r t. It follows that under this model P θ, the probability of observing r t conditional on its past is described by the elliptic density function p(r t I rt 1 ; θ) = (det H t,θ ) 1/2 g ( r t H 1 t,θ r t), (2) where the function g : R + R +. The scatter matrix H t,θ and the function g( ) are assumed to be standardized such that H t,θ is the conditional covariance 5
6 matrix of r t. While there exist elliptic distributions without finite second moment, we consider here only distributions for which the covariance matrix exists. A common choice of g( ) is the standard normal density function which we denote by φ( ). Another popular choice for g( ) is the standardized Student t density function with ν > 2 degrees of freedom: t ν (z) = Γ ( ) ν+n 2 Γ ( ) N ν 2 [π(ν 2)] 2 [ 1 + z ] N+ν 2, (3) ν 2 where Γ( ) is the gamma function. A useful result is that under (2), the squared Mahalanobis distance is i.i.d. distributed with density function d 2 t,θ = r t H 1 t,θ r t = u t u t (4) f(z) = 2πN/2 Γ(N/2) z N 1 2 g(z). (5) In the special case where g( ) is the Gaussian density function, f( ) is the chisquared density function with N degrees of freedom. We denote F 1 ( ) the quantile function associated to f( ). Under P θ, the conditional covariance matrix H t,θ is parameterized by a MGARCH model as a measurable function of the past innovation terms u t. More precisely, we assume that H t,θ = H θ (u t 1, u t 2,...), (6) with H θ ( ) a N N positive definite, symmetric matrix function of a sequence 6
7 of N-dimensional vectors. As an illustration, consider the BEKK parameterization for H t,θ (Engle and Kroner, 1995): H t,θ = C C + A H 1/2 t 1,θ u t 1u t 1 H1/2 t 1,θ A + B H t 1,θ B, (7) where A, B and C are N N parameter matrices and C is upper triangular. Denote θ the true, unknown parameter vector belonging to the parameter space Θ and g ( ) the true density function in (2). It is common to estimate θ by the value of θ Θ that maximizes the joint density of the sample S T = {r 1,..., r T } constructed under the nominal assumption that the density function of the data is g( ). This estimation approach yields the QML estimator ˆθ QML = argmax θ Θ ave t [ log det Ht,θ 2 log g ( r t H 1 t,θ r t)], (8) where ave t stands for the arithmetic average over t = 1,..., T. It coincides with the ML estimator when g( ) = g ( ). The ML estimator belongs to the broader class of M-estimators defined as follows. Definition 1 The M-estimate based on a sample S T = {r 1,..., r T } is the value of θ Θ for which the M-function M (S T ; θ, ρ) = ave t [ log det Ht,θ + σ ρ ( r t H 1 t,θ r t)], (9) is minimized, with ρ( ) a positive, non-decreasing scalar function. The scalar σ is a consistency factor. 7
8 The function ρ( ) in the above definition is called the loss function associated to the M-estimator. If ρ( ) = 2 log g( ), with g( ) a specified density function, we obtain a QML estimator. The Gaussian φ and standardized Student t ν density functions are the most popular elliptical density functions used to describe financial return series. The corresponding Gaussian and Student t ν (ν > 2) loss functions are given by [ ρ φ (z) = z ; ρ tν (z) = (N + ν) log 1 + z ]. (10) ν 2 The consistency factor σ depends on the true density function g ( ) and on ψ( ) = ρ ( ), the first derivative of the function ρ( ). Denote E Pθ the expectation under P θ. By the consistency proof in Appendix A and by formula (5), we have that σ equals { 2π σ = N/E Pθ ψ(u t u t)u N/2 t u t = N Γ(N/2) 0 ψ(t 2 )t N+1 g (t 2 )dt} 1. (11) For the M-estimator with Gaussian loss function (being the Gaussian QML estimator) the consistency factor is one for all elliptical density functions g, whatever be the dimension of the series (N). For all loss functions, we have that the consistency factor is one when the loss function is proportional to log g ( ). In this case, the M-estimator coincides with the ML estimator. In Figure 1 we illustrate for two different loss functions ρ tν with ν = 20 and ν = 4, the dependence of the consistency factor on the number of degrees of freedom ν for bivariate Student t ν distributed innovations. Note that 8
9 Figure 1: Sensitivity of the consistency factor σ for two Student t ν loss functions (ν = 20 and ν = 10) to the number of degrees of freedom ν for bivariate Student t ν distributed innovations. σ * ρ t20 ρ t ν * the consistency factor σ is most sensitive to ν when ν is small. When the econometrician believes the innovations in the MGARCH model are Student t distributed but has no idea on the value of ν, then, as a practical guideline, we recommend to calculate the consistency factor for an intermediate value of ν such as ν =. The Monte Carlo study of Subsection 5.3 documents the good estimation properties of this approach. If the return process is conditionally normal, the M-estimator with Gaussian loss function is efficient. The Monte Carlo experiment in Section 5 shows that, in the presence of outliers, the M-estimator with Gaussian loss function 9
10 is no longer reliable and that, the smaller ν is in the loss function ρ tν ( ), the more robust the M-estimator with Student t ν loss function will be. First, we give a general description of the asymptotic properties of M-estimators. 2.2 Asymptotic properties In Appendix A we prove consistency for a wide range of M-estimators, defined above as the value of θ Θ that minimizes the time series average of m (r t ; θ, ρ) = log det H t,θ + σ ρ ( r t H 1 t,θ r t). The consistency result is obtained under the following assumptions on the loss function ρ( ), on the functional form H θ ( ) of the conditional covariance matrix and on the density function g ( ) of the innovations. In parenthesis, we indicate what the assumptions are needed for. A1 (Existence and uniqueness) The loss function ρ( ) has the following properties. 1. Its derivative ψ( ) = ρ ( ) is nonnegative, nonincreasing and continuous. 2. The function ψ(z)z is bounded. Let K = sup z 0 ψ(z)z. 3. The function ψ(z)z is nondecreasing, and is strictly increasing in the interval where ψ(z)z < K. 4. There exists z 0 such that ψ(z0 2)z2 0 > N and that ψ(z)z > 0 for z z There exists a > 0 such that for every hyperplane H, the probability under P θ for an observation to lie in H is at most 1 N/K a. 10
11 A2 (Consistency) p(r t I rt 1 ; θ ) = (det H t,θ ) 1/2 g ( r t H 1 t,θ r t ). A3 (Identifiability) H t,θ = H t,θ θ = θ. A4 (Uniform convergence) The parameter space Θ is compact with E Pθ sup θ Θ m (r t ; θ, ρ) / θ <. A5 (Uniqueness) The function H θ ( ) is such that, with probability one, the gradient of the M-function can only be zero if H t,θ = E Pθ [ σ ψ(r t H 1 t,θ r t)r t r t I rt 1 ]. Proposition 1 Under A1-A5 the M-estimator in Definition 1 is a consistent estimator for θ. The key idea behind this result is that, under these assumptions, the conditional covariance matrix evaluated at the M-estimate for θ will, asymptotically, be equivalent to Maronna (1976) s M-estimator for the scale matrix, conditional on the past realizations. Consistency of the M-estimator for MGARCH models follows from the consistency of Maronna (1976) s M-estimator for the unconditional scale matrix. Assumption A1 and A2 are needed to ensure consistency of the unconditional M-estimator of scale. Note that the Gaussian and Student t loss functions satisfy assumption A1. Consistency is not guaranteed for bounded loss functions, because for bounded ρ( ), it is impossible that ψ(z)z is nondecreasing and strictly increasing in the interval where 11
12 ψ(z)z < K. The assumption A2 of elliptical symmetry guarantees the estimator to be asymptotically unbiased. When the u i s are i.i.d. distributed, but not elliptically symmetric, the M-estimator will still converge (Chapter 8 in Huber, 1981). Finally note that, under regularity conditions, M-estimators are asymptotically normal with asymptotic covariance matrix (Chapter 7 in Hayashi, 2000): { [ ]} 2 1 m(r t ; θ, ρ) E Pθ Σ θ θ { [ ]} 2 1 m(r t ; θ, ρ) E Pθ, (12) θ θ where Σ is the long run covariance matrix of m(r t ; θ ; ρ)/ θ. 3 Robust M-estimation For many applications, a fully specified MGARCH time series model with elliptical innovations is at most a very good approximation of the true data generating process. It may explain well the bulk of the observations, but realworld time series will almost always contain some outlying observations or discrepant substructures. In a conditionally heteroscedastic time series setting, outliers can be defined as observations that are extremely unusual given the past of the series. It is common to do as if there were no outliers in the data and to estimate θ by the value of θ Θ that minimizes the M-function (9). In the Monte Carlo experiment of Section 5, we show that such an approach still yields reliable results in the presence of outliers if (i) a loss function with a derivative 12
13 that decreases sufficiently fast towards zero, is used; (ii) a specification for H t,θ is taken under which innovations have a limited impact on the conditional covariance matrix process. 3.1 The loss function An intuitive view on M-estimation of MGARCH models is that it consists in finding the value of θ Θ for which the estimated conditional covariance matrix equals the conditional expectation of a weighted covariance matrix, with weights proportional to ψ(r t H 1 t,θ r t). Indeed, under the assumptions listed above, the M-estimator is asymptotically equal to the value of θ Θ which solves the following set of equations H t,θ = E Pθ [ σ ψ(r t H 1 t,θ r t)r t r t I r t 1 ]. Note how important the shape of ψ( ), the first derivative of the loss function ρ( ), is. It determines the weight given to each observation as a function of the squared Mahalanobis distance d 2 t,θ in (4). Since outliers are characterized by a large value of d 2 t,θ = r t H 1 t,θ r t, robust M-estimators are characterized by a loss function whose derivative, for large values of its argument, decreases sufficiently fast and tends to zero at infinity. The derivatives of the Gaussian and Student t ν loss functions equal ρ φ (z) = 1 ; ρ t ν (z) = N + ν ν 2 + z. (13) 13
14 We see that the M-estimator with Gaussian loss function is not robust because it gives the same weight to all realizations irrespective of their degree of outlyingness, as measured by d 2 t,θ. However, the derivative of the Student t ν loss functon is decreasing and tends to zero at infinity. The smaller ν is, the more robust the M-estimator with Student t ν loss function will be. 3.2 Impact of news on subsequent volatility predictions Most MGARCH time series models specify the elements of the conditional covariance matrix as an unbounded function of a distributed lag of the squares and cross-products of the components of the innovation vector u t. It follows that the news impact curve, which measures the effect of the innovations on the one-step ahead conditional covariance matrix (Engle and Ng, 1993), is unbounded. Moreover, because of the autoregressive structure of the process, the impact of news propagates to future values of the conditional covariance matrix process. Muler and Yohai (2002,2007) stress that because of this property of unbounded innovation propagation, neglecting outliers renders any inference on autoregressive time series processes extremely sensitive to outliers. For GARCH processes this is aggravated by the fact that the news impact curve is typically quadratic. Hence, neglected outliers may have a very adverse effect on subsequent values of the estimated volatility process. It is possible to bound the impact of outliers on the conditional covariance 14
15 matrix by specifying it as a function H θ ( ) of the weighted innovations ũ t = w(d 2 t,θ )u t instead of the raw innovations u t, where the weight function w( ) downweights extreme observations, that is observations with a large value of d 2 t,θ = r t H 1 t,θ r t = u t u t. When H θ ( ) has the BEKK specification (7), we obtain the following econometric model for conditional covariance matrix: H t,θ = H θ (ũ t 1, ũ t 2,...) = C C + w 2 (d 2 t 1,θ )A H1/2 t 1,θ u t 1u 1/2 t 1 H t 1,θ A + B Ht 1,θ B. (14) For numerical reasons, the weight function w( ) needs to be smooth. In the application we will use Tukey s bisquare weight function: w(z) = min { 1, 1 ( 1 (c/z) 2) 3 }, (15) with c the 95% quantile of u t u t, as defined in (5). Note that only the observations with squared distance above its 95% quantile are downweighted. By construction the weight function tends to zero at infinity such that the impact of outliers on the conditional covariance matrix is limited. The idea of robustifying the conditional covariance matrix against outliers in the data by downweighting past observations of r t in the GARCH specification has already been pursued by Muler and Yohai (2002,2007) in the univariate case. We follow them in calling this the Bounded Innovation Propagation (BIP) GARCH model. 15
16 Muler and Yohai (2002,2007) estimate the parameters of the univariate GARCH model by minimizing the M-function corresponding to the BIP-GARCH model and call the newly obtained estimator the BIP M-estimator. This approach readily extends to the multivariate case in which the BIP M-estimator is defined as the value of θ Θ that minimizes the BIP M-function [ M (S T ; θ, ρ) = ave log det H t,θ t ( + σ ρ r t )] 1 H t,θ r t. (16) 4 Robust diagnostic checks In the absence of outliers and under correct model specification, the squared Mahalanobis distances of the observations in terms of the M-estimate for the conditional covariance matrix, are asymptotically i.i.d. distributed with density function f( ) given in (5). In the same spirit as Rousseeuw and van Zomeren (1990) in the unconditional case, one can use this result to identify MGARCH outliers as these observations for which the estimated squared distance d 2 t,θ exceeds an extreme quantile of f( ). The above result can also be used for verifying correct model specification, using the following test statistic κ = T 1 T T t=1 w(r t H 1 t,θ r t)r t H 1 t,θ r t m, (17) s 16
17 with w( ) a weight function and m = 2πN/2 Γ(N/2) [ 2π N/2 s = Γ(N/2) 0 0 w(t 2 )t N+1 g(t 2 )dt, (18) [ ψ(t 2 )t 2 m ]2 t N 1 g(t )dt] (19) Since under the model assumptions made previously, the distances are conditionally independent and the M-estimator is consistent, one can invoke a central limit theorem to state that this statistic is asymptotically standard normal distributed under correct model specification. If a weight function is used that downweights extreme observations, then this result will be robust to a small proportion of outliers in the data. Its robustness can be increased by using the BIP conditional covariance matrix specification. The null of correct model specification can be tested by comparing κ with the critical values of the standard normal distribution. 5 Application This section illustrates the usefulness of robust estimators and robust MGARCH volatility forecasting models. We do this analysis for the BEKK model, which is a popular and representative member of the family of MGARCH models. We first report the results of a Monte Carlo study on the relative estimation performance of the Gaussian QML estimator, the M-estimator with Student t 4 loss function and the BIP counterparts of these two estimators. We find 17
18 that the BIP M-estimator with Student t 4 loss function calculated under the assumption of normality is robust to additive outliers and misspecification of the consistency factor. We then compare the forecasting performance of the BEKK and BIP-BEKK model for the daily return series of the Nasdaq and NYSE stock market indices. 5.1 A specific MGARCH Model We consider the bivariate symmetric BEKK volatility model of order one, with Student t innovations, proposed by Engle and Kroner (1995). Under this model, the data generating process is given by the following set of equations: i.i.d. r t I rt 1 t ν (0, H t,θ ) ; ν > 2 (20) H t,θ = C C + A r t 1 r t 1A + B H t 1,θ B. The parameter matrices C, A and B all denote 2 2 coefficient matrices. The matrix C is upper triangular and the matrices A and B are symmetric. Let θ be the vector of length 9 that stacks the upper diagonal elements of the matrices C, A and B one on top of the other. Observe that each element of the conditional covariance matrix is specified as a linear function of the lagged squares and cross-products of the elements of r t as well as of the lags of the elements of the conditional covariance matrix. As such, the BEKK model can replicate the stylized fact that volatility clusters over time and co-moves across assets. 18
19 The way in which the conditional covariance matrix process evolves in time depends on the parameter values and the realizations of the stochastic process {u t }. In general, one wants the process H t,θ to be positive definite and covariance stationary whatever be the realizations of {u t }. Engle and Kroner (1995) and Altay-Salih et al. (2003) show that this is guaranteed under the following constraints on the parameter space: C 11 > 0, C 22 > 0, A 11 > 0, B 11 > 0 (21) I 4 A A B B is positive definite. (22) We find the solution of this estimation problem using the DONLP2 sequential quadratic programming algorithm. 5.2 Robustness of M-estimators to additive outliers This section aims to compare the behavior of the (BIP) M-estimators with Gaussian and Student t 4 loss function for simulated data with different levels of outlier contamination. The artificial time series are generated as follows. First outlier-free time series {y t } of length 1000 are generated from model (20), 19
20 with parameter matrices 1 C C = A = B = We then add outliers to the clean series using the following probabilistic setup:. r t = y t + 5ξ t d t. (23) Under this model the occurrence of an outlier in any of the two series is governed by the random process {d t }, modeled as a sequence of i.i.d. draws from a Bernoulli distribution where the occurrence of an outlier (d t = 1) has probability ε. The transmission and the size of the outliers stems from a 2-dimensional i.i.d. vector process {ξ t } constructed such that, when an outlier occurs, it will be either in the first component, second component or both with probability 0.3, 0.3 and 0.4, respectively. The magnitude of the outlier is 5 times the conditional standard deviation of the corresponding element of y t. Since M-estimators are defined by the extremum of the M-function, it is important that the location of the extremum as a function of the parameter vector is robust to outliers in the data. In Figure 2 the M-function (9) with Gaussian and Student t 4 loss function are plotted as a function of the parameter B 22, once in the absence and once in the presence of 5% of outliers. 2 We see 1 The dependence parameter matrices A and B are similar to those in Chapter 18 of Lütkepohl (2005), while the matrix C is chosen such that it has a similar scale as A and B. 2 To save space, we do not report this figure for the BIP M-functions. The main result 20
21 that in the absence of outliers, all M-functions reach an extremum close to the true parameter value, which is 0.8. Because of the quadratic forms in the BEKK specification, the M-function reaches also a local minimum around θ 9 = 0.8. Hence, the importance of having good starting values and of imposing the bound constraints (21) in the optimization process. Note that contaminating 5% of the data with additive outliers increases strongly the variability of the M-function. The magnitude of this outlier induced variability depends on the loss function. For the Gaussian loss function, the variability of the M-function is so large that there is no more visible evidence that the M-function reaches an extremum close to the true parameter value. The outlier induced variability is much smaller for the M-function with Student t 4 loss function and it does not obscure the location of the extremum. is that the M and BIP M-functions that share the same loss function are very similar, even when 5% of the observations are outlier-contaminated. 21
22 Figure 2: Sensitivity of (BIP) M-functions with Gaussian and Student t 4 loss function to the levels of contamination ε. The M-function is taken as a function of B 22 only; the other parameters being fixed at their true values, which is 0.8 for B 22. Gaussian loss: ρ φ Student t 4 loss: ρ t No outliers No outliers 5% outliers 5% outliers M3theta9[, 3] B B 22 A current practice in MGARCH modelling of financial time series is to assume that the innovations are Student t ν distributed and to treat the number of degrees of freedom ν as an additional parameter (Pesaran and Pesaran, 2007). Figure 2 shows that M-estimation with a Student t ν loss function, where the parameter ν in the loss function is an additional estimation parameter, is no longer reliable in the presence of outliers. In Figure 2 the M-functions are 22
23 computed for Gaussian innovations (ν = ). For joint M-estimation of θ and ν to be a reliable estimation procedure, we would need that ρ t (θ) < ρ t4 (θ) for θ close to the true parameter value. We see that this is the case in the absence of outliers, but because outliers blow up the M-function with loss function ρ t ( ) much more than when the loss function ρ t4 ( ) is used, this is no longer true in the presence of outliers. It seems thus that M-estimation of θ and ν jointly is no longer reliable in the presence of outliers. Of course this is only a partial analysis. To study the effect of outlier contamination on the bias and efficiency of the (BIP) M-estimator with (un)bounded Gaussian or Student t 4 loss function, we consider 200 replications of (23) and three levels of outlier contamination (ε = 0, ε = 0.01 and ε = 0.05). For each M-estimator and for each type of simulated series, we aggregate the parameter estimates to compute the mean error (ME) and the root mean squared error (RMSE) as follows, ME = ave ˆθj θ, RMSE = ave ˆθ j θ 2, (24) j j where ave j denotes the average across the 200 replications and is the Euclidian norm operator. These summary statistics are reported in Table 1. We find that in the absence of outliers (ε = 0), the estimation of the BEKK model on the basis of the BIP-BEKK model does not seem to have a large impact on the bias of the M-estimator. 23
24 When the series are outlier contaminated (ε = 0.01 or ε = 0.05), an overall result is that outliers adversely affect the bias as well as the efficiency of all M-estimators. Regarding the relative performance of the M-estimators with Gaussian and Student t 4 loss function, we find that when 1% of the data is outlier contaminated, the M-estimator with Student t 4 loss function is clearly more robust, but when 5% of the observations are affected by outliers, then both estimators collapse. The lack of robustness of the M-estimator with a Student t 4 loss function is because of the unbounded innovation propagation in the conditional covariance matrix process. The BIP version of this estimator should thus be used. A conclusion of this simulation experiment is that, among the considered estimators, no M-estimator yields better robustness results than the BIP M-estimator with Student t 4 loss function. The estimation of the BEKK model is a non-linearly constrained non-linear optimization problem. This is a computationally difficult problem to solve, especially in the real case where data, because of outliers or model misspecification, do not exactly satisfy the model assumptions. Consider e.g. the case where the clean series has BEKK covariance dynamics, but the econometrician assumes a BIP-BEKK model. In this case, the estimation routine does not converge in approximately 10% of the cases. In Table 1 we also see that the convergence frequency of the M-estimator with Gaussian loss function deteriorates when the data contain outliers, while the convergence behavior of the 24
25 Table 1: Mean error (ME), root mean squared error (RMSE) and convergence frequency of the (BIP) M-estimators with Gaussian and Student t 4 loss function, when a proportion ε of the Gaussian BEKK simulated data are outlier contaminated. Loss function Gaussian loss Student t 4 loss BIP no yes no yes ε = 0 ME RMSE conv. freq ε = 0.01 ME RMSE conv. freq ε = 0.05 ME RMSE conv. freq M-estimator with Student t 4 loss function seems to be more robust to outlier contamination. We find thus that weighting the observations in function of their outlyingness, leads to a more stable objective function and renders the optimization procedure feasible in instances in which Gaussian QML estimation is not. 25
26 5.3 Robustness of M-estimators to misspecification of the consistency factor An issue in applying the M-estimator in practice is that it requires knowledge of the true density function through the definition of the consistency factor σ. For the M-estimator to be truly robust, it must be that misspecification of σ should impair the performance only slightly. In Table 2 we report the ME and RMSE for the Gaussian and Student t 4 M-estimator, when the innovation vector in the bivariate BEKK model (20) is Gaussian, Student t 10 or Student t 4 distributed. For each situation, we do M-estimation with a Student t 4 loss function for a consistency factor computed under the nominal assumption that the innovations that are Gaussian, Student t 10 or Student t 4 distributed. Alike Engle and González-Rivera (1991) for the univariate GARCH model, we find that the Gaussian QML estimator has a degree of inefficiency that increases with the degree of departure from normality. For Gaussian innovations, the Gaussian QML estimator has a RMSE of 0.266, while for Student t 4 innovations the RMSE is If the M-estimator with Student t 4 loss function is correctly specified, it has a RMSE of for Student t 4 innovations and a RMSE of for Gaussian innovations. Of course, in real applications the innovation density function is unknown to the econometrician. For financial time series, a Student t density function 26
27 Table 2: Mean error (ME) and root mean squared error (RMSE) of the Gaussian QML estimator and the M-estimators with Student t 4 loss function and consistency factor computed for a Student t density function with ν degrees of freedom. Loss function Gaussian loss Student t 4 loss with σ computed for: ν = ν = 10 ν = 4 Gaussian innovations ME RMSE Student t 10 innovations ME RMSE Student t 4 innovations ME RMSE is often assumed with an unknown number of degrees of freedom ν. When the assumed number of degrees of freedom does not equal the population values, the consistency factor and hence the Student t 4 M-function are misspecified and the consistency result no longer holds. Recall that for the Gaussian QML estimator, there is no problem of correct specification, since for this estimator the consistency factor is always 1. In Table 2 it seems that for the M-estimator with Student t 4 loss function, the estimation performance is more robust to overestimation of the true number of degrees of freedom than to underestima- 27
28 tion of it. Of all M-estimators with Student t ν loss function considered, the estimator with σ calculated for ν = 10 is most robust to misspecification of σ. Its finite sample bias varies between (correct specification) and (Student t 4 innovations), while for the other M-estimators the finite-sample bias more than doubles when σ is misspecified. If the true number of degrees of freedom is unknown, we therefore recommend to calculate the consistency factor σ for an intermediate value of ν such as ν = 10. In Table 2 we do not report the convergence frequencies, because for all M-estimators with Gaussian and Student t 4 loss function, it was close to 1. A failure of convergence is thus not due to the thickness of the tails nor due to misspecification of the consistency factor. The decrease in convergence frequency when the data are contaminated with additive outliers or when the BIP M-estimators are used (Table 1), indicate that a failure to converge is a manifestation of additive outliers and/or MGARCH model misspecification. 5.4 Empirical application We apply the (BIP) M-estimators with Gaussian and Student t 4 loss function to daily returns of the Nasdaq and NYSE composite indices over the period January 2, December 29, The data were downloaded from datastream and consist of 7043 observations. The series are plotted in Figure 3. Note the large negative return corresponding to the stock market crash 28
29 Figure 3: Daily returns on Nasdaq and NYSE composite index. Nasdaq NYSE 12/80 06/86 11/91 05/97 11/02 of October 19, Interestingly, the Nasdaq index had returns of similar magnitude in the year In an unconditional framework, these returns are qualified as outliers, whereas under a conditionally heteroscedastic time series model, some of these extreme returns can be explained by the model. Observe that not only there is volatility clustering in both series, but also that there is a lot of commonality in volatility across the two series. We assume that the very large majority of these observations are generated by the conditionally Gaussian BEKK time series process in (20). Return realizations such as the 1987 stock market crash are clearly outliers under this 29
30 Figure 4: Conditional BEKK and BIP-BEKK standard deviation of the daily Nasdaq and NYSE return series BEKK, Nasdaq BIP BEKK, Nasdaq BEKK, NYSE BIP BEKK, NYSE 12/80 06/86 11/91 05/97 11/02 12/80 06/86 11/91 05/97 11/02 parametric framework. Hence a robust estimator is needed. The M-estimator with Student t 4 loss function yields the following parameter estimate Ĉ C = Â = ˆB = The sequential quadratic programming algorithm failed to converge for the (BIP) M-estimator with Gaussian loss function and for the BIP M-estimator 30.
31 with Student t 4 loss function. Since in this subsection, the focus is on comparing the BEKK and BIP-BEKK volatility forecasting models, given a reliable parameter estimate, we abstract from this problem. It is of particular interest to compare the BEKK and BIP-BEKK volatility predictions for the period of the 1987 stock market crash and the period of the burst of the internet bubble in In each of these periods, large returns of similar magnitude occur but the persistence in volatility is different. The stock market crash in 1987 was not followed by returns of similar size whereas the burst of the internet bubble corresponds to a prolonged period of high volatility. As can be seen in Figure 4, we have that due to the property of a limited effect of innovations on future volatility, the BIP-BEKK model captures well this difference in reaction to shocks. In contrast, the BEKK model overestimates the persistence in the realized return variability in the period subsequent to the 1987 crash. The same interpretation holds for the patch of extreme returns on the Nasdaq index on March 27-28, At all other times, the BEKK and BIP-BEKK conditional standard deviation forecasts are similar. 31
32 6 Concluding remarks This paper introduces a class of M-estimators for multivariate GARCH time series models. We prove that under general conditions M-estimators are consistent. We study the effect of the shape of the loss function and of bounding the news impact curve on the robustness of the M-estimator to outliers. The popular Gaussian quasi-maximum likelihood estimator is shown to be very sensitive to outliers in the data. If the data is suspected to be contaminated by outliers, we recommend the use of a volatility model that has the property of bounded innovation propagation and to use a M-estimator with Student t 4 loss function. The Monte Carlo study and empirical application document the good robustness properties of this approach. This research can be extended by looking at other classes of robust estimators and by considering other multivariate GARCH models. In recent research on these models, the tendency has been to reduce the dimensionality of the problem using principal component analysis or by specifying separate univariate GARCH models for the volatility of each series and a time-varying process for the conditional correlations. Further research could show how well robust methods work for these types of models. 32
33 7 Acknowledgment The authors are grateful to Tim Bollerslev, Geert Dhaene, Sébastien Laurent and Oliver Linton for very helpful comments and suggestions. Financial support from the FWO and Research Fund K.U.Leuven is gratefully acknowledged. The usual disclaimer applies. A Consistency of the M-estimator First we derive the asymptotic limit of the gradient of the M-function. Then we show under which assumptions θ is the only value in Θ for which this probability limit is zero. Denote m (r t ; θ, ρ) = log det H t,θ + σ ρ ( r t H 1 t,θ r t). Its first derivative ṁ i (r t ; θ, ρ) = m (r t ; θ, ρ) / θ i is given by Tr [ I N σ ψ(r t H 1 t,θ r ] t)r t r t H 1 H t,θ t,θ H 1 t,θ θ, (25) i (Comte and Lieberman, 2003). The square N-dimensional matrix H t,θ / θ i holds the i-th partial derivative of the corresponding elements of H t,θ and Tr is the trace operator. Since ṁ i (r t ; θ, ρ) is a measurable function of the strictly stationary and ergodic process {r t }, it is also strictly stationary and ergodic. By the Ergodic Theorem its time series average, which is the score of the M- function, will converge to its ensemble mean, which is the mean of all possible 33
34 realizations of (25) conditional on the past of {u t }. We thus find that the asymptotic limit of the score of the M-function equals { [IN Tr E Pθ σ ψ(r t H 1 t,θ r ] t)r t r t H t,θ H 1 t,θ θ i H 1 t,θ I u t 1 } {[ = Tr E Pθ IN σ ψ(r t H 1 t,θ r ] t)r t r t H 1 H t,θ t,θ Iut 1} H 1 t,θ θ, (26) i since H t,θ = H θ (u t 1,...) is a deterministic N N matrix function of the past realizations of the process {u t }. Under the no dominance assumption A4, convergence in probability will be uniform. Borrowing results from the theory on M-estimation of unconditional scale, it can be shown that under A1-A3, the score function is zero for θ. Indeed, Maronna (1976) proves that under A1, there exists a unique solution to the following implicit equation determining H t,θ : {[ E Pθ IN σ ψ(r t H 1 t,θ r ] } t)r t r t H 1 t,θ Iut 1 = 0. (27) The solution to this equation is called the M-estimator of scatter of r t or also the pseudo-covariance matrix of r t. Since this estimator is affine equivariant and since r t is assumed to be elliptically symmetric with conditional scatter matrix H t,θ, one can show that this solution will be proportional to the true conditional scatter matrix H t,θ, with proportionality factor c equal to c = N σ { 2π N/2 Γ(N/2) 0 ψ(t 2 )t N+1 g (t 2 )dt} 1, (28) (Chapter 13 in Bilodeau and Brenner, 1999, Chapter 6 in Maronna et al., 2006). It follows that M-estimators will be consistent since σ in (11) is such 34
35 that c = 1. By the no observational equivalence assumption A3, it is only for θ that the conditional covariance matrix equals H t,θ. We may thus conclude that for loss functions satisfying A1, the limit of the M-function will have an extremum in θ. This extremum is unique under the assumption A5 that there is no θ Θ for which the score of the M-function is zero while the conditional expectation under P θ of the weighted covariance matrix is not the true conditional covariance matrix H t,θ. A further sufficient condition to ensure that the extremum of the limit is the limit of the extremum of the M-estimator is that the parameter space Θ is compact. 35
36 References [1] Altay-Salih, A., Pinar, M.C. and S. Leyffer, 2003, Constrained nonlinear programming for volatility estimation with GARCH models. SIAM review 45, [2] Bauwens, L., Laurent, S. and J.V.K. Rombouts, 2006, Multivariate GARCH models: a survey. Journal of Applied Econometrics 21, [3] Bilodeau, M. and D. Brenner, 1999, Theory of Multivariate Statistics. Springer, New York. [4] Comte, F. and O. Lieberman, 2003, Asymptotic theory for multivariate GARCH processes. Journal of Multivariate Analysis 84, [5] Engle, R.F. and G. González-Rivera, 1991, Semiparametric ARCH models. Journal of Business and Economic Statistics 9, [6] Engle, R.F. and F.K. Kroner, 1995, Multivariate simultaneous generalized ARCH. Econometric Theory 11, [7] Engle, R.F. and V.K. Ng, 1993, Measuring and testing the impact of news on volatility. The Journal of Finance 58, [8] Fiorentini, G., Sentana, E. and G. Calzolari, 2003, Maximum likelihood estimation and inference in multivariate conditionally heteroscedastic dynamic regression models with Student t innovations. Journal of Business and Economic Statistics 21,
37 [9] Franses, P.H. and H. Ghijsels, 1999, Additive outliers, GARCH and forecasting volatility. International Journal of Forecasting 15, 1-9. [10] Hayashi, F., 2000, Econometrics. Princeton University Press, Princeton. [11] Huber, P.J., 1981, Robust Statistics. Wiley, New York. [12] Jeantheau, T., 1998, Strong consistency of estimators for multivariate ARCH models. Econometric Theory 14, [13] Lütkepohl, H., 2005, New Introduction to Multiple Time Series Analysis. Springer-Verlag, Berlin. [14] Mancini, L., Ronchetti, E. and F. Trojani, 2005, Optimal conditionally unbiased bounded-influence inference in dynamic location and scale models. Journal of the American Statistical Association 105, [15] Maronna, R.A., 1976, Robust M-estimators of multivariate location and scatter. The Annals of Statistics 4, [16] Maronna, R.A., Martin, R.D. and V.J. Yohai, 2006, Robust Statistics. Wiley, New York. [17] Muler, N. and V.J. Yohai, 2002, Robust estimates for ARCH processes. Journal of Time Series Analysis 23, [18] Muler, N. and V.J. Yohai, 2007, Robust estimates for GARCH models. Journal of Statistical Planning and Inference, forthcoming. 37
38 [19] Newey, W.K. and D.G. Steigerwald, 1997, Asymptotic bias for quasi-maximumlikelihood estimators in conditional heteroskedasticity models. Econometrica 65, [20] Park, B., 2002, An outlier robust GARCH model and forecasting volatility of exchange rate returns. Journal of Forecasting 21, [21] Peng, L. and Q. Yao, 2003, Miscellanea. Least absolute deviations estimation for ARCH and GARCH models. Biometrika 90, [22] Pesaran, B. and M.H. Pesaran, 2007, Modelling volatilities and conditional correlations in futures markets with a multivariate t distribution. Cambridge working paper. [23] Rousseeuw, P.J. and B.C. van Zomeren, 1990, Unmasking multivariate outliers and leverage points. Journal of the American Statistical Association 85, [24] Sakata, S. and H. White, 1998, High breakdown point conditional dispersion estimation with application to S&P 500 daily returns volatility. Econometrica 66, [25] Van Dijk, D., Franses, P.H. and A. Lucas, 1999, Testing for ARCH in the presence of additive outliers. Journal of Applied Econometrics 14,
GARCH Models Estimation and Inference
GARCH Models Estimation and Inference Eduardo Rossi University of Pavia December 013 Rossi GARCH Financial Econometrics - 013 1 / 1 Likelihood function The procedure most often used in estimating θ 0 in
More informationMultivariate GARCH models.
Multivariate GARCH models. Financial market volatility moves together over time across assets and markets. Recognizing this commonality through a multivariate modeling framework leads to obvious gains
More informationBayesian Semiparametric GARCH Models
Bayesian Semiparametric GARCH Models Xibin (Bill) Zhang and Maxwell L. King Department of Econometrics and Business Statistics Faculty of Business and Economics xibin.zhang@monash.edu Quantitative Methods
More informationBayesian Semiparametric GARCH Models
Bayesian Semiparametric GARCH Models Xibin (Bill) Zhang and Maxwell L. King Department of Econometrics and Business Statistics Faculty of Business and Economics xibin.zhang@monash.edu Quantitative Methods
More informationVolatility. Gerald P. Dwyer. February Clemson University
Volatility Gerald P. Dwyer Clemson University February 2016 Outline 1 Volatility Characteristics of Time Series Heteroskedasticity Simpler Estimation Strategies Exponentially Weighted Moving Average Use
More informationAnalytical derivates of the APARCH model
Analytical derivates of the APARCH model Sébastien Laurent Forthcoming in Computational Economics October 24, 2003 Abstract his paper derives analytical expressions for the score of the APARCH model of
More informationDEPARTMENT OF ECONOMICS AND FINANCE COLLEGE OF BUSINESS AND ECONOMICS UNIVERSITY OF CANTERBURY CHRISTCHURCH, NEW ZEALAND
DEPARTMENT OF ECONOMICS AND FINANCE COLLEGE OF BUSINESS AND ECONOMICS UNIVERSITY OF CANTERBURY CHRISTCHURCH, NEW ZEALAND Discussion of Principal Volatility Component Analysis by Yu-Pin Hu and Ruey Tsay
More informationA Bootstrap Test for Causality with Endogenous Lag Length Choice. - theory and application in finance
CESIS Electronic Working Paper Series Paper No. 223 A Bootstrap Test for Causality with Endogenous Lag Length Choice - theory and application in finance R. Scott Hacker and Abdulnasser Hatemi-J April 200
More informationGeneralized Autoregressive Score Models
Generalized Autoregressive Score Models by: Drew Creal, Siem Jan Koopman, André Lucas To capture the dynamic behavior of univariate and multivariate time series processes, we can allow parameters to be
More informationM-estimators for augmented GARCH(1,1) processes
M-estimators for augmented GARCH(1,1) processes Freiburg, DAGStat 2013 Fabian Tinkl 19.03.2013 Chair of Statistics and Econometrics FAU Erlangen-Nuremberg Outline Introduction The augmented GARCH(1,1)
More informationGARCH Models Estimation and Inference. Eduardo Rossi University of Pavia
GARCH Models Estimation and Inference Eduardo Rossi University of Pavia Likelihood function The procedure most often used in estimating θ 0 in ARCH models involves the maximization of a likelihood function
More informationGARCH Models Estimation and Inference
Università di Pavia GARCH Models Estimation and Inference Eduardo Rossi Likelihood function The procedure most often used in estimating θ 0 in ARCH models involves the maximization of a likelihood function
More informationInference in VARs with Conditional Heteroskedasticity of Unknown Form
Inference in VARs with Conditional Heteroskedasticity of Unknown Form Ralf Brüggemann a Carsten Jentsch b Carsten Trenkler c University of Konstanz University of Mannheim University of Mannheim IAB Nuremberg
More informationGaussian kernel GARCH models
Gaussian kernel GARCH models Xibin (Bill) Zhang and Maxwell L. King Department of Econometrics and Business Statistics Faculty of Business and Economics 7 June 2013 Motivation A regression model is often
More informationEcon 423 Lecture Notes: Additional Topics in Time Series 1
Econ 423 Lecture Notes: Additional Topics in Time Series 1 John C. Chao April 25, 2017 1 These notes are based in large part on Chapter 16 of Stock and Watson (2011). They are for instructional purposes
More informationGARCH Models. Eduardo Rossi University of Pavia. December Rossi GARCH Financial Econometrics / 50
GARCH Models Eduardo Rossi University of Pavia December 013 Rossi GARCH Financial Econometrics - 013 1 / 50 Outline 1 Stylized Facts ARCH model: definition 3 GARCH model 4 EGARCH 5 Asymmetric Models 6
More informationA Non-Parametric Approach of Heteroskedasticity Robust Estimation of Vector-Autoregressive (VAR) Models
Journal of Finance and Investment Analysis, vol.1, no.1, 2012, 55-67 ISSN: 2241-0988 (print version), 2241-0996 (online) International Scientific Press, 2012 A Non-Parametric Approach of Heteroskedasticity
More informationAccounting for Missing Values in Score- Driven Time-Varying Parameter Models
TI 2016-067/IV Tinbergen Institute Discussion Paper Accounting for Missing Values in Score- Driven Time-Varying Parameter Models André Lucas Anne Opschoor Julia Schaumburg Faculty of Economics and Business
More informationLecture 8: Multivariate GARCH and Conditional Correlation Models
Lecture 8: Multivariate GARCH and Conditional Correlation Models Prof. Massimo Guidolin 20192 Financial Econometrics Winter/Spring 2018 Overview Three issues in multivariate modelling of CH covariances
More informationThe Size and Power of Four Tests for Detecting Autoregressive Conditional Heteroskedasticity in the Presence of Serial Correlation
The Size and Power of Four s for Detecting Conditional Heteroskedasticity in the Presence of Serial Correlation A. Stan Hurn Department of Economics Unversity of Melbourne Australia and A. David McDonald
More informationRobust Backtesting Tests for Value-at-Risk Models
Robust Backtesting Tests for Value-at-Risk Models Jose Olmo City University London (joint work with Juan Carlos Escanciano, Indiana University) Far East and South Asia Meeting of the Econometric Society
More informationStudies in Nonlinear Dynamics & Econometrics
Studies in Nonlinear Dynamics & Econometrics Volume 9, Issue 2 2005 Article 4 A Note on the Hiemstra-Jones Test for Granger Non-causality Cees Diks Valentyn Panchenko University of Amsterdam, C.G.H.Diks@uva.nl
More information1 Phelix spot and futures returns: descriptive statistics
MULTIVARIATE VOLATILITY MODELING OF ELECTRICITY FUTURES: ONLINE APPENDIX Luc Bauwens 1, Christian Hafner 2, and Diane Pierret 3 October 13, 2011 1 Phelix spot and futures returns: descriptive statistics
More informationIntro VEC and BEKK Example Factor Models Cond Var and Cor Application Ref 4. MGARCH
ntro VEC and BEKK Example Factor Models Cond Var and Cor Application Ref 4. MGARCH JEM 140: Quantitative Multivariate Finance ES, Charles University, Prague Summer 2018 JEM 140 () 4. MGARCH Summer 2018
More informationECONOMETRIC REVIEWS, 5(1), (1986) MODELING THE PERSISTENCE OF CONDITIONAL VARIANCES: A COMMENT
ECONOMETRIC REVIEWS, 5(1), 51-56 (1986) MODELING THE PERSISTENCE OF CONDITIONAL VARIANCES: A COMMENT Professors Engle and Bollerslev have delivered an excellent blend of "forest" and "trees"; their important
More informationDEPARTMENT OF ECONOMICS
ISSN 0819-2642 ISBN 0 7340 2601 3 THE UNIVERSITY OF MELBOURNE DEPARTMENT OF ECONOMICS RESEARCH PAPER NUMBER 945 AUGUST 2005 TESTING FOR ASYMMETRY IN INTEREST RATE VOLATILITY IN THE PRESENCE OF A NEGLECTED
More informationCHICAGO: A Fast and Accurate Method for Portfolio Risk Calculation
CHICAGO: A Fast and Accurate Method for Portfolio Risk Calculation University of Zürich April 28 Motivation Aim: Forecast the Value at Risk of a portfolio of d assets, i.e., the quantiles of R t = b r
More informationArma-Arch Modeling Of The Returns Of First Bank Of Nigeria
Arma-Arch Modeling Of The Returns Of First Bank Of Nigeria Emmanuel Alphonsus Akpan Imoh Udo Moffat Department of Mathematics and Statistics University of Uyo, Nigeria Ntiedo Bassey Ekpo Department of
More informationISSN Article. Selection Criteria in Regime Switching Conditional Volatility Models
Econometrics 2015, 3, 289-316; doi:10.3390/econometrics3020289 OPEN ACCESS econometrics ISSN 2225-1146 www.mdpi.com/journal/econometrics Article Selection Criteria in Regime Switching Conditional Volatility
More informationStock index returns density prediction using GARCH models: Frequentist or Bayesian estimation?
MPRA Munich Personal RePEc Archive Stock index returns density prediction using GARCH models: Frequentist or Bayesian estimation? Ardia, David; Lennart, Hoogerheide and Nienke, Corré aeris CAPITAL AG,
More informationA simple graphical method to explore tail-dependence in stock-return pairs
A simple graphical method to explore tail-dependence in stock-return pairs Klaus Abberger, University of Konstanz, Germany Abstract: For a bivariate data set the dependence structure can not only be measured
More informationLocation Multiplicative Error Model. Asymptotic Inference and Empirical Analysis
: Asymptotic Inference and Empirical Analysis Qian Li Department of Mathematics and Statistics University of Missouri-Kansas City ql35d@mail.umkc.edu October 29, 2015 Outline of Topics Introduction GARCH
More informationA note on adaptation in garch models Gloria González-Rivera a a
This article was downloaded by: [CDL Journals Account] On: 3 February 2011 Access details: Access Details: [subscription number 922973516] Publisher Taylor & Francis Informa Ltd Registered in England and
More informationTesting an Autoregressive Structure in Binary Time Series Models
ömmföäflsäafaäsflassflassflas ffffffffffffffffffffffffffffffffffff Discussion Papers Testing an Autoregressive Structure in Binary Time Series Models Henri Nyberg University of Helsinki and HECER Discussion
More informationDEPARTMENT OF ECONOMICS
ISSN 0819-64 ISBN 0 7340 616 1 THE UNIVERSITY OF MELBOURNE DEPARTMENT OF ECONOMICS RESEARCH PAPER NUMBER 959 FEBRUARY 006 TESTING FOR RATE-DEPENDENCE AND ASYMMETRY IN INFLATION UNCERTAINTY: EVIDENCE FROM
More informationAsymptotic inference for a nonstationary double ar(1) model
Asymptotic inference for a nonstationary double ar() model By SHIQING LING and DONG LI Department of Mathematics, Hong Kong University of Science and Technology, Hong Kong maling@ust.hk malidong@ust.hk
More informationFinancial Econometrics Return Predictability
Financial Econometrics Return Predictability Eric Zivot March 30, 2011 Lecture Outline Market Efficiency The Forms of the Random Walk Hypothesis Testing the Random Walk Hypothesis Reading FMUND, chapter
More informationSPECIFICATION TESTS IN PARAMETRIC VALUE-AT-RISK MODELS
SPECIFICATION TESTS IN PARAMETRIC VALUE-AT-RISK MODELS J. Carlos Escanciano Indiana University, Bloomington, IN, USA Jose Olmo City University, London, UK Abstract One of the implications of the creation
More informationGenerated Covariates in Nonparametric Estimation: A Short Review.
Generated Covariates in Nonparametric Estimation: A Short Review. Enno Mammen, Christoph Rothe, and Melanie Schienle Abstract In many applications, covariates are not observed but have to be estimated
More informationDiagnostic Test for GARCH Models Based on Absolute Residual Autocorrelations
Diagnostic Test for GARCH Models Based on Absolute Residual Autocorrelations Farhat Iqbal Department of Statistics, University of Balochistan Quetta-Pakistan farhatiqb@gmail.com Abstract In this paper
More informationA TIME SERIES PARADOX: UNIT ROOT TESTS PERFORM POORLY WHEN DATA ARE COINTEGRATED
A TIME SERIES PARADOX: UNIT ROOT TESTS PERFORM POORLY WHEN DATA ARE COINTEGRATED by W. Robert Reed Department of Economics and Finance University of Canterbury, New Zealand Email: bob.reed@canterbury.ac.nz
More informationQuaderni di Dipartimento. Small Sample Properties of Copula-GARCH Modelling: A Monte Carlo Study. Carluccio Bianchi (Università di Pavia)
Quaderni di Dipartimento Small Sample Properties of Copula-GARCH Modelling: A Monte Carlo Study Carluccio Bianchi (Università di Pavia) Maria Elena De Giuli (Università di Pavia) Dean Fantazzini (Moscow
More informationThe Influence of Additive Outliers on the Performance of Information Criteria to Detect Nonlinearity
The Influence of Additive Outliers on the Performance of Information Criteria to Detect Nonlinearity Saskia Rinke 1 Leibniz University Hannover Abstract In this paper the performance of information criteria
More informationThomas J. Fisher. Research Statement. Preliminary Results
Thomas J. Fisher Research Statement Preliminary Results Many applications of modern statistics involve a large number of measurements and can be considered in a linear algebra framework. In many of these
More informationROBUST ESTIMATION OF A CORRELATION COEFFICIENT: AN ATTEMPT OF SURVEY
ROBUST ESTIMATION OF A CORRELATION COEFFICIENT: AN ATTEMPT OF SURVEY G.L. Shevlyakov, P.O. Smirnov St. Petersburg State Polytechnic University St.Petersburg, RUSSIA E-mail: Georgy.Shevlyakov@gmail.com
More informationResearch Article Optimal Portfolio Estimation for Dependent Financial Returns with Generalized Empirical Likelihood
Advances in Decision Sciences Volume 2012, Article ID 973173, 8 pages doi:10.1155/2012/973173 Research Article Optimal Portfolio Estimation for Dependent Financial Returns with Generalized Empirical Likelihood
More informationThe Instability of Correlations: Measurement and the Implications for Market Risk
The Instability of Correlations: Measurement and the Implications for Market Risk Prof. Massimo Guidolin 20254 Advanced Quantitative Methods for Asset Pricing and Structuring Winter/Spring 2018 Threshold
More informationGMM and SMM. 1. Hansen, L Large Sample Properties of Generalized Method of Moments Estimators, Econometrica, 50, p
GMM and SMM Some useful references: 1. Hansen, L. 1982. Large Sample Properties of Generalized Method of Moments Estimators, Econometrica, 50, p. 1029-54. 2. Lee, B.S. and B. Ingram. 1991 Simulation estimation
More informationConfidence Intervals for the Autocorrelations of the Squares of GARCH Sequences
Confidence Intervals for the Autocorrelations of the Squares of GARCH Sequences Piotr Kokoszka 1, Gilles Teyssière 2, and Aonan Zhang 3 1 Mathematics and Statistics, Utah State University, 3900 Old Main
More informationBayesian semiparametric GARCH models
ISSN 1440-771X Australia Department of Econometrics and Business Statistics http://www.buseco.monash.edu.au/depts/ebs/pubs/wpapers/ Bayesian semiparametric GARCH models Xibin Zhang and Maxwell L. King
More informationModel Specification Test with Correlated but not Cointegrated Variables
Model Specification Test with Correlated but not Cointegrated Variables Li Gan Department of Economics, Texas A&M University, College Station, TX 77843-4228 Cheng Hsiao Department of Economics, University
More informationFinite Sample Performance of the MLE in GARCH(1,1): When the Parameter on the Lagged Squared Residual Is Close to Zero
Finite Sample Performance of the MLE in GARCH1,1): When the Parameter on the Lagged Squared Residual Is Close to Zero Suduk Kim Department of Economics Hoseo University Asan Si, ChungNam, Korea, 336-795
More informationHeteroskedasticity in Time Series
Heteroskedasticity in Time Series Figure: Time Series of Daily NYSE Returns. 206 / 285 Key Fact 1: Stock Returns are Approximately Serially Uncorrelated Figure: Correlogram of Daily Stock Market Returns.
More informationA Goodness-of-fit Test for Copulas
A Goodness-of-fit Test for Copulas Artem Prokhorov August 2008 Abstract A new goodness-of-fit test for copulas is proposed. It is based on restrictions on certain elements of the information matrix and
More informationIntroduction to Robust Statistics. Elvezio Ronchetti. Department of Econometrics University of Geneva Switzerland.
Introduction to Robust Statistics Elvezio Ronchetti Department of Econometrics University of Geneva Switzerland Elvezio.Ronchetti@metri.unige.ch http://www.unige.ch/ses/metri/ronchetti/ 1 Outline Introduction
More informationECON3327: Financial Econometrics, Spring 2016
ECON3327: Financial Econometrics, Spring 2016 Wooldridge, Introductory Econometrics (5th ed, 2012) Chapter 11: OLS with time series data Stationary and weakly dependent time series The notion of a stationary
More informationIntroduction Robust regression Examples Conclusion. Robust regression. Jiří Franc
Robust regression Robust estimation of regression coefficients in linear regression model Jiří Franc Czech Technical University Faculty of Nuclear Sciences and Physical Engineering Department of Mathematics
More informationBootstrap tests of multiple inequality restrictions on variance ratios
Economics Letters 91 (2006) 343 348 www.elsevier.com/locate/econbase Bootstrap tests of multiple inequality restrictions on variance ratios Jeff Fleming a, Chris Kirby b, *, Barbara Ostdiek a a Jones Graduate
More informationDYNAMIC CONDITIONAL CORRELATIONS FOR ASYMMETRIC PROCESSES
J. Japan Statist. Soc. Vol. 41 No. 2 2011 143 157 DYNAMIC CONDITIONAL CORRELATIONS FOR ASYMMETRIC PROCESSES Manabu Asai* and Michael McAleer** *** **** ***** The paper develops a new Dynamic Conditional
More informationGeneralized Method of Moments Estimation
Generalized Method of Moments Estimation Lars Peter Hansen March 0, 2007 Introduction Generalized methods of moments (GMM) refers to a class of estimators which are constructed from exploiting the sample
More informationWhen is a copula constant? A test for changing relationships
When is a copula constant? A test for changing relationships Fabio Busetti and Andrew Harvey Bank of Italy and University of Cambridge November 2007 usetti and Harvey (Bank of Italy and University of Cambridge)
More informationNew Statistical Model for the Enhancement of Noisy Speech
New Statistical Model for the Enhancement of Noisy Speech Electrical Engineering Department Technion - Israel Institute of Technology February 22, 27 Outline Problem Formulation and Motivation 1 Problem
More informationOn the Estimation and Application of Max-Stable Processes
On the Estimation and Application of Max-Stable Processes Zhengjun Zhang Department of Statistics University of Wisconsin Madison, WI 53706, USA Co-author: Richard Smith EVA 2009, Fort Collins, CO Z. Zhang
More informationLecture 6: Univariate Volatility Modelling: ARCH and GARCH Models
Lecture 6: Univariate Volatility Modelling: ARCH and GARCH Models Prof. Massimo Guidolin 019 Financial Econometrics Winter/Spring 018 Overview ARCH models and their limitations Generalized ARCH models
More informationCointegration Lecture I: Introduction
1 Cointegration Lecture I: Introduction Julia Giese Nuffield College julia.giese@economics.ox.ac.uk Hilary Term 2008 2 Outline Introduction Estimation of unrestricted VAR Non-stationarity Deterministic
More informationEstimation of Vector Error Correction Model with GARCH Errors. Koichi Maekawa Hiroshima University of Economics
Estimation of Vector Error Correction Model with GARCH Errors Koichi Maekawa Hiroshima University of Economics Kusdhianto Setiawan Hiroshima University of Economics and Gadjah Mada University Abstract
More informationON THE CALCULATION OF A ROBUST S-ESTIMATOR OF A COVARIANCE MATRIX
STATISTICS IN MEDICINE Statist. Med. 17, 2685 2695 (1998) ON THE CALCULATION OF A ROBUST S-ESTIMATOR OF A COVARIANCE MATRIX N. A. CAMPBELL *, H. P. LOPUHAA AND P. J. ROUSSEEUW CSIRO Mathematical and Information
More informationModel estimation through matrix equations in financial econometrics
Model estimation through matrix equations in financial econometrics Federico Poloni 1 Joint work with Giacomo Sbrana 2 1 Technische Universität Berlin (A. Von Humboldt postdoctoral fellow) 2 Rouen Business
More informationA Guide to Modern Econometric:
A Guide to Modern Econometric: 4th edition Marno Verbeek Rotterdam School of Management, Erasmus University, Rotterdam B 379887 )WILEY A John Wiley & Sons, Ltd., Publication Contents Preface xiii 1 Introduction
More informationEconometric Forecasting
Robert M. Kunst robert.kunst@univie.ac.at University of Vienna and Institute for Advanced Studies Vienna October 1, 2014 Outline Introduction Model-free extrapolation Univariate time-series models Trend
More informationLabor-Supply Shifts and Economic Fluctuations. Technical Appendix
Labor-Supply Shifts and Economic Fluctuations Technical Appendix Yongsung Chang Department of Economics University of Pennsylvania Frank Schorfheide Department of Economics University of Pennsylvania January
More informationECONOMICS 7200 MODERN TIME SERIES ANALYSIS Econometric Theory and Applications
ECONOMICS 7200 MODERN TIME SERIES ANALYSIS Econometric Theory and Applications Yongmiao Hong Department of Economics & Department of Statistical Sciences Cornell University Spring 2019 Time and uncertainty
More informationMultivariate ARCH Models: Finite Sample Properties of QML Estimators and an Application to an LM-Type Test
Multivariate ARCH Models: Finite Sample Properties of QML Estimators and an Application to an LM-ype est E M. I and G D.A. P University of Alicante Cardiff University April,004 Abstract his paper provides
More informationCombined Lagrange Multipier Test for ARCH in Vector Autoregressive Models
MEDDELANDEN FRÅN SVENSKA HANDELSHÖGSKOLAN HANKEN SCHOOL OF ECONOMICS WORKING PAPERS 563 Paul Catani and Niklas Ahlgren Combined Lagrange Multipier Test for ARCH in Vector Autoregressive Models 2016 Combined
More informationThe autocorrelation and autocovariance functions - helpful tools in the modelling problem
The autocorrelation and autocovariance functions - helpful tools in the modelling problem J. Nowicka-Zagrajek A. Wy lomańska Institute of Mathematics and Computer Science Wroc law University of Technology,
More informationDEPARTMENT OF ECONOMICS AND FINANCE COLLEGE OF BUSINESS AND ECONOMICS UNIVERSITY OF CANTERBURY CHRISTCHURCH, NEW ZEALAND
DEPARTMENT OF ECONOMICS AND FINANCE COLLEGE OF BUSINESS AND ECONOMICS UNIVERSITY OF CANTERBURY CHRISTCHURCH, NEW ZEALAND Do We Really Need Both BEKK and DCC? A Tale of Two Multivariate GARCH Models Massimiliano
More informationRobust Exponential Smoothing of Multivariate Time Series
Robust Exponential Smoothing of Multivariate Time Series Christophe Croux,a, Sarah Gelper b, Koen Mahieu a a Faculty of Business and Economics, K.U.Leuven, Naamsestraat 69, 3000 Leuven, Belgium b Erasmus
More informationDo Markov-Switching Models Capture Nonlinearities in the Data? Tests using Nonparametric Methods
Do Markov-Switching Models Capture Nonlinearities in the Data? Tests using Nonparametric Methods Robert V. Breunig Centre for Economic Policy Research, Research School of Social Sciences and School of
More informationTesting for Regime Switching in Singaporean Business Cycles
Testing for Regime Switching in Singaporean Business Cycles Robert Breunig School of Economics Faculty of Economics and Commerce Australian National University and Alison Stegman Research School of Pacific
More informationQuantile regression and heteroskedasticity
Quantile regression and heteroskedasticity José A. F. Machado J.M.C. Santos Silva June 18, 2013 Abstract This note introduces a wrapper for qreg which reports standard errors and t statistics that are
More informationA Heavy-Tailed Distribution for ARCH Residuals with Application to Volatility Prediction
ANNALS OF ECONOMICS AND FINANCE 5, 283 298 (2004) A Heavy-Tailed Distribution for ARCH Residuals with Application to Volatility Prediction Dimitris N. Politis * Department of Mathematics University of
More informationADVANCED FINANCIAL ECONOMETRICS PROF. MASSIMO GUIDOLIN
Massimo Guidolin Massimo.Guidolin@unibocconi.it Dept. of Finance ADVANCED FINANCIAL ECONOMETRICS PROF. MASSIMO GUIDOLIN a.a. 14/15 p. 1 LECTURE 3: REVIEW OF BASIC ESTIMATION METHODS: GMM AND OTHER EXTREMUM
More informationOn the threshold hyperbolic GARCH models. Citation Statistics And Its Interface, 2011, v. 4 n. 2, p
Title On the threshold hyperbolic GARCH models Authors) Kwan, W; Li, WK; Li, G Citation Statistics And Its Interface, 20, v. 4 n. 2, p. 59-66 Issued Date 20 URL http://hdl.hle.net/0722/35494 Rights Statistics
More informationEstimating GARCH models: when to use what?
Econometrics Journal (2008), volume, pp. 27 38. doi: 0./j.368-423X.2008.00229.x Estimating GARCH models: when to use what? DA HUANG, HANSHENG WANG AND QIWEI YAO, Guanghua School of Management, Peking University,
More informationJuly 31, 2009 / Ben Kedem Symposium
ing The s ing The Department of Statistics North Carolina State University July 31, 2009 / Ben Kedem Symposium Outline ing The s 1 2 s 3 4 5 Ben Kedem ing The s Ben has made many contributions to time
More informationFinancial Econometrics
Financial Econometrics Nonlinear time series analysis Gerald P. Dwyer Trinity College, Dublin January 2016 Outline 1 Nonlinearity Does nonlinearity matter? Nonlinear models Tests for nonlinearity Forecasting
More informationRobust forecasting of non-stationary time series DEPARTMENT OF DECISION SCIENCES AND INFORMATION MANAGEMENT (KBI)
Faculty of Business and Economics Robust forecasting of non-stationary time series C. Croux, R. Fried, I. Gijbels and K. Mahieu DEPARTMENT OF DECISION SCIENCES AND INFORMATION MANAGEMENT (KBI KBI 1018
More informationA Test of Cointegration Rank Based Title Component Analysis.
A Test of Cointegration Rank Based Title Component Analysis Author(s) Chigira, Hiroaki Citation Issue 2006-01 Date Type Technical Report Text Version publisher URL http://hdl.handle.net/10086/13683 Right
More informationROBUST TESTS BASED ON MINIMUM DENSITY POWER DIVERGENCE ESTIMATORS AND SADDLEPOINT APPROXIMATIONS
ROBUST TESTS BASED ON MINIMUM DENSITY POWER DIVERGENCE ESTIMATORS AND SADDLEPOINT APPROXIMATIONS AIDA TOMA The nonrobustness of classical tests for parametric models is a well known problem and various
More informationTIME SERIES ANALYSIS. Forecasting and Control. Wiley. Fifth Edition GWILYM M. JENKINS GEORGE E. P. BOX GREGORY C. REINSEL GRETA M.
TIME SERIES ANALYSIS Forecasting and Control Fifth Edition GEORGE E. P. BOX GWILYM M. JENKINS GREGORY C. REINSEL GRETA M. LJUNG Wiley CONTENTS PREFACE TO THE FIFTH EDITION PREFACE TO THE FOURTH EDITION
More informationRobust Standard Errors for Robust Estimators
Robust Standard Errors for Robust Estimators Christophe Croux K.U.Leuven Geert Dhaene K.U.Leuven Dirk Hoorelbeke K.U.Leuven Abstract A regression estimator is said to be robust if it is still reliable
More informationSPURIOUS AND HIDDEN VOLATILITY* M. Angeles Carnero, Daniel Peña and Esther Ruiz**
SPURIOUS AND HIDDEN VOLATILITY* M. Angeles Carnero, Daniel Peña and Esther Ruiz** WP-AD 4-4 Corresponding author: M. Angeles Carnero: Dpto. Fundamentos Análisis Económico. Universidad de Alicante, Campus
More informationSample Exam Questions for Econometrics
Sample Exam Questions for Econometrics 1 a) What is meant by marginalisation and conditioning in the process of model reduction within the dynamic modelling tradition? (30%) b) Having derived a model for
More informationStochastic Processes
Stochastic Processes Stochastic Process Non Formal Definition: Non formal: A stochastic process (random process) is the opposite of a deterministic process such as one defined by a differential equation.
More informationTeruko Takada Department of Economics, University of Illinois. Abstract
Nonparametric density estimation: A comparative study Teruko Takada Department of Economics, University of Illinois Abstract Motivated by finance applications, the objective of this paper is to assess
More informationLong Memory through Marginalization
Long Memory through Marginalization Hidden & Ignored Cross-Section Dependence Guillaume Chevillon ESSEC Business School (CREAR) and CREST joint with Alain Hecq and Sébastien Laurent Maastricht University
More informationEconometrics. Week 4. Fall Institute of Economic Studies Faculty of Social Sciences Charles University in Prague
Econometrics Week 4 Institute of Economic Studies Faculty of Social Sciences Charles University in Prague Fall 2012 1 / 23 Recommended Reading For the today Serial correlation and heteroskedasticity in
More informationDoes k-th Moment Exist?
Does k-th Moment Exist? Hitomi, K. 1 and Y. Nishiyama 2 1 Kyoto Institute of Technology, Japan 2 Institute of Economic Research, Kyoto University, Japan Email: hitomi@kit.ac.jp Keywords: Existence of moments,
More informationAn estimate of the long-run covariance matrix, Ω, is necessary to calculate asymptotic
Chapter 6 ESTIMATION OF THE LONG-RUN COVARIANCE MATRIX An estimate of the long-run covariance matrix, Ω, is necessary to calculate asymptotic standard errors for the OLS and linear IV estimators presented
More informationA radial basis function artificial neural network test for ARCH
Economics Letters 69 (000) 5 3 www.elsevier.com/ locate/ econbase A radial basis function artificial neural network test for ARCH * Andrew P. Blake, George Kapetanios National Institute of Economic and
More information