new results on stochastic comparisons of two-component series and parallel systems

8
Statistics and Probability Letters 82 (2012) 283–290 Contents lists available at SciVerse ScienceDirect Statistics and Probability Letters journal homepage: www.elsevier.com/locate/stapro New results on stochastic comparisons of two-component series and parallel systems Neeraj Misra, Amit Kumar Misra Department of Mathematics and Statistics, Indian Institute of Technology Kanpur, Kanpur-208016, India article info Article history: Received 27 June 2011 Received in revised form 14 October 2011 Accepted 14 October 2011 Available online 20 October 2011 Keywords: Hazard rate order Likelihood ratio order Mean residual life order Reversed hazard rate order Usual stochastic order abstract Let X 1 , X 2 , and X 3 be independent random variables with absolutely continuous distributions having the common support [0, ). We show that if X 1 hr[mrl,lr] X 3 and X 2 hr[mrl,lr] X 3 , then max{X 1 , X 2 }≤ hr[mrl,lr] max{X 1 , X 3 }. We also show that if X 2 rh[lr] X 1 and X 2 rh[lr] X 3 , then min{X 1 , X 2 }≤ rh[lr] min{X 1 , X 3 }. These results generalize and extend some of the results given in Shaked and Shanthikumar (2007, Example 1.C.36, p. 56), Joo and Mi (2010), and Da et al. (2010). © 2011 Elsevier B.V. All rights reserved. 1. Introduction A series (parallel) system functions if and only if each (at least one) of its components function (functions). Consider n components C 1 ,..., C n having random lifetimes Y 1 ,..., Y n , respectively. Then the lifetime of a series (parallel) system constructed from these components is given by min{Y 1 ,..., Y n } (max{Y 1 ,..., Y n }). Performance of two systems constructed from different sets of components can be compared through stochastic comparisons of the corresponding system lifetimes. Normally a stochastic comparison of lifetimes of two different systems is made with respect to one of the stochastic orders between lifetimes of systems. For an account on characterizations and properties of various stochastic orders, one may refer to Shaked and Shanthikumar (2007) and Müller and Stoyan (2002). A vast literature on stochastic comparisons of lifetimes of series and parallel systems exists. See, for example, Boland et al. (1994), Dykstra et al. (1997), Khaledi and Kochar (2000), Da et al. (2010), Joo and Mi (2010), Zhao and Balakrishnan (2011), and references cited therein. In this paper, we will derive some new results on stochastic comparisons of two-component series and parallel systems. These results generalize and extend some of the results known in the literature. First, we recall definitions of various stochastic orders relevant to the context of this paper. Let X and Y be random variables with the distribution functions F and G, the probability density functions f and g , the hazard functions r and µ, the reversed hazard functions ˜ r and ˜ µ, and the mean residual life functions m and l, respectively. Let ¯ F 1 F and ¯ G 1 G denote the corresponding survival functions. Suppose that F and G have the common support R + =[0, ) and {t : f (t )> 0}={t : g (t )> 0}= R + . When we say that a function is increasing (decreasing), it means that the function is non-decreasing (non-increasing). Unless otherwise stated, all the random variables considered in this study will be assumed to have absolutely continuous distributions. Corresponding author. Tel.: +91 9839425105; fax: +91 512 2597500. E-mail addresses: [email protected] (N. Misra), [email protected], [email protected] (A.K. Misra). 0167-7152/$ – see front matter © 2011 Elsevier B.V. All rights reserved. doi:10.1016/j.spl.2011.10.010

Upload: neeraj-misra

Post on 28-Oct-2016

215 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: New results on stochastic comparisons of two-component series and parallel systems

Statistics and Probability Letters 82 (2012) 283–290

Contents lists available at SciVerse ScienceDirect

Statistics and Probability Letters

journal homepage: www.elsevier.com/locate/stapro

New results on stochastic comparisons of two-component series andparallel systemsNeeraj Misra, Amit Kumar Misra ∗

Department of Mathematics and Statistics, Indian Institute of Technology Kanpur, Kanpur-208016, India

a r t i c l e i n f o

Article history:Received 27 June 2011Received in revised form 14 October 2011Accepted 14 October 2011Available online 20 October 2011

Keywords:Hazard rate orderLikelihood ratio orderMean residual life orderReversed hazard rate orderUsual stochastic order

a b s t r a c t

Let X1, X2, and X3 be independent random variables with absolutely continuousdistributions having the common support [0,∞). We show that if X1 ≤hr[mrl,lr] X3 andX2 ≤hr[mrl,lr] X3, then max{X1, X2} ≤hr[mrl,lr] max{X1, X3}. We also show that if X2 ≤rh[lr] X1and X2 ≤rh[lr] X3, then min{X1, X2} ≤rh[lr] min{X1, X3}. These results generalize and extendsome of the results given in Shaked and Shanthikumar (2007, Example 1.C.36, p. 56), Jooand Mi (2010), and Da et al. (2010).

© 2011 Elsevier B.V. All rights reserved.

1. Introduction

A series (parallel) system functions if and only if each (at least one) of its components function (functions). Considern components C1, . . . , Cn having random lifetimes Y1, . . . , Yn, respectively. Then the lifetime of a series (parallel)system constructed from these components is given by min{Y1, . . . , Yn} (max{Y1, . . . , Yn}). Performance of two systemsconstructed from different sets of components can be compared through stochastic comparisons of the correspondingsystem lifetimes. Normally a stochastic comparison of lifetimes of two different systems is made with respect to one ofthe stochastic orders between lifetimes of systems. For an account on characterizations and properties of various stochasticorders, one may refer to Shaked and Shanthikumar (2007) and Müller and Stoyan (2002). A vast literature on stochasticcomparisons of lifetimes of series and parallel systems exists. See, for example, Boland et al. (1994), Dykstra et al. (1997),Khaledi and Kochar (2000), Da et al. (2010), Joo andMi (2010), Zhao and Balakrishnan (2011), and references cited therein. Inthis paper, wewill derive some new results on stochastic comparisons of two-component series and parallel systems. Theseresults generalize and extend some of the results known in the literature. First, we recall definitions of various stochasticorders relevant to the context of this paper.

Let X and Y be random variables with the distribution functions F and G, the probability density functions f and g , thehazard functions r and µ, the reversed hazard functions r and µ, and the mean residual life functionsm and l, respectively.Let F ≡ 1− F and G ≡ 1− G denote the corresponding survival functions. Suppose that F and G have the common supportR+ = [0,∞) and {t : f (t) > 0} = {t : g(t) > 0} = R+. When we say that a function is increasing (decreasing), it meansthat the function is non-decreasing (non-increasing). Unless otherwise stated, all the random variables considered in thisstudy will be assumed to have absolutely continuous distributions.

∗ Corresponding author. Tel.: +91 9839425105; fax: +91 512 2597500.E-mail addresses: [email protected] (N. Misra), [email protected], [email protected] (A.K. Misra).

0167-7152/$ – see front matter© 2011 Elsevier B.V. All rights reserved.doi:10.1016/j.spl.2011.10.010

Page 2: New results on stochastic comparisons of two-component series and parallel systems

284 N. Misra, A.K. Misra / Statistics and Probability Letters 82 (2012) 283–290

Definition 1.1. X is said to be smaller than Y in the

(i) likelihood ratio order (written as X ≤lr Y ) if g(t)/f (t) is increasing in t ∈ R+;(ii) usual stochastic order (written as X ≤st Y ) if F(t) ≤ G(t), ∀t ∈ R+;(iii) hazard rate order (written as X ≤hr Y ) if G(t)/F(t) is increasing in t ∈ R+, or equivalently if r(t) ≥ µ(t), ∀t ∈ R+;(iv) reversed hazard rate order (written as X ≤rh Y ) if G(t)/F(t) is increasing in t ∈ (0,∞), or equivalently if r(t) ≤

µ(t),∀t ∈ (0,∞);(v) mean residual order (written as X ≤mrl Y ) if

t G(x) dx/

t F(x) dx is increasing in t ∈ R+, or equivalently if m(t) ≤

l(t),∀t ∈ R+;(vi) increasing convex order (written as X ≤icx Y ) if

t G(x) dx ≥

t F(x) dx for all t ∈ R+.

In this article we will mainly focus on two-component series and parallel systems. Let X1, X2, and X3 be statisticallyindependent and nonnegative random variables. Let Pi (Si) be the parallel (series) system consisting of two componentshaving random lifetimes X1 and Xi, i = 2, 3. Then the lifetime of Pi (Si) is given by max{X1, Xi} (min{X1, Xi}) , i = 2, 3. Forthe case when Xi has exponential distribution with hazard rate λi(>0), i = 1, 2, 3, Joo and Mi (2010, Theorem 2.2) provedthat

λ3 ≤ λ2 ≤ λ1 ⇒ max{X1, X2} ≤hr max{X1, X3}.

Daet al. (2010) generalized this result to the parallel systems consisting of general (not necessarily exponentially distributed)independent components by proving that

X1 ≤hr X2 ≤hr X3 ⇒ max{X1, X2} ≤hr max{X1, X3}. (1.1)

For general independent components Da et al. (2010) also showed that

X1 ≤mrl X2 ≤mrl X3 ⇒ max{X1, X2} ≤mrl max{X1, X3}. (1.2)

Shaked and Shanthikumar (2007, Example 1.C.36, p. 56) provide a similar result for the likelihood ratio order. Specifically,for general independent components, this result says that:

X2 ≤lr X1 ≤lr X3 ⇒ max{X1, X2} ≤lr max{X1, X3} and min{X1, X2} ≤lr min{X1, X3}. (1.3)

The purpose of this article is to strengthen results (1.1)–(1.3) by providingweaker conditions underwhich these results hold.We will also derive some additional results. The outline of the paper is as follows. Results on parallel systems are discussedin Section 2. In this section, we will prove that

X1 ≤hr[mrl,lr] X3 and X2 ≤hr[mrl,lr] X3 ⇒ max{X1, X2} ≤hr[mrl,lr] max{X1, X3},

thereby generalizing results (1.1) and (1.2) given by Da et al. (2010) and the result (1.3) given in Shaked and Shanthikumar(2007). Section 3 deals with series systems where we show that

X2 ≤rh[lr] X1 and X2 ≤rh[lr] X3 ⇒ min{X1, X2} ≤rh[lr] min{X1, X3},

thereby extending and generalizing the result (1.3) given in Shaked and Shanthikumar (2007).

2. Stochastic comparisons of parallel systems

Let X1, X2, and X3 denote independent random variables with absolutely continuous distributions having the commonsupport R+. Let Fi, fi, Fi, ri, ri, and mi denote, respectively, the distribution function, the probability density function, thesurvival function, the hazard function, the reversed hazard function, and the mean residual life function of Xi, i = 1, 2, 3.In this section we will compare random variables max{X1, X2} and max{X1, X3} with respect to the hazard rate order, themean residual life order, and the likelihood ratio order.

In the following theorem we generalize and extend the result (1.1) proved by Da et al. (2010).

Theorem 2.1. Suppose that one (or both) of the following conditions holds:

(a) X1 ≤hr X3 and X2 ≤hr X3;(b) X2 ≤hr X3 and r1(t)F3(t) ≥ r3(t)F1(t) for all t ∈ R+.

Then

max{X1, X2} ≤hr max{X1, X3}.

Proof. For t ∈ R+, consider

h1(t) =P (max{X1, X3} > t)P (max{X1, X2} > t)

=F1(t)+ F3(t)− F1(t)F3(t)F1(t)+ F2(t)− F1(t)F2(t)

=N1(t)D1(t)

, say.

Page 3: New results on stochastic comparisons of two-component series and parallel systems

N. Misra, A.K. Misra / Statistics and Probability Letters 82 (2012) 283–290 285

It can be easily verified that

D21(t)

ddt

h1(t) = r1(t)F1(t)[F3(t)− F2(t)] + r2(t)F2(t)[1 − F1(t)][F1(t)+ F3(t)− F1(t)F3(t)]

− r3(t)F3(t)[1 − F1(t)][F1(t)+ F2(t)− F1(t)F2(t)].

(a) Since r2(t) ≥ r3(t) for all t ∈ R+, it follows that

D21(t)

ddt

h1(t) ≥ F1(t)[F3(t)− F2(t)][r1(t)− r3(t)F1(t)] (2.1)

≥ 0 (since r1(t) ≥ r3(t), and X2 ≤hr X3(⇒ F2(t) ≤ F3(t))),

i.e. h1(t) is increasing in t ∈ R+.(b) From (2.1), we can write

D21(t)

ddt

h1(t) ≥ F1(t)[F3(t)− F2(t)][r1(t)F3(t)+ r1(t)F3(t)− r3(t)F1(t)]

≥ 0, ∀t ∈ R+,

where the last inequality follows from the assumption that r1(t)F3(t) ≥ r3(t)F1(t) for all t ∈ R+, and the fact thatX2 ≤hr X3 implies F2(t) ≤ F3(t) for all t ∈ R+. Thus, h1(t) is increasing in t ∈ R+. Hence the result follows. �

Remark 2.1. (i) From (2.1) it is clear that the conclusion of Theorem 2.1 remains true if X2 ≤hr X3 and r1(t) ≥ r3(t)F1(t) forall t ∈ R+.

(ii) The condition r1(t)F3(t) ≥ r3(t)F1(t),∀t ∈ R+, is equivalent to the condition r1(t)F3(t) ≥ r3(t)F1(t),∀t ∈ (0,∞).Therefore, the conclusion of Theorem 2.1 also holds if X2 ≤hr X3 and r1(t)F3(t) ≥ r3(t)F1(t) for all t ∈ (0,∞).

As a consequence of Theorem 2.1, we have the following corollary.

Corollary 2.1. Suppose that one (or more) of the following conditions holds:(a) X1 ≤hr X2 ≤hr X3;(b) X2 ≤hr X1 ≤hr X3;(c) X1 ≤st X3, X2 ≤hr X3, and r3(t)/r1(t) is decreasing in t ∈ R+;(d) X3 ≤st X1, X2 ≤hr X3, and r1(t)/r3(t) is decreasing in t ∈ (0,∞).Then

max{X1, X2} ≤hr max{X1, X3}.

Proof. Parts (a) and (b) are immediate from Theorem 2.1(a).(c) On following the proof of Corollary 2.3 of Misra et al. (2011), we can show that if X1 ≤st X3 and r3(t)/r1(t) is decreasing

in t ∈ R+, then r1(t)F3(t) ≥ r3(t)F1(t) for all t ∈ R+. Now, the result follows from Theorem 2.1(b).(d) Using the technique used in the proof of Corollary 3.2 ofMisra et al. (2011), we can show that if X3 ≤st X1 and r1(t)/r3(t) is

decreasing in t ∈ (0,∞), then r1(t)F3(t) ≥ r3(t)F1(t) for all t ∈ (0,∞). Now, the result follows fromRemark 2.1(ii). �

Recall that the result in Corollary 2.1(a) was proved by Da et al. (2010). The following examples illustrate a few situationson comparative usefulness of results given in Theorem 2.1 and Corollary 2.1.

Example 2.1. Let r1(t) = 36t + 10, r2(t) = 43t + 3, and r3(t) = 35t + 3, t ∈ R+. Obviously, r1(t) ≥ r3(t), r2(t) ≥ r3(t)for all t ∈ R+, and r1(t) ≥ (≤)r2(t) for all 0 ≤ t ≤ 1 (t ≥ 1). Moreover, it can be easily verified that r3(t)/r1(t) is strictlyincreasing on R+. For t ∈ R+, consider

φ1(t) = r1(t)F3(t)− r3(t)F1(t)

= (36t + 10)(1 − e−352 t2−3t)− (35t + 3)(1 − e−18t2−10t).

Then, φ1(0.01) = −0.000482848 and φ1(0.02) = 0.00227857. Thus, the conditions of Theorem 2.1(a) are satisfied but theconditions of Theorem 2.1(b) and Corollary 2.1(a)–(d) are not satisfied.

Example 2.2. Let r1(t) = 5t + 10, r2(t) = 12t + 3, and r3(t) = 3t + 3, t ∈ R+. Obviously, r1(t) ≥ r3(t), r2(t) ≥ r3(t)for all t ∈ R+, and r1(t) ≥ (≤) r2(t) for all 0 ≤ t ≤ 1 (t ≥ 1). Moreover, r3(t)/r1(t) is strictly increasing in t ∈ R+ andφ1(t) ≥ 0, ∀t ∈ R+ (see Misra et al., 2011, Example 2.2). Thus, the conditions of Theorem 2.1(a) and (b) are satisfied butthe conditions of Corollary 2.1(a)–(d) are not satisfied.

Example 2.3. Let F1(t) = (1 − e−t)2, F2(t) = 1 − e−2t , and F3(t) = 1 − e−t , t ∈ R+. Then we have r1(t) = 2/(et − 1) andr3(t) = 1/(et − 1), t ∈ (0,∞). We can see that X3 ≤st X1, X2 ≤hr X3, and r1(t)/r3(t) is decreasing in t ∈ (0,∞). Thus, theconditions of Corollary 2.1(d) are satisfied but the conditions of Theorem 2.1(a) and Corollary 2.1(a)–(c) are not satisfied.Moreover, on using Remark 2.1(ii), it is evident that the conditions of Theorem 2.1(b) are also satisfied.

Page 4: New results on stochastic comparisons of two-component series and parallel systems

286 N. Misra, A.K. Misra / Statistics and Probability Letters 82 (2012) 283–290

Remark 2.2. The referee in its report on an earlier version of this paper raised the question whether the conditions ofmonotonicity of the ratio of the hazard rates and that of the reversed hazard rates in Corollary 2.1(c)–(d) are necessaryfor the conclusion of Corollary 2.1 to hold. Example 2.1 suggests that the condition of monotonicity of the ratio of the hazardrates in Corollary 2.1(c) is not necessary for the conclusion of Corollary 2.1 to hold. The following example illustrates that thecondition of monotonicity of the ratio of the reversed hazard rates in Corollary 2.1(d) is also not necessary for the conclusionof Corollary 2.1 to hold.

Example 2.4. Let F1(t) = (1 − e−λ1t)α1 , F2(t) = F3(t) = (1 − e−λt)α, t ∈ R+, with λ > λ1 > 0 and α1 > α > 0. Thenwe have r1(t) = α1λ1/(eλ1t − 1) and r2(t) = r3(t) = αλ/(eλt − 1), t ∈ (0,∞). It can be verified that r3(t) ≤ r1(t) for allt ∈ (0,∞), and r1(t)/r3(t) is strictly increasing in t ∈ (0,∞). Moreover, X2 and X3 are identically distributed implies thatrandom variables max{X1, X2} andmax{X1, X3} are also identically distributed. Thus X3 ≤st X1, X2 ≤hr X3(r2(t) = r3(t),∀t ∈

R+), and r1(t)/r3(t) is not decreasing in t ∈ (0,∞), however, the conclusion of Corollary 2.1 holds.

The following lemma, whichmay be of independent interest, will be useful in proving the next theorem. The result givenin this lemma strengthens Lemma 7.1(a) of Barlow and Proschan (1975, p. 120).

Lemma 2.1. Let A = {(x, s) : 0 ≤ s ≤ x < ∞} and let K : A → R be a function such that

t |K(x, s)| dx < ∞, and∞

t K(x, s) dx ≥ 0 whenever 0 ≤ s ≤ t < ∞. Then, for any nonnegative and increasing function h : R+ → R+,∫∞

sh(x)K(x, s) dx ≥ 0, for all s ∈ R+.

Proof. Fix s ∈ R+. For a set B ⊆ R, let IB(·) denote its indicator function. Since h is a nonnegative and increasing functionthere exists a sequence {ψn}n≥1 of simple functions such that limn→∞ ψn(x) = h(x), ∀x ∈ R+, and

ψn(x) =

n−i=1

aiI[ti,∞)(x), x ∈ R+, n = 1, 2, . . . ,

for some sequences {an}n≥1 and {tn}n≥1 of positive real constants. Clearly, for each x ∈ R+, ψn(x) ≤ ψn+1(x), n = 1, 2, . . . .Let K+(x, s) = max{K(x, s), 0} and K−(x, s) = max{−K(x, s), 0}, x ≥ s, so that K+(x, s) ≥ 0, K−(x, s) ≥ 0 andK(x, s) = K+(x, s)− K−(x, s), x ≥ s. Then∫

sψn(x)K(x, s) dx =

n−i=1

ai

∫∞

max(ti,s)K(x, s) dx ≥ 0, ∀n ≥ 1,

and therefore

limn→∞

∫∞

sψn(x)K(x, s) dx ≥ 0. (2.2)

On using the monotone convergence theorem we obtain

limn→∞

∫∞

sψn(x)K+(x, s) dx =

∫∞

sh(x)K+(x, s) dx,

and

limn→∞

∫∞

sψn(x)K−(x, s) dx =

∫∞

sh(x)K−(x, s) dx.

Therefore,

limn→∞

∫∞

sψn(x)K(x, s) dx = lim

n→∞

∫∞

sψn(x)K+(x, s) dx − lim

n→∞

∫∞

sψn(x)K−(x, s) dx

=

∫∞

sh(x)K+(x, s) dx −

∫∞

sh(x)K−(x, s) dx

=

∫∞

sh(x)K(x, s) dx.

Now the result follows on using (2.2). �

The following theorem generalizes the result (1.2) proved by Da et al. (2010).

Theorem 2.2. Suppose that X1 ≤mrl X3 and X2 ≤mrl X3. Then,

max{X1, X2} ≤mrl max{X1, X3}.

Page 5: New results on stochastic comparisons of two-component series and parallel systems

N. Misra, A.K. Misra / Statistics and Probability Letters 82 (2012) 283–290 287

Proof. For t ∈ R+, define

h2(t) =

t P (max{X1, X3} > x) dx∞

t P (max{X1, X2} > x) dx

=

t [F1(x)+ F3(x)− F1(x)F3(x)] dx∞

t [F1(x)+ F2(x)− F1(x)F2(x)] dx

=

N2(t)D2(t)

, say.

We need to show that h2(t) is increasing in t ∈ R+. One can easily verify that

D22(t)

ddt

h2(t) = F1(t)[F2(t)− F3(t)]∫

tF1(x) dx + [F1(t)+ F1(t)F2(t)]

×

∫∞

tF1(x)F3(x) dx − [F1(t)+ F1(t)F3(t)]

∫∞

tF1(x)F2(x) dx, t ∈ R+. (2.3)

Fix t ∈ R+. Then the following two cases arise.Case I: F2(t) ≥ F3(t).We have

D22(t)

ddt

h2(t) ≥ [F1(t)+ F1(t)F3(t)]∫

tF1(x)[F3(x)− F2(x)] dx.

Since X2 ≤mrl X3 implies that X2 ≤icx X3 (see Shaked and Shanthikumar, 2007, p. 195), it follows that

y [F3(x)−F2(x)] dx ≥ 0for all y ∈ R+. Now, on using Lemma 2.1, we get

y F1(x)[F3(x)− F2(x)] dx ≥ 0, for all y ∈ R+. In particular, for y = t , wehave

t F1(x)[F3(x)− F2(x)] dx ≥ 0. Hence ddt h2(t) ≥ 0.

Case II: F2(t) ≤ F3(t).Since X2 ≤mrl X3 implies F2(s)

u F3(x) dx ≥ F3(s)

u F2(x) dx for all s ≤ u (see Shaked and Shanthikumar, 2007, p. 82), onusing Lemma 2.1, we obtain F2(s)

s F1(x)F3(x) ≥ F3(s)

s F1(x)F2(x) for all s ∈ R+. In particular, for s = t , we have

F2(t)∫

tF1(x)F3(x) ≥ F3(t)

∫∞

tF1(x)F2(x).

Using this in (2.3), we obtain

F3(t)D22(t)

ddt

h2(t) ≥ F1(t)F3(t)[F2(t)− F3(t)]∫

tF1(x) dx + F1(t)[F3(t)− F2(t)]

∫∞

tF1(x)F3(x) dx

≥ F1(t)F3(t)[F2(t)− F3(t)]∫

tF1(x) dx + F1(t)[F3(t)− F2(t)]

∫∞

tF1(t)F3(x) dx

= F1(t)[F3(t)− F2(t)]∫

t[F1(t)F3(x)− F3(t)F1(x)] dx

≥ 0 (since X1 ≤mrl X3).

Thus h2(t) is increasing in t ∈ R+. �

The following corollary is immediate from Theorem 2.2.

Corollary 2.2. Suppose that any one of the following conditions holds:

(a) X1 ≤mrl X2 ≤mrl X3;(b) X2 ≤mrl X1 ≤mrl X3.

Then

max{X1, X2} ≤mrl max{X1, X3}.

Recall that the result in Corollary 2.2(a) was proved by Da et al. (2010). The following example gives a practical case inwhich the assumptions of Corollary 2.2 are not satisfied but the result of Theorem 2.2 is valid.

Example 2.5. Let m1(t) =5t+24t+3 ,m2(t) =

6t+56(t+1) , and m3(t) =

7t+55t+4 , t ∈ R+. Then, it is easy to verify that m1(t) ≤ m3(t)

and m2(t) ≤ m3(t) for all t ∈ R+. Also, m2(0)− m1(0) = 1/6 > 0 and m2(1)− m1(1) = −1/12 < 0. Thus, X1 ≤mrl X3 andX2 ≤mrl X3, but there is no mean residual life ordering between the random variables X1 and X2.

Page 6: New results on stochastic comparisons of two-component series and parallel systems

288 N. Misra, A.K. Misra / Statistics and Probability Letters 82 (2012) 283–290

Now,we derive sufficient conditions underwhichmax{X1, X2} ≤lr max{X1, X3} holds, thereby generalizes the result (1.3),given in Shaked and Shanthikumar (2007, Example 1.C.36, p. 56).

Theorem 2.3. Suppose that X2 ≤lr X3 and F 21 (t)f3(t)/f1(t) (or F

21 (t)f2(t)/f1(t)) is increasing in t ∈ R+. Then,

max{X1, X2} ≤lr max{X1, X3}.

Proof. We will prove the result for the non-parenthetical part, the other part can be proved similarly. Let π2(t) and π3(t)denote the probability density functions of max{X1, X2} and max{X1, X3}, respectively. Define

h3(t) =π3(t)π2(t)

=f1(t)F3(t)+ f3(t)F1(t)f1(t)F2(t)+ f2(t)F1(t)

, t ∈ R+.

It is easy to verify that for t ∈ R+

π22 (t)

ddt

h3(t) = [2f 21 (t)− f ′

1(t)F1(t)][f3(t)F2(t)− f2(t)F3(t)] + f ′

3(t)F1(t)

× [f1(t)F2(t)+ f2(t)F1(t)] − f ′

2(t)F1(t)[f1(t)F3(t)+ f3(t)F1(t)],

where f ′

i (·) denotes the derivative of fi(·), i = 1, 2, 3. Since X2 ≤lr X3 implies f2(t)f ′

3(t)− f3(t)f ′

2(t) ≥ 0,∀t ∈ R+, it followsthat

π22 (t)

ddt

h3(t) ≥ f1(t)F1(t)[f3(t)F2(t)− f2(t)F3(t)][2f1(t)F1(t)

−f ′

1(t)f1(t)

+f ′

3(t)f3(t)

]= f1(t)F1(t)F2(t)F3(t)[r3(t)− r2(t)]

ddt

lnF 21 (t)f3(t)/f1(t)

≥ 0,

where the last inequality follows from the assumption that F 21 (t)f3(t)/f1(t) is increasing in t ∈ R+ and the fact that X2 ≤lr X3

implies X2 ≤rh X3. Hence h3(t) is increasing in t ∈ R+, i.e. max{X1, X2} ≤lr max{X1, X3}. �

The following corollary is a simple consequence of the above theorem.

Corollary 2.3. Suppose that one (or more) of the following conditions holds:

(a) X1 ≤lr X3 and X2 ≤lr X3;(b) X1 ≤lr X2 ≤lr X3;(c) X2 ≤lr X1 ≤lr X3;(d) X2 ≤lr X3 and f3(t)/r1(t) (or f2(t)/r1(t)) is increasing in t ∈ R+;(e) X2 ≤lr X3 and r3(t)/r1(t) (or r2(t)/r1(t)) is increasing in t ∈ R+.

Then,

max{X1, X2} ≤lr max{X1, X3}.

3. Stochastic comparisons of series systems

In this section, we deal with the stochastic comparisons of series systems having lifetimes min{X1, X2} and min{X1, X3}.In the following theorems, we will compare random variables min{X1, X2} and min{X1, X3} with respect to the reversedhazard rate order and the likelihood ratio order.

Theorem 3.1. Suppose that one (or both) of the following conditions holds:

(a) X2 ≤rh X1 and X2 ≤rh X3;(b) X2 ≤rh X3 and r1(t)F2(t) ≥ r2(t)F1(t) for all t ∈ (0,∞).

Then

min{X1, X2} ≤rh min{X1, X3}.

Proof. For t ∈ (0,∞), consider

h4(t) =P (min{X1, X3} ≤ t)P (min{X1, X2} ≤ t)

=F1(t)+ F3(t)− F1(t)F3(t)F1(t)+ F2(t)− F1(t)F2(t)

=N3(t)D3(t)

, say.

Page 7: New results on stochastic comparisons of two-component series and parallel systems

N. Misra, A.K. Misra / Statistics and Probability Letters 82 (2012) 283–290 289

It can be easily verified that

D23(t)

ddt

h4(t) = r1(t)F1(t)[F2(t)− F3(t)] − r2(t)F1(t)F2(t)[F1(t)+ F3(t)− F1(t)F3(t)]

+ r3(t)F1(t)F3(t)[F1(t)+ F2(t)− F1(t)F2(t)].

(a) Since r3(t) ≥ r2(t) for all t ∈ (0,∞), it follows that

D23(t)

ddt

h4(t) ≥ F1(t)[F2(t)− F3(t)][r1(t)− r2(t)F1(t)] (3.1)

≥ 0, ∀t ∈ (0,∞),

where the last inequality follows from the assumptions that r1(t) ≥ r2(t),∀t ∈ (0,∞), and the fact that X2 ≤rh X3implies F2(t) ≥ F3(t),∀t ∈ R+. It follows that h4(t) is increasing in t ∈ R+, i.e. min{X1, X2} ≤rh min{X1, X3}.

(b) From (3.1), we can write

D23(t)

ddt

h4(t) ≥ F1(t)[F2(t)− F3(t)][r1(t)F2(t)+ r1(t)F2(t)− r2(t)F1(t)]

≥ 0, ∀t ∈ (0,∞),

where the last inequality follows from the assumptions that r1(t)F2(t) ≥ r2(t)F1(t),∀t ∈ (0,∞), and the fact thatX2 ≤rh X3 implies F2(t) ≥ F3(t),∀t ∈ R+. Thus, h4(t) is increasing in t ∈ (0,∞). Hence the result follows. �

Remark 3.1. From (3.1) it is clear that the conclusion of Theorem 3.1 remains true if X2 ≤rh X3 and r1(t) ≥ r2(t)F1(t) for allt ∈ (0,∞).

As a consequence of Theorem 3.1, we have the following corollary. The proof of the following corollary, being similar tothat of Corollary 2.1, is omitted.

Corollary 3.1. Suppose that one (or more) of the following conditions holds:

(a) X2 ≤rh X1 ≤rh X3;(b) X2 ≤rh X3 ≤rh X1;(c) X2 ≤st X1, X2 ≤rh X3, and r1(t)/r2(t) is decreasing in t ∈ (0,∞);(d) X1 ≤st X2, X2 ≤rh X3, and r2(t)/r1(t) is decreasing in t ∈ R+.

Then

min{X1, X2} ≤rh min{X1, X3}.

The following examples illustrate a few situations on comparative usefulness of results given in Theorem 3.1 andCorollary 3.1.

Example 3.1. Let F1(t) = (1 − e−5t)4, F2(t) = (1 − e−35t)2, F3(t) = (1 − e−30t)5, t ∈ R+, and let φ2(t) = F1(t)/F3(t), t ∈

(0,∞). It can be easily verified thatX2 ≤rh X1 andX2 ≤rh X3. Also,wehaveφ2(0.01) = 0.00483737, φ2(0.02) = 0.00438607,and φ2(0.03) = 0.00511513, which implies that there is no reversed hazard rate ordering between the random variablesX1 and X3. For t ∈ (0,∞), consider

φ3(t) =r1(t)r2(t)

=2e30t(1 − e−35t)

7(1 − e−5t),

and

φ4(t) = r1(t)F2(t)− r2(t)F1(t)

=20e−5t

[1 − (1 − e−35t)2]

1 − e−5t−

70e−35t[1 − (1 − e−5t)4]

1 − e−35t.

Then, it can be easily verified that

[7(1 − e−5t)]2ddtφ3(t) = 7e−5t(60e35t − 70e30t + 10) ≥ 0, ∀t ∈ (0,∞),

i.e. φ3(t) is increasing in t ∈ (0,∞). Moreover, we have φ4(0.08) = 0.316523 and φ4(0.10) = −0.293543. Thus, theconditions of Theorem 3.1(a) are satisfied but the conditions of Theorem 3.1(b) and Corollary 3.1(a)–(d) are not satisfied.

Page 8: New results on stochastic comparisons of two-component series and parallel systems

290 N. Misra, A.K. Misra / Statistics and Probability Letters 82 (2012) 283–290

Example 3.2. Let Fi(t) = e−λit , t ∈ R+, i = 1, 2, 3, with λ1 > λ2 > λ3 > 0. Clearly, X1 ≤st X2, X2 ≤rh X3, and r2(t)/r1(t)is decreasing in t ∈ R+. Thus, the conditions of Corollary 3.1(d) are satisfied but the conditions of Theorem 3.1(a) andCorollary 3.1(a)–(c) are not satisfied. Since the conditions that r2(t)/r1(t) is decreasing in t ∈ R+ and X1 ≤st X2 implies thatr1(t)F2(t) ≥ r2(t)F1(t),∀t ∈ R+, which is equivalent to the condition r1(t)F2(t) ≥ r2(t)F1(t),∀t ∈ (0,∞), it follows thatthe result of Theorem 3.1(b) is also valid.

The following theorem generalizes the result (1.3), given in Shaked and Shanthikumar (2007, Example 1.C.36, p. 56). Theproof of the following theorem, being similar to that of Theorem 2.3, is omitted.

Theorem 3.2. Suppose that X2 ≤lr X3 and f1(t)/(f2(t)F 21 (t)) (or f1(t)/(f3(t)F

21 (t))) is increasing in t ∈ R+. Then,

min{X1, X2} ≤lr min{X1, X3}.

The following corollary immediately follows from above theorem.

Corollary 3.2. Suppose that one (or more) of the following conditions holds:

(a) X2 ≤lr X3 and X2 ≤lr X1;(b) X2 ≤lr X1 ≤lr X3;(c) X2 ≤lr X3 ≤lr X1;(d) X2 ≤lr X3 and r1(t)/f2(t) (or r1(t)/f3(t)) is increasing in t ∈ R+;(e) X2 ≤lr X3 and r1(t)/r2(t) (or r1(t)/r3(t)) is increasing in t ∈ R+.

Then,

min{X1, X2} ≤lr min{X1, X3}.

Acknowledgments

The authors are grateful to the anonymous referee for his valuable comments and suggestions that have led to animproved version of the article. The second author would like to acknowledge the financial assistance from the Councilof Scientific and Industrial Research, India, for carrying out this research work.

References

Barlow, R.E., Proschan, F., 1975. Statistical Theory of Reliability and Life Testing: Probability Models. Holt, Rinehart and Winston, New York.Boland, P.J., El-Neweihi, E., Proschan, F., 1994. Applications of hazard rate ordering in reliability and order statistics. Journal of Applied Probability 31,

180–192.Da, G., Ding, W., Li, X., 2010. On hazard rate ordering of parallel systems with two independent components. Journal of Statistical Planning and Inference

140, 2148–2154.Dykstra, R., Kochar, S.C., Rojo, J., 1997. Stochastic comparisons of parallel systems of heterogeneous exponential components. Journal of Statistical Planning

and Inference 65, 203–211.Joo, S., Mi, J., 2010. Some properties of hazard rate functions of systems with two components. Journal of Statistical Planning and Inference 140, 444–453.Khaledi, B.-E., Kochar, S.C., 2000. Some new results on stochastic comparisons of parallel systems. Journal of Applied Probability 37, 1123–1128.Misra, N., Misra, A.K., Dhariyal, I.D., 2011. Active redundancy allocations in series systems. Probability in the Engineering and Informational Sciences 25,

219–235.Müller, A., Stoyan, D., 2002. Comparison Methods for Stochastic Models and Risks. John Wiley & Sons, New York.Shaked, M., Shanthikumar, J.G., 2007. Stochastic Orders. Springer, New York.Zhao, P., Balakrishnan, N., 2011. New results on comparisons of parallel systems with heterogeneous gamma components. Statistics and Probability Letters

81, 36–44.