Bounds on the Lifetime Expectations of Series Systems with IFR Component Lifetimes

We consider series systems built of components which have independent identically distributed (iid) lifetimes with an increasing failure rate (IFR). We determine sharp upper bounds for the expectations of the system lifetimes expressed in terms of the mean, and various scale units based on absolute central moments of component lifetimes. We further establish analogous bounds under a more stringent assumption that the component lifetimes have an increasing density (ID) function. We also indicate the relationship between the IFR property of the components and the generalized cumulative residual entropy of the series system lifetime.


Introduction
The series system composed o n items has the structure function ϕ : {0, 1} n → {0, 1} defined as ϕ(x 1 , . . . , x n ) = min{x 1 , . . . , x n }. If X 1 , . . . , X n denote the random lifetimes of the components, then the system lifetime is X 1:n = min{X 1 , . . . , X n }, where X 1:n , . . . , X n:n denote the order statistics of X 1 , . . . , X n . The series systems are the most important objects of the reliability theory. They are composed of the minimal number of components necessary for system operating. The failure of any component results immediately in the failure of the whole system. Series structures occur often as submodules of complex systems and networks. Applications of the series systems go far beyond the reliability of technical devices. They are useful, e.g., in economy (see, Ge et al. [1]), chemistry (see, Showalter and Epstein [2]) and industry (see, Xie and Gong [3]). Many survival analyses were based on reliability models of competing risks such that the first failure from a list of risks (e.g., cancer, heart disease, accident, etc.) causes the system failure (see, Kvam and Singh [4]). Furthermore, series systems are utilized by genetic algorithms and evolutionary theory (see, e.g., Jäntschi and Bolboacȃ [5]). In the paper, we consider the series systems whose components have independent and identically distributed (iid, for short) lifetimes with a common increasing failure rate (IFR, for short).
An absolutely continuous nonnegative life distribution function F with the density function f is called the increasing failure rate distribution function if the respective failure rate function λ(x) = f (x) 1−F(x) , x > 0, is nondecreasing. The name is inherited by random variables which have IFR distribution functions. The increasing failure rate is the most classic notion used for describing a natural increasing fatigue tendency of appliances and their components during operation. The wearing tendency of the components (and the system, in consequence) is the effect of increasing fatigue of working components in the operating time, and increased burden caused by possible failures of other components. Another useful definition of the IFR property is based on the convex transform relation of the distribution with the standard exponential distribution Z(x) = 1 − exp(−x), x > 0. The order was introduced by van Zwet [6]. We say that F is the IFR distribution if it precedes the exponential distribution in the convex transform order, i.e., the composition F −1 • Z(x) is concave on R + . This definition of IFR distributions generalizes the the traditional one in two directions. Firstly, by the location-scale invariance of the convex transform order, the support of F does not need to be restricted to the positive half-axis. Secondly, the convex transform order definition admits distribution functions which have atoms at the right-ends of their supports. The IFR properties of the series systems, and their dual properties of decreasing failure rate (DFR, for short) were studied by many researchers, e.g., by Barlow and Proschan [7], Lai and Xie [8], and Hazra et al. [9].
A more stringent concept of aging is expressed by the increasing density property (we use further the acronym ID for brevity). It is obvious that if the density function f (x) is nondecreasing on the support of F, then so is the product λ(x) = f (x) 1 1−F(x) . The ID property may be also defined with use of the convex transform relation. A distribution function F has the ID property if it precedes the standard uniform distribution function U(x) = x, 0 < x < 1, in the convex transform order, which means that is concave on (0, 1). This defines the family of distributions with nondecreasing density functions on some interval, and possible atom at its right end. This definition is also more restrictive than the IFR property because Z precedes U in the convex transform order, and the order is transitive.
Recently, various versions of entropy were applied for measuring uncertainty of random life data, and order statistics and reliability system lifetimes in particular. Recall that for lifetime random variable X supported on [0, ∞) the classic entropy of Shannon [10] is defined as Rao et al. [11] defined and studied the cumulative residual entropy and the generalized cumulative residual entropy was proposed in Psarrakos and Navarro ( [12]. For instance, Baratpour [13] provided conditions for characterization of the component lifetime distribution by the cumulative residual entropy of the series system lifetime. Sunoj et al. [14] analyzed the quantile-based entropy of order statistics. Toomaj and Di Crescenzo [15] presented the relationships between the generalized cumulative residual entropy and monotone failure rate properties. They proved, e.g., that if components of the series system are iid and belong to the IFR family, then the generalized cumulative residual entropy of the conditional residual system lifetime E n (X 1:n − t|X 1:n > t), t > 0, is decreasing. Other connections between Shannon entropy (see Wang et al. [16], Jäntschi and Bolboacȃ [17], Tȃnȃsescu and Popescu [18]) and order statistics were presented, e.g., by Koné [19], Jäntschi [20], Courtier et al. [21] and Cofré et al. [22]. We assume here that the common marginal component lifetime distribution function F belongs either to the IFR or ID family and determine sharp upper bounds on the standardized expectations EX 1:n −µ σ p of series system lifetime, where µ = EX 1 is the expectation of the parent distribution function F, and σ p = (E|X − µ| p ) 1 p , p ≥ 1, is the pth root of the respective pth absolute central moment. Writing σ p , we tacitly assume that this is positive and finite, i.e., F is nondegenerate and it has a finite pth moment. Observe that our bounds are necessarily nonpositive, because EX 1:n < µ for every iid sample with a nondegenerate baseline distribution function F. In this paper, we show that for every parent distribution with increasing density and failure rate functions, and every p > 1 the bounds for standardized expectations of series system lifetimes trivially amount to 0. Far-ther, we determine strictly negative bounds expressed the scale units σ 1 , being the absolute mean deviation from the population mean. Obviously, the bounds for the IFR family occur greater than those in the ID case, because the IFR family is substantially greater than the class of increasing density distributions. We also describe the parent distributions which attain the bounds. They are elements of the families preceding the exponential and uniform distributions in the convex transform order that, in some cases, possibly have atoms at the right end of their supports. However, these extreme distributions can be approximated by sequences of absolutely continuous members of the IFR and ID classes, having even strictly increasing failure rate and density functions, respectively. Hence the bounds are attained in the limit by these sequences.
Obviously, the series systems are the particular cases (with k = n) of the k-out-of-n systems, k = 1, . . . , n, that work as long as do so at least k of its n components. The meanstandard deviation bounds on the expectations of k-out-of n systems with k = 1, n − 1, n (called the parallel, fail-safe, and series systems, respectively) when the component lifetimes are iid and belong to the ID and IFR families, were determined by Rychlik [23]. However, in the case of series systems only the trivial zero bounds were established there. Analogous results for the other k-out-of-n systems were presented by Goroncy and Rychlik in [24] for the ID case, and in [25] for the IFR family. The evaluations for the dual DD and DFR families were determined several years earlier. In particular, Danielak [26] presented positive mean-standard deviations bounds on the expectations of order statistics X k:n with ranks k relatively large with respect sample sizes n coming from populations with decreasing density and failure rate functions. These correspond to the k-out-of-n systems with k relatively far from n. Rychlik [27,28] established negative sharp bounds on the differences E(X k:n − X 1 ) for small rank order statistics from the DD and DFR families, respectively, expressed in various scale units. As we can conclude from the above short review, there are known the sharp bounds on the expected lifetimes of all k-our-of-n systems under various assumptions on monotonicity of the failure rate and density function of lifetimes of their identical components except for the series systems under the IFR and ID assumptions. In the paper, we fill this important gap.
The rest of our paper is organized as follows. In Section 2 we show that the obvious zero upper bounds for the differences between the expectations of series system and its single component lifetimes cannot be improved in the ID and IFR cases if the differences are gauged in the scale units based on central absolute moments of orders greater than 1. In Section 3 we obtain refined nontrivial strictly negative evaluations of the differences expressed in the mean absolute deviation units. It occurs that the bounds under the more restrictive assumption that the distributions of component lifetimes have a common increasing density function are substantially tighter than in the case that we only assume that the distributions have merely an increasing failure rate. Section 4 contains a brief summary of the results established in the paper.

Zero Bounds
In contrast with the DD and DFR populations, the nonpositive upper bounds on the differences between the expectations of order statistics and population means in the ID and IFR cases are possible only for the sample minima which in reliability problems represent the lifetimes of series systems. Therefore it is of a vital interest to evaluate the expectations of the system lifetimes, especially under the practically important assumption that the component lifetime distributions have the IFR property.
Theorem 1 (ID and IFR distributions). Let X 1 , . . . , X n be iid random variables with a convex baseline distribution function and finite pth absolute central moment σ p p for some p > 1. Then the following bound is sharp, and the equality is attained there in the limit by the sequences of mixtures of uniform distributions with atoms at the right ends of their support, when the contributions of the absolute continuous parts vanish.
Obviously, the mixtures described in the theorem have convex distribution functions. Since each convex distribution function has a convex cumulative failure rate, and X 1:n ≤ X 1 for all parent distributions, then the bound (1) is valid and sharp for the family of IFR distributions as well. The statement of Theorem 1 was determined in Rychlik [23] in the special case p = 2.
Proof. Consider first the family of mixtures of uniform distributions on [0, 1], and atoms at 1 with probabilities 0 < α < 1 and 1 − α, respectively. They have distribution functions which are convex on their common support [0, 1], and concave quantile functions We easily check that for iid Y 1 , . . . , Y n with the above parent distributions we have The respective pth absolute central moments are For p > 1 we have Combining (5) with (7) we obtain Obviously, replacing our Y i by . . , n, which have arbitrarily chosen means µ ∈ R, pth absolute central moments σ p p > 0, and convex distribution functions, we attain the zero bound on when α → 0 as well.

Negative Bounds
The results of the previous section show that we get noninformative zero bounds on the differences between the expectations of the system and component lifetimes if we measure the differences in non-sufficiently subtle scale units, including the standard deviation of the component lifetime in particular. The conclusions essentially change if we replace the scale units σ p , p > 1, by the mean absolute deviation from the mean σ 1 = E|X 1 − µ|. We first consider the series systems whose components have IFR lifetimes.
Theorem 2 (IFR distributions). For the minimum of n iid random variables with a common distribution that has the convex cumulative failure rate and finite mean µ and first absolute central moment σ 1 , the following inequality Proof. The idea of proof is based on the norm maximization method proposed by Goroncy and Rychlik [29] (see also Rychlik [27,28], and Goroncy and Rychlik [25]). Let f 1:n (x) = n(1 − x) n−1 , 0 < x < 1, denote the density function of the first order statistic of the standard uniform sample of size n. Let Z(x) = 1 − e −x , x > 0, be the standard exponential distribution function. We first represent as a linear functional on the elements of the Banach space L 1 (R + , e −x dx). Here is an arbitrary element of the family where Indeed, each g ∈ S 1 is nondecreasing, because this is an increasing linear transformation of a quantile function composed with the exponential distribution function, and concave, by the IFR property of F. Vanishing of (13) for g from (12) follows the representation (11) and the fact that Similarly, (11) and because EX 1:n < EX 1 for all nondegenerate iid random variables with a finite first moment (equivalently, with positive finite σ 1 ). This allows us to use the norm maximization method.
Noting that the problem of maximizing (10) is location and scale invariant, we can replace the original problem by the dual one of maximizing ||g|| 1 over the set Indeed, replacing g ∈ S 1 byg = g −T h (g) , we get T h (g) = −1 and ||g|| 1 = 1 |T h (g)| = − 1 T h (g) . It means that maximizing the norm over (14) we also maximize T h (g) over the original set (12).
The next step of our reasoning is based on the observations that (14) is a convex set, and the norm functional is convex. It follows that for every convex combination of elements of (14) we have and in consequence, as k → ∞, Now we notice that every g ∈ S 2 can be approximated by convex combinations of the elements of a particular subset of S 2 consisting of simple broken lines The slopes b(α) > 0 and intercepts a(α) are uniquely determined by the equations T 1 (g α ) = 0 and T h (g α ) = −1, but we do not need here their precise forms. The approximation procedure is standard, we restrict ourselves to outlining the main steps, and do not perform all the rigorous calculations. We first observe that each piecewise linear element of (14) is a convex combinations of some broken lines (15). The breaking points uniquely determine the elements of the combination, and the values of consecutive slopes allow us to calculate the respective mixture weights. Then we observe that each g ∈ S 2 can be linearly interpolated by sequences of piecewise linear, continuous, nondecreasing, and concave functions g k , k → ∞, with the increasing sets of knots which in the limit become dense in R + . These approximations g k do not necessarily belong to S 2 , but they tend monotonously to g on the whole positive half-axis, which implies that as k → ∞. However, the following modifications do belong to (14), and they satisfy because ||g k || 1 ≤ ||g|| 1 < ∞, the last norm amounts to 1, and the remaining three terms in the last line vanish as k → ∞. It follows that ||g k || 1 → ||g|| 1 by continuity of the norm functional. Therefore and finally ||g|| 1 ≤ sup 0<α<∞ ||g α || 1 , g ∈ S 2 , which implies that in our dual maximization problem (and the original one as well) we can restrict ourselves to considering first increasing and then constant broken lines (15) with varying breaking points.
Referring again to location-scale invariance of the problem, we consider the simplest representatives of the location-scale families. Precisely, we confine our attention to the parent distribution functions F α , α > 0, satisfying (2)). The explicit forms of the distribution functions are which means that they have absolutely continuous parts with the density functions f α (x) = α exp(−αx) on [0, 1], and jumps of heights e −α at 1. If Y 1 , . . . , Y n are iid with the distribution function F α , then We also determine (see (9)). Accordingly, say. Our goal is to maximize each B n (α), n ≥ 2, with respect to 0 < α < ∞. Fix first n = 2 and change the variable for x = 1 − e −α ∈ (0, 1) in (17). Then we obtain We develop the denominator into the power series The derivative has the following expansion where a 0 = − 1 6 < 0, a 1 = 1 3 > 0, and a i > 0, i = 2, 3, . . ., as well, because 2(−1) i (i+1)! < 1 for all i ≥ 2. This means that g 2 (0) < 0 and g 2 increases on (0, 1), where it is well defined. It follows that g 2 either always decreases on (0, 1) or it first decreases and then increases there. Anyway, the global maximum is attained at either of the end-points. Consequently, attains its global supremum either at 0 or at +∞, because functions (0, ∞) α → x = 1 − e −α ∈ (0, 1) and (0, 1) x → −1 2x ∈ −∞, − 1 2 are increasing. We easily check that Using the transformation x = 1 − e −α for n = 3 we get We verify variability of g 3 (x) = g 2 (x) 3−x analyzing the sign of Applying (18) and (19), by elementary calculations we obtain the Taylor expansion of are positive, and change the sign. We easily check that which implies that b i + c i > 0, i = 1, 2, . . ., and so the derivative of g 3 is positive. Accordingly, is increasing.
We apply the claim for establishing the same conclusion for n ≥ 4. Consider the functions with the derivatives M n (x) = n(x n−1 − 1), and analyze variability of Since we can rewrite the expression in the brackets of (20) as follows This shows that the ratios M n (x) , 0 < x < 1, are increasing for all n ≥ 3, and M n+1 (e −α ) M n (e −α ) , n ≥ 3, are decreasing with respect to α > 0. Moreover, We finally observe that is negative increasing, because so is B 3 (α), and each M i (e −α ) , i = 3, . . . , n − 1, is positive decreasing.
Summing up, for every n ≥ 2 we concluded that Distribution functions (16) which attain the bound in the limit represent the mixtures of the exponential distributions with intensity α truncated right to the interval [0, 1] with probability 1 − e −α , and atom at 1 with probability e −α . Making a standard location-scale transformation we obtain the distributions with arbitrarily chosen real µ and positive σ 1 which attain the bound in the limit. They are precisely described in the statement of Theorem 2.
Theorem 2 provides a reliable estimate of the smallest possible time distance between the lifetimes of the single component and series system under the increasing failure rate assumption. Below we also determine more stringent bounds on the expectations of series system lifetimes under the more restrictive aging condition that the component lifetimes have increasing density functions.

Theorem 3 (ID distributions).
Suppose that the assumptions of Theorem 1 hold with p = 1. If n = 2, then the bound is sharp , and it is attained in the limit by the mixtures of uniform distributions on the interval 2 with the atoms at µ + 2σ 1 (2−α) 2 when the probabilities of the atoms tend to 1. For n ≥ 3 we have If n = 3, then the equality is attained by the mixture distributions described above for any fixed α ∈ (0, 1]. For n ≥ 4, the equality in (22) is attained by the uniform distribution on [µ − 2σ 1 , µ + 2σ 1 ] only.
Bound (21) is identical with sharp bound for the lifetime of two-component system with arbitrary distributions of component lifetimes (see Goroncy [30], Corollaries 3.1 and 3.2). Surprisingly, bounds (22) coincide with the optimal evaluations for the expectation of the first order statistic based on an iid sample with symmetric unimodal parent distribution (see, Rychlik [31]).

Proof.
We again refer to the norm maximization method. We can mimic the arguments of the proof of Theorem 2, with the only difference that instead of considering the ele-ments of the space L 1 (R + , e −x dx) we take into account the functions from L 1 ([0, 1], dx). Consequently, we conclude that the possible candidates for quantile functions maximizing EX 1:n −µ σ 1 belong to the family of broken lines (2). The respective expectations of the single component and system lifetimes are presented in (3) and (4). Fixing p = 1 in (6) we obtain the respective mean absolute deviation from the mean of the form This together with (5) give say. The upper bound we look for amounts to the maximal value of the RHS of (23) with respect to 0 < α ≤ 1.
We immediately check that C 3 (α) = −1, 0 < α ≤ 1. This means that for n = 3, the upper bound is equal to −1, and it is attained by the arbitrary mixtures of standard uniform distributions with the atoms at 1. For other n we use the following result.
We consider the ratios L n (x) L n−1 (x) for n ≥ 3. Their derivatives satisfy The expression in the brackets can be rewritten as We concluded every fraction L n (x) L n−1 (x) , n ≥ 3, is increasing on (0, 1), and, in consequence, all L n (1−x) L n−1 (1−x) , n ≥ 3, are decreasing there. Note that Obviously N 1 (α) > 0. The same applies to all functions L n , because − L n (1 − α) 2(n + 1)α = EY 1:n − EY 1 , (comp. the first line of (5)), and the RHS is negative for all nondegenerate iid random variables. Also is decreasing, because the latter factor is positive increasing. Therefore (see (23)). This means that the bound for n = 2 is equal to − 1 2 , and it is attained in the limit by the convex combinations of the uniform distribution on [0, 1] and atom at 1 as the contribution of the continuous part decreases to 0.
For n ≥ 4 we have The function is increasing because each factor of the product is positive decreasing, and attains its maximum −2 n−1 n+1 at 1. This is the bound value and this is attained by the standard uniform distribution.
As we can predict, in the case of ID component lifetimes the negative bounds on the standardized expectations of series system lifetimes turn out to be more stringent than in the case of IFR components. Indeed, the absolute value of the right-hand side of (22) is 4 e n n+1 ≈ 1.47152 n n+1 times greater than its counterpart of (8). This shows that constructing systems of the component with either ID or merely IFR lifetime property significantly affects the system performance. Note that the bounds (8) and (22) as well as their ratio slightly depend on the size n of the series system.

Conclusions
In this paper we considered series systems whose components have independent and identically distributed lifetimes with either increasing failure rate and or increasing density function. The series systems are the simplest and most popular system used in technical applications, and the increase of failure rate or the density function of component lifetimes are natural analytic descriptions of the natural presumption that the components are wearing out during the system operation. It is obvious that the system lifetime is not longer than the lifetime of its any component. Here we studied the problem how much is the expectation of the system lifetime shorter than the expected lifetime of a single component under the assumption that the components have IFR and ID lifetimes. It occurs that the answer strongly depends on the units in which the difference is measured. In Section 2 we showed that the scale units based on the central absolute moments of component lifetimes with orders p > 1 are too rough for detecting the differences, and the respective bounds expressed in these units provide only the trivial zero bounds. For detecting substantial differences, in Section 3 we used the mean absolute deviations from the mean units, and we determined sharp attainable strictly negative bounds on the expected differences between the system and component lifetimes under the IFR and ID assumptions. Theorems 2 and 3 show that that bounds under the more restrictive ID assumption are essentially less than those the the IFR case.