Simulation of the continuous time random walk of the space-fractional diffusion equations

In this article, we discuss the solution of the space-fractional diffusion equation with and without central linear drift in the Fourier domain and show the strong connection between it and the α$\alpha$-stable Lévy distribution, 0<α<2$0<\alpha <2$. We use some relevant transformations of the independent variables x$x$ and t$t$, to find the solution of the space-fractional diffusion equation with central linear drift which is a special form of the space-fractional Fokker–Planck equation which is useful in studying the dynamic behaviour of stochastic differential equations driven by the non-Gaussian (Lévy) noises. We simulate the continuous time random walk of these models by using the Monte Carlo method.     

On functional limits of short- and long-memory linear processes with GARCH(1,1) noises

This paper considers the short- and long-memory linear processes with GARCH (1,1) noises. The functional limit distributions of the partial sum and the sample autocovariances are derived when the tail index α$\alpha$ is in (0,2)$(0,2)$, equal to 2, and in (2,∞)$(2,\infty )$, respectively. The partial sum weakly converges to a functional of α$\alpha$-stable process when α<2$\alpha <2$ and converges to a functional of Brownian motion when α≥2$\alpha \ge 2$. When the process is of short-memory and α<4$\alpha <4$, the autocovariances converge to functionals of α/2$\alpha /2$-stable processes; and if α≥4$\alpha \ge 4$, they converge to functionals of Brownian motions. In contrast, when the process is of long-memory, depending on α$\alpha$ and β$\beta$ (the parameter that characterizes the long-memory), the autocovariances converge to either (i) functionals of α/2$\alpha /2$-stable processes; (ii) Rosenblatt processes (indexed by β$\beta$, 1/2<β<3/4$1/2<\beta <3/4$); or (iii) functionals of Brownian motions. The rates of convergence in these limits depend on both the tail index α$\alpha$ and whether or not the linear process is short- or long-memory. Our weak convergence is established on the space of càdlàg functions on [0,1]$[0,1]$ with either (i) the _J1$J_1$ or the _M1$M_1$ topology (Skorokhod, 1956); or (ii) the weaker form S$S$ topology (Jakubowski, 1997). Some statistical applications are also discussed.     

Tempering stable processes

A tempered stable Lévy process combines both the α$\alpha$-stable and Gaussian trends. In a short time frame it is close to an α$\alpha$-stable process while in a long time frame it approximates a Brownian motion. In this paper we consider a general and robust class of multivariate tempered stable distributions and establish their identifiable parametrization. We prove short and long time behavior of tempered stable Lévy processes and investigate their absolute continuity with respect to the underlying α$\alpha$-stable processes. We find probabilistic representations of tempered stable processes which specifically show how such processes are obtained by cutting (tempering) jumps of stable processes. These representations exhibit α$\alpha$-stable and Gaussian tendencies in tempered stable processes and thus give probabilistic intuition for their study. Such representations can also be used for simulation. We also develop the corresponding representations for Ornstein–Uhlenbeck-type processes.     

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.     

Markov processes invariant under a Lie group action

We show that a Markov process in a manifold invariant under the action of a compact Lie group K$K$ induces a Lévy process in each K$K$-orbit by “forcing” it to run in the orbit.     

Comparing the minimal Hellinger martingale measure of order q to the q-optimal martingale measure

This paper investigates the relationship between the minimal Hellinger martingale measure of order q$q$ (MHM measure hereafter) and the q$q$-optimal martingale measure for any q≠1$q\ne 1$. First, we provide more results for the MHM measure; in particular we establish its complete characterization in two manners. Then we derive two equivalent conditions for both martingale measures to coincide. These conditions are in particular fulfilled in the case of markets driven by Lévy processes. Finally, we analyze the MHM measure as well as its relationship to the q$q$-optimal martingale measure for the case of a discrete-time market model.     

Multivariate CARMA processes

A multivariate Lévy-driven continuous time autoregressive moving average (CARMA) model of order (p,q$p,q$), q<p$q<p$, is introduced. It extends the well-known univariate CARMA and multivariate discrete time ARMA models. We give an explicit construction using a state space representation and a spectral representation of the driving Lévy process. Furthermore, various probabilistic properties of the state space model and the multivariate CARMA process itself are discussed in detail.     

Randomly weighted self-normalized Lévy processes

Let (_Ut,_Vt)$(U_t,V_t)$ be a bivariate Lévy process, where _Vt$V_t$ is a subordinator and _Ut$U_t$ is a Lévy process formed by randomly weighting each jump of _Vt$V_t$ by an independent random variable _Xt$X_t$ having cdf F$F$. We investigate the asymptotic distribution of the self-normalized Lévy process _Ut/_Vt$U_t/V_t$ at 0 and at ∞$\infty$. We show that all subsequential limits of this ratio at 0 (∞$\infty$) are continuous for any nondegenerate F$F$ with finite expectation if and only if _Vt$V_t$ belongs to the centered Feller class at 0 (∞$\infty$). We also characterize when _Ut/_Vt$U_t/V_t$ has a non-degenerate limit distribution at 0 and ∞$\infty$.     

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.     

Oscillation of harmonic functions for subordinate Brownian motion and its applications

In this paper, we establish an oscillation estimate of nonnegative harmonic functions for a pure-jump subordinate Brownian motion. The infinitesimal generator of such subordinate Brownian motion is an integro-differential operator. As an application, we give a probabilistic proof of the following form of relative Fatou theorem for such subordinate Brownian motion X$X$ in a bounded κ$\kappa$-fat open set; if u$u$ is a positive harmonic function with respect to X$X$ in a bounded κ$\kappa$-fat open set D$D$ and h$h$ is a positive harmonic function in D$D$ vanishing on D^c$D^c$, then the non-tangential limit of u/h$u/h$ exists almost everywhere with respect to the Martin-representing measure of h$h$.     

Measurability of semimartingale characteristics with respect to the probability law

Given a càdlàg process X$X$ on a filtered measurable space, we construct a version of its semimartingale characteristics which is measurable with respect to the underlying probability law. More precisely, let _Psem${\mathfrak{P}}_{sem}$ be the set of all probability measures P$P$ under which X$X$ is a semimartingale. We construct processes (B^P,C,ν^P)$(B^P,C,{\nu }^P)$ which are jointly measurable in time, space, and the probability law P$P$, and are versions of the semimartingale characteristics of X$X$ under P$P$ for each P∈_Psem$P\in {\mathfrak{P}}_{sem}$. This result gives a general and unifying answer to measurability questions that arise in the context of quasi-sure analysis and stochastic control under the weak formulation.     

On the Wiener–Hopf factorization for Lévy processes with bounded positive jumps

We study the Wiener–Hopf factorization for Lévy processes with bounded positive jumps and arbitrary negative jumps. We prove that the positive Wiener–Hopf factor can be expressed as an infinite product involving solutions to the equation ψ(z)=q$\psi \left(z\right)=q$, where ψ$\psi$ is the Laplace exponent. Under additional regularity assumptions on the Lévy measure we obtain an asymptotic expression for these solutions. When the process is spectrally negative with bounded jumps, we derive a series representation for the scale function. In order to illustrate possible applications, we discuss the implementation of numerical algorithms and present the results of several numerical experiments.     

Potential theory of subordinate Brownian motions with Gaussian components

In this paper we study a subordinate Brownian motion with a Gaussian component and a rather general discontinuous part. The assumption on the subordinator is that its Laplace exponent is a complete Bernstein function with a Lévy density satisfying a certain growth condition near zero. The main result is a boundary Harnack principle with explicit boundary decay rate for non-negative harmonic functions of the process in C^1,1$C^{1,1}$ open sets. As a consequence of the boundary Harnack principle, we establish sharp two-sided estimates on the Green function of the subordinate Brownian motion in any bounded C^1,1$C^{1,1}$ open set D$D$ and identify the Martin boundary of D$D$ with respect to the subordinate Brownian motion with the Euclidean boundary.     

Approximate martingale estimating functions for stochastic differential equations with small noises

An approximate martingale estimating function with an eigenfunction is proposed for an estimation problem about an unknown drift parameter for a one-dimensional diffusion process with small perturbed parameter ε$\epsilon$ from discrete time observations at n$n$ regularly spaced time points k/n$k/n$, k=0,1,…,n$k=0,1,\dots ,n$. We show asymptotic efficiency of an M$M$-estimator derived from the approximate martingale estimating function as ε→0$\epsilon \to 0$ and n→∞$n\to \infty$ simultaneously.     

Nonparametric estimation for pure jump Lévy processes based on high frequency data

In this paper, we study nonparametric estimation of the Lévy density for pure jump Lévy processes. We consider n$n$ discrete time observations with step Δ$\Delta$. The asymptotic framework is: n$n$ tends to infinity, Δ=_Δn$\Delta ={\Delta }_n$ tends to zero while n_Δn$n{\Delta }_n$ tends to infinity. First, we use a Fourier approach (“frequency domain”): this allows us to construct an adaptive nonparametric estimator and to provide a bound for the global L²${\mathbb{L}}^2$-risk. Second, we use a direct approach (“time domain”) which allows us to construct an estimator on a given compact interval. We provide a bound for L²${\mathbb{L}}^2$-risk restricted to the compact interval. We discuss rates of convergence and give examples and simulation results for processes fitting in our framework.     

A strong uniform approximation of fractional Brownian motion by means of transport processes

We construct a sequence of processes that converges strongly to fractional Brownian motion uniformly on bounded intervals for any Hurst parameter H$H$, and we derive a rate of convergence, which becomes better when H$H$ approaches 1/2$1/2$. The construction is based on the Mandelbrot–van Ness stochastic integral representation of fractional Brownian motion and on a strong transport process approximation of Brownian motion. The objective of this method is to facilitate simulation.     

Large deviations of infinite intersections of events in Gaussian processes

Consider events of the form {_Zs≥ζ(s),s∈S}$\{Z_s\ge \zeta \left(s\right),s\in S\}$, where Z$Z$ is a continuous Gaussian process with stationary increments, ζ$\zeta$ is a function that belongs to the reproducing kernel Hilbert space R$R$ of process Z$Z$, and S⊂R$S\subset \mathbb{R}$ is compact. The main problem considered in this paper is identifying the function β^∗∈R${\beta }^{\ast }\in R$ satisfying β^∗(s)≥ζ(s)${\beta }^{\ast }\left(s\right)\ge \zeta \left(s\right)$ on S$S$ and having minimal R$R$-norm. The smoothness (mean square differentiability) of Z$Z$ turns out to have a crucial impact on the structure of the solution. As examples, we obtain the explicit solutions when ζ(s)=s$\zeta \left(s\right)=s$ for s∈[0,1]$s\in [0,1]$ and Z$Z$ is either a fractional Brownian motion or an integrated Ornstein–Uhlenbeck process.