We study a class of tenable, irreducible, nondegenerate zero-balanced Pólya urn schemes. We give a full characterization of the class by sufficient and necessary conditions. Only forms with a certain cyclic structure in their replacement matrix are admissible. The scheme has a steady state into proportions governed by the principal (left) eigenvector of the average replacement matrix. We study the gradual change for any such urn containing n → ∞ balls from the initial condition to the steady state. We look at the status of an urn starting with an asymptotically positive proportion of each color after j n draws. We identify three phases of j n : the growing sublinear, the linear, and the superlinear. In the growing sublinear phase the number of balls of different colors has an asymptotic joint multivariate normal distribution, with mean and covariance structure that are influenced by the initial conditions. In the linear phase a different multivariate normal distribution kicks in, in which the influence of the initial conditions is attenuated. The steady state is not a good approximation until a certain superlinear amount of time has elapsed. We give interpretations for how the results in different phases conjoin at the ‘seam lines’. In fact, these Gaussian phases are all manifestations of one master theorem. The results are obtained via multivariate martingale theory. We conclude with some illustrating examples.

Nearly fifty years after the introduction of skew Brownian motion by Itô and McKean (1963), the first passage time distribution remains unknown. In this paper we first generalize results of Pitman and Yor (2011) and Csáki and Hu (2004) to derive formulae for the distribution of ranked excursion heights of skew Brownian motion, and then use these results to derive the first passage time distribution.

The goal is to identify the class of distributions to which the distribution of the maximum of a Lévy process with no negative jumps and negative mean (equivalently, the stationary distribution of the reflected process) belongs. An explicit new distributional identity is obtained for the case where the Lévy process is an independent sum of a Brownian motion and a general subordinator (nondecreasing Lévy process) in terms of a geometrically distributed sum of independent random variables. This generalizes both the distributional form of the standard Pollaczek-Khinchine formula for the stationary workload distribution in the M/G/1 queue and the exponential stationary distribution of a reflected Brownian motion.

Nonnegative infinitely divisible (i.d.) random variables form an important class of random variables. However, when this type of random variable is specified via Lévy densities that have infinite integrals on (0, ∞), except for some special cases, exact sampling is unknown. We present a method that can sample a rather wide range of such i.d. random variables. A basic result is that, for any nonnegative i.d. random variable X with its Lévy density explicitly specified, if its distribution conditional on X ≤ r can be sampled exactly, where r > 0 is any fixed number, then X can be sampled exactly using rejection sampling, without knowing the explicit expression of the density of X. We show that variations of the result can be used to sample various nonnegative i.d. random variables.

In this paper we show that fractional Brownian motion with H < ½ can arise as a limit of a simple class of traffic processes that we call ‘scheduled traffic models’. To our knowledge, this paper provides the first simple traffic model leading to fractional Brownnian motion with H < ½. We also discuss some immediate implications of this result for queues fed by scheduled traffic, including a heavy-traffic limit theorem.

In this paper we consider three-dimensional random tessellations that are stable under iteration (STIT tessellations). STIT tessellations arise as a result of subsequent cell division, which implies that their cells are not face-to-face. The edges of the cell-dividing polygons are the so-called I-segments of the tessellation. The main result is an explicit formula for the distribution of the number of vertices in the relative interior of the typical I-segment. In preparation for its proof, we obtain other distributional identities for the typical I-segment and the length-weighted typical I-segment, which provide new insight into the spatiotemporal construction process.

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 consider the statistical concept of causality in continuous time between filtered probability spaces, based on Granger’s definitions of causality. Then we consider some stable subspaces of $H^p$ which contain right continuous modifications of martingales $P(A \mid {\mathcal {G}}_t)$. We give necessary and sufficient conditions, in terms of statistical causality, for these spaces to coincide with $H^p$. These results can be applied to extremal measures and regular weak solutions of stochastic differential equations.

This paper is devoted to the perfect simulation of a stationary process with an at most countable state space. The process is specified through a kernel, prescribing the probability of the next state conditional to the whole past history. We follow the seminal work of Comets, Fernández and Ferrari (2002), who gave sufficient conditions for the construction of a perfect simulation algorithm. We define backward coalescence times for these kind of processes, which allow us to construct perfect simulation algorithms under weaker conditions than in Comets, Fernández and Ferrari (2002). We discuss how to construct backward coalescence times (i) by means of information depths, taking into account some a priori knowledge about the histories that occur; and (ii) by identifying suitable coalescing events.

In this paper we consider optimal stopping problems for a general class of reward functions under matrix-exponential jump-diffusion processes. Given an American call-type reward function in this class, following the averaging problem approach (see, for example, Alili and Kyprianou (2005), Kyprianou and Surya (2005), Novikov and Shiryaev (2007), and Surya (2007)), we give an explicit formula for solutions of the corresponding averaging problem. Based on this explicit formula, we obtain the optimal level and the value function for American call-type optimal stopping problems.

Dependence between coordinate extremes is a key factor in any multivariate risk assessment. Hence, it is of interest to know whether the components of a given multivariate random vector exhibit asymptotic independence or asymptotic dependence. In the latter case the structure of the asymptotic dependence has to be clarified. In the multivariate setting it is common to have an explicit form of the density rather than the distribution function. In this paper we therefore give criteria for asymptotic dependence in terms of the density. We consider distributions with light tails and restrict attention to continuous unimodal densities defined on the whole space or on an open convex cone. For simplicity, the density is assumed to be homothetic: all level sets have the same shape. Balkema and Nolde (2010) contains conditions on the shape which guarantee asymptotic independence. The situation for asymptotic dependence, treated in the present paper, is more delicate.

We present the construction of a continuous-time stochastic process which has moments that satisfy an exact scaling relation, including odd-order moments. It is based on a natural extension of the multifractal random walk construction described in Bacry and Muzy (2003). This allows us to propose a continuous-time model for the price of a financial asset that reflects most major stylized facts observed on real data, including asymmetry and multifractal scaling.

In this paper we analyze different forms of fractional relaxation equations of order ν ∈ (0, 1), and we derive their solutions in both analytical and probabilistic forms. In particular, we show that these solutions can be expressed as random boundary crossing probabilities of various types of stochastic process, which are all related to the Brownian motion B. In the special case ν = ½, the fractional relaxation is shown to coincide with Pr{sup0≤s≤t B(s) < U} for an exponential boundary U. When we generalize the distributions of the random boundary, passing from the exponential to the gamma density, we obtain more and more complicated fractional equations.

Let X = {X t : t ≥ 0} be a stationary piecewise continuous R d -valued process that moves between jumps along the integral curves of a given continuous vector field, and let S ⊂ R d be a smooth surface. The aim of this paper is to derive a multivariate version of Rice's formula, relating the intensity of the point process of (localized) continuous crossings of S by X to the distribution of X 0. Our result is illustrated by examples relating to queueing networks and stress release network models.

The ‘coupon collection problem’ refers to a class of occupancy problems in which j identical items are distributed, independently and at random, to n cells, with no restrictions on multiple occupancy. Identifying the cells as coupons, a coupon is ‘collected’ if the cell is occupied by one or more of the distributed items; thus, some coupons may never be collected, whereas others may be collected once or twice or more. We call the number of coupons collected exactly r times coupons of type r. The coupon collection model we consider is general, in that a random number of purchases occurs at each stage of collecting a large number of coupons; the sample sizes at each stage are independent and identically distributed according to a sampling distribution. The joint behavior of the various types is an intricate problem. In fact, there is a variety of joint central limit theorems (and other limit laws) that arise according to the interrelation between the mean, variance, and range of the sampling distribution, and of course the phase (how far we are in the collection processes). According to an appropriate combination of the mean of the sampling distribution and the number of available coupons, the phase is sublinear, linear, or superlinear. In the sublinear phase, the normalization that produces a Gaussian limit law for uncollected coupons can be used to obtain a multivariate central limit law for at most two other types — depending on the rates of growth of the mean and variance of the sampling distribution, we may have a joint central limit theorem between types 0 and 1, or between types 0, 1, and 2. In the linear phase we have a multivariate central limit theorem among the types 0, 1,…, k for any fixed k.

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.

We construct a process with inverse gamma increments and an asymptotically self-similar limit. This construction supports the use of long-range-dependent t subordinator models for actual financial data as advocated in Heyde and Leonenko (2005), in that it allows for noninteger-valued model parameters, as is found empirically in model estimation from data.

We study the dynamic assignment of cross-trained workers in serial production lines characterized by stochastic process times and inventory buffers between stations. Throughput maximization is the objective. Each worker is trained for a subset of tasks, where emphasis is placed on systems with each worker trained for a zone of stations with stations near the zone boundaries being served (shared) by one or more other workers as well. Using sample path comparisons, we identify structural properties of optimal worker allocation policies. We identify when (i) a worker can prioritize the job in the most downstream station (last-buffer-first-served), and (ii) only the downstream (as opposed to upstream) server should serve a single task.

We study the smooth-fit property of the American put price with finite maturity in an exponential Lévy model when the underlying stock pays dividends at a continuous rate. As in the perpetual case, a regularity property is sufficient for smooth fit to occur. We also derive conditions on the Lévy measure under which smooth fit fails.

We give a simple method to approximate multidimensional exponentially tempered stable processes and show that the approximating process converges in the Skorokhod topology to the tempered process. The approximation is based on the generation of a random angle and a random variable with a lower-dimensional Lévy measure. We then show that if an arbitrarily small normal random variable is added, the marginal densities converge uniformly at an almost linear rate, providing a critical tool to assess the performance of existing particle tracking codes.