arxiv: v1 [stat.me] 5 Apr 2012

Transcription

1 The compound class of etended Weibull power series distributions Rodrigo B. Silva a,1, Marcelo B. Pereira a,2, Cícero R. B. Dias a,3, Gauss M. Cordeiro a,4 a Universidade Federal de Pernambuco Departamento de Estatística, Cidade Universitária, Recife, PE, Brazil arxiv: v1 stat.me] 5 Apr 212 Abstract In this paper, we introduce a new class of distributions which is obtained by compounding the etended Weibull and power series distributions. The compounding procedure follows the same set-up carried out by Adamidis and Loukas (1998) and defines at least new 68 sub-models. This class includes some well-known miing distributions, such as the Weibull power series (Morais and Barreto-Souza, 21) and eponential power series (Chahkandi and Ganjali, 29) distributions. Some mathematical properties of the new class are studied including moments and generating function. We provide the density function of the order statistics and obtain their moments. The method of maimum likelihood is used for estimating the model parameters and an EM algorithm is proposed for computing the estimates. Special distributions are investigated in some detail. An application to a real data set is given to show the fleibility and potentiality of the new class of distributions. Keywords: EM algorithm, Etended Weibull distribution, Etended Weibull power series distribution, Order statistic, Power series distribution. 1. Introduction The modeling and analysis of lifetimes is an important aspect of statistical work in a wide variety of scientific and technological fields. Several distributions have been proposed in the literature to model lifetime data by compounding some useful lifetime distributions. Adamidis and Loukas (1998) introduced a two-parameter eponential-geometric (EG) distribution by compounding an eponential distribution with a geometric distribution. In the same way, the eponential Poisson (EP) and eponential logarithmic (EL) distributions were introduced and studied by Kus (27) and Tahmasbi and Rezaei (28), respectively. Recently, Chahkandi and Ganjali (29) proposed the eponential power series (EPS) family of distributions, which contains as special cases these distributions. Barreto-Souza et al. (21) and Lu and Shi (211) introduced the Weibull-geometric (WG) and Weibull-Poisson (WP) distributions which naturally etend the EG and EP distributions, respectively. In a very recent paper, Morais and Barreto-Souza (211) defined the Weibull power series (WPS) class of distributions which contains the EPS distributions as sub-models. The WPS distributions can have an increasing, decreasing and upside down bathtub failure rate function. Now, consider the class of etended Weibull (EW) distributions, as proposed by Gurvich et al. (1997), having the cumulative distribution function (cdf) G(; α, ξ) = 1 e α H(; ξ), >, α >, (1) where H(; ξ) is a non-negative monotonically increasing function which depends on a parameter vector ξ. The corresponding probability density function (pdf) is given by where h(; ξ) is the derivative of H(; ξ). g(; α, ξ) = α h(; ξ) e α H(; ξ), >, α >, (2) 1 rodrigobs29@gmail.com 2 m.p.bourguignon@gmail.com 3 cicerorafael@gmail.com 4 Corresponding author. gauss@de.ufpe.br Preprint submitted to Computational Statistics & Data Analysis June 5, 218

2 Note that many well-known models are special cases of equation (1) such as: (i) H(; ξ) = gives the eponential distribution; (ii) H(; ξ) = 2 yields the Rayleigh distribution (Burr type-x distribution); (iii) H(; ξ) = log(/k) leads to the Pareto distribution; (iv) H(; ξ) = β 1 ep(β) 1] gives the Gompertz distribution. In this article, we define the etended Weibull power series (EWPS) class of univariate distributions obtained by compounding the etended Weibull and power series distributions. The compounding procedure follows the key idea of Adamidis and Loukas (1998) or, more generally, by Chahkandi and Ganjali (29) and Morais and Barreto-Souza et al. (211). The new class of distributions contains as special models the WPS distributions, which in turn etends the EPS distributions and defines at least new 68 (17 4) sub-models as special cases. The hazard function of our class can be decreasing, increasing, bathtub and upside down bathtub. We are motivated to introduce the EWPS distributions because of the wide usage of the general class of Weibull distributions and the fact that the current generalization provides means of its continuous etension to still more comple situations. This paper is organized as follows. In Section 2, we define the EWPS class of distributions and demonstrate that there are many eisting models which can be deduced as special cases of the proposed unified model. In Section 3, we provide the density, survival and hazard rate functions and derive some useful epansions. In Section 4, we obtain its quantiles, ordinary and incomplete moments. Further, the order statistics are discussed and their moments are determined. Section 5 deals with reliability and average lifetime. Estimation of the parameters by maimum likelihood using an EM algorithm and large sample inference are investigated in Section 6. In Section 7, we present suitable constraints leading to the maimum entropy characterization of the new class. Three special cases of the proposed class are studied in Section 8. In Section 9, we provide an application to a real data set. The paper is concluded in Section The new class Our class can be derived as follows. Given N, let X 1,..., X N be independent and identically distributed (iid) random variables following (1). Here, N is a discrete random variable following a power series distribution (truncated at zero) with probability mass function p n = P (N = n) = a n θ n, n = 1, 2,..., (3) where a n depends only on n, = a n θ n and θ > is such that is finite. Table 1 summarizes some power series distributions (truncated at zero) defined according to (3) such as the Poisson, logarithmic, geometric and binomial distributions. Let X (1) = min {X i } N. The conditional cumulative distribution of X (1) N = n is given by G X(1) N=n() = 1 e nαh(;ξ), i.e., X (1) N = n follows a general class of distributions (1) with parameters nα and ξ based on the same H(; ξ) function. Hence, we obtain P (X (1), N = n) = a n θ n 1 e nαh(;ξ)], >, n 1. The EWPS class of distributions can then be defined by the marginal cdf of X (1) : F (; θ, α, ξ) = 1 C(θ e αh(;ξ) ), >. (4) 2

3 Distribution a n C (θ) C (θ) 1 Θ Poisson n! 1 e θ 1 e θ e θ log(θ + 1) θ (, ) Logarithmic n 1 log(1 θ) (1 θ) 1 (1 θ) 2 1 e θ θ (, 1) Geometric 1 θ(1 θ) 1 (1 θ) 2 2(1 θ) 3 θ(θ + 1) 1 θ (, 1) ( Binomial m n) (θ + 1) m 1 m(θ + 1) m 1 m(m 1) (θ+1) 2 m (θ 1) 1/m 1 θ (, 1) Table 1: Useful quantities for some power series distributions. The random variable X following (4) with parameters θ and α and the vector ξ of parameters is denoted by X EWPS(θ, α, ξ). Equation (4) etends several distributions which have been studied in the literature. The EG distribution (Adamidis and Loukas, 1998) is obtained by taking H(; ξ) = and = θ(1 θ) 1 with θ (, 1). Further, for H(; ξ) =, we obtain the EP (Kus, 27) and EL (Tahmasbi and Rezaei, 28) distributions by taking = e θ 1, θ >, and = log(1 θ), θ (, 1), respectively. In the same way, for H(; ξ) = γ, we obtain the WG (Barreto-Souza et al., 29) and WP (Lu and Shi, 211) distributions. The EPS distributions are obtained from (4) by miing H(; ξ) = with any listed in Table 1 (see Chahkandi and Ganjali, 29). Finally, we obtain the WPS distributions from (4) by compounding H(; ξ) = γ with any in Table 1 (see Morais and Barreto-Souza, 211). Table 2 displays some useful quantities and respective parameter vectors for each particular distribution. 3. Density, survival and hazard functions The density function associated to (4) is given by f(; θ, α, ξ) = θ α h(; ξ) e αh(; ξ) C (θ e αh(; ξ) ), >. (5) Proposition 1. The EW class of distributions with parameters cα and ξ is a limiting special case of the EWPS class of distributions when θ +, where c = min {n N : a n > }. Proof. This proof uses a similar argument to that found in Morais and Barreto-Souza (211). Define c = min {n N : a n > }. We have for >. lim F () = 1 lim θ + θ + = 1 lim θ + n=c = 1 e cαh(;ξ), a n ( θ e αh(;ξ)) n a n θ n n=c e cαh(;ξ) + a 1 c 1 + a 1 c n=c+1 n=c+1 a n θ n c e nαh(;ξ) a n θ n c We now provide an interesting epansion for (5). We have C (θ) = n a n θ n 1. By using this result in (5), it follows that f(; θ, α, ξ) = p n g(; nα, ξ), (6) where g(; nα, ξ) is given by (2). Based on equation (6), we obtain F (; θ, α, ξ) = 1 3 p n e nαh(; ξ).

4 Distribution H(; ξ) h(; ξ) α ξ References Eponential ( ) 1 α Johnson et al. (1994) Pareto ( k) log(/k) 1/ α k Johnson et al. (1994) Rayleigh ( ) 2 2 α Rayleigh (188) Weibull ( ) γ γ γ 1 α γ Johnson et al. (1994) Modified Weibull ( ) γ ep(λ) γ 1 ep(λ)(γ + λ) α γ, λ] Lai et al. (23) Weibull etension ( ) λep(/λ) β 1] β ep(/λ) β (/λ) β 1 α γ, λ, β] Xie et al. (22) Log-Weibull ( < < ) ep( µ)/σ] (1/σ) ep( µ)/σ] 1 µ, σ] White (1969) Phani ( < µ < < σ < ) ( µ)/(σ )] β β( µ)/(σ )] β 1 (σ µ)/(σ t) 2 ] α µ, σ, β] Phani (1987) Weibull Kies ( < µ < < σ < ) ( µ) β 1/(σ ) β 2 ( µ) β 1 1 (σ ) β2 1 β 1(σ ) + β2( µ)] α µ, σ, β1, β2] Kies (1958) Additive Weibull ( ) (/β1) α 1 + (/β2) α 2(α1/β1)(/β1) α (α 2/β2)(/β2) α α1, α2, β1, β2] Xie and Lai (1995) Traditional Weibull ( ) b ep(c d 1)] b b 1 ep(c d ) 1] + cd b+d 1 ep(c d ) α b, c, d] Nadarajah and Kotz (25) Gen. power Weibull ( ) 1 + (/β) α 1] θ 1 (θα1/β)1 + (/β) α 1] θ 1 (/β) α 1 1 α1, β, θ] Nikulin and Haghighi (26) Fleible Weibull etension( ) ep(α1 β/) ep(α1 β/)(α1 + β/ 2 ) 1 α1, β] Bebbington et al. (27) Gompertz ( ) β 1 ep(β) 1] ep(β) α β Gompertz (1825) Eponential power ( ) ep(λ) β ] 1 βλ ep(λ) β ](λ) β 1 1 λ, β] Smith and Bain (1975) Chen ( ) ep( b ) 1 b b 1 ep( b ) α b Chen (2) Pham ( ) (a ) β 1 β(a ) β log(a) 1 a, β] Pham (22) Table 2: Special distributions and related H(; ξ) and h(; ξ) functions. 4

5 Hence, the EWPS density function is an infinite miture of EW densities. So, some mathematical quantities (such as ordinary and incomplete moments, generating function and mean deviations) of the EWPS distributions can be obtained by knowing those quantities for the baseline density function g(; nα, ξ). The EWPS survival function is given by S(; θ, α, ξ) = C(θ e αh(; ξ) ) (7) and the corresponding hazard rate function becomes τ(; θ, α, ξ) = θα h(; ξ) e nαh(;ξ) C (θ e αh(; ξ) ) C(θ e αh(; ξ) ). (8) 4. Quantiles, moments and order statistics The EWPS distributions are easily simulated from (4) as follows: if U has a uniform U(, 1) distribution, then the solution of the nonlinear equation { X = H 1 1 C 1 ]} α log ((1 U)) θ has the EWPS(θ, α, ξ) distribution, where H 1 ( ) and C 1 ( ) are the inverse functions of H( ) and C( ), respectively. To simulate data from this nonlinear equation, we can use the matri programming language O through SolveNLE subroutine (see Doornik, 27). We now derive a general epression for the rth raw moment of X, which may be determined by using (6) and the monotone convergence theorem. So, for r N, we obtain E(X r ) = p n E(Z r ), (9) where Z is a random variable with pdf g(z; nα, ξ). The incomplete moments and moment generating function (mgf) follow by using (6) and the monotone convergence theorem: and I X (y) = M X (t) = y r f()d = p n E ( e tz). p n I Z (y) where Z is defined as before. Order statistics are among the most fundamental tools in non-parametric statistics and inference. They enter in the problems of estimation and hypothesis tests in a variety of ways. Therefore, we now discuss some properties of the order statistics for the proposed class of distributions. The pdf f i:m () of the ith order statistic for a random sample X 1,..., X m from the EWPS distribution is given by m! f i:m () = f(; θ, α, ξ) 1 C(θ e αh(; ξ) ] i 1 ) C(θe αh(; ξ) ] m i ), >, (1) (i 1)!(m i)! where f(; θ, α, ξ) is the pdf given by (5). By using the binomial epansion, we can write (1) as i 1 m! f i:m () = (i 1)!(m i)! f(; θ, α, ξ) ( ) i 1 ( 1) j S(; θ, α, ξ) m+j i, (11) j 5 j=

6 where S(; θ, α, ξ) is given by (7). The corresponding cumulative function is F i:m () = j= k=i m ( k ( 1) j j An alternative form for (1) can be obtained from (6) as f i:m () = m! (i 1)!(m i)! ) ( ) m S(; θ, α, ξ) m+j k. k i 1 ω j p n g(; nα, ξ)s(; θ, α, ξ) m+j 1, (12) j= where ω j = ( 1) j( ) i 1 j. So, the sth raw moment Xi:m comes immediately from the above equation E (X s i:m) = where Z EW(nα, ξ) is defined before. m! (i 1)!(m i)! i 1 ω j p n E Z s S(Z) m+j i], (13) j= 5. Reliability and average lifetime In the contet of reliability, the stress-strength model describes the life of a component which has a random strength X subjected to a random stress Y. The component fails at the instant that the stress applied to it eceeds the strength, and the component will function satisfactorily whenever X > Y. Hence, R = P(X > Y ) is a measure of component reliability. It has many applications, especially in engineering concepts. The algebraic form for R has been worked out for the majority of the well-known distributions. Here, we obtain the form for the reliability R when X and Y are independent random variables having the same EWPS distribution. The quantity R can be epressed as R = Substituting (4) and (5) into equation (14), we obtain R = θ α h(; ξ) e αh(; ξ) C (θe αh(; ξ) ) = 1 p n g(; nα, ξ)s(; θ, α, ξ)d, where the integral can be calculated from the baseline EW distribution. The average lifetime is given by t m = f(; θ, α, ξ)f (; θ, α, ξ)d. (14) p n e nαh(; ξ) d. ] 1 C(θe αh(;ξ) ) d Given that there was no failure prior to, the residual life is the period from time until the time of failure. The mean residual lifetime can be epressed as m( ; θ, α, ξ) = Pr(X > )] 1 = S( )] 1 p n y f( + y; θ, α, ξ)dy y g( + y; nα, ξ)dy. The last integral can be computed from the baseline EW distribution. Furthermore, m( ; θ, α, ξ) E(X) as. 6

7 6. Maimum likelihood estimation 6.1. Preliminaries Here, we determine the maimum likelihood estimates (MLEs) of the parameters of the EWPS class of distributions from complete samples only. Let X 1,..., X n be a random sample with observed values 1,..., n from an EWPS distribution with parameters θ, α and ξ. Let Θ = (θ, α, ξ) be the p 1 parameter vector. The total log-likelihood function is given by l n = l n (; Θ) = n log θ + log α log ] α + H( i ; ξ) + log h( i ; ξ) log C (θ e αh(i; ξ) ). (15) The log-likelihood can be maimized either directly by using the SAS (PROC NLMIXED) or the O program (sub-routine MaBFGS) (see Doornik, 27) or by solving the nonlinear likelihood equations obtained by differentiating (15). The components of the score function U n (Θ) = ( l n / θ, l n / α, l n / ξ) are l n α = n n α H( i ; ξ) θ l n θ = n θ nc (θ) + n H( i ; ξ)e αh(i; ξ) C (θ e αh(i; ξ) ) C (θ e αh(i; ξ) ), e αh(i; ξ) C (θ e αh(i; ξ) ) C (θ e αh(i; ξ) ) and l n ξ k = log h( i ; ξ) ξ k α H( i ; ξ) ξ k 1 + θe αh(i; ξ) C (θ e αh(i; ξ) ] ). C (θ e αh(i; ξ) ) For interval estimation on the model parameters, we require the observed information matri U θθ U θα U θξ J n (Θ) = U αθ U αα Uαξ, U θξ U αξ U ξξ whose elements are listed in Appendi A. Let Θ be the MLE of Θ. Under standard regular conditions stated in Co and Hinkley (1974) that are fulfilled for our model whenever the parameters are in the interior of the parameter space, we have that the asymptotic distribution of ) n ( Θ Θ is multivariate normal N p (, K(Θ) 1 ), where K(Θ) = lim n J n (Θ) is the unit information matri and p is the number of parameters of the compounded distribution The EM algorithm Here, we propose an EM algorithm (Dempster et al., 1977) to estimate Θ. The EM algorithm is a recurrent method such that each step consists of an estimate of the epected value of a hypothetical random variable and then maimizes the log-likelihood for the complete data. Let the complete-data be X 1,..., X n with observed values 1,..., n and the hypothetical random variables Z 1,..., Z n. The joint probability function is such that the marginal density of X 1,..., X n is the likelihood of interest. Then, we define a hypothetical complete-data distribution for each (X i, Z i ), i = 1,..., n, with a joint probability function in the form g(, z; Θ) = α z a z θ z h(; ξ) e αzh(; ξ), 7

8 where θ and α are positive, > and z N. Under this formulation, the E-step of an EM cycle requires the epectation of Z X; Θ (r) = (θ (r), α (r), ξ (r) ) as the current estimate (in the rth iteration) of Θ. The probability function of Z given X, say g(z ), is given by and its epected value is g(z ) = z a z θ θ 1 ξ) C (θe αh(i; ξ) e α(z 1)H(i; ) E(Z X) = 1 + θe αh(; ξ) C (θ e αh(; ξ) ) C (θ e αh(; ξ) ). The EM cycle is completed with the M-step by using the maimum likelihood estimation over Θ, where the missing Z s are replaced by their conditional epectations given before. The log-likelihood for the complete-data is l n( 1,..., n ; z 1,..., z n ; α, θ, ξ) n log α + log θ α z i + log h( i ; ξ) z i H( i ; ξ) n log. So, the components of the score function U n(θ) = ( l n/ θ, l n/ α, l n/ ξ) are ln θ = n n θ z i n C (θ), ln α = n n α z i H( i ; ξ) and l n ξ k = log h( i ; ξ) ξ k α z i H( i ; ξ) ξ k. From a nonlinear system of equations U n( Θ) =, we obtain the iterative procedure of the EM algorithm and ˆα (t+1) n = n z(t) i H( i ; ξ (t) ), ˆθ(t+1) = C(ˆθ(t+1) ) 1 C (ˆθ (t+1) ) n log h( i ; ˆξ (t+1) ) ξ k ˆα (t) z (t) H( i ; ˆξ (t+1) ) i =, ξ k where ˆθ (t+1), ˆα (t+1) and ˆξ (t+1) are obtained numerically. Here, for i = 1,..., n, we have z (t) i = 1 + ˆθ (t) e ˆα(t) H( i; ˆξ (t) ) C (ˆθ (t) e ˆα (t) H( i; ˆξ (t)) ) C (ˆθ (t) e ˆα(t) H( i; ˆξ (t) ) ). Note that, in each step, θ, α and ξ are estimated independently. The EWPS distributions can be very useful in modeling lifetime data and practitioners may be interested in fitting one of our models. z (t) i 7. Maimum entropy identification Shannon (1948) introduced the probabilistic definition of entropy which is closely connected with the definition of entropy in statistical mechanics. Let X be a random variable of a continuous distribution with density f. Then, the Shannon entropy of X is defined by H Sh (f) = f(; θ, α, ξ) log f(; θ, α, ξ)] d. (16) R 8

9 Jaynes (1957) introduced one of the most powerful techniques employed in the field of probability and statistics called the maimum entropy method. This method is closely related to the Shannon entropy and considers a class of density functions F = {f(; θ, α, ξ) : E f (T i (X)) = α i, i =,..., m}, (17) where T i (X), i = 1,..., m, are absolutely integrable functions with respect to f, and T (X) = a = 1. In the continuous case, the maimum entropy principle suggests deriving the unknown density function of the random variable X by the model that maimizes the Shannon entropy in (16), subject to the information constraints defined in the class F. Shore and Johnson (198) treated aiomatically the maimum entropy method. This method has been successfully applied in a wide variety of fields and has also been used for the characterization of several standard probability distributions; see, for eample, Kapur (1989), Soofi (2) and Zografos and Balakrishnan (29). The maimum entropy distribution is the density of the class F, denoted by f ME, which is obtained as the solution of the optimization problem f ME (; θ, α, ξ) = arg ma f F H Sh. Jaynes (1957, p. 623) states that the maimum entropy distribution f ME, obtained by the constrained maimization problem described above, is the only unbiased assignment we can make; to use any other would amount to arbitrary assumption of information which by hypothesis we do not have. It is the distribution which should not incorporate additional eterior information other than which is specified by the constraints. We now derive suitable constraints in order to provide a maimum entropy characterization for our class of distributions defined by (4). For this purpose, the net result plays an important role. Proposition 2. Let X be a random variable with pdf given by (5). Then, C1. E log(c (θ e αh(x; ξ) )) ] = θ E C (θ e αh(y ; ξ) ) log(c (θ e αh(y ; ξ) )) ] ; C2. E log(h(x; ξ))] = θ E C (θ e αh(y ; ξ) ) log(h(y ; ξ)) ] ; C3. E H(X; ξ)] = θ E C (θ e αh(y ; ξ) )H(Y ; ξ) ], where Y follows the EW distribution with density function (2). Proof. The constraints C1, C2 and C3 are easily obtained and therefore their demonstrations are omitted. The net proposition reveals that the EWPS distribution has maimum entropy in the class of all probability distributions specified by the constraints stated in the previous proposition. Proposition 3. The pdf f of a random variable X, given by (5), is the unique solution of the optimization problem f(; θ, α, ξ) = arg ma H Sh, h under the constraints C1, C2 and C3 presented in the Proposition 2. Proof. Let τ be a pdf which satisfies the constraints C1, C2 and C3. The Kullback-Leibler divergence between τ and f is ( ) τ(; θ, α, ξ) D(τ, f) = τ(; θ, α, ξ) log d. f(; θ, α, ξ) R Following Cover and Thomas (1991), we obtain D(τ, f) = τ(; θ, α, ξ) log τ(; θ, α, ξ)] d τ(; θ, α, ξ) log f(; θ, α, ξ)] d R R = H Sh (τ; θ, α, ξ) τ(; θ, α, ξ) log f(; θ, α, ξ)] d. R 9

10 From the definition of f and based on the constraints C1, C2 and C3, it follows that τ() log f()] d = log(θα) + θ { } R E C (θ e αh(y ; ξ) ) log h(y ; ξ)] log ] α θ ] E C (θ e αh(y ; ξ) )H(Y ; ξ) θ { ] } + E log C (θ e αh(y ; ξ) ) C (θ e αh(y ; ξ) ) = f(; θ, α, ξ) log f(; θ, α, ξ)] d = H Sh (f), R where Y is defined as before. So, we have H Sh (τ) H Sh (f) with equality if and only if τ(; θ, α, ξ) = f(; θ, α, ξ) for all, ecept for a set of measure, thus proving the uniqueness. The intermediate steps in the above proof in fact provide the following eplicit epression for the Shannon entropy of the EWPS distribution H Sh (f) = log(θα) +α θ E C (θ e αh(y ; ξ) )H(Y ; ξ) θ { } E C (θ e αh(y ; ξ) ) log h(y ; ξ)] + log ] ] θ { ]} E C (θ e αh(y ; ξ) ) log C (θ e αh(y ; ξ) ). (18) For some EWPS distributions, the above results can only be obtained numerically. 8. Special models In this section, we investigate some special cases of the EWPS class of distributions. We offer some epressions for moments and moments of the order statistics. To illustrate the fleibility of these distributions, we provide plots of the density and hazard rate functions for selected parameter values Modified Weibull geometric distribution The modified Weibull geometric (MWG) distribution is defined by the cdf (4) with H(; ξ) = γ and = θ(1 θ) 1 leading to F (; θ, α, γ, λ) = 1 (1 θ) ep ( α γ e λ) 1 θ ep ( α γ e λ, >, (19) ) where θ (, 1). The associated pdf and hazard rate function are and f(; θ, α, γ, λ) = α(1 θ)(γ + λ) γ 1 ep ( λ α γ e λ) 1 θ ep ( α γ e λ )] 2 τ(; θ, α, γ, λ) = α(γ + λ) γ 1 ep (λ) 1 θ ep ( α γ e λ ) for >, respectively. The MWG distribution contains the WG distribution (Barreto-Souza et al. (21)) as the particular choice λ =. Further, for λ = and α = 1, we obtain the EG distribution (Adamidis and Loukas (1998)). Figures 1 and 2 display the density and hazard functions of the MWG distribution for selected parameter values. The rth raw moment of the random variable X having the MWG distribution has closed-form. It is calculated from (6) as E(X r ) = p n µ r (n), (2) 1

11 Density Density Density (a) α = 2, γ = 1.5 and λ = (b) α = 2, γ =.8 and λ = (c) α =.1, γ = 6 and λ =.5 Figure 1: Plots of the MWG density functions for θ =.1 (solid line), θ =.2 (dashed line), θ =.5 (dotted line) and θ =.9 (dotdash line). where µ r (n) = r g(; nα, γ, λ)d denotes the rth raw moment of the MW distribution with parameters nα, γ and λ. Here p n corresponds to the probability function of the geometric distribution. Carrasco et al. (28) determined an infinite representation for the rth raw moment of the MW distribution with these parameters epressed as where and µ r (n) = i 1,...,i r=1 A i1,...,i r Γ(s r /γ + 1) (nα) sr/γ, (21) A i1,...,i r = a i1,..., a ir and s r = i 1,..., i r, a i = ( 1)i+1 i i 2 (i 1)! ( ) i 1 λ. γ Hence, the moments of the MWG distribution can be obtained directly from equations (2) and (21). Hazard Hazard Hazard (a) α = 2, γ = 1.5 and λ = (b) α = 2, γ =.8 and λ = (c) α =.1, γ = 6 and λ =.5 Figure 2: Plots of the MWG hazard rate functions for θ =.1 (solid line), θ =.2 (dashed line), θ =.5 (dotted line) and θ =.9 (dotdash line). The density of the ith order statistic X i:m in a random sample of size m from the MWG distribution 11

12 is given by (for i = 1,..., m) f i:m () = m! (i 1)!(m i)! i 1 j= ( (1 θ) ep α γ e λ) ] m+j i ω j p n 1 θ ep ( α γ e λ g(; nα, γ, λ), ) where g(; nα, γ, λ) denotes the MW density function with parameters nα, γ and λ. From (13), we obtain E (Xi:m) s m! i 1 ( (1 θ) ep αx γ e λx) ] m+j i = ω j p n E (i 1)!(m i)! Xs 1 θ ep ( αx γ e λx ) Pareto Poisson distribution j= The Pareto Poisson (PP) distribution is defined by taking H(; ξ) = log(/k) and = e θ 1 in (4), which yields F (; θ, α, k) = 1 ep θ (k/)α ] 1 e θ, k. 1 and The pdf and hazard functions of the PP distribution are f(; θ, α, k) = θ α kα ep θ (k/) α ] (e θ 1) α+1 τ(; θ, α, k) = θ α kα ep θ (k/) α ] α+1 {ep θ (k/) α ] 1}. We obtain the Pareto distribution as a sub-model when θ. The rth moment of the random variable X following the PP distribution becomes E(X r ) = αkr (e θ 1) In particular, setting r = 1 in (22), the mean of X reduces to µ = αk e θ 1 θ n, nα > r. (22) (n 1)! (nα r) θ n, nα > 1. (n 1)! (nα 1) Density Density Density (a) α = 2, γ = 1.5 and λ = (b) α = 2, γ =.8 and λ = (c) α =.1, γ = 6 and λ =.5 Figure 3: Plots of the PP density functions for θ =.1 (solid line), θ =.2 (dashed line), θ =.5 (dotted line) and θ =.9 (dotdash line). 12

13 Hazard Hazard Hazard (a) α = k = (b) α = 2 and k = (c) α = 7 and k = 2 Figure 4: Plots of the PP hazard functions for θ =.1 (solid line), θ =.2 (dashed line), θ =.5 (dotted line) and θ =.9 (dotdash line). From equation (13), the sth moment of the ith order statistic, for i = 1,..., m, is given by E (Xi:m) s m! i 1 ( ep(θ (k/x) α ) m+j i ] = ω j p n E X s ) 1 (i 1)!(m i)! e θ, 1 j= where p n denotes the Poisson probability function. Furthermore, after some algebra, the Shannon entropy for the PP distribution reduces to ( e θ ) 1 H Sh (f) = log θ θα e θ 1 (µ 1 αµ 2 + µ 3 ), where µ 1 = E { ( k ep θ X { µ 2 = E ep θ and ( k X ) α } ( 1 log X ) α } log ( X k )] = )] = { } 1 Chi(2θ) log(2θ) + Shi(2θ) γ 2(e θ (e 2θ 1) log k, 1) α Chi(2θ) log(2θ) + Shi(2θ) γ 2α(e θ 1) { ( ) α } ( ) α ] k k µ 3 = E θ ep θ = α θ k2α { 1 (2θ + 1)e 2θ } X X 4(e θ, 1) where is the hyperbolic cosine integral, z Chi(z) = γ + log z + Shi(z) = z cosh(t) 1 dt t sinh(t) 1 dt t is the hyperbolic sine integral and γ is the Euler-Mascheroni constant Chen logarithmic distribution The Chen logarithmic (CL) distribution is defined by the cdf (4) with H(; ξ) = ep( β ) 1 and = log(1 θ), leading to F () = 1 log { 1 θ ep α(ep( β ) 1) ]}, >, log(1 θ) 13

14 where θ (, 1). The associated pdf and hazard rate function (for > ) are f() = θαbb 1 ep { b α ep( b ) 1 ]} log(1 θ) {θ ep α(ep( b ) 1)] 1} and respectively. τ() = θαb b 1 ep b α(ep( b ) 1) ] {θ ep α(ep( b ) 1)] 1} log {1 θ ep α(ep( b ) 1)]}, Density Density Density (a) α = b = (b) α = b = (c) α = 2.5 and b = 3 Figure 5: Plots of the CL density functions for θ =.1 (solid line), θ =.2 (dashed line), θ =.5 (dotted line) and θ =.9 (dotdash line). As epected by proposition 1, we obtain the Chen distribution as a limiting special case when θ +. Hazard Hazard Hazard (a) α = 2 and b = (b) α = 3.5 and b = (c) α = 3.5 and b = 1.5 Figure 6: Plots of the CL hazard rate functions for θ =.1 (solid line), θ =.2 (dashed line), θ =.5 (dotted line) and θ =.9 (dotdash line). The density of the ith order statistic X i:m in a random sample of size m from the CL distribution is given by (for i = 1,..., m) f i:m () = m! (i 1)!(m i)! i 1 j= { m+j 1 ωj p n g(; nα, b) log 1 θ ep(α α e )]} b, where g(; nα, b) is the pdf of the Chen distribution with parameters nα and b and p n denotes the 14

15 logarithmic probability mass function and ( ) ] m+j 1 i 1 ωj = ( 1) j 1. j log(1 θ) In the same way, the sth raw moment of X i:m is obtained directly from E (X s i:m) = where Z Chen(nα, b). m! (i 1)!(m i)! i 1 ω j p n E { Z s ep nα(m + j 1)(1 ep(z b )) ]}, j= 9. Application Fonseca and França (27) studied the soil fertility influence and the characterization of the biologic fiation of N 2 for the Dimorphandra wilsonii rizz growth. For 128 plants, they made measures of the phosphorus concentration in the leaves. The data are listed in Table 3. We fit the MWG, Gompertz Poisson (GP), PP, Chen Poisson (CP) and CL models to these data. We also fit the three-parameter WG distribution introduced by Barreto-Souza et al. (21). The required numerical evaluations are implemented using the SAS (PROCNLMIXED) and R softwares Table 3: Phosphorus concentration in leaves data set. Min. Q 1 Q 2 Mean Q 3 Ma. Var Table 4: Descriptive statistics. Tables 4 and 5 display some descriptive statistics and the MLEs (with corresponding standard errors in parentheses) of the model parameters. Since the values of the Akaike information criterion (AIC), Bayesian information criterion (BIC) and consistent Akaike information criterion (CAIC) are smaller for the CL distribution compared with those values of the other models, this new distribution seems to be a very competitive model for these data. Plots of the estimated pdf and cdf of the MWG, WG, GP, PP, CP and CL models fitted to these data are displayed in Figure 7. They indicate that the CL distribution is superior to the other distributions in terms of model fitting. Table 6 lists the values of the Kolmogorov-Smirnov (K-S) statistic and the values of 2l( Θ). From these figures, we conclude that the CL distribution provides a better fit to these data than the MWG, WG, GP, PP and CP models. 1. Concluding remarks We define a new lifetime class of distributions, called the etended Weibull power series (EWPS), which generalizes the Weibull power series class of distributions proposed by Morais and Barreto-Souza (211), which in turn etends the eponential power series class of distributions (Chahkandi and Ganjali, 29). We provide a mathematical treatment of the new distribution including epansions for the density function, moments, generating function and incomplete moments. Further, eplicit epressions for the 15

16 Estimates Statistic Model θ α γ λ AIC BIC AICC MWG (.2418) ( ) (.821) (6.6673) WG (.17) (8.759) (.322) θ α β GP (1.2627) (.1473) (1.6135) θ α k PP (69.777) (.115) θ α b CP ( ) ( ) (.2634) CL (.1) (.) (.239) Table 5: MLEs of the model parameters, the corresponding SEs (given in parentheses) and the statistics AIC, BIC and AICC. Density CL MWG PP WG GP CP cdf CL MWG PP WG GP CP (a) (b) Figure 7: Estimated (a) pdf and (b) cdf for the CL, MWG, PP, WG, GP and CP models to the percentage of Phosphorus concentration in leaves data. 16

17 Model K S 2l( Θ) MWG WG Gompertz Poisson Pareto Poisson Chen Poisson Chen Logarithmic Table 6: K-S statistics and 2l( Θ) for the eceedances of phosphorus concentration in leaves data set. order statistics and Shannon entropy are derived. The EWPS density function can be epressed as a miture of EW density functions. This property is important to obtain several other results. Our formulas related with the EWPS model are manageable, and with the use of modern computer resources with analytic and numerical capabilities, they may turn into adequate tools comprising the arsenal of applied statisticians. The estimation of the model parameters is approached by the method of maimum likelihood using the EM algorithm. The observed information matri is derived. Further, maimum entropy identification for the EWPS distributions was discussed and some special models are studied in some detail. Finally, we fit the EWPS model to a real data set to show the usefulness of the proposed class. We hope that this generalization may attract wider applications in the literature of the fatigue life distributions. Acknowledgements We also gratefully acknowledge financial support from CAPES and CNPq. Appendi A. The elements of the p p information matri J n (Θ) are J θθ = n θ 2 n C ( (θ) C ) ] 2 (θ) ( ) 2 z2i 2αH(i; ξ) + θ H( i ; ξ)e z 1i z 3i 2αH(i; ξ) θ H( i ; ξ)e z 1i J αα = n n α 2 + θ z 2i H 2 ( i ; ξ)e αh(i; ξ) + θ 2 z 1i (z2i ) ] 2 J αθ = θ z 3i H 2 ( i ; ξ)e 2αH(i; ξ) J αξk J θξk J ξk ξ l z 1i z 1i (z 3i z2i 2 ) H 2 2αH(i; ξ) ( i ; ξ)e z 1i z 2i H 2 αh(i; ξ) ( i ; ξ)e z 1i H( i ; ξ) z 2i H( i ; ξ) = θ e αh(i; ξ) 1 αh( i ; ξ)] ξ k z 1i ξ k n ( ) ] 2 + αθ 2 z 3i z2i H( i ; ξ) 2αH(i; ξ) H( i ; ξ)e z 1i z 1i ξ k (z2i ) ] 2 = θα z 3i H( i ; ξ) e 2αH(i; ξ) z 2i H( i ; ξ) αh(i; ξ) α e z 1i z 1i ξ k z 1i ξ k 2 H( i ; ξ) 1 H( i ; ξ) H( i ; ξ) = α ξ k ξ l H( i ; ξ) 2 + ξ k ξ l n (z2i ) ] 2 + (αθ) 2 + z 3i H( i ; ξ) H( i ; ξ) 2αH(i; ξ) e ξ k ξ l αθ z 1i z 1i z 2i 2 H( i ; ξ) e αh(i; ξ) + α 2 θ z 1i ξ k ξ l 17 z 2i H( i ; ξ) H( i ; ξ) z 1i ξ k ξ l 1 2 H( i ; ξ) H( i ; ξ) ξ k ξ l e αh(i;ξ)

18 where z 1i = C (θe αh(i; ξ) ), z 2i = C (θe αh(i; ξ) ) and z 3i = C (θe αh(i; ξ) ), for i = 1,..., n. References References Adamidis K., Dimitrakopoulou, T., Loukas, S. (25). On a generalization of the eponential-geometric distribution. Statistics & Probability Letters, 73, Adamidis K., Loukas, S. (1998). A lifetime distribution with decreasing failure rate. Statistics & Probability Letters, 39, Barreto-Souza, W., Cribari-Neto, F. (29). A generalization of the eponential-poisson distribution. Statistics & Probability Letters, 79, Barreto-Souza, W., Morais, A.L., Cordeiro, G.M. (21). The Weibull-geometric distribution. Journal of Statistical Computation and Simulation, 81, Bebbington, M., Lai, C. D. and Zitikis, R. (27). A fleible Weibull etension. Reliability Engineering and System Safety, 92, Carrasco J.M.F., Ortega, E.M.M., Cordeiro, G.M. (28). A generalized modified Weibull distribution for lifetime modeling. Computational Statistics & Data Analysis, 53, Chahkandi, M., Ganjali, M. (29). On some lifetime distributions with decreasing failure rate. Computational Statistics & Data Analysis, 53, Chen, Z. (2). A new two-parameter lifetime distribution with bathtub shape or increasing failure rate function. Statistics and Probability Letters, 49, Dempster, A.P., Laird, N.M., Rubim, D.B. (1977). Maimum likelihood from incomplete data via the EM algorithm (with discussion). Journal of the Royal Statistical Society. Series B, 39, Doornik, J. (27) O 5: object-oriented matri programming language, 5th ed. Timberlake Consultants, London. Fonseca, M.B. and França, M.G.C. (27). A influência da fertilidade do solo e caracterização da fiação biológica de N 2 para o crescimento de Dimorphandra wilsonii rizz. Master s thesis, Universidade Federal de Minas Gerais. Ghitany, M.E., Al-Jarallah, R.A., Balakrishnan, N. (211): On the eistence and uniqueness of the MLEs of the parameters of a general class of eponentiated distributions. Statistics, DOI:1.18/ Gompertz, B. (1825). On the nature of the function epressive of the law of human mortality and on the new model of determining the value of life contingencies. Philosophical Trans. Royal Society of London, 115, Gupta, R.D. and Kundu, D. (1999). Generalized eponential distributions. Austral. NZ J. Statist., 41, Gupta, R.D. and Kundu, D. (21). Eponentiated eponential distribution: An alternative to gamma and Weibull distributions. Biometrical Journal, 43, Gupta, R.D. and Kundu, D. (21). Generalized eponential distributions: Different methods of estimations. Journal of Statistical Computation and Simulation, 69, Gupta, R.D. and Kundu, D. (27). Generalized eponential distributions: Eisting results and some recent developments. Journal of Statistical Planning and Inference, 137, Gupta, R.C., Gupta, R.D. and Gupta, P.L Modeling failure time data by Lehman alternatives. Communications in Statistics: Theory and Methods, 27, Gurvich, M., DiBenedetto, A., Ranade, S. (1997). A new statistical distribution for characterizing the random strength of brittle materials. Journal of Materials Science, 32, Jaynes, E.T. (1957). Information theory and statistical mechanics. Physical Reviews, 16, Johnson, N.L, Kotz, S. and Balakrishnan, N. (1994). Continuous Univariate Distributions volume 1. John Wiley & Sons, New York. Kapur, J.N. (1989). Maimum Entropy Models in Science and Engineering. John Wiley & Sons, New York. Kies, J. A. (1958). The strength of glass. Washington D.C. Naval Research Lab, 593. Kundu, D. and Raqab, M.Z. (25). Generalized Rayleigh distribution: Different methods of estimation. Computational Statistics & Data Analysis, 49, Kus, C. (27). A new lifetime distribution. Computational Statistics & Data Analysis, 51, Lai, C.D., Xie, M. and D. N. P. Murthy. (23). A modified weibull distribution. Transactions on Reliability, 52, Lu, W., Shi, D. (211). A new compounding life distribution: the Weibull-Poisson distribution. Journal of Applied Statistics, DOI:1.18/ Morais, A.L., Barreto-Souza, W. (211). A Compound Class of Weibull and Power Series Distributions. Computational Statistics & Data Analysis, 55, Nadarajah, S. and Kotz, S. (25). On some recent modifications of Weibull distribution. IEEE Trans. Reliability, 54, Nikulin, M. and Haghighi, F. (26). A Chi-squared test for the generalized power Weibull family for the head-and-neck cancer censored data. Journal of Mathematical Sciences, 133, Pham, H. (22). A vtub-shaped hazard rate function with applications to sys- tem safety. International Journal of Reliability and Applications, 3, Phani, K.K. (1987). A new modifiedweibull distribution function. Communications of the American Ceramic Society, 7, Rayleigh, J. W. S. (188). On the resultant of a large number of vibrations of the same pitch and of arbitrary phase. Philosophical Magazine, 1, Shannon, C.E., (1948). A mathematical theory of communication. Bell System Technical Journal, 27, Shore, J.E., Johnson, R.W. (198). Aiomatic derivation of the principle of maimum entropy and the principle of minimum cross-entropy. IEEE Transactions on Information Theory, 26, Smith, R.M. and Bain, L.J. (1975). An eponential power life-testing distribution. Communications Statistics, 4, Soofi, E.S. (2). Principal information theoretic approaches. Journal of the American Statistical Association, 95,

19 Surles, J.G. and Padgett, W.J. (21). Inference for reliability and stress-strength for a scaled Burr type X distribution. Lifetime Data Analysis., 7, Silva, R.B., Barreto-Souza, W., Cordeiro, G.M. (21). A new distribution with decreasing, increasing and upside-down bathtub failure rate. Computational Statistics & Data Analysis, 54, Tahmasbi, R., Rezaei, S. (28). A two-parameter lifetime distribution with decreasing failure rate. Computational Statistics & Data Analysis, 52, Xie, M. and Lai, C.D. (1995). Reliability analysis using additive Weibull model with bathtub-shaped failure rate function. Reliability Engineering and System Safety, 52, Xie, M., Tang, Y. and Goh, T.N. (22). A modified Weibull etension with bathtub-shaped failure rate function. Reliability Engineering & System Safety, 76, Zografos, K., Balakrishnan, N. (29). On families of beta-and generalized gamma-generated distributions and associated inference. Statistical Methodology, 6, White, J. S. (1969). The moments of log-weibull order statistics. Technometrics, 11,