We consider a nearly unstable, or near unit root, AR(1) process with regularly varying innovations. Two different approximations for the stationary distribution of such processes exist: a Gaussian approximation arising from the nearly unstable nature of the process and a heavy-tail approximation related to the tail asymptotics of the innovations. We combine these two approximations to obtain a new uniform approximation that is valid on the entire real line. As a corollary, we obtain a precise description of the regions where each of the Gaussian and heavy-tail approximations should be used.

A new approach to the solution of optimal stopping problems for one-dimensional diffusions is developed. It arises by imbedding the stochastic problem in a linear programming problem over a space of measures. Optimizing over a smaller class of stopping rules provides a lower bound on the value of the original problem. Then the weak duality of a restricted form of the dual linear program provides an upper bound on the value. An explicit formula for the reward earned using a two-point hitting time stopping rule allows us to prove strong duality between these problems and, therefore, allows us to either optimize over these simpler stopping rules or to solve the restricted dual program. Each optimization problem is parameterized by the initial value of the diffusion and, thus, we are able to construct the value function by solving the family of optimization problems. This methodology requires little regularity of the terminal reward function. When the reward function is smooth, the optimal stopping locations are shown to satisfy the smooth pasting principle. The procedure is illustrated using two examples.

Consider independent fair coin flips at each site of the lattice ℤd . A translation-equivariant matching rule is a perfect matching of heads to tails that commutes with translations of ℤd and is given by a deterministic function of the coin flips. Let Z Φ be the distance from the origin to its partner, under the translation-equivariant matching rule Φ. Holroyd and Peres (2005) asked, what is the optimal tail behaviour of Z Φ for translation-equivariant perfect matching rules? We prove that, for every d ≥ 2, there exists a translation-equivariant perfect matching rule Φ such that EZ Φ 2/3-ε < ∞ for every ε > 0.

Reflected Brownian motion is used in areas such as physiology, electrochemistry and nuclear magnetic resonance. We study the first-passage-time problem of this process which is relevant in applications; specifically, we find a Volterra integral equation for the distribution of the first time that a reflected Brownian motion reaches a nondecreasing barrier. Additionally, we note how a numerical procedure can be used to solve the integral equation.

In this paper we present an explicit solution to the infinite-horizon optimal stopping problem for processes with stationary independent increments, where reward functions admit a certain representation in terms of the process at a random time. It is shown that it is optimal to stop at the first time the process crosses a level defined as the root of an equation obtained from the representation of the reward function. We obtain an explicit formula for the value function in terms of the infimum and supremum of the process, by making use of the Wiener–Hopf factorization. The main results are applied to several problems considered in the literature, to give a unified approach, and to new optimization problems from the finance industry.

For a family of random walks {S (a)} satisfying E S 1 (a)=-a<0, we consider ladder epochs τ(a)=min {k≥1: S k (a)<0}. We study the asymptotic behaviour, as a⇒0, of P (τ(a)>n) in the case when n=n(a)→∞. As a consequence, we also obtain the growth rates of the moments of τ(a).

We consider a stationary Poisson process X of k-flats in ℝd with intensity measure Θ and a measurable set S of k-flats depending on F 1,…,F n ∈ X, x∈ℝd, and X in a specific equivariant way. If (F 1,…,F n ,x) is properly sampled (in a ‘typical way’) then Θ(S) has a gamma distribution. This result generalizes and unifies earlier work by Miles (1971), Møller and Zuyev (1996), and Zuyev (1999). As a new example, we will show that the volume of the fundamental region of a typical j-face of a stationary Poisson–Voronoi tessellation is conditionally gamma distributed. This is true in the area-biased and the area-debiased cases. In the first case the shape parameter is not integer valued. As another new example, we will show that the generalized integral-geometric contents of the (area-biased and area-debiased) typical j-face of a Poisson hyperplane tessellation are conditionally gamma distributed. In the isotropic case the contents boil down to the mean breadth of the face.

We approximate the distribution of the sum of independent but not necessarily identically distributed Bernoulli random variables using a shifted binomial distribution, where the three parameters (the number of trials, the probability of success, and the shift amount) are chosen to match the first three moments of the two distributions. We give a bound on the approximation error in terms of the total variation metric using Stein's method. A numerical study is discussed that shows shifted binomial approximations are typically more accurate than Poisson or standard binomial approximations. The application of the approximation to solving a problem arising in Bayesian hierarchical modeling is also discussed.

We consider the following two definitions of discounting: (i) multiplicative coefficient in front of the rewards, and (ii) probability that the process has not been stopped if the stopping time has an exponential distribution independent of the process. It is well known that the expected total discounted rewards corresponding to these definitions are the same. In this note we show that, the variance of the total discounted rewards is smaller for the first definition than for the second definition.

The problem we consider here is a full-information best-choice problem in which n applicants appear sequentially, but each applicant refuses an offer independently of other applicants with known fixed probability 0≤q<1. The objective is to maximize the probability of choosing the best available applicant. Two models are distinguished according to when the availability can be ascertained; the availability is ascertained just after the arrival of the applicant (Model 1), whereas the availability can be ascertained only when an offer is made (Model 2). For Model 1, we can obtain the explicit expressions for the optimal stopping rule and the optimal probability for a given n. A remarkable feature of this model is that, asymptotically (i.e. n→∞), the optimal probability becomes insensitive to q and approaches 0.580 164. The planar Poisson process (PPP) model provides more insight into this phenomenon. For Model 2, the optimal stopping rule depends on the past history in a complicated way and seems to be intractable. We have not solved this model for a finite n but derive, via the PPP approach, a lower bound on the asymptotically optimal probability.

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 Michel and Paroux (2003) the authors proposed a new proof of a well-known convergence result for the scaled elementary connected vacant component in the high intensity Boolean model towards the Crofton cell of the Poisson hyperplane process (see, e.g. Hall (1985)). In this paper we investigate the second-order term in this convergence when the two-dimensional Boolean model and the Poisson line process are coupled on the same probability space. We consider the particular case where the grains are discs with random radii. A precise coupling between the Boolean model and the Poisson line process is first established. A result of directional convergence in distribution for the difference of the two sets involved is then derived. Eventually, we show the convergence of the process, measuring the difference between the two random sets, once rescaled, as a function of the direction.

Using key tools such as Itô's formula for general semimartingales, Kunita's moment estimates for Lévy-type stochastic integrals, and the exponential martingale inequality, we find conditions under which the solutions to the stochastic differential equations (SDEs) driven by Lévy noise are stable in probability, almost surely and moment exponentially stable.

Zolotarev (1961) proved a duality result that relates stable densities with different indices. In this paper we show how Zolotarev's duality leads to some interesting results on fractional diffusion. Fractional diffusion equations employ fractional derivatives in place of the usual integer-order derivatives. They govern scaling limits of random walk models, with power-law jumps leading to fractional derivatives in space, and power-law waiting times between the jumps leading to fractional derivatives in time. The limit process is a stable Lévy motion that models the jumps, subordinated to an inverse stable process that models the waiting times. Using duality, we relate the density of a spectrally negative stable process with index 1<α<2 to the density of the hitting time of a stable subordinator with index 1/α, and thereby unify some recent results in the literature. These results provide a concrete interpretation of Zolotarev's duality in terms of the fractional diffusion model. They also illuminate a current controversy in hydrology, regarding the appropriate use of space- and time-fractional derivatives to model contaminant transport in river flows.

We derive a moderate deviation principle for word counts (which is extended to counts of multiple patterns) in biological sequences under different models: independent and identically distributed letters, homogeneous Markov chains of order 1 and m, and, in view of the codon structure of DNA sequences, Markov chains with three different transition matrices. This enables us to approximate P-values for the number of word occurrences in DNA and protein sequences in a new manner.

The goal of this paper is to provide formulae for the expectation and variance of the height and length of the ancestral recombination graph (ARG). While the formula for the expectation of the height is known (see, e.g. Krone and Neuhauser (1997)), the other formulae seem to be new. We obtain in particular (see Theorem 4.1) a very simple formula which expresses the expectation of the length of the ARG as a linear combination of the expectations of both the length of the coalescent tree and the height of the ARG. Finally, we study the speed at which the ARG comes down from infinity.

Let X be a pure-jump subordinator (i.e. nondecreasing Lévy process with no drift) with infinite Lévy measure, let X ε be the sum of jumps not exceeding ε, and let µ(ε)=E[X ε(1)]. We study the question of weak convergence of X ε/µ(ε) as ε ↓0, in terms of the limit behavior of µ(ε)/ε. The most interesting case reduces to the weak convergence of X ε/ε to a subordinator whose marginals are generalized Dickman distributions; we give some necessary and sufficient conditions for this to hold. For a certain significant class of subordinators for which the latter convergence holds, and whose most prominent representative is the gamma process, we give some detailed analysis regarding the convergence quality (in particular, in the context of approximating X itself). This paper completes, in some respects, the study made by Asmussen and Rosiński (2001).

In this paper we deal with generalized fractional kinetic equations driven by a Gaussian noise, white in time and correlated in space, and where the diffusion operator is the composition of the Bessel and Riesz potentials for any fractional parameters. We give results on the existence and uniqueness of solutions by means of a weak formulation and study the Hölder continuity. Moreover, we prove the existence of a smooth density associated to the solution process and study the asymptotics of this density. Finally, when the diffusion coefficient is constant, we look for its Gaussian index.

We study the asymptotic behavior of the tail probability of integrated stable processes exceeding power barriers. In the first part of the paper the limiting behavior of the integrals of stable processes generated by ergodic dissipative flows is established. In the second part an example with the integral of a stable process generated by a conservative flow is analyzed. Finally, the difference in the order of magnitude of the exceedance probability in the two cases is related to the dependence structure of the underlying stable process.

We introduce a class of stock models that interpolates between exponential Lévy models based on Brownian subordination and certain stochastic volatility models with Lévy-driven volatility, such as the Barndorff-Nielsen–Shephard model. The driving process in our model is a Brownian motion subordinated to a business time which is obtained by convolution of a Lévy subordinator with a deterministic kernel. We motivate several choices of the kernel that lead to volatility clusters while maintaining the sudden extreme movements of the stock. Moreover, we discuss some statistical and path properties of the models, prove absence of arbitrage and incompleteness, and explain how to price vanilla options by simulation and fast Fourier transform methods.