We give an extension of Lê’s stochastic sewing lemma. The stochastic sewing lemma proves convergence in $L_m$ of Riemann type sums $\sum _{[s,t] \in \pi } A_{s,t}$ for an adapted two-parameter stochastic process A, under certain conditions on the moments of $A_{s,t}$ and of conditional expectations of $A_{s,t}$ given $\mathcal F_s$. Our extension replaces the conditional expectation given $\mathcal F_s$ by that given $\mathcal F_v$ for $v<s$, and it allows to make use of asymptotic decorrelation properties between $A_{s,t}$ and $\mathcal F_v$ by including a singularity in $(s-v)$. We provide three applications for which Lê’s stochastic sewing lemma seems to be insufficient. The first is to prove the convergence of Itô or Stratonovich approximations of stochastic integrals along fractional Brownian motions under low regularity assumptions. The second is to obtain new representations of local times of fractional Brownian motions via discretization. The third is to improve a regularity assumption on the diffusion coefficient of a stochastic differential equation driven by a fractional Brownian motion for pathwise uniqueness and strong existence.

This paper deals with ergodic theorems for particular time-inhomogeneous Markov processes, whose time-inhomogeneity is asymptotically periodic. Under a Lyapunov/minorization condition, it is shown that, for any measurable bounded function f, the time average $\frac{1}{t} \int_0^t f(X_s)ds$ converges in $\mathbb{L}^2$ towards a limiting distribution, starting from any initial distribution for the process $(X_t)_{t \geq 0}$ . This convergence can be improved to an almost sure convergence under an additional assumption on the initial measure. This result is then applied to show the existence of a quasi-ergodic distribution for processes absorbed by an asymptotically periodic moving boundary, satisfying a conditional Doeblin condition.

In the classical framework, a random walk on a group is a Markov chain with independent and identically distributed increments. In some sense, random walks are time and space homogeneous. This paper is devoted to a class of inhomogeneous random walks on $\mathbb{Z}^d$ termed ‘Markov additive processes’ (also known as Markov random walks, random walks with internal degrees of freedom, or semi-Markov processes). In this model, the increments of the walk are still independent but their distributions are dictated by a Markov chain, termed the internal Markov chain. While this model is largely studied in the literature, most of the results involve internal Markov chains whose operator is quasi-compact. This paper extends two results for more general internal operators: a local limit theorem and a sufficient criterion for their transience. These results are thereafter applied to a new family of models of drifted random walks on the lattice $\mathbb{Z}^d$ .

For a continuous-time random walk X = {Xt, t ⩾ 0} (in general non-Markov), we study the asymptotic behaviour, as t → ∞, of the normalized additive functional $c_t\int _0^{t} f(X_s)\,{\rm d}s$, t⩾ 0. Similarly to the Markov situation, assuming that the distribution of jumps of X belongs to the domain of attraction to α-stable law with α > 1, we establish the convergence to the local time at zero of an α-stable Lévy motion. We further study a situation where X is delayed by a random environment given by the Poisson shot-noise potential: $\Lambda (x,\gamma )= {\rm e}^{-\sum _{y\in \gamma } \phi (x-y)},$ where $\phi \colon \mathbb R\to [0,\infty )$ is a bounded function decaying sufficiently fast, and γ is a homogeneous Poisson point process, independent of X. We find that in this case the weak limit has both ‘quenched’ component depending on Λ, and a component, where Λ is ‘averaged’.

A new approach to the problem of finding the distribution of integral functionals under the excursion measure is presented. It is based on the technique of excursion straddling a time, stochastic analysis, and calculus on local time, and it is done for Brownian motion with drift reflecting at 0, and under some additional assumptions for some class of Itó diffusions. The new method is an alternative to the classical potential-theoretic approach and gives new specific formulas for distributions under the excursion measure.

Let X = (Xt)t≥0 be a stochastic process issued from $x \in \mathbb{R}$ that admits a marginal stationary measure v, i.e. vPtf = vf for all t ≥ 0, where $\textbf{P}_t\,f(x)= \mathbb{E}_x[f(\textbf{X}_t)]$ . In this paper, we introduce the (resp. biorthogonal) spectral projections correlation functions which are expressed in terms of projections.” Also, update first published online date, if available. into the eigenspaces of Pt (resp. and of its adjoint in the weighted Hilbert space L2 (v)). We obtain closed-form expressions involving eigenvalues, the condition number and/or the angle between the projections in the following different situations: when X = X with X = (Xt)t ≥ 0 being a Markov process, X is the subordination of X in the sense of Bochner, and X is a non-Markovian process which is obtained by time-changing X with an inverse of a subordinator. It turns out that these spectral projections correlation functions have different expressions with respect to these classes of processes which enables to identify substantial and deep properties about their dynamics. This interesting fact can be used to design original statistical tests to make inferences, for example, about the path properties of the process (presence of jumps), distance from symmetry (self-adjoint or non-self-adjoint) and short-to-long-range dependence. To reveal the usefulness of our results, we apply them to a class of non-self-adjoint Markov semigroups studied in Patie and Savov (to appear, Mem. Amer. Math. Soc., 179p), and then time-change by subordinators and their inverses.

In this paper we study nonparametric estimation problems for a class of piecewise-deterministic Markov processes (PDMPs). Borovkov and Last (2008) proved a version of Rice's formula for PDMPs, which explains the relation between the stationary density and the level crossing intensity. From a statistical point of view, their result suggests a methodology for estimating the stationary density from observations of a sample path of PDMPs. First, we introduce the local time related to the level crossings and construct the local-time estimator for the stationary density, which is unbiased and uniformly consistent. Secondly, we investigate other estimation problems for the jump intensity and the conditional jump size distribution.

We study a risk process with dividend barrier b where the claims arrive according to a Markovian additive process (MAP). For spectrally negative MAPs, we present linear equations for the expected discounted dividends and the expected discounted penalty function. We apply results for the first exit times of spectrally negative Lévy processes and change-of-measure techniques. Explicit expressions are given when there are positive and negative claims, with phase-type distribution.

We consider a Markov-modulated Brownian motion (MMBM) with phase-dependent termination rates, i.e. while in a phase i the process terminates with a constant hazard rate ri ≥ 0. For such a process, we determine the matrix of expected local times (at zero) before termination and hence the resolvent. The results are applied to some recent questions arising in the framework of insurance risk. We further provide expressions for the resolvent and the local times at zero of an MMBM reflected at its infimum.

Let (X, J) denote a Markov-modulated Brownian motion (MMBM) and denote its supremum process by S. For some a > 0, let σ(a) denote the time when the reflected process Y := S - X first surpasses the level a. Furthermore, let σ−(a) denote the last time before σ(a) when X attains its current supremum. In this paper we shall derive the joint distribution of Sσ(a), σ−(a), and σ(a), where the latter two will be given in terms of their Laplace transforms. We also provide some remarks on scale matrices for MMBMs with strictly positive variation parameters. This extends recent results for spectrally negative Lévy processes to MMBMs. Due to well-known fluid embedding and state-dependent killing techniques, the analysis applies to Markov additive processes with phase-type jumps as well. The result is of interest to applications such as the dividend problem in insurance mathematics and the buffer overflow problem in queueing theory. Examples will be given for the former.

In this paper we determine the distributions of occupation times of a Markov-modulated Brownian motion (MMBM) in separate intervals before a first passage time or an exit from an interval. We derive the distributions in terms of their Laplace transforms, and we also distinguish between occupation times in different phases. For MMBMs with strictly positive variation parameters, we further propose scale functions.

In this paper we are interested in bounding or calculating the additive functionals of the first return time on a set for discrete-time Markov chains on a countable state space, which is motivated by investigating ergodic theory and central limit theorems. To do so, we introduce the theory of the minimal nonnegative solution. This theory combined with some other techniques is proved useful for investigating the additive functionals. This method is used to study the functionals for discrete-time birth-death processes, and the polynomial convergence and a central limit theorem are derived.

We describe the processes obtained by time reversal of a class of stationary jump diffusion processes that model the dynamics of genetic variation in populations subject to repeated bottlenecks. Assuming that only one lineage survives each bottleneck, the forward process is a diffusion on [0,1] that jumps to the boundary before diffusing back into the interior. We show that the behavior of the time-reversed process depends on whether the boundaries are accessible to the diffusive motion of the forward process. If a boundary point is inaccessible to the forward diffusion then time reversal leads to a jump diffusion that jumps immediately into the interior whenever it arrives at that point. If, instead, a boundary point is accessible then the jumps off of that point are governed by a weighted local time of the time-reversed process.

We consider a Markov-modulated Brownian motion reflected to stay in a strip [0, B]. The stationary distribution of this process is known to have a simple form under some assumptions. We provide a short probabilistic argument leading to this result and explain its simplicity. Moreover, this argument allows for generalizations including the distribution of the reflected process at an independent, exponentially distributed epoch. Our second contribution concerns transient behavior of the model. We identify the joint law of the processes defining the model at inverse local times.

In the spirit of Albrecher and Hipp (2007), and Albrecher, Renaud, and Zhou (2008) we consider a Lévy insurance risk model with tax payments of a more general structure than in the aforementioned papers, which was also considered in Albrecher, Borst, Boxma, and Resing (2009). In terms of scale functions, we establish three fundamental identities of interest which have stimulated a large volume of actuarial research in recent years. That is to say, the two-sided exit problem, the net present value of tax paid until ruin, as well as a generalized version of the Gerber–Shiu function. The method we appeal to differs from Albrecher and Hipp (2007), and Albrecher, Renaud, and Zhou (2008) in that we appeal predominantly to excursion theory.

In this paper we study the integral of the supremum process of standard Brownian motion. We present an explicit formula for the moments of the integral (or area) (T) covered by the process in the time interval [0,T]. The Laplace transform of (T) follows as a consequence. The main proof involves a double Laplace transform of (T) and is based on excursion theory and local time for Brownian motion.

The present paper generalises some results for spectrally negative Lévy processes to the setting of Markov additive processes (MAPs). A prominent role is assumed by the first passage times, which will be determined in terms of their Laplace transforms. These have the form of a phase-type distribution, with a rate matrix that can be regarded as an inverse function of the cumulant matrix. A numerically stable iteration to compute this matrix is given. The theory is first developed for MAPs without positive jumps and then extended to include positive jumps having phase-type distributions. Numerical and analytical examples show agreement with existing results in special cases.

In this paper we consider reflected diffusions with positive and negative jumps, constrained to lie in the nonnegative orthant of ℝn. We allow for the drift and diffusion coefficients, as well as for the directions of reflection, to be random fields over time and space. We provide a boundary behavior characterization, generalizing known results in the nonrandom coefficients and constant directions of the reflection case. In particular, the regulator processes are related to semimartingale local times at the boundaries, and they are shown not to charge the times the process expends at the intersection of boundary faces. Using the boundary results, we extend the conditions for product-form distributions in the stationary regime to the case when the drift and diffusion coefficients, as well as the directions of reflection, are random fields over space.

We provide a distributional study of the solution to the classical control problem due to De Finetti (1957), Gerber (1969), Azcue and Muler (2005), and Avram et al. (2007), which concerns the optimal payment of dividends from an insurance risk process prior to ruin. Specifically, we build on recent work in the actuarial literature concerning calculations of the nth moment of the net present value of dividends paid out in the optimal strategy as well as the moments of the deficit at ruin and the Laplace transform of the red period. The calculations we present go much further than the existing literature, in that our calculations are valid for a general spectrally negative Lévy process as opposed to the classical Cramér–Lundberg process with exponentially distributed jumps. Moreover, the technique we use appeals principally to excursion theory rather than integro-differential equations and, for the case of the nth moment of the net present value of dividends, makes a new link with the distribution of integrated exponential subordinators.

In this short note we show how new fluctuation identities and their associated asymptotics, given in Vigon (2002), Klüppelberg et al. (2004) and Doney and Kyprianou (2006), provide the basis for establishing, in an elementary way, asymptotic overshoot and undershoot distribitions for a general class of Lévy insurance risk processes. The results bring the earlier conclusions of Asmussen and Klüppelberg (1996) for the Cramér-Lundberg process into greater generality.