The class of life distributions for which , where , and , is studied. We prove that this class is larger than the HNBUE (HNWUE) class (consisting of those life distributions for which for x ≧ 0) and present results concerning closure properties under some usual reliability operations. We also study some shock models and a certain cumulative damage model. The class of discrete life distributions for which for 0 ≦ p ≦ 1, where , is also studied.

For the distribution functions of the stationary actual waiting time and of the stationary virtual waiting time of the GI/G/l queueing system it is shown that the tails vary regularly at infinity if and only if the tail of the service time distribution varies regularly at infinity.

For sn the sum of n i.i.d. variables xi, i = 1, …, n it is shown that if E {x1} < 0 then the distribution of sup, s1s2, …] has a regularly varying tail at + ∞ if the tail of the distribution of x1 varies regularly at infinity and conversely, moreover varies regularly at + ∞.

In the appendix a lemma and its proof are given providing necessary and sufficient conditions for regular variation of the tail of a compound Poisson distribution.

For the queueing system GI/G/1 with both waiting-line and service-line customer impatience an integral equation for the limiting waiting-time distribution function W(x) is derived and the existence of W(x) proved when the system satisfies certain conditions. Two counter examples are constructed showing that W(x) does not necessarily exist. A sufficient condition for its existence in GI/G/1 with deterministic waiting-line impatience is given. Formal solutions for M/G/l are derived for deterministic impatience and negative exponential impatience.

The conditional tail expectation in risk analysis describes the expected amount of risk that can be experienced given that a potential risk exceeds a threshold value, and provides an important measure of right-tail risk. In this paper, we study the convolution and extreme values of dependent risks that follow a multivariate phase-type distribution, and derive explicit formulae for several conditional tail expectations of the convolution and extreme values for such dependent risks. Utilizing the underlying Markovian property of these distributions, our method not only provides structural insight, but also yields some new distributional properties of multivariate phase-type distributions.

An explicit formula for the probability that a Brownian motion crosses a piecewise linear boundary in a finite time interval is derived. This formula is used to obtain approximations to the crossing probabilities for general boundaries which are the uniform limits of piecewise linear functions. The rules for assessing the accuracies of the approximations are given. The calculations of the crossing probabilities are easily carried out through Monte Carlo methods. Some numerical examples are provided.

A construction is given for a stationary sequence of random variables {Xi} which have exponential marginal distributions and are random linear combinations of order one of an i.i.d. exponential sequence {εi}. The joint and trivariate exponential distributions of Xi−1, Xi and Xi+ 1 are studied, as well as the intensity function, point spectrum and variance time curve for the point process which has the {Xi} sequence for successive times between events. Initial conditions to make the point process count stationary are given, and extensions to higher-order moving averages and Gamma point processes are discussed.

For the sample functions of the stationary virtual waiting-time process vt of the GI/G/1 queueing system some properties of the number of up- and downcrossings of level v by the vt-process during a busy cycle are investigated. It turns out that the simple fact that this number of upcrossings is equal to that of downcrossings leads in a rather easy way to basic relations for the waiting-time distributions. This approach to the study of the vt-process of the GI/G/1 system seems to be applicable to many other types of stochastic processes. As another example of this approach the infinite dam with non-constant release rate is briefly discussed.

In this paper we introduce a quantile dispersion measure. We use it to characterize different classes of ageing distributions. Based on the quantile dispersion measure, we propose a new partial ordering for comparing the spread or dispersion in two probability distributions. This new partial ordering is weaker than the well known dispersive ordering and it retains most of its interesting properties.

Let X1,…,Xn be independent exponential random variables with Xi having hazard rate . Let Y1,…,Yn be a random sample of size n from an exponential distribution with common hazard rate ̃λ = (∏i=1nλi)1/n, the geometric mean of the λis. Let Xn:n = max{X1,…,Xn}. It is shown that Xn:n is greater than Yn:n according to dispersive as well as hazard rate orderings. These results lead to a lower bound for the variance of Xn:n and an upper bound on the hazard rate function of Xn:n in terms of . These bounds are sharper than those obtained by Dykstra et al. ((1997), J. Statist. Plann. Inference65, 203–211), which are in terms of the arithmetic mean of the λis. Furthermore, let X1*,…,Xn∗ be another set of independent exponential random variables with Xi∗ having hazard rate λi∗, i = 1,…,n. It is proved that if (logλ1,…,logλn) weakly majorizes (logλ1∗,…,logλn∗, then Xn:n is stochastically greater than Xn:n∗.

Necessary and sufficient conditions are given for a distribution function in ℝ2 to be max-infinitely divisible. The d.f. F is max i.d. if Ft is a d.f. for every t > 0. This property is essential in defining multivariate extremal processes and arises in an approach to the study of the range of an i.i.d. sample.

Two continuous distributions, G, H so related that any two quantiles of H are more widely separated than the corresponding quantiles of G may be said to be ‘ordered in dispersion'; Saunders and Moran have given examples. It is shown here that distributions F (called ‘dispersive' distributions) exist, e.g. the exponential, such that if G, H are ordered in dispersion then so also are the convolutions F ∗G, F ∗H. The class of dispersive distributions is determined, and shown to coincide with the class of strongly unimodal distributions.

The classical proofs for the existence of a stationary (s, S) inventory policy that minimizes the total discounted or average cost over an infinite horizon are lengthy because they depend heavily on the optimality results for corresponding finite-horizon models. This note presents a simpler alternative. Since optimal stationary (s, S) policies are relatively simple to characterize, it is easy to construct a solution to the optimality equation which is satisfied by an (s, S) policy or an equivalent variant thereof. For the discounted model, the proof characterizes an (s, S) policy that is optimal for all initial inventory positions. This policy can be generated by a simple existing algorithm. For the average-cost model, the optimality proof is completed with some additional arguments, which are simple but novel, to overcome the normal difficulties encountered in models with unbounded one-step expected costs.

We give conditions under which the stationary distribution π of a Markov chain admits moments of the general form ∫ f(x)π(dx), where f is a general function; specific examples include f(x) = xr and f(x) = esx. In general the time-dependent moments of the chain then converge to the stationary moments. We show that in special cases this convergence of moments occurs at a geometric rate. The results are applied to random walk on [0, ∞).

A correspondence between record values and independent increment point processes is established. The asymptotic behaviour of record value sequences is studied, and results on the asymptotic behaviour of record times (for continuous F) are obtained as special cases. The joint law of the kth record value and the kth record time is also derived.

A unified way of obtaining stationary time series models with the univariate margins in the convolution-closed infinitely divisible class is presented. Special cases include gamma, inverse Gaussian, Poisson, negative binomial, and generalized Poisson margins. ARMA time series models obtain in the special case of normal margins, sometimes in a different stochastic representation. For the gamma and Poisson margins, some previously defined time series models are included, but for the negative binomial margin, the time series models are different and, in several ways, better than previously defined time series models. The models are related to multivariate distributions that extend a univariate distribution in the convolution-closed infinitely divisible class. Extensions to the non-stationary case and possible applications to modelling longitudinal data are mentioned.

A finite-state Markov process is aggregated into several groups. What can be learned about the underlying process from the aggregated one? We provide some partial answers to this question.

We consider an Mx/G/1 queueing system with N-policy and multiple vacations. As soon as the system empties, the server leaves for a vacation of random length V. When he returns, if the queue length is greater than or equal to a predetermined value N(threshold), the server immediately begins to serve the customers. If he finds less than N customers, he leaves for another vacation and so on until he finally finds at least N customers. We obtain the system size distribution and show that the system size decomposes into three random variables one of which is the system size of ordinary Mx/G/1 queue. The interpretation of the other random variables will be provided. We also derive the queue waiting time distribution and other performance measures. Finally we derive a condition under which the optimal stationary operating policy is achieved under a linear cost structure.

Let {X k , k ≥ 1} be a sequence of independent, identically distributed nonnegative random variables with common distribution function F and finite expectation μ > 0. Under the assumption that the tail probability is consistently varying as x tends to infinity, this paper investigates precise large deviations for both the partial sums S n and the random sums S N(t), where N(·) is a counting process independent of the sequence {X k , k ≥ 1}. The obtained results improve some related classical ones. Applications to a risk model with negatively associated claim occurrences and to a risk model with a doubly stochastic arrival process (extended Cox process) are proposed.

Stein's method of obtaining rates of convergence, well known in normal and Poisson approximation, is considered here in the context of approximation by Poisson point processes, rather than their one-dimensional distributions. A general technique is sketched, whereby the basic ingredients necessary for the application of Stein's method may be derived, and this is applied to a simple problem in Poisson point process approximation.

This paper extends results of Esary, Marshall and Proschan (1973) and A-Hameed and Proschan (1973). We consider the life distribution of a device subject to a sequence of shocks occurring randomly in time according to a nonstationary pure birth process: given k shocks have occurred in [0, t], the probability of a shock occurring in (t, t + Δ] is λ kλ (t)Δ + o (Δ). We show that various fundamental classes of life distributions (such as those with increasing failure rate, or those with the ‘new better than used' property, etc.) are obtained under appropriate assumptions on λ k, λ (t), and on the probability of surviving a given number of shocks.