On functional central limit theorems for dependent,heterogeneous arrays with applications to tail index and tail dependence estimation

We establish invariance principles for a large class of dependent, heterogeneous arrays. The theory equally covers conventional arrays, and inherently degenerate tail arrays popularly encountered in the extreme value theory literature including sample means and covariances of tail events and exceedances. For tail arrays we trim dependence assumptions down to a minimum leaving non-extremes and joint distributions unrestricted, covering geometrically ergodic, mixing, and mixingale processes, in particular linear and nonlinear distributed lags with long or short memory, linear and nonlinear GARCH, and stochastic volatility.     

Monte Carlo Maximum Likelihood Estimation for Generalized Long-Memory Time Series Models

An exact maximum likelihood method is developed for the estimation of parameters in a non-Gaussian nonlinear density function that depends on a latent Gaussian dynamic process with long-memory properties. Our method relies on the method of importance sampling and on a linear Gaussian approximating model from which the latent process can be simulated. Given the presence of a latent long-memory process, we require a modification of the importance sampling technique. In particular, the long-memory process needs to be approximated by a finite dynamic linear process. Two possible approximations are discussed and are compared with each other. We show that an autoregression obtained from minimizing mean squared prediction errors leads to an effective and feasible method. In our empirical study, we analyze ten daily log-return series from the S&P 500 stock index by univariate and multivariate long-memory stochastic volatility models. We compare the in-sample and out-of-sample performance of a number of models within the class of long-memory stochastic volatility models.     

A new class of models for bivariate joint tails

Summary.  A fundamental issue in applied multivariate extreme value analysis is modelling dependence within joint tail regions. The primary focus of this work is to extend the classical pseudopolar treatment of multivariate extremes to develop an asymptotically motivated representation of extremal dependence that also encompasses asymptotic independence. Starting with the usual mild bivariate regular variation assumptions that underpin the coefficient of tail dependence as a measure of extremal dependence, our main result is a characterization of the limiting structure of the joint survivor function in terms of an essentially arbitrary non-negative measure that must satisfy some mild constraints. We then construct parametric models from this new class and study in detail one example that accommodates asymptotic dependence, asymptotic independence and asymmetry within a straightforward parsimonious parameterization. We provide a fast simulation algorithm for this example and detail likelihood-based inference including tests for asymptotic dependence and symmetry which are useful for submodel selection. We illustrate this model by application to both simulated and real data. In contrast with the classical multivariate extreme value approach, which concentrates on the limiting distribution of normalized componentwise maxima, our framework focuses directly on the structure of the limiting joint survivor function and provides significant extensions of both the theoretical and the practical tools that are available for joint tail modelling.     

Stochastic volatility modelling in continuous time with general marginal distributions: Inference,prediction and model selection

We compare results for stochastic volatility models where the underlying volatility process having generalized inverse Gaussian (GIG) and tempered stable marginal laws. We use a continuous time stochastic volatility model where the volatility follows an Ornstein–Uhlenbeck stochastic differential equation driven by a Lévy process. A model for long-range dependence is also considered, its merit and practical relevance discussed. We find that the full GIG and a special case, the inverse gamma, marginal distributions accurately fit real data. Inference is carried out in a Bayesian framework, with computation using Markov chain Monte Carlo (MCMC). We develop an MCMC algorithm that can be used for a general marginal model.     

A Class of Non-Gaussian State Space Models With Exact Likelihood Inference

The likelihood function of a general nonlinear, non-Gaussian state space model is a high-dimensional integral with no closed-form solution. In this article, I show how to calculate the likelihood function exactly for a large class of non-Gaussian state space models that include stochastic intensity, stochastic volatility, and stochastic duration models among others. The state variables in this class follow a nonnegative stochastic process that is popular in econometrics for modeling volatility and intensities. In addition to calculating the likelihood, I also show how to perform filtering and smoothing to estimate the latent variables in the model. The procedures in this article can be used for either Bayesian or frequentist estimation of the model’s unknown parameters as well as the latent state variables. Supplementary materials for this article are available online.     

Estimating and Testing Nonlinear Local Dependence Between Two Time Series

ABSTRACT

The most common measure of dependence between two time series is the cross-correlation function. This measure gives a complete characterization of dependence for two linear and jointly Gaussian time series, but it often fails for nonlinear and non-Gaussian time series models, such as the ARCH-type models used in finance. The cross-correlation function is a global measure of dependence. In this article, we apply to bivariate time series the nonlinear local measure of dependence called local Gaussian correlation. It generally works well also for nonlinear models, and it can distinguish between positive and negative local dependence. We construct confidence intervals for the local Gaussian correlation and develop a test based on this measure of dependence. Asymptotic properties are derived for the parameter estimates, for the test functional and for a block bootstrap procedure. For both simulated and financial index data, we construct confidence intervals and we compare the proposed test with one based on the ordinary correlation and with one based on the Brownian distance correlation. Financial indexes are examined over a long time period and their local joint behavior, including tail behavior, is analyzed prior to, during and after the financial crisis. Supplementary material for this article is available online.     

Asymptotically Unbiased Estimation of the Coefficient of Tail Dependence

Abstract. We introduce and study a class of weighted functional estimators for the coefficient of tail dependence in bivariate extreme value statistics. Asymptotic normality of these estimators is established under a second‐order condition on the joint tail behaviour, some conditions on the weight function and for appropriately chosen sequences of intermediate order statistics. Asymptotically unbiased estimators are constructed by judiciously chosen linear combinations of weighted functional estimators, and variance optimality within this class of asymptotically unbiased estimators is discussed. The finite sample performance of some specific examples from our class of estimators and some alternatives from the recent literature are evaluated with a small simulation experiment.     

Modelling Dependence within Joint Tail Regions

Standard approaches for modelling dependence within joint tail regions are based on extreme value methods which assume max-stability, a particular form of joint tail dependence. We develop joint tail models based on a broader class of dependence structure which provides a natural link between max-stable models and weaker forms of dependence including independence and negative association. This approach overcomes many of the problems that are encountered with standard methods and is the basis for a Poisson process representation that generalizes existing bivariate results. We apply the new techniques to simulated and environmental data, and demonstrate the marked advantage that the new approach offers for joint tail extrapolation.     

Bayesian analysis of stochastic volatility models with flexible tails

An alternative distributional assumption is proposed for the stochastic volatility model. This results in extremely flexible tail behaviour of the sampling distribution for the observables, as well as in the availability of a simple Markov Chain Monte Carlo strategy for posterior analysis. By allowing the tail behaviour to be determined by a separate parameter, we reserve the parameters of the volatility process to dictate the degree of volatility clustering. Treatment of a mean function is formally integrated in the analysis.

Some empirical examples on both stock prices and exchange rates clearly indicate the presence of fat tails, in combination with high levels of volatility clustering. In addition, predictive distributions indicate a good fit with these typical financial data sets.     

Volatility model selection for extremes of financial time series

Although both widely used in the financial industry, there is quite often very little justification why GARCH or stochastic volatility is preferred over the other in practice. Most of the relevant literature focuses on the comparison of the fit of various volatility models to a particular data set, which sometimes may be inconclusive due to the statistical similarities of both processes. With an ever growing interest among the financial industry in the risk of extreme price movements, it is natural to consider the selection between both models from an extreme value perspective. By studying the dependence structure of the extreme values of a given series, we are able to clearly distinguish GARCH and stochastic volatility models and to test statistically which one better captures the observed tail behaviour. We illustrate the performance of the method using some stock market returns and find that different volatility models may give a better fit to the upper or lower tails.     

Integrated OU Processes and Non-Gaussian OU-based Stochastic Volatility Models

Abstract. In this paper, we study the detailed distributional properties of integrated non-Gaussian Ornstein–Uhlenbeck (intOU) processes. Both exact and approximate results are given. We emphasize the study of the tail behaviour of the intOU process. Our results have many potential applications in financial economics, as OU processes are used as models of instantaneous variance in stochastic volatility (SV) models. In this case, an intOU process can be regarded as a model of integrated variance. Hence, the tail behaviour of the intOU process will determine the tail behaviour of returns generated by SV models.     

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.     

Central Limit Theorems of Local Polynomial Threshold Estimator for Diffusion Processes with Jumps

Central limit theorems play an important role in the study of statistical inference for stochastic processes. However, when the non‐parametric local polynomial threshold estimator, especially local linear case, is employed to estimate the diffusion coefficients of diffusion processes, the adaptive and predictable structure of the estimator conditionally on the σ ‐field generated by diffusion processes is destroyed, so the classical central limit theorem for martingale difference sequences cannot work. In high‐frequency data, we proved the central limit theorems of local polynomial threshold estimators for the volatility function in diffusion processes with jumps by Jacod's stable convergence theorem. We believe that our proof procedure for local polynomial threshold estimators provides a new method in this field, especially in the local linear case.     

Option pricing for a stochastic volatility Lévy model with stochastic interest rates

An alternative option pricing model under a forward measure is proposed, in which asset prices follow a stochastic volatility Lévy model with stochastic interest rate. The stochastic interest rate is driven by the Hull–White process. By using an approximate method, we find a formulation for the European option in term of the characteristic function of the tail probabilities.     

Bayesian analysis of stochastic volatility models with flexible tails

An alternative distributional assumption is proposed for the stochastic volatility model. This results in extremely flexible tail behaviour of the sampling distribution for the observables, as well as in the availability of a simple Markov Chain Monte Carlo strategy for posterior analysis. By allowing the tail behaviour to be determined by a separate parameter, we reserve the parameters of the volatility process to dictate the degree of volatility clustering. Treatment of a mean function is formally integrated in the analysis.

Some empirical examples on both stock prices and exchange rates clearly indicate the presence of fat tails, in combination with high levels of volatility clustering. In addition, predictive distributions indicate a good fit with these typical financial data sets.     

The generalized linear chirp process

The linear chirp process is an important class of time series for which the instantaneous frequency changes linearly in time. Linear chirps have been used extensively to model a variety of physical signals such as radar, sonar, and whale clicks (see 1, 5 and 6). We introduce the stochastic linear chirp model and then define the generalized linear chirp (GLC) process as a special case of the G-stationary process studied by Jiang et al. (2006) to model data with time-varying frequencies. We then define GLC(p,q) processes and show that the relationship between stochastic linear chirp processes and GLC(p,q) processes is analogous to that between harmonic and ARMA models. The new methods are then applied to both simulated and actual data sets.     

Empirical Performance and Asset Pricing in Hidden Markov Models

Abstract

To improve the empirical performance of the Black-Scholes model, many alternative models have been proposed to address leptokurtic feature, volatility smile, and volatility clustering effects of the asset return distributions. However, analytical tractability remains a problem for most alternative models. In this article, we study a class of hidden Markov models including Markov switching models and stochastic volatility models, that can incorporate leptokurtic feature, volatility clustering effects, as well as provide analytical solutions to option pricing. We show that these models can generate long memory phenomena when the transition probabilities depend on the time scale. We also provide an explicit analytic formula for the arbitrage-free price of the European options under these models. The issues of statistical estimation and errors in option pricing are also discussed in the Markov switching models.     

On conditional variance estimation in nonparametric regression

In this paper we consider a nonparametric regression model in which the conditional variance function is assumed to vary smoothly with the predictor. We offer an easily implemented and fully Bayesian approach that involves the Markov chain Monte Carlo sampling of standard distributions. This method is based on a technique utilized by Kim, Shephard, and Chib (in Rev. Econ. Stud. 65:361–393, 1998) for the stochastic volatility model. Although the (parametric or nonparametric) heteroscedastic regression and stochastic volatility models are quite different, they share the same structure as far as the estimation of the conditional variance function is concerned, a point that has been previously overlooked. Our method can be employed in the frequentist context and in Bayesian models more general than those considered in this paper. Illustrations of the method are provided.     

Semi‐Parametric Models for the Multivariate Tail Dependence Function – the Asymptotically Dependent Case

Abstract. In general, the risk of joint extreme outcomes in financial markets can be expressed as a function of the tail dependence function of a high‐dimensional vector after standardizing marginals. Hence, it is of importance to model and estimate tail dependence functions. Even for moderate dimension, non‐parametrically estimating a tail dependence function is very inefficient and fitting a parametric model to tail dependence functions is not robust. In this paper, we propose a semi‐parametric model for (asymptotically dependent) tail dependence functions via an elliptical copula. Under this model assumption, we propose a novel estimator for the tail dependence function, which proves favourable compared to the empirical tail dependence function estimator, both theoretically and empirically.     

Realized Volatility and Long Memory: An Overview

The challenge of modeling, estimating, testing, and forecasting financial volatility is both intellectually worthwhile and also central to the successful analysis of financial returns and optimal investment strategies. In each of the three primary areas of volatility modeling, namely, conditional (or generalized autoregressive conditional heteroskedasticity) volatility, stochastic volatility and realized volatility (RV), numerous univariate volatility models of individual financial assets and multivariate volatility models of portfolios of assets have been established. This special issue has eleven innovative articles, eight of which are focused directly on RV and three on long memory, while two are concerned with both RV and long memory.