A Note on the Characterization of Positive Linnik Laws

Recently, Jayakumar & Pillai (1996) gave an interesting characterization of the positive Linnik laws in terms of the spectrum function of an infinitely divisible law. This paper improves their result and simplifies their proof. It proves another characterization result in terms of the Pareto law. Further, it represents the positive Linnik random variable as a function of independent gamma random variables.     

On a class of distributions on the simplex

In the present paper we define and investigate a novel class of distributions on the simplex, termed normalized infinitely divisible distributions, which includes the Dirichlet distribution. Distributional properties and general moment formulae are derived. Particular attention is devoted to special cases of normalized infinitely divisible distributions which lead to explicit expressions. As a by-product new distributions over the unit interval and a generalization of the Bessel function distribution are obtained.     

Max‐infinitely divisible models and inference for spatial extremes

For many environmental processes, recent studies have shown that the dependence strength is decreasing when quantile levels increase. This implies that the popular max‐stable models are inadequate to capture the rate of joint tail decay, and to estimate joint extremal probabilities beyond observed levels. We here develop a more flexible modeling framework based on the class of max‐infinitely divisible processes, which extend max‐stable processes while retaining dependence properties that are natural for maxima. We propose two parametric constructions for max‐infinitely divisible models, which relax the max‐stability property but remain close to some popular max‐stable models obtained as special cases. The first model considers maxima over a finite, random number of independent observations, while the second model generalizes the spectral representation of max‐stable processes. Inference is performed using a pairwise likelihood. We illustrate the benefits of our new modeling framework on Dutch wind gust maxima calculated over different time units. Results strongly suggest that our proposed models outperform other natural models, such as the Student‐t copula process and its max‐stable limit, even for large block sizes.     

Tempered stable Ornstein– Uhlenbeck processes: A practical view

We study the one-dimensional Ornstein–Uhlenbeck (OU) processes with marginal law given by tempered stable and tempered infinitely divisible distributions. We investigate the transition law between consecutive observations of these processes and evaluate the characteristic function of integrated tempered OU processes with a view toward practical applications. We then analyze how to draw a random sample from this class of processes by considering both the classical inverse transform algorithm and an acceptance–rejection method based on simulating a stable random sample. Using a maximum likelihood estimation method based on the fast Fourier transform, we empirically assess the simulation algorithm performance.     

On rates of convergence to infinitely divisible laws for dependent random variables

Large O and small o approximations of the expected value of a class of functions (modified K-functional and Lipschitz class) of the normalized partial sums of dependent random variables by the expectation of the corresponding functions of infinitely divisible random variables have been established. As a special case, we have obtained rates of convergence to the Stable Limit Laws and to the Weak Laws of Large Numbers. The technique used is the conditional version of the operator method of Trotter and the Taylor expansion.     

On some dependence structures for multidimensional Lévy driven moving averages

The Lévy copula can describe the dependence structure of a multidimensional Lévy process or a multivariate infinitely divisible random variable. Suppose the Lévy copula of a multidimensional Lévy process is known. We present the Lévy copula of the Lévy measure of the moving average driven by the multidimensional Lévy process. If there exist some special dependence structures among the components of the Lévy process, we give some dependence invariance properties after the transform of the moving average.     

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.     

The Lindsay transform of natural exponential families

Let μ be an infinitely divisible positive measure on R. If the measure ρ_μ is such that x^-2[ρ_μ(dx)—ρ_μ({0})δ₀(dx)] is the Lévy measure associated with μ and is infinitely divisible, we consider for all positive reals α and β the measure T_α,β(μ) which is the convolution of μ*^α and ρ_μ*β. For example, if μ is the inverse Gaussian law, then ρ_μ is a gamma law with paramter 3/2. Then T_α,β(μ) is an extension of the Lindsay transform of the first order, restricted to the distributions which are infinitely divisible. The main aim of this paper is to point out that it is possible to apply this transformation to all natural exponential families (NEF) with strictly cubic variance functions P. We then obtain NEF with variance functions of the form □ΔP(□Δ), where A is an affine function of the mean of the NEF. Some of these latter types appear scattered in the literature.     

Characterizations of discrete compound Poisson distributions

ABSTRACT

The aim of this paper is to give some new characterizations of discrete compound Poisson distributions. Firstly, we give a characterization by the Lévy–Khintchine formula of infinitely divisible distributions under some conditions. The second characterization need to present by row sum of random triangular arrays converges in distribution. And we give an application in probabilistic number theory, the strongly additive function converging to a discrete compound Poisson in distribution. The next characterization, is an extension of Watanabe’s theorem of characterization of homogeneous Poisson process. The last characterization will be illustrated by waiting time distributions, especially the matrix-exponential representation.     

On the two-parameter Bell–Touchard discrete distribution

Abstract

In this paper we introduce a new two-parameter discrete distribution which may be useful for modeling count data. Additionally, the probability mass function is very simple and it may have a zero vertex. We show that the new discrete distribution is a particular solution of a multiple Poisson process, and that it is infinitely divisible. Additionally, various structural properties of the new discrete distribution are derived. We also discuss two methods (moments and maximum likelihood) to estimate the model parameters. The usefulness of the proposed distribution is illustrated by means of real data sets to prove its versatility in practical applications.     

Some families of group divisible designs

Generalizing methods of constructions of Hadamard group divisible designs due to Bush (1979), some new families of semi-regular or regular group divisible designs are produced. Furthermore, new nonisomorphic solutions for some known group divisible designs are given, together with useful group divisible designs not listed in Clatworthy (1973).     

Accurate and efficient numerical calculation of stable densities via optimized quadrature and asymptotics

Stable distributions are an important class of infinitely divisible probability distributions, of which two special cases are the Cauchy distribution and the normal distribution. Aside from a few special cases, the density function for stable distributions has no known analytic form and is expressible only through the variate’s characteristic function or other integral forms. In this paper, we present numerical schemes for evaluating the density function for stable distributions, its gradient, and distribution function in various parameter regimes of interest, some of which had no preexisting efficient method for their computation. The novel evaluation schemes consist of optimized generalized Gaussian quadrature rules for integral representations of the density function, complemented by asymptotic expansions near various values of the shape and argument parameters. We report several numerical examples illustrating the efficiency of our methods. The resulting code has been made available online.     

A new infinitely divisible discrete distribution with applications to count data modeling

A new discrete distribution involving geometric and discrete Pareto as special cases is introduced. The distribution possesses many interesting properties like decreasing hazard rate, zero vertex uni-modality, over-dispersion, infinite divisibility and compound Poisson representation, which makes the proposed distribution well suited for count data modeling. Other issues including closure property under minima, comparison of its distribution tail with other distributions via actuarial indices are discussed. The method of proportion and maximum likelihood method are presented for parameter estimation. Finally the performance of the proposed distribution over other classical and newly proposed infinitely divisible distributions are discussed.     

Integer‐valued Trawl Processes: A Class of Stationary Infinitely Divisible Processes

This paper introduces a new continuous‐time framework for modelling serially correlated count and integer‐valued data. The key component in our new model is the class of integer‐valued trawl processes, which are serially correlated, stationary, infinitely divisible processes. We analyse the probabilistic properties of such processes in detail and, in addition, study volatility modulation and multivariate extensions within the new modelling framework. Moreover, we describe how the parameters of a trawl process can be estimated and obtain promising estimation results in our simulation study. Finally, we apply our new modelling framework to high‐frequency financial data.     

On Infinitely Divisible Exponential Dispersion Model Related to Poisson-Exponential Distribution

We construct a univariate exponential dispersion model comprised of discrete infinitely divisible distributions. This model emerges in the theory of branching processes. We obtain a representation for the Lévy measure of relevant distributions and characterize their laws as Poisson mixtures and/or compound Poisson distributions. The regularity of the unit variance function of this model is employed for the derivation of approximations by the Poisson-exponential model. We emphasize the role of the latter class. We construct local approximations relating them to properties of special functions and branching diffusions.     

Characterizations,length-biasing,and infinite divisibility

SupposeL(X) is the law of a positive random variableX, andZ is positive and independent ofX. Admissible solution pairs (L(X),L(Z)) are sought for the in-law equation $\hat X \cong X o Z$ °Z, where $L\left( {\hat X} \right)$ is a weighted law constructed fromL(X), and ° is a binary operation which in some sense is increasing. The class of weights includes length biasing of arbitrary order. When ° is addition and the weighting is ordinary length biasing, the class of admissibleL(X) comprises the positive infinitely divisible laws. Examples are given subsuming all known specific cases. Some extensions for general order of length-biasing are discussed.     

Characterization by invariance under length-biasing and random scaling

A positive random variable X with law L(X) and finite moment of order r > 0 has an induced length-biased law of order r, denoted by L(X_r). Let V ⩾ 0 be independent of X_r. A characterization problem seeks solution pairs (L(X), L(V)) for the “in-law” equation X ≅ VX_r, where ≅ denotes equality in law. A renewal process interpretation asks when is the random rescaling of the stationary total lifetime VX₁ equal in law to a typical lifetime X? Solutions are known in special cases.A comprehensive existence/uniqueness theory is presented, and many consequences are explored. Unique solutions occur when − log X and − log V have spectrally positive infinitely divisible laws. Particular cases are explored.Connections with the stationary lifetime law of renewal theory also are investigated.     

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.     

On the optimalizes of the mbgdd's under mixed effects model

In the present paper, under the assumption of a mixed effects model with random block effects, the type 1 optimality of the most balanced group divisible designs of type 1 has been established within the general class of all proper and connected block designs with k<v.