Normal Inverse Gaussian Distributions and Stochastic Volatility Modelling

The normal inverse Gaussian distribution is defined as a variance-mean mixture of a normal distribution with the inverse Gaussian as the mixing distribution. The distribution determines an homogeneous Lévy process, and this process is representable through subordination of Brownian motion by the inverse Gaussian process. The canonical, Lévy type, decomposition of the process is determined. As a preparation for developments in the latter part of the paper the connection of the normal inverse Gaussian distribution to the classes of generalized hyperbolic and inverse Gaussian distributions is briefly reviewed. Then a discussion is begun of the potential of the normal inverse Gaussian distribution and Lévy process for modelling and analysing statistical data, with particular reference to extensive sets of observations from turbulence and from finance. These areas of application imply a need for extending the inverse Gaussian Lévy process so as to accommodate certain, frequently observed, temporal dependence structures. Some extensions, of the stochastic volatility type, are constructed via an observation-driven approach to state space modelling. At the end of the paper generalizations to multivariate settings are indicated.     

Pricing longevity-linked derivatives using a stochastic mortality model

Abstract

We propose a 2-factor MBMM model with exponential Lévy process to develop a stochastic mortality process. The two components are fitted by two independent NIG distributions. Compared to Lee–Carter model or 1-factor MBMM model, our mortality model explains more variation and improves the goodness of fit by including the second time component. Based on the improved model, we price three longevity-linked financial instruments, namely the longevity bond, q-forward and s-forward. The pricing is demonstrated on English and Welsh males aged 65 in 2013. Results indicate that the 2-factor MBMM model gives the highest price for mortality-related type of contract.     

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.     

Lévy‐based Modelling in Brain Imaging

Abstract. A substantive problem in neuroscience is the lack of valid statistical methods for non‐Gaussian random fields. In the present study, we develop a flexible, yet tractable model for a random field based on kernel smoothing of a so‐called Lévy basis. The resulting field may be Gaussian, but there are many other possibilities, including random fields based on Gamma, inverse Gaussian and normal inverse Gaussian (NIG) Lévy bases. It is easy to estimate the parameters of the model and accordingly to assess by simulation the quantiles of test statistics commonly used in neuroscience. We give a concrete example of magnetic resonance imaging scans that are non‐Gaussian. For these data, simulations under the fitted models show that traditional methods based on Gaussian random field theory may leave small, but significant changes in signal level undetected, while these changes are detectable under a non‐Gaussian Lévy model.     

Prior Distributions on Measure Space

A measure is the formal representation of the non-negative additive functions that abound in science. We review and develop the art of assigning Bayesian priors to measures. Where necessary, spatial correlation is delegated to correlating kernels imposed on otherwise uncorrelated priors. The latter must be infinitely divisible (ID) and hence described by the Lévy–Khinchin representation. Thus the fundamental object is the Lévy measure, the choice of which corresponds to different ID process priors. The general case of a Lévy measure comprising a mixture of assigned base measures leads to a prior process comprising a convolution of corresponding processes. Examples involving a single base measure are the gamma process, the Dirichlet process (for the normalized case) and the Poisson process. We also discuss processes that we call the supergamma and super-Dirichlet processes, which are double base measure generalizations of the gamma and Dirichlet processes. Examples of multiple and continuum base measures are also discussed. We conclude with numerical examples of density estimation.     

Lévy Processes,Phase-Type Distributions,and Martingales

Lévy processes are defined as processes with stationary independent increments and have become increasingly popular as models in queueing, finance, etc.; apart from Brownian motion and compound Poisson processes, some popular examples are stable processes, variance gamma processes, CGMY Lévy processes (tempered stable processes), NIG (normal inverse Gaussian) Lévy processes, and hyperbolic Lévy processes. We consider here a dense class of Lévy processes, compound Poisson processes with phase-type jumps in both directions and an added Brownian component. Within this class, we survey how to explicitly compute a number of quantities that are traditionally studied in the area of Lévy processes, in particular two-sided exit probabilities and associated Laplace transforms, the closely related scale function, one-sided exit probabilities and associated Laplace transforms coming up in queueing problems, and similar quantities for a Lévy process with reflection in 0. The solutions are in terms of roots to polynomials, and the basic equations are derived by purely probabilistic arguments using martingale optional stopping; a particularly useful martingale is the so-called Kella-Whitt martingale. Also, the relation to fluid models with a Brownian component is discussed.     

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.     

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.     

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.     

Bayesian Conditional Mean Estimation in Log‐Normal Linear Regression Models with Finite Quadratic Expected Loss

Log‐normal linear regression models are popular in many fields of research. Bayesian estimation of the conditional mean of the dependent variable is problematic as many choices of the prior for the variance (on the log‐scale) lead to posterior distributions with no finite moments. We propose a generalized inverse Gaussian prior for this variance and derive the conditions on the prior parameters that yield posterior distributions of the conditional mean of the dependent variable with finite moments up to a pre‐specified order. The conditions depend on one of the three parameters of the suggested prior; the other two have an influence on inferences for small and medium sample sizes. A second goal of this paper is to discuss how to choose these parameters according to different criteria including the optimization of frequentist properties of posterior means.     

On the Normal Inverse Gaussian Stochastic Volatility Model

In this article, the normal inverse Gaussian stochastic volatility model of Barndorff-Nielsen is extended. The resulting model has a more flexible lag structure than the original one. In addition, the second-and fourth-order moments, important properties of a volatility model, are derived. The model can be considered either as a generalized autoregressive conditional heteroscedasticity model with nonnormal errors or as a stochastic volatility model with an inverse Gaussian distributed conditional variance. A simulation study is made to investigate the performance of the maximum likelihood estimator of the model. Finally, the model is applied to stock returns and exchange-rate movements. Its fit to two stylized facts and its forecasting performance is compared with two other volatility models.     

An asymmetric multivariate weibull distribution

Abstract

A class of multivariate laws as an extension of univariate Weibull distribution is presented. A well known representation of the asymmetric univariate Laplace distribution is used as the starting point. This new family of distributions exhibits some similarities to the multivariate normal distribution. Properties of this class of distributions are explored including moments, correlations, densities and simulation algorithms. The distribution is applied to model bivariate exchange rate data. The fit of the proposed model seems remarkably good. Parameters are estimated and a bootstrap study performed to assess the accuracy of the estimators.     

Moments of nonparametric probability density functions of entropy estimators applied to testing the inverse Gaussian distribution

The inverse Gaussian (IG) distribution is widely used to model data and then it is important to develop efficient goodness of fit tests for this distribution. In this article, we introduce some new test statistics for examining the IG goodness of fit based on correcting moments of nonparametric probability density functions of entropy estimators. These tests are consistent against all alternatives. Critical points and power of the tests are explored by simulation. We show that the proposed tests are more powerful than competitor tests. Finally, the proposed tests are illustrated by real data examples.     

Bayesian Survival Analysis in Proportional Hazard Models with Logistic Relative Risk

Abstract.  The traditional Cox proportional hazards regression model uses an exponential relative risk function. We argue that under various plausible scenarios, the relative risk part of the model should be bounded, suggesting also that the traditional model often might overdramatize the hazard rate assessment for individuals with unusual covariates. This motivates our working with proportional hazards models where the relative risk function takes a logistic form. We provide frequentist methods, based on the partial likelihood, and then go on to semiparametric Bayesian constructions. These involve a Beta process for the cumulative baseline hazard function and any prior with a density, for example that dictated by a Jeffreys-type argument, for the regression coefficients. The posterior is derived using machinery for Lévy processes, and a simulation recipe is devised for sampling from the posterior distribution of any quantity. Our methods are illustrated on real data. A Bernshtĕn–von Mises theorem is reached for our class of semiparametric priors, guaranteeing asymptotic normality of the posterior processes.     

Generalized Inverse Gaussian Distributions and their Wishart Connections

The matrix generalized inverse Gaussian distribution (MGIG) is shown to arise as a conditional distribution of components of a Wishart distributio n. In the special scalar case, the characterization refers to members of the class of generalized inverse Gaussian distributions (GIGs) and includes the inverse Gaussian distribution among others     

Multilevel particle filters for Lévy-driven stochastic differential equations

Statistics and Computing - We develop algorithms for computing expectations with respect to the laws of models associated to stochastic differential equations driven by pure Lévy processes. We...     

Regression Properties for Asymmetric Generalized Scale Mixtures of Multivariate Gaussian Variables

Analytical properties of regression and the variance–covariance matrix of asymmetric generalized scale mixture of multivariate Gaussian variables are presented. The analysis includes an in-depth analytical investigation of the first two conditional moments of the mixing variable. Exact computable expressions for the prediction and the conditional variance are presented for the generalized hyperbolic distribution using the inversion theorem for Fourier transforms. An application to financial log returns is demonstrated via the classical Euler approximation. The methodology is illustrated by analyzing the regression of intraday log returns for CISCO against the corresponding data from S&P 500.     

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.     

Kolmogorov-Smirnov Type Test-Statistics For The Gamma,Erlang-2 And The Inverse Gaussian Distributions When The Parameters Are Unknown.

Critical values are presented for the Kolmogorov-Smirnov type test statistics for the following three cases: (i) the gamma distribution when both the scale and the shape parameters are not known, (ii) the scale parameter of the gamma distribution is not known and (iii) the inverse Gaussian distribution when both the parameters are unknown. This study was motivated by the necessity to fit the gamma, the Erlang-2 and the inverse Gaussian distributions to the interpurchase times of individuals for coffee in marketing research.