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.

The inverse stochastic dominance of degree r is a stochastic order of interest in several branches of economics. We discuss it in depth, the central point being the characterization in terms of the weak r-majorization of the vectors of expected order statistics. The weak r-majorization (a notion introduced in the paper) is a natural extension of the classical (reverse) weak majorization of Hardy, Littlewood and Pòlya. This work also shows the equivalence between the continuous majorization (of higher order) and the discrete r-majorization. In particular, our results make it clear that the cases r = 1, 2 differ substantially from those with r ≥ 3, a fact observed earlier by Muliere and Scarsini (1989), among other authors. Motivated by this fact, we introduce new stochastic orderings, as well as new social inequality indices to compare the distribution of the wealth in two populations, which could be considered as natural extensions of the first two dominance rules and the S-Gini indices, respectively.

We approximate the distribution of the sum of independent but not necessarily identically distributed Bernoulli random variables using a shifted binomial distribution, where the three parameters (the number of trials, the probability of success, and the shift amount) are chosen to match the first three moments of the two distributions. We give a bound on the approximation error in terms of the total variation metric using Stein's method. A numerical study is discussed that shows shifted binomial approximations are typically more accurate than Poisson or standard binomial approximations. The application of the approximation to solving a problem arising in Bayesian hierarchical modeling is also discussed.

A necessary condition for a rational Laplace–Stieltjes transform to correspond to a matrix exponential distribution is that the pole of maximal real part is real and negative. Given a rational Laplace–Stieltjes transform with such a pole, we present a method to determine whether or not the numerator polynomial admits a transform that corresponds to a matrix exponential distribution. The method relies on the minimization of a continuous function of one variable over the nonnegative real numbers. Using this approach, we give an alternative characterization for all matrix exponential distributions of order three.

A model for describing the lifetimes of coherent systems, where the failures of components may have an impact on the lifetimes of the remaining components, is proposed. The model is motivated by the definition of sequential order statistics (cf. Kamps (1995)). Sequential order statistics describe the successive failure times in a sequential k-out-of-n system, where the distribution of the remaining components' lifetimes is allowed to change after every failure of a component. In the present paper, general component lifetimes which can be influenced by failures are considered. The ordered failure times of these components can be used to extend the concept of sequential order statistics. In particular, a definition of sequential order statistics based on exchangeable components is proposed. By utilizing the system signature (cf. Samaniego (2007)), the distribution of the lifetime of a coherent system with failure-dependent exchangeable component lifetimes is shown to be given by a mixture of the distributions of sequential order statistics. Furthermore, some results on the joint distribution of sequential order statistics based on exchangeable components are given.

We consider coherent and mixed reliability systems composed of elements with independent and identically distributed lifetimes. We present upper bounds on variances of system lifetimes, expressed in terms of variances of single components. We also discuss attainability conditions and some special cases and examples.

Directionally convex ordering is a useful tool for comparing the dependence structure of random vectors, which also takes into account the variability of the marginal distributions. It can be extended to random fields by comparing all finite-dimensional distributions. Viewing locally finite measures as nonnegative fields of measure values indexed by the bounded Borel subsets of the space, in this paper we formulate and study directionally convex ordering of random measures on locally compact spaces. We show that the directionally convex order is preserved under some of the natural operations considered on random measures and point processes, such as deterministic displacement of points, independent superposition, and thinning, as well as independent, identically distributed marking. Further operations on Cox point processes such as position-dependent marking and displacement of points are shown to preserve the order. We also examine the impact of the directionally convex order on the second moment properties, in particular on clustering and on Palm distributions. Comparisons of Ripley's functions and pair correlation functions, as well as examples, seem to indicate that point processes higher in the directionally convex order cluster more. In our main result we show that nonnegative integral shot noise fields with respect to the directionally convex ordered random measures inherit this ordering from the measures. Numerous applications of this result are shown, in particular to comparison of various Cox processes and some performance measures of wireless networks, in both of which shot noise fields appear as key ingredients. We also mention a few pertinent open questions.