1. Introduction
The idea of using Laplace transforms to characterize ageing classes was introduced in reliability literature by Vinogradov [Reference Vinogradov40], who obtained the first result of this kind in the context of the fundamental class of life distributions having increasing failure rate (IFR). Subsequently, Block and Savits, in a seminal paper [Reference Block and Savits11], went on to develop Laplace transform characterization results in the spirit of Vinogradov for other well-known ageing classes, namely increasing failure rate average (IFRA), decreasing mean residual life (DMRL), new better than used (NBU), and new better than used in expectation (NBUE). Analogous results for the harmonic new better than used in expectation (HNBUE) class were obtained in [Reference Klefsjö25, Reference Klefsjö26]. Further, a new notion of life distribution, new better than used in convex ordering (NBUC) was introduced and characterized in [Reference Cao and Wang12]. In the same vein, the present authors attempt to develop new Block–Savits type characterization theorems based on Laplace transforms.
For a non-negative random variable X with distribution function (DF) F, the Laplace–Stieltjes transform of F is given by
$\phi_F(s) = \int_0 ^\infty {\text{e}}^{-sx}\, {\text{d}} F(x)$
,
$s>0$
. Define
and
\begin{equation} \begin{cases} \alpha_{n+1}(s) = s^{n+1}a_n(s), \quad & n\geq0,\,s\geq 0, \\ \alpha_0(s) = 1, & \text{for all }\, s \geq 0. \end{cases}\end{equation}
Then we have
\begin{alignat}{2} a_{n}(s) & = \int_0^\infty \frac{x^n}{n!}{\text{e}}^{-sx}\bar{F}(x)\,{\text{d}} x, & & n\geq0,\,s\geq 0, \nonumber \\ \alpha_{n+1}(s) & = s\int_0^\infty\frac{(sx)^n}{n!}{\text{e}}^{-sx}\bar{F}(x)\,{\text{d}} x, \quad & & n\geq0,\,s\geq 0, \end{alignat}
where
$\bar F=1-F$
is the survival function of X.
Among various possible ways of characterizing a life distribution, the one based on
$\alpha_n(s)$
has special significance due to the specific forms that arise in this approach. More specifically, conditions on the sequence
$\{\alpha_n(s),{n\geq1}\}$
are the exact analogs of properties possessed by the discrete survival probabilities
$\bar{P}_k$
in the context of a homogeneous Poisson shock model. However, unlike
$\alpha_n(s)$
,
$\bar{P}_k$
fails to characterize the ageing classes using shock models as in this case the results obtained are only sufficient but not necessary. Furthermore,
$\alpha_n(s)$
can be viewed as the discrete survival probabilities for a cumulative damage shock model with random threshold X independent of the damages
$Y_1, Y_2, \ldots$
coming from successive shocks, where
$Y_1, Y_2, \ldots$
are also mutually independent and each follows an exponential distribution with rate parameter
$s>0$
(see [Reference Block and Savits11] for details). Considering a random variable
$X_n(s)$
that follows a gamma distribution with shape and scale parameters n and s respectively, an interpretation of
$\alpha_n(s)$
in a probabilistic sense is that it can be looked upon as the expectation of
$\bar{F}(X_n(s))$
. In view of the manifold interpretations discussed above, characterization results based on
$\alpha_n(s)$
have a very definite appeal. Such results will henceforth be referred to as Block–Savits (BS) type characterization results. The present authors noticed that there are some important reliability classes such as decreasing mean time to failure (DMTTF), increasing mean inactivity time (IMIT), decreasing mean residual life in harmonic average (DMRLHA), new better than used in failure rate (NBUFR), and new better than used average in failure rate (NBAFR), as well as several renewal classes, namely new better than used renewal failure rate (NBURFR), new better than average used renewal in failure rate (NBARFR), new better than renewal used in expectation (NBRUE), and harmonic new better than renewal used in expectation (HNBRUE) for which BS type characterization results are as yet absent in the literature. The present article seeks to provide concrete answers to these problems. For a quick glance at these classes as well as the interrelationships hold between them, see [Reference Abouammoh and Ahmed1–Reference Abouammoh, Ahmed and Barry3, Reference Bhattacharjee, Abouammoh, Ahmed and Barry9, Reference Chandra and Roy13, Reference Deshpande, Kochar and Singh14, Reference Klefsjö27, Reference Li and Xu30, Reference Loh32]. Figure 1 presents a schematic representation depicting connections between the various ageing classes that are the focus of this article.
Interrelationships between ageing classes.

We have managed to obtain the desired characterization results for some of the above classes by exploiting the techniques of [Reference Block and Savits11]. However, for some other families it appears that this method of proof requires the imposition of certain restrictions on the underlying DF, leading to comparatively weaker results, while in other scenarios the calculations turn out to be extremely involved. To find a way out, we introduce a novel technique based on an important connection between F and its equilibrium (or first derived) distribution. Our method is described in Section 2.3, and applications thereof can be found in Sections 2.4–2.6 and 3.2. The simplicity of the proofs demonstrates the effectiveness of the approach. Besides characterizing the new classes, this connection has also provided several useful insights. Exploiting this connection and the theory of total positivity, we have managed to prove that, irrespective of the presence of any ageing criterion,
$\alpha_n(s)$
inherits a certain property of
$\bar{F}$
for any continuous DF F. Interestingly enough, the result also goes the other way. This is analogous to [Reference Block and Savits11, Corollary 3.7], with the density replaced by the survival function in our case. Further, a particular case of our result yields a necessary and sufficient condition for a continuous DF F to be IFR. Our approach leads to a simple proof of Vinogradov’s characterization of the IFR class for the continuous case. We have also furnished a short and elegant proof for the result of Block and Savits in the context of the DMRL, NBUE, and HNBUE classes. In these cases, however, no restrictions are necessary. Although our approach is much easier to adopt to prove the characterization results for certain ageing classes, it has certain limitations as well (see Section 5).
While our main goal has been to characterize the classes of life distributions in terms of the Laplace transform, our investigations have necessitated the development of some new auxiliary results and even a new definition. To the best of our knowledge, the reliability literature so far does not contain any definition of the NBUFR class except in the absolutely continuous case. In order to retain the somewhat general flavour of our results, we have introduced a general definition of the NBUFR class without this restriction. Furthermore, we have managed to relax a condition of [Reference Shokry, Ahmed, Rakha and Hewedi39, Theorem 2.1] for the said class. Motivated by [Reference Shokry, Ahmed, Rakha and Hewedi39], a moment inequality for the NBAFR class has been derived and, further utilizing this result, the characterization for this class has been established. For F belonging to the NBAFR class, an upper bound for the quantity
$r_F(0)$
(the definition of
$r_F(0)$
is given in the appropriate subsection) is also obtained. Moving on to the renewal classes, namely NBRUE, HNBRUE, NBURFR, and NBARFR, we have established that they possess finite moments of all orders.
We therefore organize this article as follows. In Section 2 we characterize the ageing classes namely DMTTF, IMIT, DMRLHA, NBUFR, and NBAFR. Some useful results related to the ageing class under consideration, along with a brief discussion, are also presented in the appropriate subsections before establishing our main theorems. We also derive the relation between
$\alpha_n(s)$
and its analogue for the equilibrium distribution; this is exploited in some of the proofs. Section 3 is devoted to characterizing the renewal classes, NBURFR, NBARFR, NBRUE, and HNBRUE. Finally, in Section 4 we focus on the issue of convolutions. We wind up with a brief discussion in Section 5.
Throughout the paper, F is considered to be a life distribution with
$F(0) = 0$
, and the term ‘increasing’ (decreasing) means ‘non-decreasing’ (non-increasing).
2. Reliability classes and their characterizations
We start out by proving the following proposition, which will be extensively used in many of the theorems in this paper. Note that the result holds not only for a DF but also for a broader class.
Proposition 2.1. Let
$g\colon\mathbb{R^+}\to\mathbb{R^+} $
be any measurable function. Then
Proof. From the survival function of an Erlang distribution with parameters n and
$\lambda$
, we have the identity
\begin{equation} \int_t^\infty\frac{\lambda^n}{\Gamma(n)}{\text{e}}^{-\lambda u}u^{n-1}\,{\text{d}} u = \sum_{k=0}^{n-1}\frac{{\text{e}}^{-\lambda t}(\lambda t)^k}{k!} \quad \text{for all } t > 0 \text{ and } n\in\mathbb N,\,\lambda>0. \end{equation}
Now, applying Tonelli’s theorem [Reference Axler6, p. 129] and the non-negativity of the integrand yield, for all
$n\geq0$
and
$s>0$
,
\begin{align*} s\int_0^\infty\frac{{\text{e}}^{-st}(st)^n}{n!}\bigg(\int_0^tg(y)\,{\text{d}} y\bigg)\,{\text{d}} t & = \int_0^\infty\int_y^\infty\frac{s{\text{e}}^{-st}(st)^n}{n!}g(y)\,{\text{d}} t\,{\text{d}} y \\ & = \sum_{i=0}^n\int_0^\infty\frac{{\text{e}}^{-sy}(sy)^i}{i!}g(y)\,{\text{d}} y. \end{align*}
Lemma 2.1. Let
$s=s(n,x)>0$
and
$\lim_{n\to\infty}{n}/{s}=x>0$
. Then, as
$n\to \infty$
,
$G_n$
converges in distribution to G where
\begin{equation} G_n(y) = s\int _0^y\frac{{\text{e}}^{-su}(su)^{n}}{n!}\,{\text{d}} u, \qquad G (y) = \begin{cases} 1, & y\geq x, \\ 0, & y < x. \end{cases} \end{equation}
Proof. The result is stated without proof in [Reference Block and Savits11]. For convenience, we include a formal proof here. Let
$Y_1,Y_2,\ldots$
be mutually independent exponential random variables with parameter s. Then
$G_n(y)$
is nothing but the DF of
$\sum_{i=1}^{n+1}Y_i$
. Now, applying central limit theorem we obtain that
$\sum_{i=1}^{n+1}Y_i$
converges in law to a normal distribution with mean x and variance
$x^2/(n+1)$
, which tends to 0 as
$n\to\infty$
; in other words, the DF
$G_n(y)$
converges weakly to a one-point distribution concentrated at the point x.
Remark 2.1. Using identity (2.1),
$\bar{G}_n(y)$
can be expressed as
$\bar{G}_n(y) = \sum_{k=0}^n{{\text{e}}^{-sy}(sy)^k}/{k!}$
.
Lemma 2.2. [[Reference Block and Savits11]] Let
$x > 0$
be a continuity point of the DF F. Let
$s = s(n,x)$
and
$\lim_{n\to\infty}{n}/{s}=x$
. Then
$\lim_{n\to\infty}\alpha_{n+1}(s) = \bar{F}(x)$
, where
$\alpha_{n+1}(s)$
is as defined in (1.2).
Before proceeding further, we recall an important definition and a crucial result.
Definition 2.1. [[Reference Karlin19]] Let A and B be two subsets of
$\mathbb R$
. A function
$K\colon A\times B\to\mathbb R$
is said to be totally positive of order r (TP
$_r$
) if
\begin{equation*} \begin{vmatrix} K(x_1,y_1) &\quad K(x_1,y_2) & \quad\cdots & \quad K(x_1,y_m) \\ K(x_2,y_1) & \quad K(x_2,y_2) & \quad \cdots & \quad K(x_2,y_m) \\ \quad \vdots & \quad \vdots & \quad & \quad \vdots \\ K(x_m,y_1) & \quad K(x_m,y_2) & \quad \cdots & \quad K(x_m,y_m) \end{vmatrix} \geq 0 \end{equation*}
for all
$x_1 < x_2 < \cdots < x_m \in A$
and
$y_1 < y_2 < \cdots < y_m\in B$
,
$m=1,2,\ldots, r$
.
Now consider a function
$f\colon I\to\mathbb R$
, where I is an ordered set of
$\mathbb R$
. Let
where m is arbitrary but finite and
$S(x_1,x_2,\ldots,x_m)$
denotes the number of sign changes of the sequence
$x_1,x_2,\ldots,x_m$
, disregarding zero terms [Reference Barlow and Proschan8, p. 93]. We now state the variation diminishing property (VDP) of a totally positive function.
Lemma 2.3. ([Reference Barlow and Proschan8].) Let K(x,y) be TP
$_r$
on
$A\times B$
, where
$A,B\subset\mathbb R$
, and let f be a bounded and Borel-measurable function defined on B. Let the transformation
$g(x) = \int_BK(x,y)f(y)\,{\text{d}} y$
be finite for each
$x\in A$
. Then
$S(g)\leq S(f)$
provided
$S(f) \leq r-1$
. Furthermore, if
$S(g)= S(f) \leq r-1$
, then the functions f and g exhibit the same sequence of signs in their respective arguments as they traverse the domain of definition from left to right.
Having completed the necessary prerequisites, we now move on to the main results.
2.1. The decreasing mean time to failure class
The notion of a mean time to failure (MTTF) function was first introduced in [Reference Barlow and Proschan7] in the context of an age replacement policy where the replacement of an item takes place at the failure of the item or at a predetermined time t, whichever comes first. It is a standard practice to adopt an age replacement policy in order to control in-service failure of systems. Let
$X_{[t]}$
denote the time to the first in-service failure of an item under an age replacement policy with the age replacement time t. Then the MTTF function of F is defined as
This function itself has a special importance in the literature [Reference Kayid, Ahmad, Izadkhah and Abouammoh21, Reference Klefsjö27, Reference Li and Xu31]. We now have the following formal definition of the DMTTF and increasing mean time to failure (IMTTF) classes.
Definition 2.2.
F is said to belong to a DMTTF (IMTTF) class if
$M_F(t)$
is decreasing (increasing) in t.
Another ageing family, NBRU
$_{\text{rh}}$
(new better than renewal used in the reversed hazard rate order), which was introduced in [Reference Li and Xu31], turns out to be identical to the DMTTF class. The preservation of the DMTTF class under some basic reliability operations was investigated in [Reference Knopik28], while the weak convergence issues as well as closure under convolution for this class were studied in [Reference Knopik29]. An extensive study on the monotonic as well as non-monotonic classes based on the MTTF function can be found in [Reference Izadi, Sharafi and Khaledi18, Reference Khan, Bhattacharyya and Mitra23]. However, characterization results for this class are so far absent in the literature. The next theorem addresses this. We first prove a lemma that is essential for the next theorem.
Lemma 2.4. Let
$s>0$
. Suppose that, for each
$c>0$
,
$({1}/{s})\sum_{i=1}^{n+1}\alpha_i(s)-c(1-\alpha_{n+1}(s))$
has at most one change of sign from
$+$
to
$-$
. Then
$({1}/{s})\sum_{i=1}^{n+1}\alpha_i(s)/(1-\alpha_{n+1}(s))$
is decreasing in n, where
$\alpha_{n}(s)$
is defined as in (1.2).
Proof. We need to show that, for all
$0\leq n_1 < n_2$
and
$s>0$
,
$$ \frac{1}{s}\sum_{i=1}^{n_1+1}\alpha_i(s)/(1-\alpha_{n_1+1}(s)) \geq \frac{1}{s}\sum_{i=1}^{n_2+1}\alpha_i(s)/(1-\alpha_{n_2+1}(s)). $$
Define
$g(k)= ({1}/{s})\sum_{i=1}^{k+1}\alpha_i(s)$
and
Now,
$({1}/{s})\sum_{i=1}^{k+1}\alpha_i(s)-c(1-\alpha_{k+1}(s))$
has at most one change of sign from
$+$
to
$-$
for each
$c>0$
. So this happens for
$c_1$
and
$c_2$
in particular. Thus,
$g(k)-c_j(1-\alpha_{k+1}(s))$
has at most one change of sign from
$+$
to
$-$
,
$j=1,2$
. Consequently, by choice of
$c_j$
,
$j=1,2$
, we have
Now, if possible, let
Then, taking
$n_1<k<n_2$
where
$k\in \mathbb N$
, we have
\begin{alignat*}{2} g(k) & \leq c_1(1-\alpha_{k+1}(s)), \qquad & & (k > n_1) \\ & < c_2(1-\alpha_{k+1}(s)), & & (\text{by ({2.3})}) \\ & \leq g(k), & & (k < n_2), \end{alignat*}
which is a contradiction. Therefore, we have
${g(n_1)}/({1-\alpha_{n_1+1}(s)}) \geq {g(n_2)}/({1-\alpha_{n_2+1}(s)})$
for all
$0\leq n_1 < n_2$
and
$s>0$
, i.e.
$({1}/{s})\sum_{i=1}^{n+1}\alpha_i(s)/(1-\alpha_{n+1}(s))$
is decreasing in
$n\geq 0$
for all
$s>0$
.
Theorem 2.1. F is DMTTF if and only if
$\sum_{k=0}^{n-1} \alpha_k(s)/(1-\alpha_n(s))$
is decreasing in
$n\geq 1 $
and for every
$s>0$
.
Proof. Necessity. For any life DF F we have, by Proposition 2.1,
Therefore, for
$n\geq 0$
,
$s >0$
, and for any
$c> 0$
,
\begin{equation} s\int_0^\infty\frac{{\text{e}}^{-st}(st)^n}{n!}\bigg[\int_0^t\bar{F}(y)\,{\text{d}} y - cF(t)\bigg]\,{\text{d}} t = \frac{1}{s}\sum_{i=1}^{n+1}\alpha_i(s) - c(1-\alpha_{n+1}(s)). \end{equation}
Observe that the kernel
$K(n,t)= {\text{e}}^{-st}(st)^n/n!$
is TP
$_2$
in
$t\in(0,\infty$
) for
$n\in \{0,1,2,\ldots\}$
. Since F is DMTTF, the variation diminishing property (VDP) of K(n,t) together with (2.4) imply that, for any
$c> 0$
, the expression
$({1}/{s})\sum_{i=1}^{n+1}\alpha_i(s)-c(1-\alpha_{n+1}(s))$
has at most one sign change from
$+$
to
$-$
. Therefore, by virtue of Lemma 2.4,
$({1}/{s})\sum_{i=1}^{n+1}\alpha_i(s)/(1-\alpha_{n+1}(s))$
is decreasing in
$n\geq 0$
for all
$s>0$
. This further implies that
$\sum_{i=1}^{n}\alpha_i(s)/(1-\alpha_{n}(s))$
is decreasing in
$n> 0$
for all
$s>0$
. Now adding 1 to this expression and using the condition
$\alpha_0(s)=1$
, the necessary part is proved.
Sufficiency. It is given that
$\sum_{k=0}^{n-1} \alpha_k(s)/ (1-\alpha_n(s))$
is decreasing in
$n> 0$
for all
$s>0$
. Therefore, for all
$n,k \geq 0$
,
\begin{equation} (1-\alpha_{n+k+1}(s))\sum_{i=1}^{n+1}\alpha_i(s) \geq (1-\alpha_{n+1}(s))\sum_{i=1}^{n+k+1}\alpha_i(s). \end{equation}
Now,
\begin{align*} \sum_{i=1}^{n+1}\alpha_i(s) = s\sum_{i=0}^{n}\int_0^\infty\frac{{\text{e}}^{-su}(su)^{i}}{i!}\bar{F}(u)\,{\text{d}} u & = s\int_0^\infty\Bigg(\sum_{i=0}^n\frac{{\text{e}}^{-su}(su)^{i}}{i!}\Bigg)\bar{F}(u)\,{\text{d}} u \\ & = s\int_0^\infty\bar{G}_n(u)\bar{F}(u)\,{\text{d}} u. \end{align*}
Let
$\lim_{n \to \infty}({n}/{s}) = x > 0$
. Then, using Lemma 2.1 and appealing to the dominated convergence theorem (DCT), we get
\begin{equation} \lim_{n\to\infty}\frac{1}{s}\sum_{i=1}^{n+1}\alpha_i(s) = \lim_{n\to\infty}\int_0^\infty\bar{G}_n(u)\bar{F}(u)\,{\text{d}} u = \int_0^x\bar{F}(u)\,{\text{d}} u. \end{equation}
Now let x and
$x+y$
be two continuity points of F, where
$x,y>0$
. Also let
$s=s(n,x)={n}/{x}$
and
$k= [{ny}/{x}]$
, where [w] indicates the greatest integer less than or equal to w. Then
$\lim_{n\to\infty}({n+k})/{s} = x+y$
. Now, taking the limit as
$n\to\infty$
on both sides of (2.5) and using Lemma 2.2 and (2.6), we finally obtain
$F(x+y)\int_0^x\bar{F}(u)\,{\text{d}} u \geq F(x)\int_0^{x+y}\bar{F}(u)\,{\text{d}} u$
. This inequality holds for all
$x,y \geq 0$
since the continuity points of F are dense in
$[0,\infty)$
. Therefore, F is DMTTF.
2.2. The increasing mean inactivity time class
For any fixed
$t>0$
, the inactivity time (IT) [Reference Ruiz and Navarro38] defined by
$X_{(t)}\;{:}{=}\;[t-X\mid X\leq t]$
is the random time that has elapsed until time t after the failure occurs. Then the mean inactivity time (MIT) function is given by
This function has widespread applications in diverse fields and has also been called the mean waiting time (MWT) function [Reference Finkelstein15] or mean past life (MPL) function [Reference Asadi5, Reference Goliforushani and Asadi17]. It was pointed out in [Reference Finkelstein15] that the MIT function can be treated as the dual of the mean residual life (MRL) function. In recent work, preservation of MIT ordering for coherent systems was studied in [Reference Rao and Naqvi37]. The monotonic behaviour of
$m_F(t)$
induces the following ageing class.
Definition 2.3. A distribution function F is said to belong to the increasing mean inactivity time (IMIT) class if
$m_F(t)$
is increasing in t.
Note that the concept of a decreasing mean inactivity time class does not hold for non-negative distributions [Reference Chandra and Roy13]. The behaviour of the excess lifetime of a renewal process with IMIT interarrivals was investigated in [Reference Li and Xu30]. Characterization results for the IMIT class using some ordering relations such as right spread order, increasing convex order, etc. can be found in [Reference Ahmad, Kayid and Pellerey4]. Closure under weak convergence of the IMIT class was established in [Reference Khan, Bhattacharyya and Mitra24]. In order to obtain a BS type characterization result, we proceed as follows. We consider a quantity analogue of
$\alpha_n(s)$
in this situation. This was first considered in [Reference Kebir22] while characterizing various probabilistic orderings using the Laplace transform. Define
and
$\bar{\alpha}_{n+1}(s)=s^{n+1}\bar{a}_n(s)$
for
$n\geq 0$
,
$s>0$
. Therefore,
Recalling (1.2), it is clear that
\begin{equation} \begin{cases} \bar{\alpha}_{n+1}(s) = 1-\alpha_{n+1}(s), & n\geq0,\,s>0, \\ \bar{\alpha}_0(s) = 0, & s > 0. \end{cases}\end{equation}
Observe that, under the setup of the cumulative damage shock model with random threshold as described in the introduction,
$\bar{\alpha}_{n}(s)$
denotes the discrete probability distribution for that specific model.
$\bar{\alpha}_{n}(s)$
can also be viewed as
$\bar{\alpha}_{n}(s) = \mathbb{E}[F(X_n(s))]$
, where the random variable
$X_n(s)$
has a Gamma(n, s) distribution. Next, we establish the characterization result for the IMIT class.
Theorem 2.2. F is IMIT if and only if
$\sum_{k=0}^{n-1}\bar{\alpha}_k(s)/\bar{\alpha}_n(s)$
is increasing in
$n\geq 1$
and for all
$s>0$
.
Proof.
Necessity. Using Proposition 2.1, it is easy to check that for
$n\geq 0$
,
$s>0$
, and for any
$c\geq0$
,
\begin{equation*} s\int_0^\infty\frac{{\text{e}}^{-st}(st)^n}{n!}\bigg[\int_0^tF(y)\,{\text{d}} y - cF(t)\bigg]\,{\text{d}} t = \frac{1}{s}\sum_{i=1}^{n+1}\bar{\alpha}_i(s) - c\bar{\alpha}_{n+1}(s). \end{equation*}
Since F is IMIT, the function
$\int_0^tF(y)\,{\text{d}} y-cF(t)$
can undergo at most one change of sign from
$-$
to
$+$
. Thus, using the VDP of the kernel
$K(n,t) = {\text{e}}^{-st}(st)^n/n!$
, we obtain that
$\sum_{i=1}^{n+1}\bar{\alpha}_i(s)/\bar{\alpha}_{n+1}(s)$
is increasing in
$n\geq 0$
and for all
$s>0$
. Observing that
$\bar{\alpha}_0(s) = 0$
from (2.7), we get that
$\sum_{i=0}^{n}\bar{\alpha}_i(s)/\bar{\alpha}_{n}(s)$
is increasing in
$n> 0$
for all
$s>0$
. We complete the proof by adding
$-1$
to this expression.
Sufficiency. From the given condition we have
\begin{equation} \bar{\alpha}_{n+k+1}(s)\sum_{i=1}^{n+1}\bar{\alpha}_i(s) \leq \bar{\alpha}_{n+1}(s)\sum_{i=1}^{n+k+1}\bar{\alpha}_i(s) \quad \text{for all } n,k \geq 0,\,s>0. \end{equation}
A calculation analogous to the DMTTF case yields
$({1}/{s})\sum_{i=1}^{n+1}\bar{\alpha}_i(s) = \int_0^\infty\bar{G}_n(u)F(u)\,{\text{d}} u$
. Changing the order of integration and substituting
$\bar{G}_n(u)$
from (2.2), we have
where
$h(y)=\int_0^yF(u)\,{\text{d}} u$
. Let
$\lim_{n\to\infty}({n}/{s}) = x > 0$
; applying [Reference Billingsley10, Theorem 5.2] for
$n\to\infty$
we obtain
\begin{equation} \frac{1}{s}\sum_{i=1}^{n+1}\bar{\alpha}_i(s) \to h(x) = \int_0^xF(u)\,{\text{d}} u. \end{equation}
Consider x and
$x+y$
to be two continuity points of F, where
$x,y>0$
. As before, let
$s={n}/{x}$
and
$k=[{ny}/{x}]$
. Letting
$n\to\infty$
in (2.8) and using Lemma 2.2 in conjunction with (2.9), we obtain
$F(x+y)\int_0^xF(u)\,{\text{d}} u \leq F(x)\int_0^{x+y}F(u)\,{\text{d}} u$
. The proof is now completed by noting that the set of continuity points of F is dense in
$[0,\infty)$
.
2.3. Some important connections
The techniques employed by Block and Savits serve quite well in establishing characterization results for a number of ageing classes. However, in some cases either the method fails or the calculations become intractable. Searching for an alternative route, we have developed a new approach involving the equilibrium distribution of F. Let the DF of the equilibrium (or first derived) distribution of F(t) be denoted by
$F_{(1)}(t)\;{:}{=}\; ({1}/{\mu})\int_0^t\bar{F}(u)\,{\text{d}} u$
for all
$t\geq0$
, where
$\mu$
is the mean of F. For a life-distribution function F, it is natural to assume that the expected lifetime is positive and finite. We start out by looking at an analogue of
$\alpha_n(s)$
for the equilibrium distribution. Let
$\phi_{F_{(1)}}(s) = \int_0^\infty{\text{e}}^{-sx}\,{\text{d}} F_{(1)}(x)$
,
$s>0$
, so that
and
\begin{equation} \begin{cases} \alpha'_{n+1}(s) =s^{n+1}a^{\prime}_n(s), & n\geq0,\,s\geq 0, \\ \alpha'_0(s) = 1 & \text{for all } s \geq 0. \end{cases}\end{equation}
Note that we use a
$'$
to denote the quantities corresponding to (1.1) for the equilibrium distribution. We now have the following important relation, which is the key to our approach.
Proposition 2.2. For
$\alpha_{n}(s)$
and
$\alpha'_{n}(s)$
as defined in (1.1) and (2.10),
Proof. We have, for all
$n\geq0$
and
$s\geq 0$
,
\begin{align*} \alpha'_{n+1}(s) & = s\int_0^\infty\frac{(st)^n}{n!}{\text{e}}^{-st}\bar{F}_{(1)}(t)\,{\text{d}} t \\ & = s\int_0^\infty\frac{(st)^n}{n!}{\text{e}}^{-st}\bigg(\frac{1}{\mu}\int_t^\infty\bar{F}(u)\,{\text{d}} u\bigg)\,{\text{d}} t \\ & = \frac{s}{\mu}\int_0^\infty\int_0^u\frac{(st)^n}{n!}{\text{e}}^{-st}\bar{F}(u)\,{\text{d}} t\,{\text{d}} u \\ & = \frac{1}{\mu s}\sum_{k=n+1}^\infty s\int_0^\infty\frac{{\text{e}}^{-su}(su)^k}{k!}\bar{F}(u)\,{\text{d}} u \qquad \text{(using Proposition 2.1)} \\ & = \frac{1}{\mu s}\sum_{k=n+2}^\infty\alpha_{k}(s). \end{align*}
The proof is completed by noting that
From Proposition 2.2, it is easy to observe that
At this juncture, it would be appropriate to highlight a few comments on
$\alpha'_{n}(s)$
. To this end, we have the following lemma.
Lemma 2.5. Let
$s=s(n,x)>0$
and
$\lim_{n\to\infty}({n}/{s})=x>0$
. Then
$\lim_{n\to\infty}\alpha'_{n}(s) = \bar{F}_{(1)}(x)$
.
Proof. Using Proposition 2.2, the proof follows directly from the facts that
as
$n\to\infty$
and
$\mu s = \sum_{k=1}^\infty\alpha_{k}(s)$
.
Again,
$\alpha'_{n}(s)$
can be written as
where
$X_n(s)$
is a random variable having a gamma distribution with parameters n and s. Thus, by virtue of Lemma 2.5 and (2.12) it is easy to observe that properties of
$\alpha'_{n}(s)$
can be translated into properties of
$\bar{F}_{(1)}$
, and vice versa. The same was shown to be true in [Reference Block and Savits11] for the functions
$\alpha_{n}(s)$
and
$\bar F$
. Now consider a cumulative damage shock model with random threshold
$X_e$
having DF
$F_{(1)}$
. Let the damages
$Y_1, Y_2, \ldots$
caused by the successive shocks be mutually independent, each having exponential distribution with parameter
$s>0$
. Also,
$Y_1, Y_2, \ldots$
are independent of the threshold
$X_{e}$
. Then
$\alpha'_{n}(s)$
can be considered as the discrete survival probabilities for that particular model.
Let us now consider a homogeneous Poisson process
$\{N_s(t)\colon t\geq 0\}$
with intensity
$s>0$
. Then
As in [Reference Block and Savits11], we define
\begin{equation} \begin{cases} \beta_{n}(s) = \mathbb{P}(N_s(X)=n), & n\geq0,\,s>0, \\ \beta_n(0) = 0, & n > 0, \\ \beta_0(0) = 1. & \end{cases}\end{equation}
Then, for
$n\geq0$
,
$s\geq 0$
, we obtain
Therefore,
$\beta_{n}(s) = \alpha_n(s) - \alpha_{n+1}(s)$
. Considering the above quantities for the equilibrium distribution
$F_{(1)}$
, this leads to
where
From (2.11) and (2.14) we now note that
Observe that (2.15) is another connection between the DF and its equilibrium distribution by means of the Laplace transform. Let
$f_{(1)}(x) = {\bar{F}(x)}/{\mu}$
be the density of
$F_{(1)}$
. Then, in the presence of the relation in (2.15), we can restate [Reference Block and Savits11, Corollary 3.7] for the first derived distribution. Before going to the lemma, we first recall the definition of PF
$_r$
.
Definition 2.4. A function
$f(x)$
defined on
$({-}\infty,\infty)$
is said to be a P
$\acute{\text{o}}$
lya frequency function of order r (PF
$_r$
) if
$K(x,y) = f(x - y)$
is TP
$_r$
for all
$-\infty < x,y < \infty$
.
Lemma 2.6. Let
$f_{(1)}$
be a continuous density on
$[0,\infty)$
. Then
$f_{(1)}$
is PF
$_{r}$
in
$x>0$
if and only if
$\beta'_{n}(s)$
is PF
$_{r}$
in
$n \geq 0$
for any
$s>0$
.
Note that ‘
$f_{(1)}$
is PF
$_{r}$
in
$x>0$
’ is equivalent to ‘
${\bar{F}(x)}/{\mu}$
is PF
$_{r}$
in
$x>0$
’, i.e.
$\bar{F}(x)$
is PF
$_{r}$
in
$x>0$
, as
$\mu >0$
. Again using (2.15), ‘
$\beta'_{n}(s)$
is PF
$_{r}$
in
$n \geq 0$
for any
$s>0$
’ is equivalent to saying ‘
$({1}/{\mu s})\alpha_{n+1}(s)$
is PF
$_{r}$
in
$n \geq 0$
for any
$s>0$
’, i.e. if and only if, for any
$s>0$
,
$\alpha_{n+1}(s)$
is PF
$_{r}$
in
$n \geq 0$
. Thus we have the following proposition.
Proposition 2.3. Let
$\bar{F}$
be continuous on
$[0,\infty)$
. Then
$\bar{F}(x)$
is PF
$_r$
in
$x>0$
if and only if
$\alpha_{n+1}(s)$
is PF
$_{r}$
in
$n \geq 0$
for any
$s>0$
.
The particular case
$r=2$
yields the following corollary, which is nothing but Vinogradov’s result in the continuous case.
Corollary 2.1. A continuous life DF F is IFR if and only if
$\alpha_{n+1}(s)$
is PF
$_{2}$
in
$n \geq 1$
for any
$s>0$
.
Therefore, exploiting the theory of total positivity in conjunction with our results involving the equilibrium distribution, we arrive at the idea of characterizing the IFR class in the special case when F is continuous.
Recall the theorem in [Reference Vinogradov40] that states that a necessary and sufficient condition for a DF F to be an IFR class is
$a^{2}_n(s) \geq a^{}_{n-1}(s)a^{}_{n+1}(s)$
for all
$n \geq 0$
, i.e., using (1.1),
This is equivalent to
${\alpha^{}_{n+1}(s)}/{\alpha^{}_{n}(s)}$
being decreasing in
$n=0,1,2,\ldots$
for all
$s>0$
. Now it is easy to prove that F has decreasing mean residual life if and only if
$F_{(1)}$
is IFR. Using this together with Proposition 2.2, we can provide an alternative proof of Block and Savits’ characterization for the DMRL class. For completeness, we first state the theorem as given in [Reference Block and Savits11] and then present our proof.
Theorem 2.3 ([Reference Block and Savits11].) For non-negative distributions F, F is DMRL if and only if
is decreasing in
$n = 1,2,\ldots$
for all
$s > 0$
.
Proof.
$F_{(1)}$
is IFR if and only if
${\alpha'_{n+1}(s)}/{\alpha'_{n}(s)}$
is decreasing in
$n=0,1,2,\ldots$
and
$s>0$
. Replacing
$\alpha'_{n}(s)$
by using Proposition 2.2, we have that F is DMRL if and only if
\begin{alignat*}{2} & \sum_{k=n+1}^\infty\alpha_{k}(s)/\sum_{k=n}^\infty\alpha_{k}(s) & & \text{is decreasing in } n>0 \text{ and } s>0 \\ \Leftrightarrow\ & \Bigg[\alpha_{n}(s)+\sum_{k=n+1}^\infty\alpha_{k}(s)\Bigg]/\sum_{k=n+1}^\infty\alpha_{k}(s) \quad & & \text{is increasing in } n>0 \text{ and } s>0 \\ \Leftrightarrow\ & \sum_{k=n+1}^\infty\alpha_{k}(s)/\alpha_{n}(s) & & \text{is decreasing in } n>0 \text{ and } s>0 \\ \Leftrightarrow\ & \sum_{k=n}^\infty\alpha_{k}(s)/\alpha_{n}(s) & & \text{is decreasing in } n>0 \text{ and } s>0. \end{alignat*}
This illustrates our new approach of obtaining BS type characterizations using the machinery developed here. We shall now apply this technique to prove further interesting results.
2.4. The decreasing mean residual life in harmonic average class
In [Reference Deshpande, Kochar and Singh14], the authors illustrated various new concepts of ageing criteria that arise as a consequence of positive ageing of equipment, and the notion of this particular ageing class is one such. Let
$e_F(t)=({1}/{\bar{F}(t)})\int_t^\infty\bar{F}(u)\,{\text{d}} u$
be the mean residual life (MRL) function of F. We have the following definition.
Definition 2.5.
F is said to be decreasing mean residual life in harmonic average (DMRLHA) (increasing mean residual life in harmonic average (IMRLHA)) if
$\big[({1}/{t})\int_0^t{1}/{e_F(x)}\,{\text{d}} x\big]^{-1}$
is a decreasing (increasing) function in t.
The interrelationships of the DMRLHA class with the other ageing classes are also discussed in [Reference Deshpande, Kochar and Singh14]. There are well-established BS type characterization results for various ageing classes such as IFRA, NBUE, HNBUE, etc. in the existing literature [Reference Block and Savits11, Reference Klefsjö26, Reference Klefsjö27, Reference Vinogradov40]. We have already demonstrated the possibility of obtaining such characterizations using the appropriate ageing properties of the corresponding equilibrium distribution in the context of the DMRL class. Motivated by this, we now try to exploit the same technique to obtain new BS type characterizations. In particular, our new approach yields a simple proof of the characterization result for the DMRLHA class, whereas the standard techniques, outlined in [Reference Block and Savits11], require an involved calculation. We thus adopt this alternative method of proof.
Theorem 2.4. F is DMRLHA if and only if
$\big[\sum_{k=n+1}^\infty\alpha_{k}(s)/\sum_{k=1}^\infty\alpha_{k}(s)\big]^{1/n}$
is decreasing in
$n\geq 1$
for all
$s> 0$
.
Proof. According to [Reference Deshpande, Kochar and Singh14, Proposition 5.2], F is DMRLHA
$\Leftrightarrow$
$F_{(1)}$
is IFRA. From [Reference Block and Savits11, Theorem 1.1] we have that
$F_{(1)}$
is IFRA
$\Leftrightarrow$
$[\alpha_{n}'(s)]^{1/n}$
is decreasing in
$n\geq1$
for all
$s>0$
$\Leftrightarrow$
$\big[\sum_{k=n+1}^\infty\alpha_{k}(s)/\sum_{k=1}^\infty\alpha_{k}(s)\big]^{1/n}$
is decreasing in
$n\geq 1$
for all
$s>0$
, by using Proposition 2.2. Combining these two observations, F is DMRLHA
$\Leftrightarrow$
$\big[\sum_{k=n+1}^\infty\alpha_{k}(s)/\sum_{k=1}^\infty\alpha_{k}(s)\big]^{1/n}$
is decreasing in
$n\geq 1$
for all
$s>0$
.
As discussed earlier, an interesting observation in the context of BS type characterization results is that we can often anticipate the exact forms that would appear in the characterization results of the classes. However, for the DMRLHA class this is not at all obvious. Next we present an intuitive justification behind the specific form that arises in Theorem 2.4.
Utilizing the inversion formula of the MRL function [Reference Meilijson35],
\begin{equation*} \bar{F}(t) = \begin{cases} \dfrac{e_F(0)}{e_F(t)}\exp\bigg[{-}\displaystyle\int_0^t\dfrac{{\text{d}} u}{e_F(u)}\bigg], & 0\leq t < {F}^{-1}(1), \\ 0, & t\geq {F}^{-1}(1), \end{cases}\end{equation*}
where
${F}^{-1}(1)\;{:}{=}\; \sup\{t\mid F(t)<1\}$
, we have
where
$\mu = e_F(0)$
is the mean of the DF F. Therefore, F is DMRLHA if and only if
\begin{align*} \bigg[\frac{1}{t}\int_0^t\frac{1}{e_F(x)}\,{\text{d}} x\bigg]^{-1} \text{ is decreasing in } t & \Leftrightarrow \ln\bigg(\frac{\mu}{e_F(t)\bar{F}(t)}\bigg)^{1/t} \text{ is increasing in } t \\ & \Leftrightarrow \bigg(\frac{\mu}{e_F(t)\bar{F}(t)}\bigg)^{1/t} \text{ is increasing in } t.\end{align*}
Observing that
$e_F(t)\bar{F}(t)$
corresponds to
$({1}/{s})\sum_{k=n+1}^\infty\alpha_{k}(s)$
and
$\mu s = \sum_{k=1}^\infty \alpha_{k} (s)$
, the expression for the DMRLHA class no longer comes as a surprise.
2.5. The new better than used in failure rate class
The new better than used in failure rate (NBUFR) ageing class was introduced in [Reference Abouammoh and Ahmed1]. The definition given there for this class is reproduced here.
Definition 2.6. An absolutely continuous DF F is said to be NBUFR (NWUFR) with failure rate
$r_F(t)$
if
$r_F(0) \leq (\geq)\, r_F(t)$
for all
$t>0$
, where
The following equivalent form of this definition was also suggested in [Reference Abouammoh and Ahmed1].
Definition 2.7. For an absolutely continuous DF F such that
$\lim_{x\to 0^+}{F(x)}/{x}$
exists, is said to have the NBUFR (NWUFR) property if and only if
$\bar{F}(x+t) \leq (\geq)\, \bar{F}(x){\text{e}}^{-r_F(0)t}$
for all
$x,t\geq 0$
, where
$r_F(0)$
is defined by (2.16).
Preservation properties under various reliability operations as well as shock model theory for both NBUFR and NWUFR have been discussed in [Reference Abouammoh and Ahmed1, Reference Gohout and Kuhnert16]. The non-monotonic analogue of the NBUFR class, namely new worse then better than used in failure rate (NWBUFR), was introduced in [Reference Pandey and Mitra36]. For the NWBUFR class, weak convergence issues, reliability bounds, and interrelationships with some other non-monotonic ageing families are well-studied in [Reference Majumder and Mitra33]. However, the literature so far does not contain any definition of the NBUFR class without the assumption of the existence of the density. We try to redress this. Our definition for this class is motivated by the definition of the IFR class; the notion is broadened by dropping the requirement for absolute continuity. Taking our cue from this, we adopt the definition of
$r_F(0)$
as in [Reference Loh32]:
Henceforth,
$r_F(0)$
is defined to be the quantity given in (2.17). We thus propose the following definition.
Definition 2.8. A DF F is NBUFR (NWUFR) if
$\bar{F}(x+t) \leq (\geq)\, \bar{F}(x){\text{e}}^{-r_F(0)t}$
for all
$x,t\geq 0$
.
Lemma 2.7. Let F be NBUFR such that
$\mu _{r+s+2} = \mathbb{E}(X^{r+s+2}) < \infty$
for integers
$r,s\geq 0$
. Define
$\lambda_r\;{:}{=}\; \mu_r/r!$
for all
$r\geq 0$
. Then
$\lambda_{r+s+2}\leq \lambda_{r+1}/[r_F(0)]^{s+1}$
.
Proof. [Reference Shokry, Ahmed, Rakha and Hewedi39, Theorem 2.1] proves the same result by assuming absolute continuity of F. We note that even if we consider the notion of the NBUFR class as given in Definition 2.8, the proof goes through verbatim.
In fact, we shall show later that the condition
$\mathbb{E}(X^{r+s+2}) < \infty$
can be dropped (see Remark 2.3). Next we state the main result.
Theorem 2.5. F is NBUFR if and only if
$\alpha_{n+m}(s) \leq {1}/{[1+r_F(0)/s]^n}\alpha_m(s)$
for all
$n, m >0$
and
$s>0$
.
Proof.
Necessity. Fix
$s>0$
. From the definition of
$a_{n}(s)$
, we have
where
$\bar{G}(u)= {\text{e}}^{-su}\bar{F}(u)$
. Note that
$\bar{G}(u)$
is NBUFR with
$r_G(0)= (r_F(0)+s)$
if and only if
$\bar{F}(u)$
is NBUFR with
$r_F(0)$
as defined by (2.17). Note that
where
$\mu_{n+1}$
is the
$(n+1)$
th moment of G. Now, using Lemma 2.7 for the DF G, we have
$\lambda_{n+m+2}\leq \lambda_{m+1}/[r_G(0)]^{n+1}$
for all
$n, m \geq 0$
. Utilizing (2.18) gives
$a_{n+m+2}\leq a_{m+1}/[r_G(0)]^{n+1}$
for all
$n, m \geq0$
,
$s>0$
. Now, considering (1.1), we further obtain
$\alpha_{n+m+2}(s) \leq s^{n+1}/[r_G(0)]^{n+1} \alpha_{m+1}(s)$
. Finally, substituting the value of
$r_G(0)$
, the necessary part is proved.
Sufficiency. We have
Now let x and
$x+y$
be two continuity points of F, where
$x,y>0$
. For
$n=1,2,\ldots$
, define
$s={n}/{x}$
and
$m=[{ny}/{x}]$
. Now letting
$n\to\infty$
and using Lemma 2.2, we have, from (2.19),
$\bar{F}(x+y) \leq {\text{e}}^{-r_F(0)x}\bar{F}(y)$
. The continuity points of F are dense in
$[0,\infty)$
, hence the above inequality holds for all
$x,y \geq 0$
. Therefore, F is NBUFR.
Remark 2.2. Note that the condition of Theorem 2.5 is equivalent to
Thus, (2.20) can also be treated as the necessary and sufficient condition for characterizing the NBUFR class.
At this juncture, it is appropriate to highlight an interesting application of our results. Theorem 2.5 in conjunction with Proposition 2.2 can be used to provide a short and elegant alternative proof of Block and Savits’ type characterization result for the NBUE class. The simplicity of this proof makes it rather attractive.
Theorem 2.6 ([Reference Block and Savits11].) For non-negative distributions F, F is NBUE if and only if
for all
$n>0$
and all
$s>0$
.
Proof. From [Reference Deshpande, Kochar and Singh14, Proposition 6.3],
$F_{(1)}$
is NBUFR
$\Leftrightarrow$
F is NBUE. Also, from (2.20) we have that
$F_{(1)}$
is NBUFR if and only if
Now it is easy to observe that
$r_{F_{(1)}}(0)={1}/{\mu}$
, where
$\mu$
is the mean of F. Proposition 2.2 together with (2.21) shows that, for all
$n>0$
,
$s>0$
,
i.e.
$[1+\mu s]\sum_{k=n+2}^\infty\alpha_{k}(s) \leq \mu s\sum_{k=n+1}^\infty\alpha_{k}(s)$
, i.e.
$\sum_{k=n+1}^\infty\alpha_{k}(s) \leq \alpha_{n}(s)\sum_{k=1}^\infty\alpha_{k}(s)$
. As
$\alpha_{0}(s)=1$
, we finally obtain that F is NBUE if and only if
$\sum_{k=n}^\infty\alpha_{k}(s) \leq \alpha_{n}(s)\sum_{k=0}^\infty\alpha_{k}(s)$
for all
$n>0$
and
$s>0$
.
2.6. The new better than used average in failure rate class
The new better than used average in failure rate (NBAFR) class was introduced in [Reference Loh32] as a generalization of the NBU class. A noticeable property of the NBAFR class is that it neither contains nor is contained in the NBUE family, which is an immediate generalization of the NBU class. First we look into the definition of the NBAFR class as proposed in [Reference Loh32].
Definition 2.9. A life distribution function F is said to belong to the NBAFR (NWAFR) class if
$t^{-1}\log\bar{F}(t) \leq (\geq)\,\lim_{s\to 0^+}s^{-1}\log\bar{F}(s)$
for all
$t>0$
.
The preservation properties of the NBAFR class under well-known reliability operations were also investigated in [Reference Loh32]. A non-monotonic generalization of this class was introduced in [Reference Majumder and Mitra34], which also investigated a moment bound, weak convergence, and related issues for this class. Motivated by [Reference Shokry, Ahmed, Rakha and Hewedi39, Theorem 2.1], we now give the following proposition that gives a moment inequality for the NBAFR class.
Proposition 2.4. Let F be NBAFR with
$\mu_{n}= \mathbb{E}(X^{n})$
for
$n>0$
. Define
$\lambda_n\;{:}{=}\;\mu_n/n!$
,
$n\geq0$
. Then
$\lambda_{n}\leq 1/[r_F(0)]^{n}$
for all
$n>0$
, where
$r_F(0)$
is defined in (2.17).
Proof. Note that a DF
$F\in \text{NBAFR}$
has finite moments of all orders [Reference Majumder and Mitra34, Corollary 3.1]. Therefore,
$\mu_{n}=\mathbb{E}(X^{n})$
is finite for all
$n>0$
. Again using the definition of
$r_F(0)$
given in (2.17) we obtain from Definition 2.9 that F is NBAFR if and only if
$\bar{F}(t) \leq {\text{e}}^{-r_F(0)t}$
for all
$t>0$
. Now, for all
$n\geq 0$
, we have
$\int_0^\infty x^n\bar{F}(x)\,{\text{d}} x \leq \int_0^\infty x^n{\text{e}}^{-r_F(0)x}\,{\text{d}} x$
. Some elementary algebra yields
$\lambda_{n+1} \leq 1/[r_F(0)]^{n+1}$
for all
$n\geq 0$
.
Remark 2.3 Note that
$\text{NBUFR} \subset \text{NBAFR}$
. Thus the condition
$\mu_{r+s+2}=\mathbb{E}(X^{r+s+2})<\infty$
of Lemma 2.7 (as well as in [Reference Shokry, Ahmed, Rakha and Hewedi39, Theorem 2.1]) can be dropped.
Now it would be appropriate to propose the characterization result for the NBAFR class.
Theorem 2.7. F is NBAFR if and only if
$\alpha_{n}(s) \leq {1}/{[1+r_F(0)/s]^n}$
for all
$n>0$
,
$s>0$
, where
$r_F(0)$
is defined in (2.17).
Proof.
Necessity. Fix
$s>0$
. Consider the DF G defined by
$\bar{G}(u) = {\text{e}}^{-su}\bar{F}(u)$
. Then G is NBAFR with
$r_G(0) = (r_F(0)+s)$
if and only if F is NBAFR with
$r_F(0)$
as defined in (2.17). As before, we have
Using Proposition 2.4 for the DF G we have
$\lambda_{n+1} \leq {1}/{[r_G(0)]^{n+1}}$
for all
$n\geq 0$
,
$s>0$
. Replacing
$r_G(0)$
, we get that if F is NBAFR then
$a_{n}(s) \leq {1}/{[r_F(0)+s]^{n+1}}$
for all
$n\geq0$
,
$s>0$
. Finally, using (1.1), the proof of the necessary part is completed.
Sufficiency. Now we have
Let
$x>0$
be a continuity point of F and define
$s=s(n,x)$
such that
$s={n}/{x}$
. Then letting
$n\to\infty$
in (2.22) and using Lemma 2.2 we obtain
$\bar{F}(x) \leq {\text{e}}^{-r_F(0)x}$
. Since the continuity points of F are dense, the above inequality holds for all
$x>0$
. Hence F is NBAFR.
Using Theorem 2.7, we now furnish an alternative proof for the BS type characterization for the HNBUE class as given in [Reference Klefsjö26].
Theorem 2.8. F is HNBUE if and only if
$\sum_{k=n+1}^\infty\alpha_k(s) \leq (\alpha(s)-1)(1-{1}/{\alpha(s)})^n$
for all
$n > 0$
and
$s>0$
, where
$\alpha(s) = \sum_{k=0}^\infty\alpha_k(s) = 1 + \mu s$
.
Proof. From [Reference Deshpande, Kochar and Singh14, Proposition 6.4] we have that
$F_{(1)}$
is NBAFR
$\Leftrightarrow F$
is HNBUE. Again, Theorem 2.7 gives that
$F_{(1)}$
is NBAFR if and only if
$\alpha'_n(s) \leq {1}/{[1+r_{F_{(1)}}(0)/s]^n}$
for all
$n > 0$
and
$s > 0$
. Note that
$r_{F_{(1)}}(0) = 1/\mu$
. Therefore, the proof is obtained simply by substituting
$\alpha'_n(s)$
from Proposition 2.2.
For the NBAFR class, another significant insight from considering the equilibrium distribution is that it enables us to obtain an upper bound for
$r_F(0)$
.
Proposition 2.5. If F is NBAFR, then the reciprocal of the mean of F gives an upper bound for
$r_F(0)$
.
Proof. We have
\begin{align*} \bar{F}_{(1)} (t) = \frac{1}{\mu}\int_t^\infty\bar{F}(u)\,{\text{d}} u & \leq \frac{1}{\mu}\int_t^\infty{\text{e}}^{-r_F(0)u}\,{\text{d}} u \qquad \text{(using the definition of NBAFR)} \\ & = \frac{1}{\mu}\frac{{\text{e}}^{-r_F(0)t}}{r_F(0)}. \end{align*}
Taking the
$\log$
of both sides, we have
$\log\bar{F}_{(1)}(t) \leq -r_F(0)t - \log\mu r_F(0)$
, i.e.
$-{\log\bar{F}_{(1)}(t)}/{t} \geq r_F(0) + {\log\mu r_F(0)}/{t}$
. Taking the limit
$t\to 0^+$
on both sides, we obtain
Suppose
$\mu r_F(0) > 1$
. In this case the right-hand side of (2.23) tends to
$+\infty$
as
$t\to 0^+$
, which implies that
$r_{F_{(1)}}(0) = +\infty$
; this is a contradiction since
$r_{F_{(1)}}(0)={1}/{\mu}$
. Therefore,
$\mu r_F(0) \not> 1$
. This observation leads to an upper bound of
$r_F(0)$
, which is
$ r_F(0)\leq{1}/{\mu}$
when F is NBAFR.
Note that our conclusion agrees with [Reference Majumder and Mitra34, Theorem 3.2(iii)] for
$k=1$
as we can look upon the NBAFR class as a special case of the NWBAFR class with change point
$t_0 = 0$
.
3. Characterizations of renewal classes
Consider a renewal process
$\{N(t)\colon t\geq0\}$
with interarrival times
$X_1,X_2,\ldots$
that are independent and with a common DF F and finite mean
$\mu>0$
. Here, N(t) is the number of failures by time t under an instantaneous replacement or perfect repair policy. Let the remaining life
$L(t)= \sum_{i=1}^{N(t)+1}X_i-t$
of a unit in use converge in distribution to a random variable
$X_{e}$
as
$t\to \infty$
. Then the distribution function of
$X_{e}$
, called the asymptotic remaining life, is given by
$F_{(1)}(t)= ({1}/{\mu})\int_0^t \bar{F}(u)\,\text{d}u$
for all
$t\geq0$
, which is nothing more than the equilibrium distribution of F. In reliability theory, it is standard practice to compare the lifetime of a new unit with the asymptotic remaining life of a used unit in terms of various reliability functions. In fact, such comparisons give rise to various ageing classes, termed renewal classes in the literature. These kind of ageing notions naturally arise in actuarial risk theory, maintenance models, replacement policies, and so on. In this section, we focus on those renewal classes that arise from two important reliability functions: the failure rate function and the MRL function.
3.1. The NBURFR and NBARFR classes
Proposed in [Reference Abouammoh and Ahmed2], the notions of new better (worse) than used renewal failure rate, NBURFR (NWURFR), and new better (worse) than average used renewal in failure rate, NBARFR (NWARFR), classes arise from the comparison between the failure rate function of a new unit and that of the asymptotic remaining life of a used unit. Let
$r_{F_{(1)}}(t)$
be the failure rate functions of
$F_{(1)}$
and
$r_F(0)$
as defined in (2.17). Then the formal definitions of these classes are as follows.
Definition 3.1. A life DF F is said to be NBURFR (NWURFR) if
$r_F(0) \leq (\geq)\, r_{F_{(1)}}(t)$
for all
$t \geq 0$
. Alternatively,
$\int_t^\infty\bar{F}(x)\,{\text{d}} x \leq (\geq)\, \bar{F}(t)/r_F(0)$
for all
$t \geq 0$
.
Definition 3.2 A life DF F is said to be NBARFR (NWARFR) if
$r_F(0) \leq (\geq)\, ({1}/{t}){}\int_0^t r_{F_{(1)}}(x)\,{\text{d}} x$
for all
$t \geq 0$
; equivalently,
$\int_t^\infty\bar{F}(x)\,{\text{d}} x \leq (\geq)\, \mu{\text{e}}^{-r_F(0)t}$
for all
$t \geq 0$
.
Theorem 3.1. F is NBURFR if and only if
$\sum_{k=n+1}^\infty\alpha_k(s) \leq ({1}/({r_F(0)/s}))\alpha_n(s)$
for all
$n\geq0$
and
$s>0$
, where
$r_F(0)$
is given by (2.17).
Proof. For the necessary part, we have, for
$n\geq0$
,
$s\geq 0$
,
\begin{align*} \sum_{i=n+1}^\infty\alpha_i(s) & = \sum_{i=n}^\infty s\int_0^\infty\frac{(sy)^i}{i!}{\text{e}}^{-sy}\bar{F}(y)\,{\text{d}} y \\ & = s^2\int_0^\infty\frac{{\text{e}}^{-st}(st)^n}{n!}\bigg(\int_t^\infty\bar{F}(y)\,{\text{d}} y\bigg)\,{\text{d}} t \qquad \text{(by virtue of Proposition 2.1)} \\ & \leq s^2\int_0^\infty\frac{{\text{e}}^{-st}(st)^{n-1}}{\Gamma(n)}\frac{\bar{F}(t)}{r_F(0)}\,{\text{d}} t = \frac{1}{r_F(0)/s}\alpha_n(s). \end{align*}
For the other part, given that
$n\geq0$
and
$s > 0$
,
Now let
$x>0$
be a continuity point of F and
$\lim_{n\to \infty}({n}/{s})=x$
. Then, appealing to the DCT and arguing as in the derivation of (2.6), we get
Therefore, letting
$n\to\infty$
in condition (3.1), we obtain
$\int_x^\infty\bar{F}(u)\,{\text{d}} u \leq \bar{F}(x)/r_F(0)$
. Note that this inequality holds for all
$x>0$
. Thus, F is NBURFR.
Theorem 3.2. F is NBARFR if and only if
$\sum_{k=n+1}^\infty\alpha_k(s)\leq\sum_{k=1}^\infty\alpha_k(s)/{[1+r_F(0)/s]^n}$
for all
$n\geq0$
and
$s>0$
, where
$r_F(0)$
is given by (2.17).
Proof. We omit the proof as it is quite similar to those of the previous theorems.
3.2. The NBRUE and HNBRUE classes
Exploiting the concept of renewal theory, two more classes, new better (worse) than renewal used in expectation, NBRUE (NWRUE), and harmonic new better (worse) than renewal used in expectation, HNBRUE (HNWRUE), were introduced in [Reference Abouammoh, Ahmed and Barry3], where the classes are analogues of NBUE and HNBUE in the sense that here the MRL function of a new unit is compared with that of the asymptotic remaining life of a used unit under a repeated renewal policy. Considering
$e_{F_{(1)}}(t)$
to be the MRL function of
$F_{(1)}$
, we have the following definitions for these classes.
Definition 3.3. A life DF F is said to be NBRUE (NWRUE) if
$e_{F_{(1)}}(t)\leq (\geq)\,\mu$
for
$t \geq 0$
, or, equivalently,
$\int_t^\infty\int_x^\infty\bar{F}(u)\,{\text{d}} u\,{\text{d}} x \leq (\geq)\,\mu\int_t^\infty\bar{F}(u)\,{\text{d}} u$
for all
$t \geq 0$
.
Definition 3.4. A life DF F is said to be HNBRUE (HNWRUE) if
for
$t \geq 0$
, which is equivalent to
$\int_t^\infty\int_x^\infty\bar{F}(u)\,{\text{d}} u\,{\text{d}} x \leq (\geq)\,{\text{e}}^{-t/\mu}\int_0^\infty\int_x^\infty \bar{F}(u)\,{\text{d}} u\,{\text{d}} x$
for all
$t \geq 0$
.
Several important issues such as the preservation properties of these classes under various reliability operations, interrelationships with the other classes, etc. are discussed in [Reference Bhattacharjee, Abouammoh, Ahmed and Barry9].
However, establishing BS type characterizations for these classes seems impossible using the techniques of [Reference Block and Savits11]. Even if some tedious manipulation yields the necessary part, the proof for the sufficiency part requires the boundedness of the function
$h(y)= \int_y^\infty\int_v^\infty\bar{F}(u)\,{\text{d}} u\,{\text{d}} v$
in order to apply the DCT or even [Reference Billingsley10, Theorem 5.2], which, in turn, results in imposing a strong assumption on the underlying DF F. We face a similar problem when dealing with the HNBRUE class.
Under the circumstances, we are forced to seek an alternative route for obtaining BS type characterization results in the context of these two classes. We noticed that these classes have connections with some other renewal classes in terms of the equilibrium distribution. The findings are recorded in the following lemma.
Lemma 3.1.
-
(i) F is NBRUE if and only if
$F_{(1)}$
is NBURFR. -
(ii) F is HNBRUE if and only if
$F_{(1)}$
is NBARFR.
Proof. Trivial from [Reference Abouammoh and Ahmed2, Theorem 2.3].
Lemma 3.1 makes a couple of useful additions to the list compiled in [Reference Deshpande, Kochar and Singh14, p. 756]; furthermore, it also enables us to adopt our new approach in proving the characterization results. We now have the following two theorems.
Theorem 3.3. F is NBRUE if and only if
$\sum_{i=n+1}^\infty\sum_{j=i+1}^\infty\alpha_j(s)\leq\sum_{k=1}^\infty\alpha_k(s)\sum_{i=n+1}^\infty\alpha_i(s)$
for all
$n \geq 0$
,
$s>0$
.
Proof. From Theorem 3.1 we have that
$F_{(1)}$
is NBURFR if and only if
Using this with Proposition 2.2 we obtain
\begin{equation*} \frac{1}{\mu s}\sum_{i=n+1}^\infty\sum_{j=i+1}^\infty\alpha_j(s) \leq \frac{1}{r_{F_{(1)}}(0)/s}\frac{1}{\mu s}\sum_{i=n+1}^\infty\alpha_i(s) \quad \text{for all } n\geq0,\,s>0. \end{equation*}
Recalling that
$r_{F_{(1)}}(0)={1}/{\mu}$
and using Lemma 3.1(i) we finally obtain that F is NBRUE if and only if
\begin{equation*} \sum_{i=n+1}^\infty\sum_{j=i+1}^\infty\alpha_j(s) \leq \sum_{k=1}^\infty\alpha_k(s)\sum_{i=n+1}^\infty\alpha_i(s) \quad \text{for all } n \geq 0, \,s>0.\end{equation*}
Theorem 3.4. F is HNBRUE if and only if
$$ \sum_{i=n+1}^\infty\sum_{j=i+1}^\infty\alpha_j(s) \leq \bigg(1-\frac{1}{\alpha(s)}\bigg)^n\sum_{i=0}^\infty\sum_{j=i+2}^\infty\alpha_j(s) $$
for all
$n \geq 0$
and
$s>0$
, where
$\alpha(s) = \sum_{k=0}^\infty \alpha_k(s)=1+\mu s$
.
Proof. For the sake of brevity, the proof is omitted.
As the HNBRUE family is analogous to HNBUE, it is natural to ask whether this class also possesses finite moments of all order. To answer this, we first look into the following proposition.
Proposition 3.1. For any life DF F with equilibrium DF
$F_{(1)}$
, the rth moment of
$F_{(1)}$
exists if and only if the
$(r+1)$
th moment of F exists.
Proof. Let
$\mu_{r}$
be the rth moment of F. Then we have
\begin{align*} \mu_{r+1} & = (r+1)\int_0^\infty x^r\bar{F}(x)\,{\text{d}} x \\ & = (r+1)\int_0^\infty\int_0^x ry^{r-1}\bar{F}(x)\,{\text{d}} y\,{\text{d}} x \\ & = (r+1)\int_0^\infty ry^{r-1}\int_y^\infty\bar{F}(x)\,{\text{d}} x\,{\text{d}} y \\ & = \mu(r+1)\int_0^\infty ry^{r-1}\bar{F}_{(1)}(y)\,{\text{d}} y = \mu(r+1)\mu^{(1)}_{r}, \end{align*}
where
$\mu^{(1)}_{r}$
is the rth moment of
$F_{(1)}$
. The result follows.
Corollary 3.1. Suppose
$F\in\mathcal{C}$
, where
$\mathcal{C}$
is any one of the classes {NBRUE}, {HNBRUE}, {NBURFR}, {NBARFR}. Then F has finite moments of all orders.
Proof. We first demonstrate a proof for the HNBRUE class. From the definition, it is clear that F is HNBRUE if and only if
where
$\bar{F}_{(2)}(t)$
is the survival function of the second derived distribution of
$F(t)$
. Now, this inequality shows that the second derived distribution of F possesses finite moments of all orders. Therefore, successive applications of Proposition 3.1 yield the desired result.
The proof for the NBARFR class is evident from Proposition 3.1.
Also, since NBRUE
$\subset$
HNBRUE and NBURFR
$\subset$
NBARFR, the results for the NBRUE and NBURFR classes are immediate.
4. Convolution
As pointed out in [Reference Block and Savits11], a nice application of the characterization results for the ageing classes using the Laplace transform is that the convolution of the continuous random variable can be derived from their discrete versions, simply by using the characterization result for the class and the convolution formula for
$\alpha_n(s)$
as suggested in [Reference Block and Savits11, Section 4], whenever the class is given to be closed under convolution. However, note that for the IMIT class, the convolution formula for
$\alpha_n(s)$
is not directly applicable. Therefore, as pointed out in [Reference Block and Savits11], by exploiting the properties of the Poisson process we can also obtain a convolution formula for
$\bar{\alpha}_n(s)$
. Considering X and Y to be two independent non-negative random variables, we have
and a similar expression can also be obtained for Y. Then, from (2.13), we have
Hence, using the fact that a Poisson process has stationary and independent increments, we obtain the following convolution formula for
$\bar{\alpha}_{n+1}(s)$
:
It is worth noting that if the convolution of the discrete IMIT class can be proved then, by virtue of Theorem 2.2 and this formula, we can directly conclude that closure under convolution also holds for this class in the continuous case.
5. Discussion
The characterizations for the dual classes such as IMTTF, NWUFR, NWAFR, IMRLHA, etc. can be obtained by some suitable and obvious modifications of the theorems and their corresponding proofs. Note that the finiteness of the moments of all orders for the NBUFR and NBAFR classes can also be shown from Proposition 3.1 and the equivalence relation of the NBAFR class with HNBUE via the equilibrium distribution. Further, it is worth noting that, for a delayed renewal process
$\{N_D(t)\colon t\geq 0\}$
with independent interarrival times
$X_1', X_2', X_3',\ldots$
where
$X_2', X_3',\ldots$
follow some common DF F and
$X_1'$
follows
$F_{(1)}$
(the first derived distribution of F), the sum
$\sum_{k=1}^\infty \alpha_k(s)$
is, in fact, the reciprocal of the Laplace–Stieltjes transform of the expected number of renewals in time [0, t) [Reference Karlin and Taylor20].
An important direction of this work is to uncover the significance of the equilibrium distribution in the context of characterizing the reliability classes. A few important points should be noted in this regard. Adopting our suggested approach is more sensible when the characterization of the equilibrium distribution is given and the concern is to characterize the original DF F. Apparently, it is also not always possible to find a related equivalence with the derived distribution for each and every reliability class. Despite these limitations, we recommend the use of this methodology whenever possible, as it typically results in short and elegant proofs.
Acknowledgements
The authors would like to convey their sincere thanks to the anonymous referees and the Editor for their careful reading of the manuscript and valuable comments and suggestions which have led to a substantial improvement of an earlier version of this paper. The second author is grateful to Dr. Ruhul Ali Khan for drawing his attention to [Reference Block and Savits11].
Funding information
The first author is grateful to the Council of Scientific & Industrial Research (CSIR), Ministry of Science & Technology, Govt. of India for financial support.
Competing interests
There were no competing interests to declare which arose during the preparation or publication process of this article.
