Estimation of a tail index based on minimum density power divergence

In this paper, we consider the minimum density power divergence estimator for the tail index of heavy tailed distributions in strong mixing processes. It is shown that the estimator is consistent and asymptotically normal under regularity conditions. The simulation results demonstrate that the estimator is robust in the presence of outliers.     

Kernel density estimation for spatial processes: the L1 theory

The purpose of this paper is to investigate kernel density estimators for spatial processes with linear or nonlinear structures. Sufficient conditions for such estimators to converge in L₁ are obtained under extremely general, verifiable conditions. The results hold for mixing as well as for nonmixing processes. Potential applications include testing for spatial interaction, the spatial analysis of causality structures, the definition of leading/lagging sites, the construction of clusters of comoving sites, etc.     

On the chernoff-savage theorem for dependent sequences

Summary Given a sequence of ϕ-mixing random variables not necessarily stationary, a Chernoff-Savage theorem for two-sample linear rank statistics is proved using the Pyke-Shorack [5] approach based on weak convergence properties of empirical processes in an extended metric. This result is a generalization of Fears and Mehra [4] in that the stationarity is not required and that the condition imposed on the mixing numbers is substantially relaxed. A similar result is shown to hold for strong mixing sequences under slightly stronger conditions on the mixing numbers. Research partially supported by the National Research Council of Canada under Grant No. A-3954.     

Weak convergence of non-stationary multivariate marked processes with applications to martingale testing

This paper establishes the weak convergence of a class of marked empirical processes of possibly non-stationary and/or non-ergodic multivariate time series sequences under martingale conditions. The assumptions involved are similar to those in Brown's martingale central limit theorem. In particular, no mixing conditions are imposed. As an application, we propose a test statistic for the martingale hypothesis and we derive its asymptotic null distribution. Finally, a Monte Carlo study shows that the asymptotic results provide good approximations for small and moderate sample sizes. An application to the S&P 500 is also considered.     

Sequences of elliptical distributions and mixtures of normal distributions

Two conditions are shown under which elliptical distributions are scale mixtures of normal distributions with respect to probability distributions. The issue of finding the mixing distribution function is also considered. As a unified theoretical framework, it is also shown that any scale mixture of normal distributions is always a term of a sequence of elliptical distributions, increasing in dimension, and that all the terms of this sequence are also scale mixtures of normal distributions sharing the same mixing distribution function. Some examples are shown as applications of these concepts, showing the way of finding the mixing distribution function.     

Estimation for stochastic damping hamiltonian systems under partial observation—I. Invariant density

In this paper, we study the non-parametric estimation of the invariant density of some ergodic hamiltonian systems, using kernel estimators. The main result is a central limit theorem for such estimators under partial observation (only the positions are observed). The main tools are mixing estimates and refined covariance inequalities, the main difficulty being the strong degeneracy of such processes. This is the first paper of a series of at least two, devoted to the estimation of the characteristics of such processes: invariant density, drift term, volatility.     

Testing independence of variates in an infinitely divisible random vector

A method is given for testing the independence of variates in an infinitely divisible random vector and for testing the independence of several subsets of the variates. Applications to stochastic processes are indicated.     

A note on testing hypotheses for stationary processes in the frequency domain

In a recent paper, Eichler (2008) [11] considered a class of non- and semiparametric hypotheses in multivariate stationary processes, which are characterized by a functional of the spectral density matrix. The corresponding statistics are obtained using kernel estimates for the spectral distribution and are asymptotically normally distributed under the null hypothesis and local alternatives. In this paper, we derive the asymptotic properties of these test statistics under fixed alternatives. In particular, we also show weak convergence but with a different rate compared to the null hypothesis. We also discuss potential statistical applications of the asymptotic theory by means of a small simulation study.     

A note on weak convergence of mean residual life of stationary mixing random variables

For a sequence of strictly stationary uniform or strong mixing we estimate the mean residual time of the marginal distribution from the first n observations. Under appropriate conditions it is shown that the estimate converges weakly to a well-defined Gaussian process even when the sample size is random.     

Some geometric mixed integer-valued autoregressive (INAR) models

In this paper, we introduce some mixed integer-valued autoregressive models of orders 1 and 2 with geometric marginal distributions, denoted by MGINAR(1) and MGINAR(2), using a mixture of the well-known binomial and the negative binomial thinning. The distributions of the innovation processes are derived and several properties of the model are discussed. Conditional least squares and Yule-Walker estimators are obtained, and some numerical results of the estimations are presented. A real-life data example is investigated to assess the performance of the models.     

Cramér–Karhunen–Loève representation and harmonic principal component analysis of functional time series

We develop a doubly spectral representation of a stationary functional time series, and study the properties of its empirical version. The representation decomposes the time series into an integral of uncorrelated frequency components (Cramér representation), each of which is in turn expanded in a Karhunen–Loève series. The construction is based on the spectral density operator, the functional analogue of the spectral density matrix, whose eigenvalues and eigenfunctions at different frequencies provide the building blocks of the representation. By truncating the representation at a finite level, we obtain a harmonic principal component analysis of the time series, an optimal finite dimensional reduction of the time series that captures both the temporal dynamics of the process, as well as the within-curve dynamics. Empirical versions of the decompositions are introduced, and a rigorous analysis of their large-sample behaviour is provided, that does not require any prior structural assumptions such as linearity or Gaussianity of the functional time series, but rather hinges on Brillinger-type mixing conditions involving cumulants.     

On the mode-change problem for random measures

The classical change-point problem in modern terms, i.e., the mode-change problem, is stated for sufficiently general set-indexed random processes, namely for random measures. A method is shown for solving this problem both in the general form and for the intensity of compound Poisson random measures. The results obtained are novel for the change-point problem, too.     

Estimating structural VARMA models with uncorrelated but non-independent error terms

The asymptotic properties of the quasi-maximum likelihood estimator (QMLE) of vector autoregressive moving-average (VARMA) models are derived under the assumption that the errors are uncorrelated but not necessarily independent nor martingale differences. Relaxing the martingale difference assumption on the errors considerably extends the range of application of the VARMA models, and allows one to cover linear representations of general nonlinear processes. Conditions are given for the asymptotic normality of the QMLE. Particular attention is given to the estimation of the asymptotic variance matrix, which may be very different from that obtained in the standard framework.     

A nonparametric regression cross spectrum for multivariate time series

We consider dependence structures in multivariate time series that are characterized by deterministic trends. Results from spectral analysis for stationary processes are extended to deterministic trend functions. A regression cross covariance and spectrum are defined. Estimation of these quantities is based on wavelet thresholding. The method is illustrated by a simulated example and a three-dimensional time series consisting of ECG, blood pressure and cardiac stroke volume measurements.     

Generation of the Exact Distribution and Simulation of Matched Nucleotide Sequences on a Phylogenetic Tree

Nucleotide sequences are often generated by Monte Carlo simulations to address complex evolutionary or analytic questions but the simulations are rarely described in sufficient detail to allow the research to be replicated. Here we briefly review the Markov processes of substitution in a pair of matching (homologous) nucleotide sequences and then extend it to k matching nucleotide sequences. We describe calculation of the joint distribution of nucleotides of two matching sequences. Based on this distribution, we give a method for simulation of the divergence matrix for n sites using the multinomial distribution. This is then extended to the joint distribution for k nucleotide sequences and the corresponding 4 ^k divergence array, generalizing Felsenstein (Journal of Molecular Evolution, 17, 368–376, 1981), who considered stationary, homogeneous and reversible processes on trees. We give a second method to generate matched sequences that begins with a random ancestral sequence and applies a continuous Markov process to each nucleotide site as in Rambaut and Grassly (Computer Applications in the Biosciences, 13, 235–238, 1997); further, we relate this to an equivalent approach based on an embedded Markov chain. Finally, we describe an approximate method that was recently implemented in a program developed by Jermiin et al. (Applied Bioinformatics, 2, 159–163, 2003). The three methods presented here cater for different computational and mathematical limitations and are shown in an example to produce results close to those expected on theoretical grounds. All methods are implemented using functions in the S-plus or R languages.     

Quasi-maximum likelihood estimation in GARCH processes when some coefficients are equal to zero

The asymptotic distribution of the quasi-maximum likelihood (QML) estimator is established for generalized autoregressive conditional heteroskedastic (GARCH) processes, when the true parameter may have zero coefficients. This asymptotic distribution is the projection of a normal vector distribution onto a convex cone. The results are derived under mild conditions. For an important subclass of models, no moment condition is imposed on the GARCH process. The main practical implication of these results concerns the estimation of overidentified GARCH models.     

Asymmetric GARCH processes featuring both threshold effect and bilinear structure

A class of asymmetric GARCH models is proposed by combining threshold effect and bilinear structure. The class is referred to as threshold-bilinear GARCH processes. A simulation study demonstrates that the class exhibits diverse asymmetries in volatilities, accommodating existing asymmetric models. Stationarity and existence of moments are discussed. Applications to Korean stock prices are illustrated.     

Moderate deviations for quadratic forms in Gaussian stationary processes

Moderate deviations limit theorem is proved for quadratic forms in zero-mean Gaussian stationary processes. Two particular cases are the cumulative periodogram and the kernel spectral density estimator. We also derive the exponential decay of moderate deviation probabilities of goodness-of-fit tests for the spectral density and then discuss intermediate asymptotic efficiencies of tests.     

State space models on special manifolds

This paper concerns modeling time series observations in state space forms considered on the Stiefel and Grassmann manifolds. We develop a state space model relating the time series observations to a sequence of unobserved state or parameter matrices assuming the matrix Langevin noise processes on the Stiefel manifolds. We show a Bayes method for estimating the state matrices by the posterior modes. We consider a further extended state space model where two sequences of unobserved state matrices are involved. A simple state space model on the Grassmann manifolds with matrix Langevin noise processes is also investigated.     

On same-realization prediction in an infinite-order autoregressive process

Let observations come from an infinite-order autoregressive (AR) process. For predicting the future of the observed time series (referred to as the same-realization prediction), we use the least-squares predictor obtained by fitting a finite-order AR model. We also allow the order to become infinite as the number of observations does in order to obtain a better approximation. Moment bounds for the inverse sample covariance matrix with an increasing dimension are established under various conditions. We then apply these results to obtain an asymptotic expression for the mean-squared prediction error of the least-squares predictor in same-realization and increasing-order settings. The second-order term of this expression is the sum of two terms which measure both the goodness of fit and model complexity. It forms the foundation for a companion paper by Ing and Wei (Order selection for same-realization predictions in autoregressive processes, Technical report C-00-09, Institute of Statistical Science, Academia Sinica, Taipei, Taiwan, ROC, 2000) which provides the first theoretical verification that AIC is asymptotically efficient for same-realization predictions. Finally, some comparisons between the least-squares predictor and the ridge regression predictor are also given.