We consider the problem of allocating k active spares to n components of a series system in order to optimize its lifetime. Under the hypotheses that lifetimes of n components are identically distributed with distribution function F(⋅), lifetimes of k spares are identically distributed with distribution function G(⋅), lifetimes of components and spares are independently distributed, and that ln(G(x))/ln(F(x)) is increasing in x, we show that the strategy of balanced allocation of spares optimizes the failure rate function of the system. Furthermore, under the hypotheses that lifetimes of n components are stochastically ordered, lifetimes of k spares are identically distributed, and that lifetimes of components and spares are independently distributed, we show that the strategy of balanced allocation of spares is superior to the strategy of allocating a larger number of components to stronger components. For coherent systems consisting of n identical components with n identical redundant (spare) components, we compare strategies of component and system redundancies under the criteria of reversed failure rate and likelihood ratio orderings. When spares and original components do not necessarily match in their life distributions, we provide a sufficient condition, on the structure of the coherent system, for the strategy of component redundancy to be superior to the strategy of system redundancy under reversed failure rate ordering. As a consequence, we show that, for r-out-of-n systems, the strategy of component redundancy is superior to the strategy of system redundancy under the criterion of reversed failure rate ordering. When spares and original components match in their life distributions, we provide a necessary and sufficient condition, on the structure of the coherent system, for the strategy of component redundancy to be superior to the strategy of system redundancy under the likelihood ratio ordering. As a consequence, we show that, for r-out-of-n systems, with spares and original components matching in their life distributions, the strategy of component redundancy is superior to the strategy of system redundancy under the likelihood ratio ordering.

We investigate stochastic comparisons between exponential family distributions and their mixtures with respect to the usual stochastic order, the hazard rate order, the reversed hazard rate order, and the likelihood ratio order. A general theorem based on the notion of relative log-concavity is shown to unify various specific results for the Poisson, binomial, negative binomial, and gamma distributions in recent literature. By expressing a convolution of gamma distributions with arbitrary scale and shape parameters as a scale mixture of gamma distributions, we obtain comparison theorems concerning such convolutions that generalize some known results. Analogous results on convolutions of negative binomial distributions are also discussed.

In this paper, a new approach is proposed to investigate Blackwell-type renewal theorems for weighted renewal functions systematically according to which of the tails of weighted renewal constants or the underlying distribution is asymptotically heavier. Some classical results are improved considerably.

The upper tail behaviour is explored for a stopped random product ∏j=1 N X j , where the factors are positive and independent and identically distributed, and N is the first time one of the factors occupies a subset of the positive reals. This structure is motivated by a heavy-tailed analogue of the factorial n!, called the factoid of n. Properties of the factoid suggested by computer explorations are shown to be valid. Two topics about the determination of the Zipf exponent in the rank-size law for city sizes are discussed.

The representation of the reliability function of the lifetime of a coherent system as a mixture of the reliability function of order statistics associated with the lifetimes of its components is a very useful tool to study the ordering and the limiting behaviour of coherent systems. In this paper, we obtain several representations of the reliability functions of residual lifetimes of used coherent systems under two particular conditions on the status of the components or the system in terms of the reliability functions of residual lifetimes of order statistics.

We present a method to obtain both exact values and sharp estimates for the total variation distance between binomial and Poisson distributions with the same mean λ. We give a simple efficient algorithm, whose complexity order is to compute exact values. Such an algorithm can be further simplified for moderate sample sizes n, provided that λ is neither close to from the left nor close to from the right. Sharp estimates, better than other known estimates in the literature, are also provided. The 0s of the second Krawtchouk and Charlier polynomials play a fundamental role.

In this paper we study the properties of total time on test transforms of order n and examine their applications in reliability analysis. It is shown that the successive transforms produce either distributions with increasing or bathtub-shaped failure rates or distributions with decreasing or upside bathtub-shaped failure rates. The ageing properties of the baseline distribution is compared with those of transformed distributions, and a partial order based on nth-order transforms and their implications are discussed.

By considering k-out-of-n systems with independent and nonidentically distributed components, we discuss stochastic monotone properties of the residual life and the inactivity time. We then present some stochastic comparisons of two systems based on the residual life and inactivity time.

We study a particular example of a recursive distributional equation (RDE) on the unit interval. We identify all invariant distributions, the corresponding ‘basins of attraction’, and address the issue of endogeny for the associated tree-indexed problem, making use of an extension of a recent result of Warren.

We describe the genealogy of individuals with infinite descent in a supercritical continuous-state branching process.

The distributions of the run occurrences for a sequence of independent and identically distributed (i.i.d.) experiments are usually obtained by combinatorial methods (see Balakrishnan and Koutras (2002, Chapter 5)) and the resulting formulae are often very tedious, while the distributions for non i.i.d. experiments are generally intractable. It is therefore of practical interest to find a suitable approximate model with reasonable approximation accuracy. In this paper we demonstrate that the negative binomial distribution is the most suitable approximate model for the number of k-runs: it outperforms the Poisson approximation, the general compound Poisson approximation as observed in Eichelsbacher and Roos (1999), and the translated Poisson approximation in Rollin (2005). In particular, its accuracy of approximation in terms of the total variation distance improves when the number of experiments increases, in the same way as the normal approximation improves in the Berry-Esseen theorem.

For an order statistic (X 1:n ,…,X n:n ) of a collection of independent but not necessarily identically distributed random variables, and any i ∈ {1,…,n}, the conditional distribution of (X i+1:n ,…,X n:n ) given X i:n > s is shown to be stochastically increasing in s. This answers a question by Hu and Xie (2006).

A comparison theorem is stated for Markov processes in Polish state spaces. We consider a general class of stochastic orderings induced by a cone of real functions. The main result states that stochastic monotonicity of one process and comparability of the infinitesimal generators imply ordering of the processes. Several applications to convex type and to dependence orderings are given. In particular, Liggett's theorem on the association of Markov processes is a consequence of this comparison result.

The tail of risk neutral returns can be related explicitly with the wing behaviour of the Black-Scholes implied volatility smile. In situations where precise tail asymptotics are unknown but a moment generating function is available we establish, under easy-to-check Tauberian conditions, tail asymptotics on logarithmic scales. Such asymptotics are enough to make the tail-wing formula (see Benaim and Friz (2008)) work and so we obtain, under generic conditions, a limiting slope when plotting the square of the implied volatility against the log strike, improving a lim sup statement obtained earlier by Lee (2004). We apply these results to time-changed exponential Lévy models and examine several popular models in more detail, both analytically and numerically.

Let S 0 := 0 and Sk := ξ 1 + ··· + ξk for k ∈ ℕ := {1, 2, …}, where {ξk : k ∈ ℕ} are independent copies of a random variable ξ with values in ℕ and distribution pk := P{ξ = k}, k ∈ ℕ. We interpret the random walk {Sk : k = 0, 1, 2, …} as a particle jumping to the right through integer positions. Fix n ∈ ℕ and modify the process by requiring that the particle is bumped back to its current state each time a jump would bring the particle to a state larger than or equal to n. This constraint defines an increasing Markov chain {Rk (n) : k = 0, 1, 2, …} which never reaches the state n. We call this process a random walk with barrier n. Let Mn denote the number of jumps of the random walk with barrier n. This paper focuses on the asymptotics of Mn as n tends to ∞. A key observation is that, under p 1 > 0, {Mn : n ∈ ℕ} satisfies the distributional recursion M 1 = 0 and for n = 2, 3, …, where In is independent of M 2, …, M n−1 with distribution P{In = k} = pk / (p 1 + ··· + p n−1), k ∈ {1, …, n − 1}. Depending on the tail behavior of the distribution of ξ, several scalings for Mn and corresponding limiting distributions come into play, including stable distributions and distributions of exponential integrals of subordinators. The methods used in this paper are mainly probabilistic. The key tool is to compare (couple) the number of jumps, Mn , with the first time, Nn , when the unrestricted random walk {Sk : k = 0, 1, …} reaches a state larger than or equal to n. The results are applied to derive the asymptotics of the number of collision events (that take place until there is just a single block) for β(a, b)-coalescent processes with parameters 0 < a < 2 and b = 1.

In this paper we study a single-period optimal portfolio problem in which the aim of the investor is to maximize the expected utility. We assume that the return of every security in the market is a mixture of some common underlying source of risks. A sufficient condition to order the optimal allocations is obtained, and it is shown that several models studied in the literature before are special cases of the proposed model. In the course of the analysis concepts in stochastic orders are employed, and a new characterization of the likelihood ratio order is obtained.

The copula of a multivariate distribution is the distribution transformed so that one-dimensional marginal distributions are uniform. We review a different transformation of a multivariate distribution which yields standard Pareto for the marginal distributions, and we call the resulting distribution the Pareto copula. Use of the Pareto copula has a certain claim to naturalness when considering asymptotic limit distributions for sums, maxima, and empirical processes. We discuss implications for aggregation of risk and offer some examples.

We make a correction to an important result by Cline [D. B. H. Cline, ‘Convolutions of distributions with exponential tails’, J. Austral. Math. Soc. (Series A)43 (1987), 347–365; D. B. H. Cline, ‘Convolutions of distributions with exponential tails: corrigendum’, J. Austral. Math. Soc. (Series A)48 (1990), 152–153] on the closure of the exponential class under convolution power mixtures (random summation).

Recently Li and Yam (2005) studied which ageing properties for series and parallel systems are inherited for the components. In this paper we provide new results for the increasing convex and concave orders, the increasing mean residual life (IMRL), decreasing failure rate (DFR), the new worse than used in expectation (NWUE), the increasing failure rate in average (IFRA), the decreasing failure rate in average (DFRA), and the new better than used in the convex order (NBUC) ageing classes.

We prove that the Bartlett spectrum of a stationary, infinitely divisible (ID) random measure determines ergodicity, weak mixing, and mixing. In this context, the Bartlett spectrum plays the same role as the spectral measure of a stationary Gaussian process.