Estimation de la matrice de variance asymptotique des estimateurs du QMV de modèles ARMA faibles multivariés

In this Note, we consider the problems of estimating the asymptotic variance of the quasi-maximum likelihood estimator (QMLE) of vector autoregressive moving-average (VARMA) models under the assumption that the errors are uncorrelated but not necessarily independent (i.e. weak VARMA). We first give expressions for the derivatives of the VARMA residuals in terms of the parameters of the models. Secondly we give an explicit expression of the asymptotic variance of the QMLE, in terms of the VAR and MA polynomials, and of the second- and fourth-order structure of the noise. We deduce a consistent estimator of the asymptotic variance of the QMLE.     

Tests portmanteau multivariés d'adéquation de modèles VARMA faibles

In this Note, we consider portmanteau tests for testing the adequacy of vector autoregressive moving-average (VARMA) models under the assumption that the errors are uncorrelated but not necessarily independent. We relax the standard independence assumption to extend the range of application of the VARMA models, allowing us to treat linear representations of general nonlinear processes. We first study the joint distribution of the quasi-maximum likelihood estimator (QMLE) and the noise empirical autocovariances. We thus obtain the asymptotic distribution of residual empirical autocovariances and autocorrelations under weak assumptions on the noise. We deduce the asymptotic distribution of the Ljung–Box (or Box–Pierce) portmanteau statistics for VARMA models with nonindependent innovations. We propose a method to adjust the critical values of the portmanteau tests.     

Modified Schwarz and Hannan–Quinn information criteria for weak VARMA models

Numerous multivariate time series admit weak vector autoregressive moving-average (VARMA) representations, in which the errors are uncorrelated but not necessarily independent nor martingale differences. These models are called weak VARMA by opposition to the standard VARMA models, also called strong VARMA models, in which the error terms are supposed to be independent and identically distributed (iid). This article considers the problem of order selection of the weak VARMA models by using the information criteria. It is shown that the use of the standard information criteria are often not justified when the iid assumption on the noise is relaxed. As a consequence, we propose the modified versions of the Schwarz or Bayesian information criterion and of the Hannan and Quinn criterion for identifying the orders of weak VARMA models. Monte Carlo experiments show that the proposed modified criteria estimate the model orders more accurately than the standard ones. An illustrative application using the squared daily returns of financial series is presented.     

Estimation des ordres de modèles ARMA faibles multivariés

In this Note, we consider the problem of order selection of vector autoregressive moving-average (VARMA) models under the assumption that the errors are uncorrelated, but not necessarily independent. These models are called weak VARMA by opposition to the standard VARMA models, also called strong VARMA models, in which the error terms are supposed to be iid. This selection is based on minimizing an information criterion, especially that introduced by Akaike. The theoretical foundations of the Akaike information criterion (AIC) are not more established when the iid assumption on the noise is relaxed. We propose a modified AIC criterion, and which may be very different from the standard AIC criterion.     

Diagnostic checking of multivariate nonlinear time series models with martingale difference errors

In this article, we derive the asymptotic distribution of residual autocovariance and autocorrelation matrices for a general class of multivariate nonlinear time series models by assuming only that the error term is a martingale difference sequence. Two types of applications are developed: global test statistics of the portmanteau type and one-lag test statistics, which describe the residual correlation at individual lags. To illustrate the proposed methodology, simulation results are reported for diagnosing multivariate threshold time series models. The following test statistics are compared: the classical test statistics presuming independent errors and the proposed methodology which supposes only martingale difference errors.     

Multivariate Local Polynomial Fitting for Martingale Nonlinear Regression Models

Local polynomial modelling is a useful tool for nonlinear time series analysis. For nonlinear regression models with martingale difference errors, this paper presents a simple proof of local linear and local quadratic fittings under apparently minimal short-range dependence condition. Explicit formulae for the asymptotic bias and asymptotic variance are given, which facilitate numerical evaluations of these important quantities. The general theory is applied to nonparametric partial derivative estimation in nonlinear time series. A bias-adjusted method for constructing confidence intervals for first-order partial derivatives is described. Two examples, including the sunspots data, are used to demonstrate the use of local quadratic fitting for modelling and characterizing nonlinearity in time series data.     

Affine-invariant aligned rank tests for the multivariate general linear model with VARMA errors

We develop optimal rank-based procedures for testing affine-invariant linear hypotheses on the parameters of a multivariate general linear model with elliptical VARMA errors. We propose a class of optimal procedures that are based either on residual (pseudo-)Mahalanobis signs and ranks, or on absolute interdirections and lift-interdirection ranks, i.e., on hyperplane-based signs and ranks. The Mahalanobis versions of these procedures are strictly affine-invariant, while the hyperplane-based ones are asymptotically affine-invariant. Both versions generalize the univariate signed rank procedures proposed by Hallin and Puri (J. Multivar. Anal. 50 (1994) 175), and are locally asymptotically most stringent under correctly specified radial densities. Their AREs with respect to Gaussian procedures are shown to be convex linear combinations of the AREs obtained in Hallin and Paindaveine (Ann. Statist. 30 (2002) 1103; Bernoulli 8 (2002) 787) for the pure location and purely serial models, respectively. The resulting test statistics are provided under closed form for several important particular cases, including multivariate Durbin-Watson tests, VARMA order identification tests, etc. The key technical result is a multivariate asymptotic linearity result proved in Hallin and Paindaveine (Asymptotic linearity of serial and nonserial multivariate signed rank statistics, submitted).     

Asymptotic inference of unstable periodic ARCH processes

This paper studies asymptotic properties of the quasi maximum likelihood and weighted least squares estimates (QMLE and WLSE) of the conditional variance slope parameters of a strictly unstable ARCH model with periodically time varying coefficients (PARCH in short). The model is strictly unstable in the sense that its parameters lie outside the strict periodic stationarity domain and its boundary. Obtained from the regression form of the PARCH, the WLSE is a variant of the least squares method weighted by the square of the conditional variance evaluated at any fixed value in the parameter space. In calculating the QMLE and WLSE, the conditional variance intercepts are set to any arbitrary values not necessarily the true ones. The theoretical finding is that the QMLE and WLSE are consistent and asymptotically Gaussian with the same asymptotic variance irrespective of the fixed conditional variance intercepts and the weighting parameters. So because of its numerical complexity, the QMLE may be dropped in favor of the WLSE which enjoys closed form.     

ASYMPTOTIC NORMALITY OF QUASI MAXIMUM LIKELIHOOD ESTIMATE IN GENERALIZED LINEAR MODELS

For the Generalized Linear Model (GLM), under some conditions including that the specification of the expectation is correct, it is shown that the Quasi Maximum Likelihood Estimate (QMLE) of the parameter-vector is asymptotic normal. It is also shown that the asymptotic covariance matrix of the QMLE reaches its minimum (in the positive-definte sense) in case that the specification of the covariance matrix is correct.     

Multistage weighted least squares estimation of ARCH processes in the stable and unstable cases

For an ARCH model, we propose a multistage weighted least squares (WLS) estimate which consists of repeated WLS procedures until the corresponding asymptotic variance equals that of the quasi-maximum likelihood estimate (QMLE). At every stage, the current estimate is of a WLS type weighted by the squared conditional variance evaluated at the estimate of the previous stage. Initially, the weighting parameter is any fixed and known value in the parameter space. The procedure provides, without any moment requirement, an asymptotically Gaussian estimate having the same asymptotic distribution as the QMLE even in the unstable case. Apart from the initialization stage, two additional stages are required in the stable case to obtain the same asymptotic distribution as the QMLE, while in the unstable case only one stage is enough. So in all, the proposed procedure involves three stages WLS in the stable case and two stages WLS in the unstable case.     

On consistent testing for serial correlation of unknown form in vector time series models

Multivariate autoregressive models with exogenous variables (VARX) are often used in econometric applications. Many properties of the basic statistics for this class of models rely on the assumption of independent errors. Using results of Hong (Econometrica 64 (1996) 837), we propose a new test statistic for checking the hypothesis of non-correlation or independence in the Gaussian case. The test statistic is obtained by comparing the spectral density of the errors under the null hypothesis of independence with a kernel-based spectral density estimator. The asymptotic distribution of the statistic is derived under the null hypothesis. This test generalizes the portmanteau test of Hosking (J. Amer. Statist. Assoc. 75 (1980) 602). The consistency of the test is established for a general class of static regression models with autocorrelated errors. Its asymptotic slope is derived and the asymptotic relative efficiency within the class of possible kernels is also investigated. Finally, the level and power of the resulting tests are also studied by simulation.     

Inference for 2-D GARCH models

The purpose of this paper is, in the first step, to consider a class of GMM estimators with interesting asymptotic properties and a reasonable number of computations for two dimensionally indexed Generalized Autoregressive Conditional Heteroscedasticity (GARCH) model. In the second step, we use the central limit theorem of Huang (1992) for spatial martingale differences to establish the LAN property for general two-dimensional discrete models on a regular grid with Gaussian errors. We then apply this result to the spatial GARCH model and derive the limit distribution of the maximum likelihood estimators of the parameters. Results of numerical simulations are presented.     

On asymptotic theory for multivariate GARCH models

The paper investigates the asymptotic theory for a multivariate GARCH model in its general vector specification proposed by Bollerslev, Engle and Wooldridge (1988) [4], known as the VEC model. This model includes as important special cases the so-called BEKK model and many versions of factor GARCH models, which are often used in practice. We provide sufficient conditions for strict stationarity and geometric ergodicity. The strong consistency of the quasi-maximum likelihood estimator (QMLE) is proved under mild regularity conditions which allow the process to be integrated. In order to obtain asymptotic normality, the existence of sixth-order moments of the process is assumed.     

一类Kaplan-Meier估计泛函的一些大样本性质

设F,G分别表示某寿命随机变量与删失随机变量的分布函数,在不假定F、G连续的情况下该文使用点过程鞅方法证明了Kaplan－Meier估计的一类泛函的渐近正态性,并建立了一个均方误差不等式和一个概率不等式．     

Asymptotic asset pricing and bubbles

We define the concept of asymptotic superreplication, and prove a duality principle of asset pricing for sequences of financial markets (e.g., weakly converging financial markets and large financial markets) based on contiguous sequences of equivalent local martingale measures. This provides a pricing mechanism to calculate the fundamental value of a financial asset in the asymptotic market. We introduce the notion of asymptotic bubbles by showing that this fundamental value can be strictly lower than the current price of the asset. In the case of weakly converging markets, we show that this fundamental value is equal to an expectation of the terminal value of the asset in the weak-limit market. From a practical perspective, we relate the asymptotic superreplication price to a limit of quantile-hedging prices. This shows that even when a price process is a true martingale, it can have properties similar to a bubble, up to a set of small probability. For practical applications, we give examples of weakly converging discrete-time models (e.g. some GARCH models) and large financial models that present bubbles.     

The quantile process under random censoring

In this paper we discuss the asymptotic properties of quantile processes under random censoring. In contrast to most work in this area we prove weak convergence of an appropriately standardized quantile process under the assumption that the quantile regression model is only linear in the region, where the process is investigated. Additionally, we also discuss properties of the quantile process in sparse regression models including quantile processes obtained from the Lasso and adaptive Lasso. The results are derived by a combination of modern empirical process theory, classical martingale methods and a recent result of Kato (2009).     

Embedding and asymptotic expansions for martingales

The paper develops a way of embedding general martingales in continuous ones in such a way that the quadratic variation of the continuous martingale has conditional cumulants (given the original martingale) that are explicitly given in terms of optional and predictable variations of the original process. Bartlett identities for the conditional cumulants are also found. A main corollary to these results is the establishment of second (and in some cases higher) order asymptotic expansions for martingales.Research supported in part by National Science Foundation grant DMS 93-05601 and Army Research Office grant DAAH04-1-0105     

Detecting and modeling nonlinearity in the gas furnace data

Using a modification of the Hinich, J Time Ser Anal 3(3):169–176, (1982) bispectrum test for nonlinearity and Gaussianity, the residuals of the Tiao and Box, J Am Stat Assoc 76:802–816, (1981) constrained and unconstrained VAR models for the gas furnace data reject the assumption of Gaussianity and linearity over a grid of bandwidths for estimating the bispectrum. These findings call into question the specification of the linear VAR and VARMA models assumed by Tiao and Box, J Am Stat Assoc 76:802–816, (1981). Utilizing the alternative Hinich J Nonparametr Stat 6:205–221, (1996) nonlinearity test, the residuals of the VAR model were shown to exhibit episodic nonlinearity. The sensitivity of the findings to outliers is investigated by estimating and testing the residuals of L1 and MINIMAX models from 1–6 lags. Building on the linear dynamic specification, a multivariate adaptive regression splines (MARS) model is estimated, using two software implementations, and shown to remove the nonlinearity in the residuals. Leverage plots were used to illustrate the “cost” of imposing a linearity assumption. Out-of-sample forecasting tests from 1–6 periods ahead found that using the sum-of-squared errors criteria, the MARS model out performed ACE, GAM and projection pursuit models.     

一般半相依回归系统的协方差改进估计

本文讨论了由两个等阶的回归方程组成的半相依系统,运用协方差改进法获得了参数的一个迭代估计序列,并证明了它的协方差阵已知时,处处收敛到最佳线性无偏估计,同时其协方差阵在矩阵偏序意义下单调性,并且给出了当迭代次数亦趋于无穷时,保证其具有相合性的一个条件。     

Testing the equality of error distributions from k independent GARCH models

In this paper we study the problem of testing the null hypothesis that errors from $k$ independent parametrically specified generalized autoregressive conditional heteroskedasticity (GARCH) models have the same distribution versus a general alternative. First we establish the asymptotic validity of a class of linear test statistics derived from the $k$ residual-based empirical distribution functions. A distinctive feature is that the asymptotic distribution of the test statistics involves terms depending on the distributions of errors and the parameters of the models, and weight functions providing the flexibility to choose scores for investigating power performance. A Monte Carlo study assesses the asymptotic performance in terms of empirical size and power of the three-sample test based on the Wilcoxon and Van der Waerden score generating functions in finite samples. The results demonstrate that the two proposed tests have overall reasonable size and their power is particularly high when the assumption of Gaussian errors is violated. As an illustrative example, the tests are applied to daily individual stock returns of the New York Stock Exchange data.