Melamed's theorem states that, for a Jackson queueing network, the equilibrium flow along a link follows a Poisson distribution if and only if no customers can travel along the link more than once. Barbour and Brown (1996) considered the Poisson approximate version of Melamed's theorem by allowing the customers a small probability p of travelling along the link more than once. In this note, we prove that the customer flow process is a Poisson cluster process and then establish a general approximate version of Melamed's theorem that accommodates all possible cases of 0 ≤ p < 1.

In this paper we study the Wiener-Hopf factorization for a class of Lévy processes with double-sided jumps, characterized by the fact that the density of the Lévy 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évy processes, defined by the fact that the density of the Lévy 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.

The distributions that result from zero-truncating mixed Poisson (ZTMP) distributions and those obtained from mixing zero-truncated Poisson (MZTP) distributions are characterised based on their probability generating functions. One consequence is that every ZTMP distribution is an MZTP distribution, but not vice versa. These characterisations also indicate that the size-biased version of a Poisson mixture and, under certain regularity conditions, the shifted version of a Poisson mixture are neither ZTMP distributions nor MZTP distributions.

Let {X i }i=1 n be a sequence of random variables with two possible outcomes, denoted 0 and 1. Define a random variable S n,m to be the maximum number of 1s within any m consecutive trials in {X i }i=1 n . The random variable S n,m is called a discrete scan statistic and has applications in many areas. In this paper we evaluate the distribution of discrete scan statistics when {X i }i=1 n consists of exchangeable binary trials. We provide simple closed-form expressions for both conditional and unconditional distributions of S n,m for 2m ≥ n. These results are also new for independent, identically distributed Bernoulli trials, which are a special case of exchangeable trials.

We show that the conjecture of Kannan, Lovász, and Simonovits on isoperimetric properties of convex bodies and log-concave measures is true for log-concave measures of the form ρ(∣x∣B) dx on ℝn and ρ(t,∣x∣B) dx on ℝ1+n, where ∣x∣B is the norm associated to any convex body B already satisfying the conjecture. In particular, the conjecture holds for convex bodies of revolution.

In this paper we further investigate the problem considered by Mizuno (2006) in the special case of identically distributed signals. Specifically, we first propose an alternative sufficient condition of crossing type for the convex order to hold between the conditional expectations given signal. Then, we prove that the bivariate (2,1)-increasing convex order ensures that the conditional expectations are ordered in the convex sense. Finally, the L 2 distance between the quantity of interest and its conditional expectation given signal (or expected conditional variance) is shown to decrease when the strength of the dependence increases (as measured by the (2,1)-increasing convex order).

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.

We introduce and motivate the study of (n + 1) × r arrays X with Bernoulli entries X k,j and independently distributed rows. We study the distribution of which denotes the number of consecutive pairs of successes (or runs of length 2) when reading the array down the columns and across the rows. With the case r = 1 having been studied by several authors, and permitting some initial inferences for the general case r > 1, we examine various distributional properties and representations of S n for the case r = 2, and, using a more explicit analysis, the case of multinomial and identically distributed rows. Applications are also given in cases where the array X arises from a Pólya sampling scheme.

For several pairs (P, Q) of classical distributions on ℕ0, we show that their stochastic ordering P ≤stQ can be characterized by their extreme tail ordering equivalent to P({k *})/Q({k *}) ≥ 1 ≥ limk→k* P({k})/Q({k}), with k * and k * denoting the minimum and the supremum of the support of P + Q, and with the limit to be read as P({k *})/Q({k *}) for finite k *. This includes in particular all pairs where P and Q are both binomial (b n1,p1 ≤stb n2,p2 if and only if n 1 ≤ n 2 and (1 - p 1)n1 ≥ (1 - p 2)n2 , or p 1 = 0), both negative binomial (b − r1,p1 ≤stb − r2,p2 if and only if p 1 ≥ p 2 and p 1 r1 ≥ p 2 r2 ), or both hypergeometric with the same sample size parameter. The binomial case is contained in a known result about Bernoulli convolutions, the other two cases appear to be new. The emphasis of this paper is on providing a variety of different methods of proofs: (i) half monotone likelihood ratios, (ii) explicit coupling, (iii) Markov chain comparison, (iv) analytic calculation, and (v) comparison of Lévy measures. We give four proofs in the binomial case (methods (i)-(iv)) and three in the negative binomial case (methods (i), (iv), and (v)). The statement for hypergeometric distributions is proved via method (i).

In this paper we propose a general bivariate random effect model with special emphasis on frailty models and environmental effect models, and present some stochastic comparisons. The relationship between the conditional and the unconditional hazard gradients are derived and some examples are provided. We investigate how the well-known stochastic orderings between the distributions of two frailties translate into the orderings between the corresponding survival functions. These results are used to obtain the properties of the bivariate multiplicative model and the shared frailty model.

Expectations of unbounded functions of dependent nonnegative integer-valued random variables are approximated by the expectations of the functions of independent copies of these random variables. The Lindeberg method is used.

Li and Shaked (2007) introduced the family of generalized total time on test transform (TTT) stochastic orders, which is parameterized by a real function h that can be used to capture the preferences of a decision maker. It is natural to look for properties of these orders when there is an uncertainty in determining the appropriate function h. In this paper we study these orders when h is nondecreasing. We note that all these orders are location independent, and we characterize the dispersive order, and the location-independent riskier order, by means of the generalized TTT orders with nondecreasing h. Further properties, which strengthen known properties of the dispersive order, are given. A useful nontrivial closure property of the generalized TTT orders with nondecreasing h is obtained. Applications in poverty comparisons, risk management, and reliability theory are described.

Let X1 and X2 be two independent and nonnegative random variables with distributions F1 and F2, respectively. This paper proves that if both F1 and F2 are of Weibull type and fulfill certain easily verifiable conditions, then the distribution of the product X1X2, called the product convolution of F1 and F2, belongs to the class 𝒮* and, hence, is subexponential.

We consider a system with Poisson arrivals and independent and identically distributed service times, where requests in the system are served according to the state-dependent (Cohen's generalized) processor-sharing discipline, where each request receives a service capacity that depends on the actual number of requests in the system. For this system, we derive expressions as well as tight insensitive upper bounds for the moments of the conditional sojourn time of a request with given required service time. The bounds generalize and extend corresponding results, recently given for the single-server processor-sharing system in Cheung et al. (2006) and for the state-dependent processor-sharing system with exponential service times by the authors (2008). Analogous results hold for the waiting times. Numerical examples for the M/M/m-PS and M/D/m-PS systems illustrate the given bounds.

In this paper, the componentwise increasing convex order, the upper orthant order, the upper orthant convex order, and the increasing directionally convex order for random vectors are generalized to hierarchical classes of integral stochastic order relations. The elements of the generating classes of functions possess nonnegative partial derivatives up to some given degrees. Some properties of these new stochastic order relations are studied. Particular attention is paid to the comparison of weighted sums of the respective components of ordered random vectors. By providing a unified derivation of standard multivariate stochastic orderings, the present paper shows how some well-known results derive from a common principle.

Consider a sequence of exchangeable or independent binary trials ordered on a line or on a circle. The statistics denoting the number of times an F-S string of length (at least)k 1 + k 2, that is, (at least)k 1 failures followed by (at least) k 2 successes in n such trials, are studied. The associated waiting time for the rth occurrence of an F-S string of length (at least) k 1 + k 2 in linearly ordered trials is also examined. Exact formulae, lower/upper bounds and approximations are derived for their distributions. Mean values and variances of the number of occurrences of F-S strings are given in exact formulae too. Particular exchangeable and independent sequences of binary random variables, used in applied research, combined with numerical examples clarify further the theoretical results.

System signatures are useful tools in the study and comparison of coherent systems. In this paper, we define and study a similar concept, called the joint signature, for two coherent systems which share some components. Under an independent and identically distributed assumption on component lifetimes, a pseudo-mixture representation based on this joint signature is obtained for the joint distribution of the lifetimes of both systems. Sufficient conditions are given based on the respective joint signatures of two pairs of systems, each with shared components, to ensure various forms of bivariate stochastic orderings between the joint lifetimes of the two pairs of systems.

Measures of divergence or discrepancy are used either to measure mutual information concerning two variables or to construct model selection criteria. In this paper we focus on divergence measures that are based on a class of measures known as Csiszár's divergence measures. In particular, we propose a measure of divergence between residual lives of two items that have both survived up to some time t as well as a measure of divergence between past lives, both based on Csiszár's class of measures. Furthermore, we derive properties of these measures and provide examples based on the Cox model and frailty or transformation model.

In this note we deal with the allocation of independent and identical active redundancies to a k-out-of-n system with the usual stochastic order among its independent components. The optimal policy is proved both to assign more redundancies to the weaker component and to majorize all other policies. This improves the corresponding one in Hu and Wang (2009) and serves as a nice supplement to that in Misra, Dhariyal and Gupta (2009) as well.

We apply Stein's method for probabilistic approximation by a compound geometric distribution, with applications to Markov chain hitting times and sequence patterns. Bounds on our Stein operator are found using a complex analytical approach based on generating functions and Cauchy's formula.