Search

In this paper we investigate the limiting behavior of the failure rate for the convolution of two or more life distributions. In a previous paper on mixtures, Block, Mi and Savits (1993) showed that the failure rate behaves like the limiting behavior of the strongest component. We show a similar result here for convolutions. We also show by example that unlike a mixture population, the ultimate direction of monotonicity does not necessarily follow that of the strongest component.

In this paper, we describe a class of Wiener functionals that are ‘indeterminate by their moments’, that is, whose distributions are not uniquely determined by their moments. In particular, it is proved that the integral of a geometric Brownian motion is indeterminate by its moments and, moreover, shown that previous proofs of this result are incorrect. The main result of this paper is based on geometric inequalities in Gauss space and on a generalization of the Krein criterion due to H. L. Pedersen.

On the basis of the Karlin and McGregor result, which states that the transition probability functions of a birth and death process can be expressed via the introduction of an orthogonal polynomial system and a spectral measure, we investigate in this paper how the Laplace transforms and the distributions of different transient characteristics related to excursions of a birth and death process can be expressed by means of the basic orthogonal polynomial system and the spectral measure. This allows us in particular to give a probabilistic interpretation of the series introduced by Stieltjes to study the convergence of the fundamental continued fraction associated with the system. Throughout the paper, we pay special attention to the case when the birth and death process is ergodic. Under the assumption that the spectrum of the spectral measure is discrete, we show how the distributions of different random variables associated with excursions depend on the fundamental continued fraction, the orthogonal polynomial system and the spectral measure.

In this paper we analyze a single server two-queue model with Bernoulli schedules. This discipline is very flexible and contains the exhaustive and 1-limited disciplines as special cases. We formulate the queueing system as a Riemann boundary value problem with shift. The boundary value problem is solved by exploring a Fredholm integral equation around the unit circle. Some numerical examples are presented at the end of the paper.

This paper is concerned with the distibutions and the moments of the area and the local time of a random walk, called the Bernoulli meander. The limit behavior of the distributions and the moments is determined in the case where the number of steps in the random walk tends to infinity. The results of this paper yield explicit formulas for the distributions and the moments of the area and the local time for the Brownian meander.

In this paper we address the question of finding the point which maximizes the pth moment of the exit time of planar Brownian motion from a given domain. We present a geometrical method for excluding parts of the domain from consideration which makes use of a coupling argument and the conformal invariance of Brownian motion. In many cases the maximizing point can be localized to a relatively small region. Several illustrative examples are presented.

Goldman (2010) proved that the distribution of a stationary determinantal point process (DPP) Φ can be coupled with its reduced Palm version Φ0,! such that there exists a point process η where Φ=Φ0,!∪η in distribution and Φ0,!∩η=∅. The points of η characterize the repulsive nature of a typical point of Φ. In this paper we use the first-moment measure of η to study the repulsive behavior of DPPs in high dimensions. We show that many families of DPPs have the property that the total number of points in η converges in probability to 0 as the space dimension n→∞. We also prove that for some DPPs, there exists an R∗ such that the decay of the first-moment measure of η is slowest in a small annulus around the sphere of radius √nR∗. This R∗ can be interpreted as the asymptotic reach of repulsion of the DPP. Examples of classes of DPP models exhibiting this behavior are presented and an application to high-dimensional Boolean models is given.

Kac's formula for Brownian functionals and Levy's local time decomposition are shown to be useful tools in analysing Brownian excursion properties. These tools are applied to maximum, local time and area distributions. Some curious connections between some of these distributions are explained by simple The original purpose of this paper was to find a workable expression for the transform of the Brownian excursion area. We wanted to use two well-known results: Kac's formula for Brownian functionals and Levy's local time decompos- ition. We actually found that these two powerful tools lead also to alternative and sometimes simpler proofs for other results on Brownian excursion prob- abilities: maximum and local time distributions. We were also able to derive a useful general result on Brownian excursion symmetric additive functionals. Finally, we also wanted to understand some curious relations that had been observed by earlier authors: connections between hitting-time densities and maximum distributions. Using simple probabilistic arguments and a formula for Jacobi's third 0 function, we can explain these connections in direct terms. The paper is organized as follows. In Section 1, we summarize basic notations and known results. Sections 2 and 3 are short surveys of Kac's and Levy's results. These tools are applied in Sections 4 and 5 to Brownian excursion maximum and local time. Section 6 contains the main results of the paper: a useful new form for the transform of the Brownian excursion area density and a simple relation for functionals based on a symmetric positive function. Section 7 explains those curious relations alluded to earlier between Brownian excursion probabilities. 1. Basic notations and known results

This paper deals with the existence of solutions of a stochastic integral equation of the Volterra type and their asymptotic behaviour. Investigations of this paper use the concept of a measure of non-compactness in Banach space and fixed-point theorem of Darbo type. An application to a stochastic model for chemotherapy is also presented.

Consider an absolutely continuous distribution on [0, ∞) with finite mean μ and hazard rate function h(t) ≤ b for all t. For bμ close to 1, we would expect F to be approximately exponential. In this paper we obtain sharp bounds for the Kolmogorov distance between F and an exponential distribution with mean μ, as well as between F and an exponential distribution with failure rate b. We apply these bounds to several examples. Applications are presented to geometric convolutions, birth and death processes, first-passage times, and to decreasing mean residual life distributions.

In this paper we use the Siegert formula to derive alternative expressions for the moments of the first passage time of the Ornstein-Uhlenbeck process through a constant threshold. The expression for the nth moment is recursively linked to the lower-order moments and consists of only n terms. These compact expressions can substantially facilitate (numerical) applications also for higher-order moments.

In this paper three models are considered. They are the infinitely-many-neutral-alleles model of Ethier and Kurtz (1981), the two-parameter infinitely-many-alleles diffusion model of Petrov (2009), and the infinitely-many-alleles model with symmetric dominance Ethier and Kurtz (1998). New representations of the transition densities are obtained for the first two models and the ergodic inequalities are provided for all three models.

In this paper we consider a one dimensional stochastic system described by an elliptic equation. A spatially varying random coefficient is introduced to account for uncertainty or imprecise measurements. We model the logarithm of this coefficient by a Gaussian process and provide asymptotic approximations of the tail probabilities of the derivative of the solution.

In this paper we discuss an exponential integrator scheme, based on spatial discretization and time discretization, for a class of stochastic partial differential equations. We show that the scheme has a unique stationary distribution whenever the step size is sufficiently small, and that the weak limit of the stationary distribution of the scheme as the step size tends to 0 is in fact the stationary distribution of the corresponding stochastic partial differential equations.

In this paper we discuss a nondecreasing lower bound for the Poisson cumulative distribution function (CDF) at z standard deviations above the mean λ, where z and λ are parameters. This is important because the normal distribution as an approximation for the Poisson CDF may overestimate or underestimate its value. A sharp nondecreasing lower bound in the form of a step function is constructed. As a corollary of the bound's properties, for a given percent α and parameter λ, the minimal z is obtained such that, for any Poisson random variable with the mean greater or equal to λ, its αth percentile is at most z standard deviations above its mean. For Poisson distributed control parameters, the corollary allows simple policies measuring performance in terms of standard deviations from a benchmark.

This paper provides an analysis on a discrete version of the Ornstein-Uhlenbeck (OU) process which reflects the small discrete movements caused by the tick size effect. This discrete OU process is derived from matching the first two moments to those of the standard OU process in an infinitesimal sense. We discuss the distributional convergence from the discrete to the continuous processes, and show that the convergence speed is in the second order of the step (tick) size. We also provide some analytical results for the proposed discrete OU process itself, including the closed-form formula of the moment generating function and a full characterisation of the steady state distribution. These results enable us to examine the convergence order explicitly.

In this paper we explore the first passage times of constant-elasticity-of-variance (CEV) processes with two-sided reflecting barriers. The explicit Laplace transforms of the first passage times are derived. Our results can include analytic formulae concerning Laplace transforms of first passage times of reflected Ornstein–Uhlenbeck processes, reflected geometric Brownian motions, and reflected square-root processes.

In this paper we study the functional central limit theorem (CLT) for stationary Markov chains with a self-adjoint operator and general state space. We investigate the case when the variance of the partial sum is not asymptotically linear in n, and establish that conditional convergence in distribution of partial sums implies the functional CLT. The main tools are maximal inequalities that are further exploited to derive conditions for tightness and convergence to the Brownian motion.

The risk processes considered in this paper are generated by an underlying Markov process with a regenerative structure and an independent sequence of independent and identically distributed claims. Between the arrivals of claims the process increases at a rate which is a nonnegative function of the present value of the Markov process. The intensity for a claim to occur is another nonnegative function of the value of the Markov process. The claim arrival times are the regeneration times for the Markov process. Two-sided claims are allowed, but the distribution of the positive claims is assumed to have a Laplace transform that is a rational function. The main results describe the joint Laplace transform of the time at ruin and the deficit at ruin. The method used consists in finding partial eigenfunctions for the generator of the joint process consisting of the Markov process and the accumulated claims process, a joint process which is also Markov. These partial eigenfunctions are then used to find a martingale that directly leads to an expression for the desired Laplace transform. In the final section, three examples are given involving different types of the underlying Markov process.

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.

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.

The theory of L 2-spectral gaps for reversible Markov chains has been studied by many authors. In this paper we consider positive recurrent general state space Markov chains with stationary transition probabilities. Replacing the assumption of reversibility with a weaker assumption, we still obtain a simple necessary and sufficient condition for the spectral gap property of the associated Markov operator in terms of the isoperimetric constant. We show that this result can be applied to a large class of Markov chains, including those that are related to positive recurrent finite-range random walks on Z.

In this paper we consider a general Lévy process X reflected at a downward periodic barrier At and a constant upper barrier K, giving a process VKt=Xt +LAt−LKt. We find the expression for a loss rate defined by lK=ELK1 and identify its asymptotics as K→∞ when X has light-tailed jumps and EX1

A simple model for a randomly oscillating variable is suggested, which is a variant of the two-state random velocity model. As in the latter model, the variable keeps rising or falling with constant velocity for some time before randomly reversing its direction. In contrast however, its propensity to reverse depends on its current value and it is for this desirable feature that the model is proposed here. This feature has two implications: (a) neither the changing variable nor its velocity is Markovian, although the joint process is, and (b) the linear differential equations arising in the case of our model do not have constant coefficients. The results given in this paper are meant to illustrate the straightforward nature of some of the calculations involved and to highlight the relationship with one-dimensional diffusions.

In this paper we study the Wiener-Hopf factorization for a class of L\'evy processes with double-sided jumps, characterized by the fact that the density of the L\'evy measure is given by an infinite series of exponential functions with positive coefficients. We express the Wiener-Hopf factors as infinite products over roots of a certain transcendental equation, and provide a series representation for the distribution of the supremum/infimum process evaluated at an independent exponential time. We also introduce five eight-parameter families of L\'evy processes, defined by the fact that the density of the L\'evy measure is a (fractional) derivative of the theta-function, and we show that these processes can have a wide range of behavior of small jumps. These families of processes are of particular interest for applications, since the characteristic exponent has a simple expression, which allows efficient numerical computation of the Wiener-Hopf factors and distributions of various functionals of the process., Comment: 12 pages, published online at http://projecteuclid.org/euclid.jap/1294170516

A sequence of random variables is said to be extended negatively dependent (END) if the tails of its finite-dimensional distributions in the lower-left and upper-right corners are dominated by a multiple of the tails of the corresponding finite-dimensional distributions of a sequence of independent random variables with the same marginal distributions. The goal of this paper is to establish the strong law of large numbers for a sequence of END and identically distributed random variables. In doing so we derive some new inequalities of large deviation type for the sums of END and identically distributed random variables being suitably truncated. We also show applications of our main result to risk theory and renewal theory.

This paper is concerned with properties of the algebraic degree of the Laplace-Stieltjes transform of phase-type (PH) distributions. The main problem of interest is: given a PH generator, how do we find the maximum and the minimum algebraic degrees of all irreducible PH representations with that PH generator? Based on the matrix exponential (ME) order of ME distributions and the spectral polynomial algorithm, a method for computing the algebraic degree of a PH distribution is developed. The maximum algebraic degree is identified explicitly. Using Perron-Frobenius theory of nonnegative matrices, a lower bound and an upper bound on the minimum algebraic degree are found, subject to some conditions. Explicit results are obtained for special cases.

In this paper, a stochastic integral of Ornstein–Uhlenbeck type is represented to be the sum of two independent random variables: one has a tempered stable distribution and the other has a compound Poisson distribution. In distribution, the compound Poisson random variable is equal to the sum of a Poisson-distributed number of positive random variables, which are independent and identically distributed and have a common specified density function. Based on the representation of the stochastic integral, we prove that the transition distribution of the tempered stable Ornstein–Uhlenbeck process is self-decomposable and that the transition density is a C∞-function.

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 consider two classes of reflected Ornstein–Uhlenbeck (OU) processes: the reflected OU process with jumps and the Markov-modulated reflected OU process. We prove that their stationary distributions exist. Furthermore, for the jump reflected OU process, the Laplace transform (LT) of the stationary distribution is given. As for the Markov-modulated reflected OU process, we derive an equation satisfied by the LT of the stationary distribution.

In this paper we prove that a sequence of scaled generalized Jiřina processes can converge weakly to a nonlinear diffusion process with Lévy jumps under certain conditions.

In this paper we generalize existing results for the steady-state distribution of growth-collapse processes. We begin with a stationary setup with some relatively general growth process and observe that, under certain expected conditions, point- and time-stationary versions of the processes exist as well as a limiting distribution for these processes which is independent of initial conditions and necessarily has the marginal distribution of the stationary version. We then specialize to the cases where an independent and identically distributed (i.i.d.) structure holds and where the growth process is a nondecreasing Lévy process, and in particular linear, and the times between collapses form an i.i.d. sequence. Known results can be seen as special cases, for example, when the inter-collapse times form a Poisson process or when the collapse ratio is deterministic. Finally, we comment on the relation between these processes and shot-noise type processes, and observe that, under certain conditions, the steady-state distribution of one may be directly inferred from the other.

Chiu and Yin (2005) found the Laplace transform of the last time a spectrally negative Lévy process, which drifts to ∞, is below some level. The main motivation for the study of this random time stems from risk theory: what is the last time the risk process, modeled by a spectrally negative Lévy process drifting to ∞, is 0? In this paper we extend the result of Chiu and Yin, and we derive the Laplace transform of the last time, before an independent, exponentially distributed time, that a spectrally negative Lévy process (without any further conditions) exceeds (upwards or downwards) or hits a certain level. As an application, we extend a result found in Doney (1991).

In this paper we introduce the transformed two-parameter Poisson–Dirichlet distribution on the ordered infinite simplex. Furthermore, we prove the central limit theorem related to this distribution when both the mutation rate θ and the selection rate σ become large in a specified manner. As a consequence, we find that the properly scaled homozygosities have asymptotical normal behavior. In particular, there is a certain phase transition with the limit depending on the relative strength of σ and θ.

In Lyons, Pemantle and Peres (1995), a martingale change of measure method was developed in order to give an alternative proof of the Kesten–Stigum L log L theorem for single-type branching processes. Later, this method was extended to prove the L log L theorem for multiple- and general multiple-type branching processes in Biggins and Kyprianou (2004), Kurtz et al. (1997), and Lyons (1997). In this paper we extend this method to a class of superdiffusions and establish a Kesten–Stigum L log L type theorem for superdiffusions. One of our main tools is a spine decomposition of superdiffusions, which is a modification of the one in Englander and Kyprianou (2004).

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.

For a spectrally negative Lévy process X on the real line, let S denote its supremum process and let I denote its infimum process. For a > 0, let τ(a) and κ(a) denote the times when the reflected processes Ŷ := S − X and Y := X − I first exit level a, respectively; let τ−(a) and κ−(a) denote the times when X first reaches S τ(a) and I κ(a), respectively. The main results of this paper concern the distributions of (τ(a), S τ(a), τ−(a), Ŷ τ(a)) and of (κ(a), I κ(a), κ−(a)). They generalize some recent results on spectrally negative Lévy processes. Our approach relies on results concerning the solution to the two-sided exit problem for X. Such an approach is also adapted to study the excursions for the reflected processes. More explicit expressions are obtained when X is either a Brownian motion with drift or a completely asymmetric stable process.

We show in this paper that if a stationary traffic source is regulated by a leaky bucket with leak rate ρ and bucket size σ, then the amount of information generated in successive time intervals is dominated, in the increasing convex ordering sense, by that of a Poisson arrival process with rate ρ/σ, with each arrival bringing an amount of information equal to σ. By exploiting this property, we then show that the mean value in the stationary regime of the content of a buffer drained at constant rate and fed with the superposition of regulated flows is less than the mean value of the same buffer fed with an adequate Poisson process, whose characteristics depend upon the regulated input flows.

In this paper we are interested in a nonlinear parabolic evolution equation occurring in rheology. We give a probabilistic interpretation to this equation by associating a nonlinear martingale problem with it. We prove the existence of a unique solution, P, to this martingale problem. For any t, the time marginal of P at time t admits a density ρ(t,x) with respect to the Lebesgue measure, where the function ρ is the unique weak solution to the evolution equation in a well-chosen energy space. Next we introduce a simulable system of n interacting particles and prove that the empirical measure of this system converges to P as n tends to ∞. This propagation-of-chaos result ensures that the solution to the equation of interest can be approximated using a Monte Carlo method. Finally, we illustrate the convergence in some numerical experiments.

We discuss the limit behavior of the partial sums process of stationary solutions to the (autoregressive) AR(1) equation X t = a t X t−1 + ε t with random (renewal-reward) coefficient, a t , taking independent, identically distributed values A j ∈ [0,1] on consecutive intervals of a stationary renewal process with heavy-tailed interrenewal distribution, and independent, identically distributed innovations, ε t , belonging to the domain of attraction of an α-stable law (0 < α ≤ 2, α ≠ 1). Under suitable conditions on the tail parameter of the interrenewal distribution and the singularity parameter of the distribution of A j near the unit root a = 1, we show that the partial sums process of X t converges to a λ-stable Lévy process with index λ < α. The paper extends the result of Leipus and Surgailis (2003) from the case of finite-variance X t to that of infinite-variance X t .

Let {X(s, t), s = (s1, s2) ∈ ℝ2, t ∈ ℝ} be a stationary random field defined over a discrete lattice. In this paper, we consider a set of domain of attraction criteria giving the notion of extremal index for random fields. Together with the extremal-types theorem given by Leadbetter and Rootzen (1997), this will give a characterization of the limiting distribution of the maximum of such random fields.

In this paper, we present an iterative procedure to calculate explicitly the Laplace transform of the distribution of the maximum for a Lévy process with positive jumps of phase type. We derive error estimates showing that this iteration converges geometrically fast. Subsequently, we determine the Laplace transform of the law of the upcrossing ladder process and give an explicit pathwise construction of this process.

We provide a scaling for compound Poisson distributions that leads (under certain conditions on the Fourier transform) to a weak convergence result as the parameter of the distribution tends to infinity. We show that the limiting probability measure belongs to the class of stable Cauchy laws with Fourier transformt↦ exp(−c|t|− iatlog|t|). We apply this convergence result to the standard discrete Luria–Delbrück distribution and derive an integral representation for the corresponding limiting density, as an alternative to that found in a closely related paper of Kepler and Oprea. Moreover, we verify local convergence and we derive an integral representation for the distribution function of the limiting continuous Luria–Delbrück distribution.

Consider a geometric Brownian motion X t (ω) with drift. Suppose that there is an independent source that sends signals at random times τ 1 < τ 2 < ⋯. Upon receiving each signal, a decision has to be made as to whether to stop or to continue. Stopping at time τ will bring a reward S τ , where S t = max(max0≤u≤t X u , s) for some constant s ≥ X 0. The objective is to choose an optimal stopping time to maximize the discounted expected reward E[e−r τ i S τ i | X 0 = x, S 0 = s], where r is a discount factor. This problem can be viewed as a randomized version of the Bermudan look-back option pricing problem. In this paper, we derive explicit solutions to this optimal stopping problem, assuming that signal arrival is a Poisson process with parameter λ. Optimal stopping rules are differentiated by the frequency of the signal process. Specifically, there exists a threshold λ* such that if λ>λ*, the optimal stopping problem is solved via the standard formulation of a ‘free boundary’ problem and the optimal stopping time τ * is governed by a threshold a * such that τ * = inf{τ n : X τ n ≤a * S τ n }. If λ≤λ* then it is optimal to stop immediately a signal is received, i.e. at τ * = τ 1. Mathematically, it is intriguing that a smooth fit is critical in the former case while irrelevant in the latter.

In this paper, we show that the mean comparison theorem, which is valid for Brownian motion, cannot be extended to Poisson processes. A counterexample in the Poisson case for which the mean comparison theorem does not hold is provided.

A continuous-time random walk is a simple random walk subordinated to a renewal process used in physics to model anomalous diffusion. In this paper we show that, when the time between renewals has infinite mean, the scaling limit is an operator Lévy motion subordinated to the hitting time process of a classical stable subordinator. Density functions for the limit process solve a fractional Cauchy problem, the generalization of a fractional partial differential equation for Hamiltonian chaos. We also establish a functional limit theorem for random walks with jumps in the strict generalized domain of attraction of a full operator stable law, which is of some independent interest.

Classes of life distributions based on the moment-generating-function order are investigated in this paper. It is shown firstly that the class ℳ is closed under both convex linear combination and geometric compounding. Secondly, the class NBUmg (new better than used in the moment-generating-function order) is proved to be closed under increasing star-shaped transformations. Finally, the interplay between the stochastic comparison of the excess lifetime of a renewal process and the NBUmg interarrivals is studied.

This paper is motivated by the modelling of leaching of bacteria through soil. A semi-linear process X t − may be used to describe the soil-drying process between rain showers. This is a backward recurrence time process that corresponds to the renewal process of instances of rain. If a bacterium moves according to another process h, then the fact that h(t) stays above X t − means that the bacterium never hits a dry patch of soil and so survives. We describe a critical behaviour of h that separates the cases when survival is possible with a positive probability from the cases when this probability vanishes. An explicit formula for the survival probability is obtained in case h is linear and rain showers follow a Poisson process.

Continuous-time random walks incorporate a random waiting time between random jumps. They are used in physics to model particle motion. A physically realistic rescaling uses two different time scales for the mean waiting time and the deviation from the mean. This paper derives the scaling limits for such processes. These limit processes are governed by fractional partial differential equations that may be useful in physics. A transfer theorem for weak convergence of finite-dimensional distributions of stochastic processes is also obtained.

The maximum drawdown at time T of a random process on [0,T] can be defined informally as the largest drop from a peak to a trough. In this paper, we investigate the behaviour of this statistic for a Brownian motion with drift. In particular, we give an infinite series representation of its distribution and consider its expected value. When the drift is zero, we give an analytic expression for the expected value, and for nonzero drift, we give an infinite series representation. For all cases, we compute the limiting (T → ∞) behaviour, which can be logarithmic (for positive drift), square root (for zero drift) or linear (for negative drift).