Characterization of dependence of multidimensional Lévy processes using Lévy copulas

This paper suggests Lévy copulas in order to characterize the dependence among components of multidimensional Lévy processes. This concept parallels the notion of a copula on the level of Lévy measures. As for random vectors, a version of Sklar's theorem states that the law of a general multivariate Lévy process is obtained by combining arbitrary univariate Lévy processes with an arbitrary Lévy copula. We construct parametric families of Lévy copulas and prove a limit theorem, which indicates how to obtain the Lévy copula of a multivariate Lévy process X from the ordinary copula of the random vector X_t for small t.     

Exact asymptotics of supremum of a stationary Gaussian process over a random interval

Let {X(t):t∈[0,∞)} be a centered stationary Gaussian process. We study the exact asymptotics of P(sup_s∈[0,T]X(s)>u), as u→∞, where T is an independent of {X(t)} nonnegative random variable. It appears that the heaviness of T impacts the form of the asymptotics, leading to three scenarios: the case of integrable T, the case of T having regularly varying tail distribution with parameter λ∈(0,1) and the case of T having slowly varying tail distribution.     

Scaling limits of Gaussian vector fields

For Gaussian vector fields {X(t) ∈ Rⁿ:t ∈ R^d} we describe the covariance functions of all scaling limits Y(t) = $\text{L}$lim_α↓0 B^?1(α) X(αt) which can occur when B(α) is a d × d matrix function with B(α) → 0. These matrix covariance functions $\text{r(t, s) = EY(t) Y}{}^1\text{(s)}$ are found to be homogeneous in the sense that for some matrix L and each α > 0, $\text{(}{}^1\text{) r(αt, αs) = α}{}^{\mathrm{L1}}\text{r(t, s) α}{}^{\mathrm{L}}$. Processes with stationary increments satisfying (¹) are further analysed and are found to be natural generalizations of Lévy's multiparameter Brownian motion.     

The singularity spectrum of Lévy processes in multifractal time

Let X=(Xt)_t?0 be a Lévy process and μ a positive Borel measure on R₊. Suppose that the integral of μ defines a continuous increasing multifractal time . Under suitable assumptions on μ, we compute the singularity spectrum of the sample paths of the process X in time μ defined as the process (XF(t))_t?0.A fundamental example consists in taking a measure μ equal to an “independent random cascade” and (independently of μ) a suitable stable Lévy process X. Then the associated process X in time μ is naturally related to the so-called fixed points of the smoothing transformation in interacting particles systems.Our results rely on recent heterogeneous ubiquity theorems.     

Lévy-Driven Carma Processes

Properties and examples of continuous-time ARMA (CARMA) processes driven by Lévy processes are examined. By allowing Lévy processes to replace Brownian motion in the definition of a Gaussian CARMA process, we obtain a much richer class of possibly heavy-tailed continuous-time stationary processes with many potential applications in finance, where such heavy tails are frequently observed in practice. If the Lévy process has finite second moments, the correlation structure of the CARMA process is the same as that of a corresponding Gaussian CARMA process. In this paper we make use of the properties of general Lévy processes to investigate CARMA processes driven by Lévy processes {W(t)} without the restriction to finite second moments. We assume only that W (1) has finite r-th absolute moment for some strictly positive r. The processes so obtained include CARMA processes with marginal symmetric stable distributions.     

A bivariate Lévy process with negative binomial and gamma marginals

The joint distribution of X and N, where N has a geometric distribution and X is the sum of N IID exponential variables (independent of N), is infinitely divisible. This leads to a bivariate Lévy process {(X(t),N(t)),t≥0}, whose coordinates are correlated negative binomial and gamma processes. We derive basic properties of this process, including its covariance structure, representations, and stochastic self-similarity. We examine the joint distribution of (X(t),N(t)) at a fixed time t, along with the marginal and conditional distributions, joint integral transforms, moments, infinite divisibility, and stability with respect to random summation. We also discuss maximum likelihood estimation and simulation for this model.     

Uniform Hölder exponent of a stationary increments Gaussian process: Estimation starting from average values

Let {X(t)}_t∈R be a stationary increments Gaussian process satisfying some assumptions. By using the notion of generalized quadratic variation we build a strongly consistent and asymptotically normal estimator of the uniform Hölder exponent of X, over a compact interval. Our estimator is obtained starting from average values of the process over a regular grid.     

Regularity of semigroups generated by Lévy type operators via coupling

By adopting the coupling method, we obtain new verifiable sufficient conditions about the _Cb(R^d)$C_b\left({\mathbb{R}}^d\right)$-Feller continuity, the Lipschitz continuity and the strong Feller continuity of the semigroups associated with Lévy type operators. These results easily apply to jump–diffusion processes and stochastic differential equations driven by Lévy processes. Our results also yield the criterion for the e$e$-property (namely the characterization about the equi-continuity of semigroups acting on bounded Lipschitz functions) of Lévy type operators, and show that both genuine Lévy processes and the Ornstein–Uhlenbeck type processes are e$e$-processes.     

On–off fluid models in heavy traffic environment

We consider fluid models with infinite buffer size. Let {Z _N (t)} be the net input rate to the buffer, where {{Z _N (t)} is a superposition of N homogeneous alternating on–off flows. Under heavy traffic environment {{Z _N (t)} converges in distribution to a centred Gaussian process with covariance function of a single flow. The aim of this paper is to prove the convergence of the stationary buffer content process {X _N ^* (t)} in the fNth model to the buffer content process {X _N (t)} in the limiting Gaussian model. This revised version was published online in June 2006 with corrections to the Cover Date.     

Sina?ˇ's condition for real valued Lévy processes

We prove that the upward ladder height subordinator H associated to a real valued Lévy process ξ has Laplace exponent φ that varies regularly at ∞ (respectively, at 0) if and only if the underlying Lévy process ξ satisfies Sina?ˇ's condition at 0 (respectively, at ∞). Sina?ˇ's condition for real valued Lévy processes is the continuous time analogue of Sina?ˇ's condition for random walks. We provide several criteria in terms of the characteristics of ξ to determine whether or not it satisfies Sina?ˇ's condition. Some of these criteria are deduced from tail estimates of the Lévy measure of H, here obtained, and which are analogous to the estimates of the tail distribution of the ladder height random variable of a random walk which are due to Veraverbeke and Grübel.     

Median,concentration and fluctuations for Lévy processes

We estimate a median of f(_Xt)$f\left(X_t\right)$ where f$f$ is a Lipschitz function, X$X$ is a Lévy process and t$t$ is an arbitrary time. This leads to concentration inequalities for f(_Xt)$f\left(X_t\right)$. In turn, corresponding fluctuation estimates are obtained under assumptions typically satisfied if the process has a regular behavior in small time and a, possibly different, regular behavior in large time.     

Fluctuations of the empirical quantiles of independent Brownian motions

We consider iid Brownian motions, B_j(t), where B_j(0) has a rapidly decreasing, smooth density function f. The empirical quantiles, or pointwise order statistics, are denoted by B_j:n(t), and we consider a sequence Q_n(t)=B_j(n):n(t), where j(n)/n→α∈(0,1). This sequence converges in probability to q(t), the α-quantile of the law of B_j(t). We first show convergence in law in C[0,∞) of F_n=n^1/2(Q_n−q). We then investigate properties of the limit process F, including its local covariance structure, and Hölder-continuity and variations of its sample paths. In particular, we find that F has the same local properties as fBm with Hurst parameter H=1/4.     

Distribution and characteristic functions for correlated complex Wishart matrices

Let A(t) be a complex Wishart process defined in terms of the M×N complex Gaussian matrix X(t) by A(t)=X(t)X(t)^H. The covariance matrix of the columns of X(t) is Σ. If X(t), the underlying Gaussian process, is a correlated process over time, then we have dependence between samples of the Wishart process. In this paper, we study the joint statistics of the Wishart process at two points in time, t₁, t₂, where t₁<t₂. In particular, we derive the following results: the joint density of the elements of A(t₁), A(t₂), the joint density of the eigenvalues of Σ^-1A(t₁),Σ^-1A(t₂), the characteristic function of the elements of A(t₁), A(t₂), the characteristic function of the eigenvalues of Σ^-1A(t₁),Σ^-1A(t₂). In addition, we give the characteristic functions of the eigenvalues of a central and non-central complex Wishart, and some applications of the results in statistics, engineering and information theory are outlined.     

Convergence to type I distribution of the extremes of sequences defined by random difference equation

We study the extremes of a sequence of random variables (R_n) defined by the recurrence R_n=M_nR_n−1+q, n≥1, where R₀ is arbitrary, (M_n) are iid copies of a non-degenerate random variable M, 0≤M≤1, and q>0 is a constant. We show that under mild and natural conditions on M the suitably normalized extremes of (R_n) converge in distribution to a double-exponential random variable. This partially complements a result of de Haan, Resnick, Rootzén, and de Vries who considered extremes of the sequence (R_n) under the assumption that P(M>1)>0.     

On a stochastic nonlinear equation arising from 1D integro-differential scalar conservation laws

In this paper we study the initial problem for a stochastic nonlinear equation arising from 1D integro-differential scalar conservation laws. The equation is driven by Lévy space-time white noise in the following form:

(t_∂−A)u+x_∂q(u)=f(u)+g(u)Ft,x