#### Abstract

We review some of the recent developments in the area of stochastic comparisons of order statistics and sample spacings. We consider the cases when the parent observations are identically as well as nonidentically distributed. But most of the time we will be assuming that the observations are independent. The case of independent exponentials with unequal scale parameters as well as the proportional hazard rate model is discussed in detail.

#### 1. Introduction

The simplest and the most popular method of comparing the magnitudes of two random variables is through their means and medians. It may happen that in some cases the median of is larger than that of , while the mean of is smaller than the mean of . However, this confusion will not arise if the random variables are stochastically ordered. Similarly, the same may happen if one would like to compare the variability of with that of based only on numerical measures like standard deviation, and so forth. Besides, these characteristics of distributions might not exist in some cases. In most cases one can express various forms of knowledge about the underlying distributions in terms of their survival functions, quantile functions, hazard rate functions, mean residual functions, and other suitable functions of probability distributions. These methods are much more informative than those based only on few numerical characteristics of distributions. Comparisons of random variables based on such functions usually establish partial orders among them. We call them as stochastic orders.

Stochastic models are usually sufficiently complex in various fields of statistics. Obtaining bounds and approximations for their characteristics is of practical importance. That is, the approximation of a stochastic model either by a simpler model or by a model with simple constituent components might lead to convenient bounds and approximations for some particular and desired characteristics of the model. The study of changes in the properties of a model, as the constituent components vary, is also of great interest. Accordingly, since the stochastic components of models involve random variables, the topic of stochastic orders among random variables plays an important role in these areas. Books by Muller and Stoyan [1] and Shaked and Shanthikumar [2] give excellent treatment of this topic.

Order statistics and spacings are of great interest in many areas of statistics and they have received a lot of attention from many researchers. Let be random variables. The th order statistic, the th smallest of ’s, is denoted by . In reliability engineering, an component system that works if and only if at least of the components work is called a -out-of- system. The lifetime of a -out-of- system can be represented as . A parallel system is a -out-of- system while a series system is an -out-of- system. Thus, the study of lifetimes of -out-of- systems is equivalent to the study of the stochastic properties of order statistics. Spacings, the differences between successive order statistics, and their functions are also important in statistics, in general, and in particular in the context of life testing and reliability models. The books by David and Nagaraja [3], and Arnold et al. [4]; and two volumes of papers on this topic by Balakrishnan and Rao [5, 6] are excellent sources of information on this topic. But most of this work has been confined to the case when the observations are independent and identically distributed (i.i.d.). In many practical situations, like in reliability theory, the observations are not necessarily i.i.d. Only during the last two decades or so this topic has got the attention of researchers. Some important early references for this case are Sen [7], David and Nagaraja [3], Shaked and Tong [8], Bapat and Beg [9], Boland et al. [10], Kochar [11], and Nappo and Spizzichino [12], Boland et al. [13] and Balakrishnan [14], among others.

Some interesting partial ordering results on order statistics and spacings from independent but nonidentically random variables have been obtained by Pledger and Proschan [15], Proschan and Sethuraman [16], Bapat and Kochar [17], Boland et al. [18], Kochar and Kirmani [19], Kochar and Korwar [20], Kochar and Rojo [21], Dykstra et al. [22], Kochar and Ma [23], Bon and Pǎltǎnea [24], Kochar [25], Khaledi and Kochar [26–30]. The book by Shaked and Shanthikumar [2] gives an excellent description of the various results on this topic till 2007.

In this review paper, we will discuss some latest developments on the topic of stochastic comparisons of order statistics and spacings. In Section 2, we introduce the required notation and definitions. Sections 3 and 4 are devoted to stochastic comparisons of order statistics in one-sample and two-sample problems, respectively. In Section 5, we discuss the topic of stochastic orderings among spacings in one-sample problem and two sample problems. Section 6 is devoted to some applications of these results. Throughout this chapter *increasing* means *nondecreasing* and *decreasing* means *nonincreasing*; and we will be assuming that all distributions under study are absolutely continuous.

#### 2. Definitions and Some Preliminaries

In this section, we recall some basic definitions of stochastic orders and their properties. Let and be univariate random variables with distribution functions and , survival functions and , density functions and ; and hazard rates () and (), respectively. Let and be the left and the right endpoints of the support of .

##### 2.1. Magnitude Orders

First we give some definitions of stochastic orders to compare the magnitudes of two random variables.

*Definition 2.1. *is said to be *stochastically* smaller than (denoted by ) if

It is easy to see that if and only if It can be shown that this is also equivalent to for all increasing function for which expectations exist.

*Definition 2.2. * is said to be smaller than in hazard rate ordering (denoted by ) if is increasing in .

It can be shown that is equivalent to the inequalities In other words, the conditional distributions, given that the random variables are at least of a certain size, are all stochastically ordered (in the standard sense) in the same direction. Thus, if and represent the survival times of different models of an appliance that satisfy this ordering, one model is better (in the sense of stochastic ordering) when the appliances are new, the same appliance is better when both are one month old, and in fact is better no matter how much time has elapsed. It is clearly useful to know when this strong type of stochastic ordering holds since qualitative judgements are then easy to make. In case the hazard rates exist, it is easy to see that , if and only if, for every . The hazard rate ordering is also known as uniform stochastic ordering in the literature.

*Definition 2.3. * is said to be smaller than in the reverse hazard rate order, denoted by , if

It is shown in Ross [31] that if and are two independent random variables, then may not imply
However (2.5) will hold if and satisfy a stronger ordering called *likelihood ratio ordering* as defined below.

*Definition 2.4. * is said to be smaller than in likelihood ratio ordering (denoted by ) if is increasing in .

The mean residual life function of a random variable is defined as . An interesting order based on the mean residual life is defined as follows.

*Definition 2.5. * is said to be smaller than in the mean residual life order, denoted by , if

When the supports of and have a common left end-point, we have the following chain of implications among the above stochastic orders:

*Definition 2.6. * is said to be smaller than in the *increasing convex order* (denoted by ) if

Note that (2.8) holds if and only if for every convex function for which the above expectations exist. Also note that .

Another order closely related to the likelihood ratio order is the joint likelihood ordering as introduced by Shanthikumar and Yao [32].

*Definition 2.7. *For a bivariate random variable , is said to be smaller than according to joint likelihood ordering, denoted by , if and only if
where

It can be seen that where denotes the joint density of .

As pointed out by Shanthikumar and Yao [32], joint likelihood ratio ordering between the components of a bivariate random vector may not imply likelihood ratio ordering between their marginal distributions unless the random variables are independent, but it does imply stochastic ordering between them, that is,

A bivariate function is called *arrangement increasing* (AI). Hollander et al. [33] have studied many interesting properties of such functions, though, apparently, they did not relate it to the notion of likelihood ratio ordering.

The above idea can be extended to compare the components of an -dimensional vector . We define if the joint density of is an *arrangement increasing function*. (See, Marshall et al. [34] for the definition of an arrangement increasing function on .)

In a different context, Robertson and Wright [35] studied a subclass of arrangement increasing functions on , which they call as functions, as described below. Let and be two vectors on such that and . We will denote this partial ordering between the vectors by .

*Definition 2.8. * A real-valued function defined on a set is said to be on if .

As mentioned earlier, an function is arrangement increasing but the converse is not true. It is easy to see that the joint density of a bivariate random vector is if and only if the conditional density of given is monotonically increasing for each fixed .

The usual likelihood ratio order has the following multivariate version.

*Definition 2.9. * Let and be two -dimensional random vectors with absolutely continuous (or discrete) distribution functions and let and denote their continuous (or discrete) density functions, respectively. Suppose that
for every and . Then is said to be smaller than in the multivariate likelihood ratio order (denoted as ).

*Definition 2.10. *The random vector is smaller than the random vector in the multivariate stochastic order (denoted by ) if for all increasing functions .

It is known that multivariate likelihood ratio ordering implies multivariate stochastic ordering, but the converse is not true. Also if two random vectors are ordered according to multivariate likelihood ratio (stochastic) ordering, then their corresponding subsets of components are also ordered accordingly. For more details on stochastic orderings, see Chapters 1 and 4 of Shaked and Shanthikumar [2].

##### 2.2. Variability Orders

One of the basic criteria for comparing variability in probability distributions is that of dispersive ordering. Let and be the right continuous inverses (quantile functions) of and , respectively.

*Definition 2.11. * is less dispersed than (denoted by ) if

Note that (2.15) is equivalent to

When (2.16) holds, Doksum [36] called this ordering as tail-ordering and used it to find bounds on powers and efficiencies of nonparametric tests. Deshpande and Kochar [37] pointed out that tail ordering is same as dispersive ordering and obtained some new results for this partial order.

A consequence of is that and which in turn implies as well as , where are two independent copies of . For details, see Saunders [38], Saunders and Moran [39], Lewis and Thompson [40], Deshpande and Kochar [37], Bagai and Kochar [41], Bartoszewicz [42, 43], and Section 2.B of Shaked and Shanthikumar [2].

*Connection between Hazard Rate Order and Dispersive Order*

Kochar [44] observed that
where, for example, means with the notation that denotes the th order statistic of a random sample on .

Equation (2.17) indicates that hazard rate ordering not only compares the magnitudes of two random variables, but it also has perhaps some connection with the variability between the random variables. On differentiating (2.16), one can easily see that when the random variables and admit densities. This lead Bagai and Kochar [41] to prove the following connection between hazard rate ordering and dispersive ordering.

Theorem 2.12. *(a) If or have increasing failure rate (IFR), then . **(b) If or have decreasing failure rate (DFR), then . *

For a review on applications of dispersive ordering, see Joen et al. [45].

A weaker variability ordering is *right spread order* as introduced in Fernandez-Ponce et al. [46]. It was also independently studied by Shaked and Shanthikumar[47] who call it as *excess wealth order*.

*Definition 2.13. * is said to be smaller than in the *excess wealth order* (denoted by ) if

Dispersive ordering implies excess wealth ordering which in turn implies that the commonly used measures of variability like variances are ordered. One may refer to Shaked and Shanthikumar [2] for a comprehensive discussion of this order.

##### 2.3. Skewness Orderings

Skewness describes the departure of a distribution from symmetry, where one tail of the density is more “stretched out" than the other. Several partial orders have been introduced in the literature to compare the relative skewness of probability distributions. van Zwet [48] introduced the concept of *convex transform order* to compare two distributions according to skewness as defined below.

*Definition 2.14. * is said to be smaller than in *the convex transform order*, denoted by if and only if, is convex in on the support of .

If , then is more skewed than as explained in van Zwet [48] and Marshall et al. [34]. The convex transform order is also called more IFR (increasing failure rate) order in reliability theory, since when and exist, the convexity of means that is increasing in . Thus can be interpreted to mean that ages faster than in some sense. Gamma distributions are ordered according to the convex transform order in terms of their shape parameters.

Another well-known partial order to compare the skewness of two probability distributions is *star order*.

*Definition 2.15. * is said to be smaller than in *the star order*, denoted by (or ) if the function is star shaped in the sense that is increasing in on the support of .

The star order is also called more IFRA (increasing failure rate in average) order in reliability theory, since the average failure of at is Thus can be interpreted in terms of average failure rates as is increasing in . Note that has an increasing failure rate if and only if is star-ordered with respect to exponential distribution.

The function is known as Lorenz curve in the economics literature. It is often used to express inequality in incomes. Based on Lorenz curve, the Lorenz order has been proposed in Economics to compare income inequalities.

*Definition 2.16. * is said to be smaller than in the Lorenz order, denoted by , if

It is known in the literature (Marshall et al. [34] ) that, where denote the coefficient of variation of .

All the above partial orders are scale invariant. A good discussion of the star order and Lorenz order can be found in Barlow and Proschan [49], Marshall et al. [34], and Kochar [50].

##### 2.4. Dependence Concepts and Orderings

Let be a bivariate ramdom vector with joint distribution and with density function . Lehmann [51] introduced many partial orderings for dependence.

*Definition 2.17. *(a) and are *dependent* or *likelihood* dependent if their joint density is totally positive of order 2 in and , or more precisely if
whenever , .

(b) is stochastically increasing in , denoted by SI, if is increasing in for all , or equivalently,

(c) is right tail increasing in , denoted by RTI, if is increasing in for all .

(d) is left tail decreasing in , denoted by LTD, if is decreasing in for all .

Finally, random variables and are *associated* (written ) if for all pairs of increasing binary functions and . As shown in Barlow and Proschan [49, page 143] the following chain of implications holds among the above notions of positive dependence:

These concepts of bivariate dependence can be easily extended to the multivariate case. A function is said to be *multivariate total positivity of order 2* (denoted by ) if
where and . Random variables are said to be dependent if their joint density function is . It is shown in Kemperman [52] (see also Block and Ting, [53]) that if the support of a random vector is a lattice (i.e., if and are in the support of then so are and ) then is if and only if, its density function is in each pair of its variables when the other variables are held fixed. See Karlin and Rinott [54] for more details on properties of functions. Also refer to Joe [55] and Nelsen [56] for a comprehensive discussion.

Observing that when and are continuous, inequality (2.27) can be written as where stands for the th quantile of the marginal distribution of , and denotes the conditional distribution of given . Avérous et al. [57] proposed the following definition to measure the relative degree of monotone dependence between two pairs of bivariate random variables and .

*Definition 2.18. * is said to be less *stochastic increasing* in than is in , denoted by , if and only if, for , and ,
where stands for the th quantile of the marginal distribution of , and denotes the conditional distribution of given , for .

Dolati et al. [58] proposed the following weaker dependence order based on RTI criteria, called more RTI order.

*Definition 2.19. * is said to be less *right-tail increasing* (RTI) in than is in , denoted by , if and only if, for , and ,
where stands for the th quantile of the marginal distribution of , and denotes the conditional distribution of given , for .

It is easy to see that both more SI order and more RTI order are *copula-based* orders, and more SI order implies more RTI order which in turn implies more concordance ordering (i.e., the two copulas are ordered). For the concept of *copula*, please refer to Nelsen [59] for more details.

As observed in Avérous et al. [57] and Genest et al. [60], there is a close connection between the above concepts of more dependence and the notion of dispersive ordering.

##### 2.5. Notions of Majorization and Related Orderings

One of the basic tools in establishing various inequalities in statistics and probability is the notion of majorization.

Let denote the increasing arrangement of the components of the vector .

*Definition 2.20. *The vector is said to majorize the vector (written ) if for and .

Functions that preserve the majorization ordering are called Schur-convex functions. The vector is said to majorize the vector weakly (written ) if for . Marshall et al. [34] provides extensive and comprehensive details on the theory of majorization and its applications in statistics.

Bon and Pǎltǎnea [24] have considered a preorder on , which they call as a*-larger order*.

*Definition 2.21. *A vector in is said to be -larger than another vector also in (written ) if .

Let denote the vector of logarithms of the coordinates of . It is easy to verify that It is known that for all concave functions (cf. [34]). From this and (2.33), it follows that when The converse is, however, not true. For example, the vectors but majorization does not hold between these two vectors.

#### 3. Stochastic Comparisons of Order Statistics in the One-Sample Problem

In this section, we compare order statistics from a single sample according to various stochastic orders.

##### 3.1. Magnitude Orderings between Order Statistics in the One-Sample Problem

Let be a set of independent and identically distributed random variables. It is easy to see that , for all . Boland et al. [18] extended this result from usual stochastic order to hazard rate order. Using the definition of likelihood ratio ordering, it is not hard to prove that for . Shaked and Shanthikumar [2] considered the problem of comparing order statistics from samples with possibly unequal sample sizes. They showed that if random variables ’s are iid, then and . Raqab and Amin [61] strengthened this result and proved that , whenever and .

It is interesting to investigate the above stochastic relations among order statistics when the random variables are independent but not identically distributed. Boland et al. [18] showed that if random variables are independent and , , then , . They also proved that if ’s are independent and , , then , . Assuming , Bapat an Kochar [17] proved that , .

##### 3.2. Variability Orderings between Order Statistics in the One-Sample Problem

David and Groeneveld [62] proved that if ’s are iid random variables with a common decreasing failure rate (DFR) distribution, then , for . Kochar [11] strengthened this result to prove that under the same conditions, , . Khaledi and Kochar [27] further proved that if ’s are iid with a distribution, then , whenever and .

##### 3.3. Skewness Orderings between Order Statistics in the One-Sample Problem

Arnold and Villasenor [63], Arnold and Nagaraja [64], Wilfling [65], and Kleiber [66], among others, studied Lorenz order relations between order statistics from uniform and other distributions. In particular, Arnold and Villasenor [63] proved the following result on Lorenz ordering between the order statistics from uniform distributions.

Theorem 3.1. *Let denote the th order statistic of a random sample of size from a uniform distribution over , . Then *(a)*, *(b)*, *(c)*, *(d)*.
*

The last inequality may be described as “sample medians exhibit less variability as sample size increases.” Arnold and Villasenor [63] wonder about the conditions on , , , and under which holds.

Kochar [67] answered this question in Theorem 3.2 below in which sufficient conditions on the parent distribution are obtained under which holds. Many of the previously known results follow from this general result as particular cases as star ordering implies Lorenz ordering.

Theorem 3.2. *For , let denote the th order statistic of a random sample of size from a distribution with reverse hazard rate . If
**
then for and ,
*

*Remark 3.3. *Arnold and Villasenor [63] mention (3.2) as a sufficient condition for the relation
to hold. We have a more general result.

The following theorem is also proved in in Kochar [67].

Theorem 3.4. *For , let denote the th order statistic of a random sample of size from a distribution with its hazard rate satisfying the condition,
**Then
*

The above theorems immediately lead to the following result because of the relation between star ordering and Lorenz ordering.

Corollary 3.5. *If for , denotes the th order statistic of a random sample of size from a distribution satisfying *(a)*condition (3.2), then
for and . *(b)*condition (3.5), then
for and . *

*Example 3.6. *The condition (3.2) is satisfied by the power function distribution with distribution function, , , . Therefore, the conclusions of Theorem 3.2 and Corollary 3.5(a) hold for this distribution. Arnold and Villasenor [63] also conjectured that for this distribution,
Its proof immediately follows from Theorem 3.2 and Corollary 3.5.

*Example 3.7. *Let has Pareto distribution with , , . Condition (3.5) is satisfied by this distribution. Hence the conclusions of Theorem 3.4 and Corollary 3.5(b) hold for this distribution.

##### 3.4. Dependence among Order Statistics in the One-Sample Problem

Assuming that are independent and identically random variables, Bickel [68] showed that When the parent distribution has an increasing hazard rate and a decreasing reverse hazard rate, Tukey [69] showed that for either and ; or and .

Kim and David [70] proved that if both the hazard and the reverse hazard rates of the ’s are increasing, then inequality (3.11) remains valid when and ; However, the inequality (3.11) is reversed when and .

Avérous et al. [57] used the more SI concept to study the relative degree of dependence between two pairs of random variables. They proved the following result.

Theorem 3.8. *Let and be the order statistics associated with two independent random samples of sizes and from the same continuous distribution. Then, for , , and , , , , one hold that
*

As a direct consequence, we have the following result.

Corollary 3.9. *Let be order statistics from the same continuous distribution. Then, *(1)*for ; *(2)*for ; *(3)*for . *

It can be seen from the above result that the dependence between the components of a pair of order statistics decreases in the sense of SI ordering as and get further apart.

*Remark 3.10. * Since the copula of a pair of order statistics of a random sample is independent of the parent distribution and since the concept of more SI is copula based, it follows that in Theorem 3.8 and Corollary 3.9, the two samples could be from different distributions.

As explained in Avérous et al. [57], the following result follows immediately from Theorem 3.8.

Corollary 3.11. * Under the assumptions of Theorem 3.8, one has
**
where represents Spearman’s rho, Kendall’s tau, Gini’s coefficient, or indeed any other copula-based measure of concordance satisfying the axioms of Scarsini [71].*

###### 3.4.1. Kendall’s and Spearman’s for Order Statistics

Two popular nonparametric measures of association for bivariate random variables are Kendall’s and Spearman’s , which measure different aspects of the dependence structure. In terms of dependence properties, Spearman’s is a measure of average quadrant dependence, while Kendall’s is a measure of average likelihood ratio dependence (cf. [59, 72]).

Avérous et al. [57] made an important observation that in the case of a random sample from a continuous distribution with cdf , the copula of a pair of order statistics is independent of the parent distribution . As a result the value of any copula-based measure of dependence like Kendall’s tau or Spearman’s coefficient for any pair of order statistics will be the same for all continuous distributions Schmitz [73] derived the following formulas: X. Li and Z. Li [74] proved a conjecture in Schmitz [73] that

Avérous et al. [57] used a combinatorial approach to prove the following formula of Kendall’s for any pair of order statistics from the same continuous distribution:

Subsequently, Chen [75] developed three new formulas to compute : where can be computed by any one of the following formulas. (1)Formula 1: (2)Formula 2: (3)Formula 3: He further showed the following compound inequality: from which, (3.15) follows immediately.

Recently, Navarro and Balakrishnan [76] have also studied this problem and have obtained alternate expressions for computing these measures of dependence.

###### 3.4.2. Dependence Orderings among Order Statistics When the Observations Are Not Identically Distributed

The proportional hazard rates (PHRs) model is a popular model in survival analysis, reliability theory, and other branches of statistics. Let denote the lifetime of the component of a reliability system with survival function , . Then they have proportional hazard rates if there exist constants and a (cumulative hazard) function such that for . Clearly then the hazard rate of is (assuming it exists). An example of such a situation is when the components have independent exponential lifetimes with respective hazard rates .

Boland et al. [10] studied in detail the dependence properties of order statistics when the observations are not necessarily identically distributed. They proved the following dependence result for the PHR model.

Theorem 3.12. * Let be independent random variables with differentiable densities and follow the PHR model on an interval. Then for , is SI in .*

They also gave a counterexample to illustrate that, in general, is not SI in . However, they showed that for , is RTI in .

Theorem 3.13. * Let be independent random variables. Then for any , RTI and LTD. *

They showed with the help of a counter example that, in general, the relation RTI may not hold for .

This topic has been further developed by Hu and Xie [77], where they exploited the negative dependence of occupancy numbers in the balls and bins experiment. They proved the following result.

Theorem 3.14. *Let be independent random variables. For , and ,*(1)*if , then
is increasing in ; *(2)*if event is either or , then the conditional probability in (3.22) is decreasing in for each ; and *(3)*if each has a continuous distribution function, and if is either or , then (3.22) is decreasing in for each , where . *

Dubhashi and Häggström [78] further extended the above result to multivariate stochastic comparisons.

Theorem 3.15. * Let be independent random variables. Then
*

Subsequently, Theorem 3.14 was further extended by Hu and Chen [79] as follows.

Theorem 3.16. *Let be independent random variables. *(1)*If , then
is increasing in for all ; *(2)*If , then
is decreasing in for all . *

Recently, Zhuang et al. [80] discussed the dependence among order statistics in the sense of multivariate stochastic comparisons, which extends the results in Hu and Chen [79] and Dubhashi and Häggström [78].

Theorem 3.17. * Let be independent random variables. *(1)*If , then
*(2)*If , then
*

#### 4. Stochastic Comparisons of Order Statistics in the Two-Sample Problem

Let be a set of independent random variables and be another set of independent random variables. Ross [31] proved that if , , then . A consequence of this result is that . Lynch et al. [81] generalized this result from stochastic ordering to hazard rate ordering though under a stronger condition. They proved that if , , then , . A similar result for likelihood ratio ordering has been proved by Chan et al. [82]. They proved that if , , then , . Lillo et al. [83] in the following theorem strengthened this result to the case when the number of ’s and ’s is not necessarily equal.

Theorem 4.1. *Let be independent random variables and let be another set of independent random variables, all having absolutely continuous distributions. Then for all implies whenever and . *

In the next theorem we establish dispersive ordering between order statistics when the random samples are drawn from different distributions.

Theorem 4.2. *Let be a random sample of size from a continuous distribution and let be a random sample of size from another continuous distribution . If either or is , then
*

*Proof. * Let be a DFR distribution. The proof for the case when is DFR is similar. Khaledi and Kochar [27] proved that for and . Bartoszewicz [42] proved that if then . Combining these we get the required result.

Since the property together with the condition that either or is DFR implies that , we get the following result from the above theorem.

Corollary 4.3. *Let be a random sample of size from a continuous distribution and be a random sample of size from another continuous distribution . If either or is , then
*

##### 4.1. Stochastic Comparisons of Order Statistics When Observations Are Heterogeneous

Many researchers have investigated the effect on the survival function, the hazard rate function and other characteristics of the time to failure of this system when we switch the vector to another vector say in the PHR model. Pledger and Proschan [15], for the first time, studied this problem and proved the following interesting result among many other results.

Theorem 4.4. *Let and be two random vectors of independent lifetimes with proportional hazards with and as the constants of proportionality. Then for *

Proschan and Sethuraman [16] strengthened this result from component wise stochastic ordering to multivariate stochastic ordering. That is, under the assumptions of Theorem 4.4, they proved that

These results show that if the parent observations have proportional hazard rates, then their order statistics are stochastically larger when the proportionality parameters are dispersed in the sense of majorization.

Boland et al. [18] showed with the help of the following counterexample that (4.3) can not be strengthened from stochastic ordering to hazard rate ordering for when .

*Example 4.5. * Let be independent exponential random vector with hazard rate vector and be independent exponential random vector with hazard rate vector . It is easily seen that
However,
Hence,

Khaledi and Kochar [28] proved the following result for the largest order statistics under *p*-larger ordering.

Theorem 4.6. * Let and be two random vectors of independent lifetimes with proportional hazards with and as the constants of proportionality. Then
*

The following example due to Khaledi and Kochar [28] shows that Theorem 4.6 may not hold for other order statistics.

*Example 4.7. * Let be independent exponential random vector with hazard rate vector and be independent exponential random vector with hazard rate vector . It is easily seen that
However,
which implies

Khaledi and Kochar [28] also showed that in, general,

Torrado and Veerman [84] formulated a simple criteria in terms of the ’s to see whether two -out-of- systems are stochastically ordered.

###### 4.1.1. Comparison of Order Statistics from Heterogenous and Homogeneous Samples

Let be continuous independent random variables such that the c.d.f. of is , , where ’s are not necessarily all equal. Let be a random sample of size from a distribution with c.d.f . Sen [7] proved a very general result that whereas for the first order statistic, , the inequality is reversed for the largest order statistic, that is, .

Dykstra et al. [22] studied the above problem for the exponential distributions. They proved the following stronger result on hazard rate ordering for largest order statistics.

Theorem 4.8. *Let be independent exponential random variables with having hazard rate , ; and let be a random sample of size from an exponential distribution with common hazard rate . Then for
**
as well as
*

Khaledi and Kochar [28] proved that (4.14) and (4.15) continue to hold with , the geometric mean of the s. This result gives better bounds on the hazard rate and measures of variability of a parallel system with heterogeneous exponential components.

Khaledi and Kochar [28] also extended the results (4.14) and (4.15) from the exponential case to the PHR model as stated in the next theorem.

Theorem 4.9. *Let be independent random variables with having survival function , . Let be a random sample with common population survival distribution , where , then *(i)*; *(ii)* if has decreasing failure rate (DFR). *

These results give nice bounds on the hazard rate and the variance of a parallel system with components which are independent following the PHR model in terms of the case when they are independent and identically distributed.

Kochar and Xu [85] further strengthened the conclusion (4.14) of Theorem 4.8 of Dykstra et al. [22] from hazard rate ordering to likelihood ratio order for the PHR model, which includes exponential as a special case.

Theorem 4.10. *Let be independent random variables with having survival function , . Let be a random sample with common population survival distribution , where , then
*

The following example due to Kochar and Xu [85] shows that the result of Khaledi and Kochar [27] above cannot be strengthened from the hazard rate order to the likelihood ratio order under *p*-larger ordering.

*Example 4.11. *Let be independent exponential random variables with having hazard rate , , and be a random sample of size from an exponential distribution with common hazard rate . Then, the reversed hazard rate of is
Similarly, the reversed hazard rate of is
Let , , and , then,
Thus,
which implies that

Bon and Pǎltǎnea [86] proved the following result on stochastic ordering between th order statistics from heterogeneous and homogeneous samples from exponential distributions.

Theorem 4.12. *Let be independent exponential random variables with hazard rates , . Let be independent exponential random variables with a common hazard rate . Then,
**
where
**
a result which also follows from Khaledi and Kochar [27] for the case when .*

Kochar and Xu [87] proved the following results on excess wealth ordering between largest order statistics from heterogeneous and homogeneous samples.

Theorem 4.13. * Let be independent exponential random variables with having hazard rate , . Let be a random sample of size from an exponential distribution with common hazard rate . Then,
**
where
*

Theorem 4.14. * Let be independent random variables with having survival function , , and let be another random sample with the common survival function . If is DFR, then
**
where is given in (4.24).*

As a consequence, we have the following result.

Corollary 4.15. * Let be independent random variables with having survival function , , and let be another random sample with the common survival function . If is DFR, then
**
where is given in (4.24).*

*Remark 4.16. *In the special case of exponential distribution, , Khaledi and Kochar [28] gave a lower bound on the variance of in terms of the geometric mean,
Since the right spread order implies the variance being ordered (cf. Shaked and Shanthikumar [2]), Corollary 4.15 provides a better lower bound for the variance,
as .

###### 4.1.2. Skewness Comparisons of Order Statistics from Two Samples

Kochar and Xu [87] for the first time studied the problem of comparing order statistics according to shape orderings like convex ordering and star ordering. They proved that a parallel system with heterogenous exponential component lifetimes is more skewed (according to convex transform order) than the one with i.i.d. exponential components. This is stated in the following theorem.

Theorem 4.17. *Let be independent exponential random variables with having hazard rate , . Let be a random sample of size from an exponential distribution with common hazard rate . Then,
*

*Remark 4.18. *Theorem 4.17 means that a parallel system with homogeneous exponential components ages faster than than a system with heterogenous exponential components in the sense of “more IFR” property. Note that a parallel system with homogeneous exponential components is IFR (Barlow and Proschan [49]). However, a parallel system with heterogenous exponential components may not be IFR.

*Remark 4.19. *It is interesting to note that, unlike the magnitude and variability orders, no restriction on the parameters is needed for Theorem 4.17 to hold as the convex transform order is scale invariant. Intuitively, due to the heterogeneity, the largest order statistic from a heterogenous sample will be more skewed than that from a homogeneous sample. Theorem 4.17 confirms this fact for exponential samples.

The following result, which is of independent interest in economics, is a direct consequence of Theorem 4.17

Corollary 4.20. * Let be independent exponential random variables with having hazard rate , . Let be a random sample of size from an exponential distribution with common hazard rate . Then,
*

From Barlow and Proschan [49, page 60], it follows that, Using (2.25), one get the following lower bound on the coefficient of variation of ,

Kochar and Xu [88] extended Theorem 4.17 to the PHR family as follows.

Theorem 4.21. *Let be independent random variables with having survival function , , and let be a random sample from a distribution with the common survival distribution where , the geometric mean of ’s. If
**
then
**
where is the cumulative hazard rate function, and is the hazard rate function of .*

It can be seen that distributions like Weibull, Lomax, and Pareto satisfy the condition (4.34). Kochar and Xu [88] proved the following result on star ordering for th order statistics.

Theorem 4.22. * Let be i.i.d. random variables with the common survival distribution , and let be another set of i.i.d. random variables with common survival distribution , and let be a random sample from a distribution with common survival distribution where . If
**
then
**
and hence,
**
where
**
and . *

###### 4.1.3. Stochastic Comparisons of Second Order Statistics

In the next theorem, we give a collection of results on various types of stochastic orders for the second order statistics.

Theorem 4.23. *Let be independent exponential random variables with parameters , respectively; and let be independent random variables with common hazard rate . *(a)*(Pǎltǎnea [89])
*(b)* (Zhao and Balakrishnan [90])
*