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 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.

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.

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.

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 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.

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.

We consider a countable number of independent random uniform lines in the hyperbolic plane (in the sense of the theory of geometrical probability) which divide the plane into an infinite number of convex polygonal regions. The main purpose of the paper is to compute the mean number of sides, the mean perimeter, the mean area and the second order moments of these quantities of such polygonal regions. For the Euclidean plane the problem has been considered by several authors, mainly Miles [4]-[9] who has taken it as the starting point of a series of papers which are the basis of the so-called stochastic geometry. GEOMETRICAL PROBABILITY; RANDOM LINES; RANDOM POLYGONS; PLANE OF CONSTANT CURVATURE; HYPERBOLIC PLANE; CONVEX DOMAINS; POISSON LINE PROCESS

The paper represents a generalization of one of the main theoretical results of my Ph.D. thesis. The work is an outgrowth of work first begun by E. J. Hannan and a correct 'conjecture' by P. Whittle. The main theorem of this paper proves the existence of, and gives an explicit formula for, the asymptotic best linear predictor of a certain type of non-stationary multivariate time series from noise distorted observations. The non-stationarity arises from the fact that the signal satisfies a difference equation, which when considered as a polynomial, has only elementary divisors. The proof is accomplished by showing, through Hilbert space and harmonic analysis methods, that the generating function is a limit of the generating functions of the stationary analogue; that is, where the difference function has elementary divisors. Finally, it is shown that this asymptotic generating function exactly predicts null solutions to the difference equation. The proof is direct and due to E. J. Hannan.

This paper applies the invariance principle to the problem of testing a distribution on a compact homogeneous space for uniformity. The notion of using a reduction by invariance in such a situation is due to Ajne[1], who considers tests invariant under rotation on a circle. In his paper, he derives the distribution of the maximal invariant and gives the general form of the most powerful invariant test for uniformity on the circle.

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

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 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.

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.

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.

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).

This paper introduces long-range dependence for a stationary random field on a plane lattice, derives an exact power-law correlation model and other models with long-range dependence on the lattice, and explores the close connection between short-range dependent correlation functions and absolutely summable double sequences.

We analyse in this paper the fluid weighted fair queueing system with two classes of customers, who arrive according to Poisson processes and require arbitrarily distributed service times. In a first step, we express the Laplace transform of the joint distribution of the workloads in the two virtual queues of the system by means of unknown Laplace transforms. Such an unknown Laplace transform is related to the distribution of the workload in one queue provided that the other queue is empty. We explicitly compute the unknown Laplace transforms by means of a Wiener—Hopf technique. The determination of the unknown Laplace transforms can be used to compute some performance measures characterizing the system (e.g. the mean waiting time for each class) which we compute in the exponential service case.

In this paper, we apply coupling methods to study strong ergodicity for Markov processes, and sufficient conditions are presented in terms of the expectations of coupling times. In particular, explicit criteria are obtained for one-dimensional diffusions and birth-death processes to be strongly ergodic. As a by-product, strong ergodicity implies that the essential spectra of the generators for these processes are empty.

In this paper, we show that the time-average distributions of excess, age, and spread are given by the solution of first-order differential equations. These differential equations can be directly derived in a simple, unified way using a general level-crossing formula based on the balance of up and down crossings on sample paths, which may be helpful for the intuitive interpretation.

In this paper, the class of periodic Ornstein-Uhlenbeck processes is defined. It is shown that periodic Ornstein-Uhlenbeck processes are stationary Markov random fields and the class of stationary distributions is characterized. In particular, any self-decomposable distribution is the stationary distribution of some periodic Ornstein-Uhlenbeck process. As examples, gamma periodic Ornstein-Uhlenbeck processes and Gaussian periodic Ornstein-Uhlenbeck processes are considered.

In this paper we consider the stochastic wave equation in one spatial dimension driven by a two-parameter Gaussian noise which is white in time and has general spatial covariance. We give conditions on the spatial covariance of the driving noise sufficient for the string to have finite expected energy and calculate this energy as a function of time. We show that these same conditions on the spatial covariance of the driving noise are also sufficient to guarantee that the energy of the string has a version which is continuous almost surely.

In this paper, the instantaneous availability of a system maintained under periodic inspection is investigated using random walk models. Two cases are considered. In the first model, the system is repaired or modified and it is assumed to be as good as new upon periodic inspection and maintenance. In the second model, the system is not modified after the inspection if the system is still working, and the condition of the system is assumed to be the same as that before the inspection. For both models the failures only can be found through the inspection. Perfect repair or replacement of a failed system is assumed to be carried out, but the time it takes can be constant or of a random length. The relationship between this problem and the random walk model in a two-dimensional plane is described. Several new results are also shown.