Model selection for Lévy measures in diffusion processes with jumps from discrete observations

We deal with parametric inference and selection problems for jump components in discretely observed diffusion processes with jumps. We prepare several competing parametric models for the Lévy measure that might be misspecified, and select the best model from the aspect of information criteria. We construct quasi-information criteria (QIC), which are approximations of the information criteria based on continuous observations.     

Estimation of causal continuous‐time autoregressive moving average random fields

We estimate model parameters of Lévy‐driven causal continuous‐time autoregressive moving average random fields by fitting the empirical variogram to the theoretical counterpart using a weighted least squares (WLS) approach. Subsequent to deriving asymptotic results for the variogram estimator, we show strong consistency and asymptotic normality of the parameter estimator. Furthermore, we conduct a simulation study to assess the quality of the WLS estimator for finite samples. For the simulation, we utilize numerical approximation schemes based on truncation and discretization of stochastic integrals and we analyze the associated simulation errors in detail. Finally, we apply our results to real data of the cosmic microwave background.     

Local asymptotic normality for Student-Lévy processes under high-frequency sampling

There is considerable interest in parameter estimation in Lévy models. The maximum likelihood estimator is widely used because under certain conditions it enjoys asymptotic efficiency properties. The toolkit for Lévy processes is the local asymptotic normality which guarantees these conditions. Although the likelihood function is not known explicitly, we prove local asymptotic normality for the location and scale parameters of the Student-Lévy process assuming high-frequency data. In addition, we propose a numerical method to make maximum likelihood estimates feasible based on the Monte Carlo expectation-maximization algorithm. A simulation study verifies the theoretical results.     

Expansion and estimation of Lévy process functionals in nonlinear and nonstationary time series regression

In this article, we develop a series estimation method for unknown time-inhomogeneous functionals of Lévy processes involved in econometric time series models. To obtain an asymptotic distribution for the proposed estimators, we establish a general asymptotic theory for partial sums of bivariate functionals of time and nonstationary variables. These results show that the proposed estimators in different situations converge to quite different random variables. In addition, the rates of convergence depend on various factors rather than just the sample size. Finite sample simulations are provided to evaluate the finite sample performance of the proposed model and estimation method.     

Estimation of the characteristics of a Lévy process

We consider a Lévy process that is e.g. used in finance to model stock price developments. We want to estimate the characteristics of that process, based on historical data where we assume that we have discrete, high frequency observations. We introduce a threshold estimation method and show consistency and in the case of finite activity asymptotic normality of these estimators.     

Parametric Estimation for Subordinators and Induced OU Processes

Abstract.  Consider a stationary sequence of random variables with infinitely divisible marginal law, characterized by its Lévy density. We analyse the behaviour of a so-called cumulant M-estimator, in case this Lévy density is characterized by a Euclidean (finite dimensional) parameter. Under mild conditions, we prove consistency and asymptotic normality of the estimator. The estimator is considered in the situation where the data are increments of a subordinator as well as the situation where the data consist of a discretely sampled Ornstein–Uhlenbeck (OU) process induced by the subordinator. We illustrate our results for the Gamma-process and the Inverse-Gaussian OU process. For these processes we also explain how the estimator can be computed numerically.     

On the likelihood function of small time variance Gamma Lévy processes

We investigate the likelihood function of small generalized Laplace laws and variance gamma Lévy processes in the short time framework. We prove the local asymptotic normality property in statistical inference for the variance gamma Lévy process under high-frequency sampling with its associated optimal convergence rate and Fisher information matrix. The location parameter is required to be given in advance for this purpose, while the remaining three parameters are jointly well behaved with an invertible Fisher information matrix. The results are discussed with relation to equivalent formulations of the variance gamma Lévy process, that is, as a time-changed Brownian motion and as a difference of two independent gamma processes.     

Dependence Estimation for High‐frequency Sampled Multivariate CARMA Models

The paper considers high‐frequency sampled multivariate continuous‐time autoregressive moving average (MCARMA) models and derives the asymptotic behaviour of the sample autocovariance function to a normal random matrix. Moreover, we obtain the asymptotic behaviour of the cross‐covariances between different components of the model. We will see that the limit distribution of the sample autocovariance function has a similar structure in the continuous‐time and in the discrete‐time model. As a special case, we consider a CARMA (one‐dimensional MCARMA) process. For a CARMA process, we prove Bartlett's formula for the sample autocorrelation function. Bartlett's formula has the same form in both models; only the sums in the discrete‐time model are exchanged by integrals in the continuous‐time model. Finally, we present limit results for multivariate MA processes as well, which are not known in this generality in the multivariate setting yet.     

Tempered Mittag-Leffler Lévy processes

In this article, we introduce tempered Mittag-Leffler Lévy processes (TMLLP). TMLLP is represented as tempered stable subordinator delayed by a gamma process. Its probability density function and Lévy density are obtained in terms of infinite series and Mittag-Leffler function, respectively. Asymptotic forms of the tails and moments are given. A step-by-step procedure of the parameters estimation and simulation of sample paths is given. We also provide main results available for Mittag-Leffler Lévy processes (MLLP) and some extensions which are not available in a collective way in a single article. Our results generalize and complement the results available on Mittag-Leffler distribution and MLLP in several directions. Further, the asymptotic forms of the moments of the first-exit times of the TMLLP are also discussed.     

Estimation for Non-Negative Lévy-Driven CARMA Processes

Continuous-time autoregressive moving average (CARMA) processes with a nonnegative kernel and driven by a nondecreasing Lévy process constitute a useful and very general class of stationary, nonnegative continuous-time processes that have been used, in particular, for the modeling of stochastic volatility. Brockwell, Davis, and Yang (2007) derived efficient estimates of the parameters of a nonnegative Lévy-driven CAR(1) process and showed how the realization of the underlying Lévy process can be estimated from closely-spaced observations of the process itself. In this article we show how the ideas of that article can be generalized to higher order CARMA processes with nonnegative kernel, the key idea being the decomposition of the CARMA process into a sum of dependent Ornstein–Uhlenbeck processes.     

Assessing the association between two spatial or temporal sequences

This paper deals with the codispersion coefficient for spatial and temporal series. We present some results and simulations concerning the codispersion coefficient in the context of spatial models. The results obtained are immediate consequences of the asymptotic normality of the sample codispersion coefficient and show certain limitations of the coefficient. New simulation studies provide information about the performance of the coefficient with respect to other coefficients of spatial association. The behavior of the codispersion coefficient under additively contaminated processes is also studied via Monte Carlo simulations. In the context of time series, explicit expressions for the asymptotic variance of the sample version of the coefficient are given for autoregressive and moving average processes. Resampling methods are used to compute the variance of the coefficient. A real data example is presented to explore how well the codispersion coefficient captures the comovement between two time series in practice.     

Parameter inference for time series with regular and seasonal unit roots

It is common to have both regular and seasonal roots present in many time series data. It may occur that one or both of the roots are just close but not equal to unity. Parameter inference for this situation is considered both when the time series has a finite or an infinite variance. Asymptotic char-acterizations of the test statistics were obtained via functionals of Ornstein-Uhlenbeck processes and Lévy processes. Tabulations for the large sample distributions are obtained. The results will be useful in applications deciding whether both regular and seasonal differencing are needed in fitting a time series model.     

Quantile Regression for Location‐Scale Time Series Models with Conditional Heteroscedasticity

This paper considers quantile regression for a wide class of time series models including autoregressive and moving average (ARMA) models with asymmetric generalized autoregressive conditional heteroscedasticity errors. The classical mean‐variance models are reinterpreted as conditional location‐scale models so that the quantile regression method can be naturally geared into the considered models. The consistency and asymptotic normality of the quantile regression estimator is established in location‐scale time series models under mild conditions. In the application of this result to ARMA‐generalized autoregressive conditional heteroscedasticity models, more primitive conditions are deduced to obtain the asymptotic properties. For illustration, a simulation study and a real data analysis are provided.     

Fractional Lévy stable motion time-changed by gamma subordinator

Abstract

In this paper a new stochastic process is introduced by subordinating fractional Lévy stable motion (FLSM) with gamma process. This new process incorporates stochastic volatility in the parent process FLSM. Fractional order moments, tail asymptotic, codifference and persistence of signs long-range dependence of the new process are discussed. A step-by-step procedure for simulations of sample trajectories and estimation of the parameters of the introduced process are given. Our study complements and generalizes the results available on variance-gamma process and fractional Laplace motion in various directions, which are well studied processes in literature.     

A note on non-negative continuous time processes

Summary.  Recently there has been much work on developing models that are suitable for analysing the volatility of a continuous time process. One general approach is to define a volatility process as the convolution of a kernel with a non-decreasing Lévy process, which is non-negative if the kernel is non-negative. Within the framework of time continuous autoregressive moving average (CARMA) processes, we derive a necessary and sufficient condition for the kernel to be non-negative. This condition is in terms of the Laplace transform of the CARMA kernel, which has a simple form. We discuss some useful consequences of this result and delineate the parametric region of stationarity and non-negative kernel for some lower order CARMA models.     

Asymptotic Properties of QML Estimators for VARMA Models with Time‐dependent Coefficients

This paper is about vector autoregressive‐moving average models with time‐dependent coefficients to represent non‐stationary time series. Contrary to other papers in the univariate case, the coefficients depend on time but not on the series' length n. Under appropriate assumptions, it is shown that a Gaussian quasi‐maximum likelihood estimator is almost surely consistent and asymptotically normal. The theoretical results are illustrated by means of two examples of bivariate processes. It is shown that the assumptions underlying the theoretical results apply. In the second example, the innovations are marginally heteroscedastic with a correlation ranging from ?0.8 to 0.8. In the two examples, the asymptotic information matrix is obtained in the Gaussian case. Finally, the finite‐sample behaviour is checked via a Monte Carlo simulation study for n from 25 to 400. The results confirm the validity of the asymptotic properties even for short series and the asymptotic information matrix deduced from the theory.     

On Some Applications of a Symbolic Representation of Non Centered Lévy Processes

By using a symbolic technique known in the literature as the classical umbral calculus, we characterize two classes of polynomials related to Lévy processes: the Kailath-Segall and the time-space harmonic polynomials. We provide the Kailath-Segall formula in terms of cumulants and we recover simple closed-forms for several families of polynomials with respect to not centered Lévy processes, such as the Hermite polynomials with Brownian motion, Poisson-Charlier polynomials with Poisson processes, actuarial polynomials with Gamma processes, first kind Meixner polynomials with Pascal processes, and Bernoulli, Euler, and Krawtchuk polynomials with suitable random walks.