A Class of Truncated Binomial Lifetime Distributions

In this paper, a new lifetime class with decreasing failure rate is introduced by compounding truncated binomial distribution with any proper continuous lifetime distribution. The properties of the proposed class are discussed, including a formal proof of its probability density function, distribution function and explicit algebraic formulae for its reliability and failure rate functions. A simple EM-type algorithm for iteratively computing maximum likelihood estimates is presented. The Fisher information matrix is derived in order to obtain the asymptotic covariance matrix. This new class of distributions generalizes several distributions which have been introduced and studied in the literature.


Introduction
The study of life length of organisms, structures, materials, etc., is very important in the biological and engineering sciences.A substantial part of such study is devoted to modeling the lifetime data by a failure distribution.The exponential, Rayleigh, and Weibull distributions are the most commonly used distributions in reliability and life testing.These distributions have several desirable properties and nice physical interpretations.Unfortunately, however, the exponential and Rayleigh distributions have constant and increasing failure rates, respectively.The Weibull distribution generalizes both these distributions which may have increasing, constant, or decreasing failure rates.
Situations where the failure rate function decreases with time have been reported by several authors.Indicative examples are business mortality (Lomax [5]), failure in the air-conditioning equipment of a fleet of Boeing 720 aircrafts or in semiconductors from various lots combined (Proschan [6]) and the life of integrated circuit modules (Saunders and Myhre [7]).In general, a population is expected to exhibit decreasing failure rate (DFR) when its behavior over time is characterized by "work hardening" (in engineering terms) or "immunity" (in biological terms); sometimes the broader term "infant mortality" is used to denote the DFR phenomenon.The resulting improvement of reliability with time might have occurred by means of actual physical changes that caused self-improvement or simply it might have been due to population heterogeneity.Indeed, (Proschan [6]) provided that the DFR property is inherent to mixtures of distributions with constant failure rate (see also (McNolty et al. [8]) for other properties of exponential mixtures and (Gleser [9]) demonstrated the converse for any gamma distribution with shape parameter less than one.In addition, (Gurland and Sethuramm [10]) give examples illustrating that such results may hold for mixtures of distributions with rapidly increasing failure rate.A mixture of truncated geometric distribution and exponential with DFR was introduced by (Adamidis and Loukas [11]).The exponential-Poisson (EP) distribution was proposed by (Kus [2]), and generalized by (Hemmati et al. [12]) using Wiebull distribution and the exponential-logarithmic distribution was discussed by (Tahmasbi and Rezaei [3]).(Silvia et al. [13]) did a new distribution with decreasing, increasing and upside down bathtub failure rate.A twoparameter distribution family with decreasing failure rate arising by mixing power-series distribution has been introduced by (Chahkandi and Ganjali [4]).A Weibull power series class of distributions with Poisson was pre-sented by (Morais and Barreto-Souza [14]).(Morais [15]) in a master degree thesis presented a class of generalized Beta distributions, Pareto power series and Weibull power series.Lately, (Alkarni and Oraby [16]) and (Alkarni [17]) obtained a class of truncated Poisson and logarithmic distributions with any continuous lifetime distribution.
A further exponentiated type distribution has been introduced and studied in the literature.The exponential-Weibull (EW) distribution was proposed by (Mudholkar and Srivastava [18]) to extend the GE distribution.This distribution was also studied by (Mudholkar et al. [19]), (Mudholkar and Hutson [20]) and (Nassar and Eissa [21]).(Nadarajah and Kotz [22]) introduced four more exponentiated type distributions: the exponentiated gamma, exponentiated Weibull, exponentiated Gumbel and exponentiated Fréchet distributions by generalizing the gamma, Weibull, Gumbel and Fréchet distributions in the same way that the GE distribution extends the exponential distribution.(Barreto-Souza and Cribari-Neto [23]) introduced the generalized exponential-Poisson distribution which extends the exponential-Poisson distribution in the same way that the GE distribution extends the exponential distribution.
In this paper we generalize the work of (Bakouchetal.[24]) to a class of several lifetime continuous distributions and hence any mixture of continuous lifetime with truncated binomial distribution such as exponential, Weibull, pareto becomes a special case of this class.This paper is organized as follow.In Section 2, the new class of binomial lifetime distributions with its probability and distribution functions is introduced.In Section 3, the corresponding survival and hazard rate functions with some of their properties are derived.In Section 4, maximum likelihood estimate of the unknown parameters is obtained based on a random sample via EM algorithm.In Section 5, the entropy for the binomial lifetime distributions class is discussed.

The Class
, , N T T  be independent and identically distributed (iid) random variables with probability density function (pdf) given by Here, is a zero truncated binomial random variable with probability mass function given by then, the pdf of the random variable X is obtained as And hence the cumulative distribution function The proof of the results in ( 1) and ( 2) are presented in the following theorem.
and is a zero truncated binomial random variable with probability mass function and hence the joint pdf of X and is obtained as The marginal pdf and cdf of are given by We denote a random variable X with pdf and cdf (1) and (2) by . This new class of distributions generalizes several distributions which have been introduced and studied in the literature.For instance using the probability density and its distribution function of exponential distribution in (1), we obtain the binomial exponential distribution (Bakouch et al. [24]) and using Wiebull probability density and its distribution function gives binomial Wiebull distribution (Morais and Barreto-Souza ( [14]).The model is obtained under the concept of population heterogeneity (through the process of compounding).An interpretation of the proposed model is as follows: a situation where failure (of a device for example) occurs due to the presence of an unknown number, of initial defects of same kind (a number of semiconductors from a defective lot, for example).The Ts represent their lifetimes and each defect can be detected only after causing failure, in which case it is repaired perfectly (Adamidis and Loukas [1]).Then, the distributional assumptions given earlier lead to any of the BL distributions for modeling the time to the first failure N X .Table 1 shows the probability function and the distribution function for some lifetime distributions.
Some lifetime distributions are excluded from this table such as Gamma and lognormal distribution stable since they do not have closed forms.They still can be applied in this class numerically.
The quantile q th q x of the BL distribution, the inverse of the distribution function   X q F x  q is the same as the inverse of the distribution function for any continuous lifetime with distribution function   .T F .

Survival and Hazard Functions
Since the BL is not part of the exponential family, there is no simple form for moments see for instant (Kus [2]) for the exponential case.Survival function (also known reliability function) (sf) and hazard function (known as failure rate function) (hf) for the BL class are given in the following theorem.
and is a zero truncated binomial random variable with probability mass function  with where and are independent.If ,0 and Respectively.
Proof: Using ( 1) and ( 2), survival function (also known reliability function) and hazard function (known as failure rate function) for the BL class are given respectively by Table 2 summarizes the survival functions and hazard rate functions for some distributions of the class.
The hazard function for BL class is decreasing because the DFR property follows from the result of (Barlow et al. [25]) on mixtures.

Estimation
In what follows, we discuss the estimation of the BL class parameters.Let, and the associated score function is given by ; and For 1, , .i k   The maximum likelihood estimates(MLE) of  , say  , is obtained by solving the nonlinear system The solution of this nonlinear system of equations has not a closed form, but can be found numerically by using software such as MATHEMATICA, MAPLE, Ox and R.
For interval estimation and hypothesis tests on the model parameters, we require the information matrix.The where the elements of   m I  are the second partial derivatives of ( 5), ( 6) and (7).Under the regular conditions stated in (Cox and Hinkley [26]), that are fulfilled for our model whenever the parameters are in the interior of the parameter space, we have that the asymptotic distribu-  is the unit information matrix.
1 e 1 1 e ln 1 1 e Based on the underline distribution, the maximum likelihood estimation of the parameters can be found analytically using an EM algorithm.Newton-Raphson algorithm is one of the standard methods to determine the MLEs of the parameters.To employ the algorithm, second derivatives of the log-likelihood are required for all iteration.EM algorithm is a very powerful tool in handling the incomplete data problem (Dempster et al., [27]; McLachlan and Krishnan, [28]).It is an iterative method by repeatedly replacing the missing data with estimated values and updating the parameter estimates.It is especially useful if the complete data set is easy to analyze.As pointed out by (Little and Rubin [29]), the EM algorithm will converge reliably but rather slowly (as compared to the Newton-Raphson method) when the amount of information in the missing data is relatively large.Recently, EM algorithm has been used by several authors such as (Adamidis and Loukas [1]), (Adamidis [30]), (Ng et al. [31]), (Karlis [32]) and (Adamidis et al. [11]).
To estimate , EM algorithm is a recurrent method such that each step consists of an estimate of the expected value of a hypothetical random variable and later maximizes the log-likelihood of the complete data.Let the complete data be is the likelihood of interest.Then, we define a hypothetical complete-data distribution for each 1 , , , . Thus, it is straight forward to verify that the E-step of an EM cycle requires the computation of the conditional expectation . The EM cycle is completed with M-step, which is complete data maximum likelihood over

Entropy for the Class
If X is a random variable having an absolutely continuous cumulative distribution function Statistical entropy is a probabilistic measure of uncertainty or ignorance about the outcome of a random experiment, and is a measure of a reduction in that uncertainty.Since Shannon's [33] pioneering work on the mathematical theory of communication, entropy (8) has been used as a major tool in information theory and in almost every branch of science and engineering.Numerous entropy and information indices, among which there is the Renyi entropy, have been developed and used in various disciplines and contexts.Information theoretic principles and methods have become integral parts of probability and statistics and have been applied in various branches of statistics and related fields.
the basic uncertainty measure for distribution F (called the entropy of F ) is defined as