Journal of Statistical North-Holland
Planning
and Inference
On a nonparametric and its dual Murari
385
39 (1994) 385-397
family of life distributions
Mitra
Indian Statistical
Institute,
Calcutta,
India
Sujit K. Basu Indian Institute Received
of Management
Calcutta,
India
22 June 1992/ Revised manuscript
received
25 January
1993
Abstract In this paper, we discuss a nonparametric family of life distributions (called the NWBUE family) which includes the IDMRL class ofdistributions introduced by Guess, Hollander and Proschan as well as all BFR distributions. We prove two inequalities which are then used to obtain bounds for moments of a NWBUE life distribution. The bounds thus obtained are shown to be related to the moments of an appropriate negative exponential distribution, and a characterization of the exponential distribution is derived as a consequence. Closure under weak convergence and the equivalence of weak convergence and moment convergence in the NWBUE family are established under mild conditions. Similar properties in respect of the dual family comprising of NBWUE distributions have also been explored. AMS Subject Classification:
Primary
62N05, 90B25; secondary
Key words: Different nonparametric classes convergence; convex and concave functions.
60F05.
of life distributions
based
on notions
of ageing;
weak
1. Introduction The concept of ageing is of vital importance in reliability theory. ‘No ageing’ means that the age of a component has no effect on the distribution of its residual lifetime. ‘Positive ageing’ describes the situation where the residual lifetime tends to decrease, in some probabilistic sense, with increasing age of the component. This phenomenon has been described in the literature in terms of differnt ageing notions which have given rise to a whole variety of nonparametric classes of life distributions, e.g. IFR, IFRA, DMRL, NBU, NBUE, etc. (see e.g. Barlow and Proschan, 1981). We have
Correspondence to: Dr. SK. Basu, Indian Institute P.O. Box No. 16757, Calcutta 700027, India.
037%3758/94/$07.00 0 1994-Elsevier SSDI 0378-3758(93)E0048-L
of Management,
Science B.V. All rights reserved.
Joka,
Diamond
Harbour
Road,
386
M.
a corresponding concept of negative (beneficial) like DFR, DFRA, IMRL, NWU and NWUE. In many
practical
situations,
however,
ageing and the resulting
the effect of age is initially
dual classes beneficial
(a
‘burn-in’ phase where negative ageing takes place), but after a certain period, the effect of age is adverse indicating a ‘wear-out’ phase where the ageing is positive. Situations displaying infant mortality or work hardening such ageing. This kind of nonmonotonic
of certain tools are common examples of ageing is typically modelled using life
distributions displaying bathtub failure rates (BFR). Guess et al. (1986) have introduced the class of IDMRL distributions to study the above phenomenon through the mean residual life (MRL) function. Let us introduce some notations and a few definitions to make these notions more explicit. Let F be a life distribution function (d.f.) and let F( .)= 1 -F( .) denote the corresponding survival function. Let the cumulative hazard function be denoted by R(t)+ -log,F(t). The mean residual life (MRL) function at age t20 is defined as
e&):=8(X-t,x>+ whenever
m
F(x) dx.
F(t) > 0.
Definition 1.1. A life distribution F is said to be a BFR (‘upturned bathtub failure rate’ (UBFR)) distribution if there exists a to 20 such that R(t) is concave (convex) on [0, to) and convex (concave) on [to, co). Definition 1.2. A life distribution F is said to be an ‘increasing initially, then decreasing mean residual life’ (IDMRL) (‘decreasing initially, then increasing mean residual life’ (DIMRL)) distribution if there exists a co 20 such that eF(t) is nondecreasing (nonincreasing) on [0, to) and nonincreasing (nondecreasing) on [to, co); in such a case, we shall say that F is IDMRL (to) (DIMRL(t,)). It is to be noted that the above definition of a BFR distribution is extremely general and extends the idea of distributions possessing a bathtub-shaped failure rate to situations
where the failure rate itself does not exist.
In this paper, we discuss a nonparametric class of life distributions which not only provides models for the kind of nonmonotonic ageing we are talking about, but also includes the IDMRL class of Guess et al. (1986) as well as the family of BFR distributions. To this end, we have the following definition. Definition 1.3. A life distribution F having support on [0, co) (and finite mean P) is said to be ‘new worse then better than used in expectation’ (NWBUE) (‘new better then worse than used in expectation’ (NBWUE)) if there exists a point x0 BO such that eF (-4
>(<)eF(0)
for
<(>)ep(0)
for x2x0.
x
(1.1)
M. Mitra. S.K. Basu 1 Nonparametric
387
.family of l(f>distributions
We shall refer to such an x,, (which need not be unique)
as a change point of the d.f. F;
we shall write F NWBUE(xO) (NBWUE(x,)) to indicate that the life distribution F is NWBUE (NBWUE) in the above sense and x0 is a change point of F. Let %‘Fbe the collection of all change points of an NWBUE (or NBWUE) life d.f. F; it is plain to see that for a continuous NWBUE (NBWUE)
d.f. F, %YFis either a singleton or a closed interval. Note that an F is NBUE (NWUE) if 0~57~ while it is NWUE (NBUE) if
COEWF. Example
1.4. Consider F(x)=
The corresponding
eF (x)
the life distribution
(2/V + 4)’ >
Odx
{ (4x/9)exp(l-x2)/6,
lbx
MRL function
=
whose survival
function
is given by
is given by
2+x,
Odx
i 31%
x> 1.
It is clear that F is NWBUE(xO)
with x0=3/2.
In the next section, we discuss the interrelationships amongst NWBUE, BFR and IDMRL classes of life distributions. In Section 3, we obtain some useful bounds for the moments of NWBUE distributions; we also provide a characterization of the exponential distribution as a consequence of some of the results discussed in this section. In Section 4, we prove closure under weak convergence and the equivalence of weak convergence and moment convergence in the NWBUE family of life distributions and also furnish a corollary which is an interesting by-product of our results. Our observations relating to the dual class comprising of NBWUE distributions are given in Section 5.
2. Relationships In this IDMRL
amongst NWBUE,
section,
we explore
BFR and IDMRL
the interrelationships
classes among
NWBUE,BFR
and
classes of life distributions.
Proposition 2.1. If a continuous then F is NWBUE.
and strictly increasing
Proof. Since F is continuous and strictly increasing, (TTT)-transform Ic/F(.) of F as
life d.jI F is BFR with mean p,
we can write the total time on test
388
M. Mitra, S.K. Basu 1 Nonparametric
where F(x)=t. O
family
that F is NWBUE
It then follows
of life disrributions
if and only if there exists a tb,
1, such that for t-c&,
dt *F(t) =
(2.1)
for tat;.
i 3t
Theorem 3.1 of Deshpande and Suresh (1990), to be referred to as DST below, informs us that F is BFR if and only if there exists a t0, 0 < t,, < 1 such that $F(t) is convex in [0, to) and concave in [to, 11. We now discuss the proof in 4 different cases as follows: Case I: qF(t) does not intersect the line g(t) = t, 0 6 t < 1. In view of DST, then F is either NBUE or NWUE according as $F(t) 3or d t, tE[O, l] and the conclusion of the theorem follows. Case II: First crossing ofg(t)=t by the curve GF(t) isfrom above. Let this point of intersection be t*; then, $F(t)= t, t~[t*, I], for, otherwise, either the DST or the fact that (cIF(l)= 1 will be violated. As such, F is NBUE and the theorem follows. Case III: First crossing ofg(t) = t by the curve GF(t) is from below and there are no further crossings. Here also tiF(t) = t, V tE[t *, 11, and the theorem follows. Case IV: First crossing of g(t) = t by the curve $F(t) is from below and there is a second crossing. In this case, in view of the DST the second crossing must be from 0 above, and the proof will follow as in Case II. The converse example.
of the above proposition
Example 2.2. Consider
is false as can be seen from the following
the life distribution
whose survival
I
e- ‘/x2,
x/*exp(
The corresponding
l
MRL function
!
+(x)=
r(x)= 1
Odx
X,
16x<2,
4/x,
2
co,
with x0 = 4. The failure rate function
1,
Odx
2/x,
16x<2,
(x/4-
reveals that F is not BFR.
26x
is given by
1,
SO that F is NWBUE(xO),
is given by
O
[ ePx,
F(x) =
function
l/x),
2
00,
which is given by
389
M. Mitra. S.K. Basu/ Nonparametric family of lifr distributions
Proposition
is IDMRL(t,),
2.3. IfF
then F is NWBUE(tb)
We omit the proof since it is trivial. Below we provide an example to convey without
with tb> to
that a life distribution
can be NWBUE
being IDMRL.
Example 2.4. Consider
F(x) =
the life distribution
I
{e/(e+cc)}*
{ej(e+cc))2 exp
i
--
having
the following
( e+2c(-p)z, ’f ( e+2c(+)O
function:
BGX
e+*u-*p+x
e+*a--2p+y
1
(x2-y2)
8+x,
o
e+2cr-x,
cr
e+2a-2p+x,
BbX
Y(e+2a-2p+y)jx,
x3y.
I
X3Y,
I ’
2 y(e+2c(-2p+y)
where the parameters CI,fi, y and f3 are such that 0
44 =
survival
B d y, 2cr > B and 8 > 0. The
It is easy to observe that the life distribution F is NWBUE with xo=y(8+2c(-2/3+y)/8; however, F is not IDMRL whenever cc?
for a life distribution
to be NWBUE,
its MRL function
has to lie above
the mean life ,u till x0 and below it then onwards. This property seems to imply a natural nonmonotonic ageing phenomenon which cannot be captured, as illustrated in Examples 2.2 and 2.4, through either BFR or IDMRL notions.
3. Inequalities
and moment bounds
We would need the following Lemma 3.1. If F is NWBUE(x,)
basic inequality. with (jnite)
mean p, then (3.1)
390
M. Mitra, S.K. Basu / Nonparametric family of life distributions
Further, if ep(x) is bounded above by M, then m
F(u)du d pe -xlP,-xo(M-‘-K’),
Proof. The probability fci,(x)=F(x)/p,
x30
vx>xO.
(3.2)
density function of the first derived distribution is given by and let r C1J(.) denote the corresponding failure rate function.
Using the fact that F is NWBUE(xe), x < or 3 x0. Writing
it follows that r(,,(x)
> l/p according
and using the above inequalities, (3.1) follows easily. If eF(x)d M, then (3.2) follows simply by noting that rC1,(x)= l/eF(x).
as
0
Proposition 3.2. Let F be NWBUE(x& x0 < co with (finite) mean ,u and 4( .) be a nondecreasing (nonincreasing) function on [0, co). Then
som
F(x)dx
4(x)
d jXy $(x)exp (~)
(F)dx.
(3.3)
Proof. We shall prove the result for nondecreasing 4; the other case can be treated similarly. Without loss of generality, we take 4 to be nonnegative. Let Z be a random variable whose d.f. is FC1,, defined in the previous proof and U be another random variable such that U-x0 is exponentially distributed with mean p. Then, as 4 is nondecreasing, in view of (3.1), we have
s,”
4 (4 6) dx
=P
j
0
(h9(4(Z)
>x)dx
P(Z>d-‘(x))dx
s
s m
b(m)
=P
9(Z>4-‘(x))dx+p
0
0(-w)
cc d
&(x0)
+
P(U
P s
@(+I)
>
4-‘(x))dx
P(Z>d-‘(x))dx
391
M. Mitra. S.K. Basu / Nonparametric family of life distributions
Now,
s
s
Q(.w)
m
9($(U)> Q(xo)
x)dx=b&U)-
9(4(u)>
x)dx
0
4(x0) =ac$(u)-
WW)>x)dx s0 Q(m)
ccP(U> 4-‘(x))dx,
=&c#I(U)s0
the integrand
being unity,
s
since 4 is nondecreasing.
Thus,
om 6(x)~(x)dxQ&(xo)+~C~@(U)--(xo)l=@’+(U),
which completes Corollary 3.3.
the proof.
0
IfF is NWBUE(xe),
x0 < co, with (finite) mean ,u, then x*-l exp(-(x-xo)/p)dx
for r>
1.
(3.4)
(')
Proof. Take &J(X)= x r-1 in Proposition and nondecreasing for r > 1. 0
3.1 and note that 4 is nonincreasing
The above result leads to the following distribution where x0 is finite. Corollary 3.4. Zf F is NWBUE(xo),
bounds
for moments
for r < 1
of a NWBUE(xo)
x0 < co with (finite) mean p, then
a, xI-1 dFX* < reXoip e-“‘“dx, (3) s x0
for r > 1. cc)
(3.5)
Given a NWBUE(xo) distribution, x0 < co, with finite mean p, we shall now obtain bounds for the moments of the above distribution in terms of the moments of an appropriate negative exponential distribution. Consider the following negative exponential distribution defined by the survival function
392
M. Mirra. S.K. Basu 1 Nonparametric family of life disrributions
Simple calculations
show that for r> 1, 30 x’-l
BGX’=xh+rexoiP
eexiP dx.
(3.6)
s x0 If r> 1 is an integer, ‘-I (xolP)j 1) 1 jl.
c?~Xr=X;+/L’r(r+
j=O
The following
corollary
(3.7)
.
now gives the various
Corollary 3.5. Let F be NWBUE(xe), dejined above. Then
moment
bounds.
x0 < cc with &ire)
mean p, and let G be as
(i) bcX*< cofor all r>O. (ii) &FX’<&GX’f~r all r> 1. (iii) f~Y~X’Qxb+p’r(r+ 1) CJ’&r (xo/p)j/j!for (iv) dFX’
all integers r3 1.
Proof. (i) This follows from (3.6). (ii) For r> 1, by (3.5), (3.8) co
x’-1 e-“‘“dx=d,X’.
s x0 (iii) This follows from the previous (iv) For rB 1, by (3.5),
part and (3.7).
30 X’-l eexip dx s x0 00 Gre XOIP x'-l e-X’pdx=p’r(r+ s0
8FX’
l)exO’P.
Remarks. (A) Taking x0 = 0 (i.e. if F is NBUE) in Corollary 3.3, we get the usual well-known bounds for the NBUE situation. (B) Although Corollary 3.5(ii) and (iv) give the same bounds for the NBUE case, it is clear that the former, which has an interesting probabilistic significance, is sharper than the latter, which, on the other hand, has useful applications as illustrated in the next section. The next distribution.
result
gives
an
interesting
characterization
of
the
exponential
393
M. Mitra, S.K. Basu / Nonparametric family of life distributions
Theorem
3.6. Let F be NWBUE(xe),
x0 < cc with finite mean p and let G be as in
Corollary 3.4. Then F is the exponential distribution if&FX’ = bcX’ Proof.
As &FX’= gGXr for some r > 1, then equality
(3.8) so that x,, = 0 and hence F is NBUE some r> 1. Consequently, F is HNBUE Basu and Bhattacharjee
(1984, Lemma
4. Weak convergence of NWBUE We start this section
for some r > 1.
holds in the string of inequalities
the relation 6, X’= F (r + 1) $ for the above relation and hence, by 0 2.4) F is exponential. satisfying satisfying
distributions
with the following
theorem.
Theorem 4.1. Let F,, n = 1,2, . . . be a sequence ofNWBUE(xO,) means p,,. Suppose that
life distributions with
(i) F,,-+F in law, where F is a continuous d$ and (ii) the sequences {p,,} and {x0”} are bounded. Then F is NWBUE. Further
s m
lim x*dF,(x)= rl+CC 0
s m
xrdF(x)
for every r>O.
(4.1)
0
Proof. Let ,u,,d B for all n 3 1 and let p be the mean of F. An application of Fatou’s lemma together with (i) above shows that p< co. We first prove that p,,-+p as a-03; for A >O, to be chosen suitably,
sa
{F,(x)-F(x)}dx=
0
Suppose
A (F,--)(x)dx+ s 0
cc F,(x)dxs A
xoll GM, for all n 2 1. For A > M > xOn, E> 0, by Lemma
“F(x)dx s A
3.1,
s m
12n=
F,(x)dx~~L,exp(-(A-xo,)/~~)~Be-‘A-M)’B
A
which can be made smaller than ~13 by choosing a sufficiently large A. Since p< 03, 1Z31 M, sufficiently large, we have 1Z2,,( -c&/3 for all nk 1 and II,1
(4.2)
394
M. Mitra, S.K. Basu / Nonparametric family of life distribulions
Since {x0,,} is bounded, there exists a subsequence {x0,,,} such that x0.,+/I< co as k-co. For x x, Vk B kO which implies 1
i”F,,(u)duBpL,,.
F,,(x). s Taking
as k-co,
limits
convergence
theorem 1
via condition
(i) of the theorem,
(4.2) and the dominated
we get
mF(u) du b p.
(4.3)
F(x), s Similarly
for x>fi, 1
we can show that coF(u) du
(4.4)
F(x), s
The proof of the first part of the theorem now follows from (4.3) and (4.4). To prove (4.1) we first note via (4.2) that it holds for r = 1. Let X,, y1= 1,2, . . . , and X be random variables having d.f.s F, and F, respectively. Then, as sup,, 6X, < co, {X.‘} is uniformly integrable for Y< 1, this together with (i) establishes (4. l), for r < 1. On the other hand, for r> 1, by Corollary 3.5 (iv), dX,‘<~~~(r+l)exo”‘@~. As {x0,,} is bounded and {pn} converges, { xOn/p,,} is bounded Vr> 1, by the above relation. So, (XL} is uniformly integrable proves (4.1) for r> 1. 0 Lemma 4.2. A d.jYwhich is NWBUE(x,J, sequence.
so that sup” 8X,* < cc for each r> 1, and this
x0 < cc, is uniquely determined
by its moment
Proof. Let F be NWBUE(xo), x0 < co with mean p. By Corollary 3.5(iv), it follows easily that the power series clY!o (u’/r!) bFX’ has a nonnull radius of convergence. The lemma now follows from a result of Loeve (1963, p. 217). 0 In Theorem 4.1, we showed that under a couple of conditions, weak convergence of a sequence of NWBUE d.f.s implies convergence of moments of the sequence of d.f.s to the corresponding moments of the limiting d.f.; the following deals with the converse of this result. d.js with xon < CC Theorem 4.3. Let F,,, n= 1,2, . . . , be a sequence of NWBUE(xon) Vn 2 1 and suppose that F is NWBUE(xo), x0 < co such that for all integers r 2 1, m m xrdF(x). lim x’dF,(x)= (4.5) n+30 s 0 s0 Then F,+F
in law.
M. Mitra, SK.
Proof. By Lemma
4.2, F is uniquely
this and (4.5) the limiting (F,,} happens
determined
distribution
to be necessarily
If condition
395
Basu / Nonparametric ,famiiy of life distributions
(ii) in Theorem
F and this completes 4.1 is replaced CC xdF(x)
Cc lim xdF,,(x)= n-tm s 0
by its moment
of every weakly
sequence.
convergent
the proof.
In view of
subsequence
of
0
by (
(4.6)
s0
then we will have an interesting a proposition.
conclusion
which we present
below in the form of
Proposition 4.4. Suppose in Theorem 4.1, we replace condition (ii) by (4.6). lf {x0,,}
is
unbounded, then F is N WUE. Proof. Let p be the mean of F; consider
there exists
(xonk} such that xOnr >x, Vk3 1. So,
a subsequence
Taking
any x > 0. As {x0”} is unbounded,
limits as k+ CC and arguing 1
F(x).
as in Theorem
4.1, we get
s m-
F(u) du >, pu,
so that F is NWUE.
0
Remark. It is trivial to observe from Proposition 4.4 that under (4.6), the weak limit of a sequence of NWUE d.f.s is also NWUE, provided the limiting distribution is continuous. The counterpart by Basu and Bhattacharjee
of this result concerning NBUE (1984) with more generality.
distributions
was given
5. The dual class In this section, family comprising
we extend similar investigations of NBWUE distributions.
as in Sections
3 and 4 to the dual
Remark 5.1. Through the same approach as in the proof of Proposition 2.1, it can easily be seen that if a continuous and strictly increasing life distribution F is UBFR (with finite mean p), then F is NBWUE. Also, as in Proposition 2.3, if F is DIMRL(t,), then F is NBWUE(tb) with t&2to.
396
M. Mitra, S.K. Basu 1 Nonparametric family
oflife distributions
Remark 5.2. The following example confirms that unlike NWBUE life distributions, an NBWUE distribution having a finite change point may not possess finite moments of order higher than one. Example
5.3. Consider
the life distribution
whose
F(x) = a(x 2 + c()- ‘, x > 0, where a > 0 is a constant. (9+/2
< co; however,
for I> 1, &rX’=
is given
by
The mean of this distribution
survival
function
is
co.
On the other hand, its failure rate function
is of the form
-1
r(x)=2
( > x+”
X
which is strictly increasing in (0, &) and strictly decreasing in (&, UBFR and hence, from Remark 5.1, we conclude that F is NBWUE. its MRL function
we even conclude
that F has a unique
finite change
co). As such, F is In fact, analysing
point,
In view of the above example, it is evident that a result analogous to Corollary 3S(ii) is, in general, false for NBWUE life distributions having finite change points. As such, issues contained in subsequent results in the same corollary do not seem relevant in the context of NBWUE distributions. Remark 5.4. From version of Theorem
the following 4.1 is false.
example,
Example 5.5. For n = 1,2, . . . , consider defined by the survival functions
we further
the sequence
conclude
of NBWUE(0)
that the NBWUE
life distributions
i.(x)=r,exp(-i)+&exp(-t). where O 1. Now, as n+ co, F,,+F in law, F being the d.f. of the exponential distribution having mean ,K But eF,Xm+ 00 > dFXm, though, both the sequences {,u”} and {xon} here are bounded. However, the question of closure of the NBWUE class under weak convergence still remains open.
M. Mitra. S.K. Basu 1 Nonparametric family of life distributions
Remark 5.6. In view of Remark
5.2, issues dealt with in Lemma
4.2 and Theorem
391
4.3
are, in general, not relevant in the context of NBWUE distributions. However, proceeding as in the proofs of Theorem 4.1 and Proposition 4.4, we can prove the following
version
of the latter.
Theorem 5.7. Let F,, n= 1,2, . . . , be a sequence ofNBWUE(xon) means pL,. Suppose
(i) F,-+F
Zf{xon}
in law, where F is a continuous
m m lim x dF, (x) = xdF(x) n-a, s 0 s0
(ii)
is bounded
life distributions with
that
th en F is NBWUE;
From the above theorem, we observe weak convergence provided (5.1) holds.
d.5;
(< co).
(5.1)
otherwise, F is NBUE.
easily that the NBUE
class is closed under
References Barlow, R.E. and F. Proschan (1981). Statistical Theory ofReliability and Life Testing. To Begin With, Silver Spring, MD. Basu, S.K. and M.C. Bhattacharjee (1984). On weak convergence within the HNBUE family of life distributions. J. Appl. Prob. 21, 6544660. Despande, J.V. and R.P. Suresh (1990). Non-monotonic ageing. Stand. J. Statist. 17, 257-262. Guess, F., M. Hollander and F. Proschan (1986). Testing exponentially versus a trend change in mean residual life. Ann. Statist. 14, 1388-1398. Loeve, M. (1963). Probability Theory, 3rd ed. Van Nostrand, New York.