arxiv: v2 [stat.me] 6 Dec 2012

Transcription

1 A new class of faigue life disribuions Marcelo Bourguignon a,1, Rodrigo B. Silva a,2, Gauss M. Cordeiro a,3 a Universidade Federal de Pernambuco Deparameno de Esaísica, Cidade Universiária, Recife, PE, Brazil arxiv: v2 [sa.me] 6 Dec 2012 Absrac In his paper, we inroduce he Birnbaum-Saunders power series class of disribuions which is obained by compounding Birnbaum-Saunders power series disribuions. The new class of disribuions has as a paricular case he wo-parameer Birnbaum-Saunders disribuion. The hazard rae funcion of he proposed class can be increasing upside-down bahub shaped. We provide imporan mahemaical properies such as momens, order saisics, esimaion of he parameers inference for large sample. Three special cases of he new class are invesigaed wih some deails. We illusrae he usefulness of he new disribuions by means of wo applicaions o real daa ses. Keywords: Birnbaum-Saunders disribuion; Maximum likelihood esimaion; Order saisic; Power series disribuion. 1. Inroducion Birnbaum Saunders (1969 inroduced a wo-parameer disribuion in order o model he failure ime disribuion for faigue failure caused under cyclic loading, called he Birnbaum-Saunders disribuion (BS. Since hen, his disribuion has been exensively used for modelling failure imes of faiguing maerials in several fields such as environmenal medical sciences, engineering, demography, biological sudies, acuarial, economics, finance insurance. A rom variable T following he BS(α, disribuion wih shape parameer α > 0 scale parameer > 0 is defined by [ (αz 2 1/2 2 αz T = ], 2 where Z is a sard normal rom variable. The cumulaive disribuion funcion (cdf of T BS(α, has he form F BS (; α, = Φ(υ, > 0, (1 where υ = α 1 ρ(/, ρ(z = z 1/2 z 1/2 Φ( is he sard normal cumulaive funcion. The corresponding probabiliy densiy funcion (pdf is [ f BS (; α, = κ(α, 3/2 ( + exp τ(/ ] 2α 2, > 0, (2 where κ(α, = exp(α 2 /(2α 2π τ(z = z + z 1. The parameer is he median of he disribuion, i.e. F BS ( = Φ(0 = 1/2. For any k > 0, kt BS(α, k. Addiionally, he reciprocal propery holds, i.e., T 1 BS(α, 1, see Saunders (1974. The expeced value variance are E(T = (1 + α 2 /2 Var(T = (α 2 (1 + 5α 2 /4, respecively. The BS disribuion has a close relaion o he normal disribuion has been widely sudied in he saisical lieraure. For example, Ng e al. (2006 proposed a simple bias-reducion mehod o reduce he biases of he maximum likelihood esimaors (MLEs for he BS model. They also discussed 1 Corresponding auhor. m.p.bourguignon@gmail.com 2 rodrigobs29@gmail.com 3 gauss@de.ufpe.br Preprin submied o Elsevier Sepember 22, 2018

2 a Mone Carlo EM-algorihm for compuing heses esimaors. Kundu e al. (2008 invesigaed he shape of he BS hazard funcion. In he las few years, several classes of disribuions were proposed by compounding some useful lifeime power series disribuions. Chahki Ganjali (2009 sudied he exponenial power series (EPS family of disribuions, which conains as special cases he exponenial geomeric (EG, exponenial Poisson (EP exponenial logarihmic (EL disribuions, inroduced by Adamidis Loukas (1998, Kus (2007 Tahmasbi Rezaei (2008, respecively. Morais Barreo-Souza (2011 defined he Weibull power series (WPS class of disribuions which includes as sub-models he EPS disribuions. The WPS disribuions can have an increasing, decreasing upside down bahub failure rae funcion. The generalized exponenial power series (GEPS disribuions were proposed by Mahmoudi Jafari (2012 in a similar way o ha found in Morais Barreo-Souza (2011. In a very recen paper, Silva emphe al. (2012 inroduced he exended Weibull power series (EWPS class of disribuions which includes as special models he EPS WPS disribuions. In a similar manner, we inroduce he Birnbaum-Saunders power series (BSPS class of disribuions, which is defined by compounding he Birnbaum-Saunders power series disribuions. The main aim of his paper is o inroduce a new class of disribuions, which exends he Birnbaum-Saunders disribuion, wih he hope ha he new disribuion may have a beer fi in cerain pracical siuaions. Addiionally, we will provide a comprehensive accoun of he mahemaical properies of he proposed new class of disribuions. The hazard rae funcion of he proposed class can be increasing upsidedown bahub shaped. We are moivaed o inroduce he BSPS disribuions because of he wide usage of (1 he fac ha he curren generalizaion provides means of is coninuous exension o sill more complex siuaions. This paper is unfolds as follows. In Secion 2, a new class of disribuions is obained by mixing he BS zero runcaed power series disribuions. The mixing procedure was previously carried ou by Morais Barreo-Souza (2011 Silva e al. (2012. In Secion 3, some properies of he new class are discussed. In Secion 4, he esimaion of he model parameers is performed by he mehod of maximum likelihood. In Secion 5, we inroduce sudy hree special cases of he BSPS class of disribuions. In Secion 6, wo illusraive applicaions based on real daa are provided. Finally, concluding remarks are presened in Secion The new class Le T 1,..., T N be independen idenically disribued (iid BS rom variables wih cdf (1 pdf (2. We assume ha he rom variable N has a zero runcaed power series disribuion wih probabiliy mass funcion p n = P (N = n = a n θ n, n = 1, 2,..., (3 C(θ where a n depends only on n, C(θ = n=1 a n θ n θ > 0 is such ha C(θ is finie. Table 1, repored in Morais Barreo-Souza (2011, summarizes some power series disribuions (runcaed a zero defined according o (3 such as he Poisson, logarihmic, geomeric binomial disribuions. By assuming ha he rom variables N T are independen, we define X = min(t 1,..., T N. Then, P (X x, N = n = a n θ n C(θ [1 (1 Φ(υn ], x > 0, n 1. The BSPS class of disribuions is defined by he marginal cdf of X: F BSPS (x; θ, α, = 1 C [θ(1 Φ(υ], x > 0. (4 C(θ Hereafer, he rom variable X following (4 wih parameers θ, α is denoed by X BSPS(θ, α,. The pdf o (4 is where f BS ( is given in (2. f BSPS (x; θ, α, = θf BS (x C [θ(1 Φ(υ], x > 0, (5 C(θ 2

3 Disribuion a n C(θ C (θ C (θ C(θ 1 Θ Poisson n! 1 e θ 1 e θ e θ log(θ + 1 θ (0, Logarihmic n 1 log(1 θ (1 θ 1 (1 θ 2 1 e θ θ (0, 1 Geomeric 1 θ(1 θ 1 (1 θ 2 2(1 θ 3 θ(θ θ (0, 1 Binomial ( m n (θ + 1 m 1 m(θ + 1 m1 m(m 1 (θ + 1 2m (θ 1 1/m 1 θ (0, 1 Table 1: Useful quaniies for some power series disribuions. Proposiion 1. The classical BS disribuion wih parameers α is a limiing special case of he proposed model when θ 0. Proof. For x > 0, we have lim F BSPS(x = 1 lim θ 0 + θ 0 + = 1 lim θ 0 + a n [θ(1 Φ(υ] n n=1 a n θ n n=1 1 Φ(υ + a a 1 1 a n θ n1 (1 Φ(υ n n=2 a n θ n1 n=2 = Φ(υ. We now derive wo useful expansions for he densiy funcion (5. We have C (θ = n=1 n a n θ n1 hen f BSPS (; θ, α, = f BS (; α, n=1 By using he binomial heorem in equaion (6, we can wrie ( n 1 where ω n,k = (1 k n p n. k f BSPS (x; θ, α, = f BS (x; α, n a n θ n C(θ [1 Φ(υ]n1. (6 n1 ω n,k Φ(υ k, (7 n=1 k=0 Proposiion 2. The densiy funcion of X can be expressed as an infinie linear combinaion of densiies of minimum order saisics of T. Proof. We know ha C (θ = n=1 n a n θ n1. Therefore, f BSPS (x; θ, α, = p n f T(1 (x; α,, n, x > 0, where f T(1 (x; α,, n is he densiy funcion of T (1 = min(t 1,..., T n, for fixed n, given by n=1 f T(1 (; α,, n = n f BS ([1 Φ(v] n1, > 0. 3

4 Remark 1. Le Y = max(t 1,..., T N, hen he cdf pdf of Y are Noe ha F Y (y; θ, α, = C(θ Φ(v, y > 0 (8 C(θ f Y (y; θ, α, = θf BS (y C (θ Φ(v, y > 0. C(θ f Y (y; θ, α, = p n f T(n (y; α,, n, y > 0, n=1 where f T(n (y; α,, n is he densiy funcion of T (n = max(t 1,..., T n, for fixed n, given by f T(n (; α,, n = n f BS ( Φ(v n1, > 0. The disribuion wih cdf (8 is called he complemenary Birnbaum-Saunders power series (CBSPS disribuion. This class of disribuions is a suiable model in a complemenary risk problem based in he presence of laen risks which arise in several areas such as public healh, acuarial science, biomedical sudies, demography indusrial reliabiliy (Basu Klein, However, in his work, we do no focus on his alernaive class of disribuions. The BSPS survival funcion becomes S BSPS (x; θ, α, = he corresponding hazard rae funcion reduces o C [θ(1 Φ(υ], x > 0 C(θ h BSPS (x; θ, α, = θf BS (x C [θ(1 Φ(υ] C [θ(1 Φ(υ], x > Quaniles, order saisics momens The quanile funcion of X, say Q(u, is given by = Q(u = α δ 2 + [ (α δ ] 1/2 where δ = Φ 1 { C 1 [(1 uc(θ] /θ }, Φ 1 ( is he inverse cumulaive funcion of he sard normal disribuion, C 1 ( is he inverse funcion of C( u is a uniform rom number on he uni inerval (0, 1. Le X 1,..., X m be a rom sample wih densiy funcion (5 define he ih order saisic by T i:m. The densiy funcion of T i:m, say f i:m (x; θ, α,, is given by 2, f i:m (x; θ, α, = = 1 B(i, m i + 1 f BSPS(x; θ, α, F BSPS (x; θ, α, i1 [1 F BSPS (x; θ, α, ] mi f BS (x; α, i1 ( i 1 (1 k B(i, m i + 1 k k=0 S BSPS (x; θ, α, m+ki, where S BSPS ( is he BSPS survival funcion. Proposiion 3. The BSPS order saisics can be expressed as a linear combinaion of exponeniaed Birnbaum-Saunders (EBS disribuions wih parameers α, s + 1. In fac, f i:m (x; θ, α, = i1 k=0 j=0 m+ki+j 4 s=0 ω k,j,s h s+1 (x, α,,

5 where (1 k+s( i1 k ( m+k+ji s dj ω k,j,s = (s + 1B(i, m i + 1C(θ m+ki, d 0 = a m+ki 1, d r = 1 rd 0 0 r [j(m + k i r + j] a j+1 d rj j=1 h s+1 (, α, denoes he EBS densiy funcion wih parameers α, s + 1. Many of he imporan characerisics feaures of a disribuion are obained hrough he momens. The ordinary momens of X can be derived from he probabiliy weighed momens (PWMs (Greenwood e al., 1979 of he BS disribuion formally defined for p r non-negaive inegers by [ τ p,r = κ(α, p3/2 ( + exp τ(/ ] 2α 2 Φ r (υd. (9 The inegral (9 can be easily compued numerically in sofware such as MAPLE, MATHEMATICA, Ox R, see Cordeiro Lemone (2011. An alernaive expression for τ p,r can be derived from Cordeiro Lemone (2011 as τ p,r = p 2 r r j=0 ( r A(k 1,..., k j j k 1,...,k j=0 2s j+j m=0 (1 m ( 2sj + j m (2sj+j2m/2 I ( p + (2s j + j 2m/2, α. Here, s j = k k j, A(k 1,..., k j = α 2sjj a k1... a kj, a k = (1 k 2 (12k/2 { π(2k + 1k!} 1 I(p, α = K p+1/2(α 2 + K p1/2 (α 2 2K 1/2 (α 2, where he funcion K ν (z denoes he modified Bessel funcion of he hird kind wih ν represening is order z he argumen. Is inegral represenaion is K ν (z = 0.5 exp{z cosh( ν}d. Thus, he sh momen of X can be obained from equaions (7 (9 as E(X s = n1 ω n,k τ s,k. n=1 k=0 Expressions for he sh momens of he order saisics X 1:m,..., X m,m wih cumulaive funcion (4 can be obained using a resul due o Baraka Abdelkader (2004. We have for i = 1,..., m. E(X s i:m = s = s m j=mi+1 m j=mi+1 ( ( j 1 m (1 jm+i1 x s1 S BSPS (x, θ, α, j dx m i j 0 (1 jm+i1 ( ( j 1 m C(θ j x s1 C [θ(1 Φ(υ] j dx, m i j 0 4. Maximum likelihood esimaion Here, we discuss maximum likelihood esimaion inference for he BSPS disribuion. Le x 1,..., x n be a rom sample from X le Θ = (θ, α, be he vecor of he model parameers. The log-likelihood funcion for Θ reduces o l(θ, α, = n{log(θ log[c(θ] + log[κ(α, ]} α 2 τ ( xi + log(x i + log{c [θ(1 Φ(v i ]}. 5 log(x i +

6 The componens of he score vecor U(Θ are given by 1 φ(υ i C [θ(1 Φ(υ i ] U α (Θ = n (1 + 2α α n ( xi α 3 τ U θ (Θ = n θ nc (θ C(θ + n + θ α 3 n C [θ(1 Φ(υ i ], θ φ(υ i ρ(x i / C [θ(1 Φ(υ i ] C [θ(1 Φ(υ i ] α U (Θ = n α 2 ( xi + x i θ 2α φ(υ i τ( x i / C [θ(1 Φ(υ i ] C [θ(1 Φ(υ i ], where φ( is he sard normal densiy v i = α 1 ρ(x i /θ. Seing hese equaions o zero, U(Θ = 0, solving hem simulaneously yields he MLE Θ of Θ. These equaions canno be solved analyically saisical sofware can be used o compue hem numerically using ieraive echniques such as he Newon-Raphson algorihm. Ofen wih lifeime daa reliabiliy sudies, one encouners censoring. A very simple rom censoring mechanism ha is ofen realisic is one in which each individual i is assumed o have a lifeime X i a censoring ime C i, where X i C i are independen rom variables. Suppose ha he daa consis of n independen observaions x i = min(x i, C i δ i = I(X i C i is such ha δ i = 1 if X i is a ime o even δ i = 0 if i is righ censored for i = 1,..., n. The censored likelihood L(θ, α, for he model parameers is L(θ, α, n [f BSPS (x i ; θ, α, ] δi [S BSPS (x i ; θ, α, ] 1δi. In order o perform inerval esimaion hypohesis ess on he model parameers α, θ, he normal approximaion for he MLE can be applied. Noe ha under cerain condiions for he parameers, he asympoic disribuion of n( Θ Θ is mulivariae normal N 3 (0, I 1 (Θ, where J n (θ is he observed informaion marix I(Θ = lim n J n (Θ. Noe also ha an 100(1 γ% (0 < γ < 1/2 asympoic confidence inerval for he ih parameer θ i in Θ is specified by ACI i = ( θ i z 1γ/2 Ĵ θ i,θ i, ˆθi + z 1γ/2 Ĵ θ i,θ i, where Ĵ θi,θi ss for he ih diagonal elemen of he inverse of he observed informaion marix esimaed a Θ, i.e., J n ( Θ 1, for i = 1, 2, 3, z 1γ/2 is he 1 γ/2 sard normal quanile. For esing he goodness-of-fi of he BSPS disribuion for comparing i wih some sub-models, we can use he likelihood raio (LR saisics. The LR saisic for esing he null hypohesis H 0 : Θ 1 = Θ (0 1 versus he alernaive hypohesis H 1 : Θ 1 Θ (0 1 is given by w = 2{l( Θ l( Θ}, where Θ Θ are he MLEs under he alernaive null hypoheses, respecively. Under he null hypohesis, H 0, w χ d 2 k, where k is he dimension of he subse Θ 1 of ineres. 5. Special cases In his secion, we invesigae some special models of he BSPS class of disribuions. We offer some expressions for he momens momens of he order saisics. We illusrae he flexibiliy of hese disribuions provide plos of he densiy hazard rae funcions for seleced parameer values Birnbaum-Saunders geomeric disribuion The Birnbaum-Saunders geomeric (BSG disribuion is defined by he cdf (4 wih C(θ = θ(1θ 1 leading o (1 θ[1 Φ(v] F BSG (x; θ, α, = 1, x > 0, (10 1 θ[1 Φ(v] where θ (0, 1. 6

7 Proposiion 4. The disribuion of he form (10 is geomeric minimum sable. Proof. The proof follows easily of he argumens given by Marshall Olkin (1997, p We omi he deails. The associaed densiy hazard rae funcions are f BSG (x; θ, α, = (1 θf BS (x [1 θ(1 Φ(v] 2 h BSG (x; θ, α, = f BS (x [1 Φ(v]{1 θ[1 Φ(v]} = h BS (x 1 θ[1 Φ(v], (11 for x > 0, respecively. From equaion (11, we noe ha h BSG (x; θ, α, /h BS (x; α, is decreasing in x. I can also be shown ha lim h BSG(x; θ, α, = 0 lim h BSG(x; θ, α, = 1 x 0 x 2α 2. Corollary 1. Le X BSG(θ, α, T BS(α,. Then h BSG (x; θ, α, h BS (; α, for all θ (0, 1. Proof. I is sraighforward. Cancho e al. (2012 proposed sudied he geomeric Birnbaum-Saunders regression wih cure rae. However, hey considered C(θ = (1 θ 1 insead of C(θ = θ(1 θ 1. In Figure 1, we plo he densiy hazard rae funcions of he BSG disribuion for seleced parameer values. We can verify ha his disribuion has an upside-down bahub or an increasing failure rae funcion depending on he values of is parameers. Expressions for he densiy funcion he momens of he order saisics X 1:m,..., X m:m from a rom sample of he BSG disribuion are given by f i:m (x = f BS(x; α, i1 ( [ ] m+ki i 1 (1 θφ(υ (1 k B(i, m i + 1 k 1 θφ(υ E(X s i:m = s m j=mi+1 k=0 (1 jm+i1 θ j (1 θ j ( j 1 m i ( m x s1 [1 θφ(υ] j dx, j which are easily deermined numerically. Le x 1,, x n be a rom sample of size n from X BSG(δ, α,. The log-likelihood funcion for he vecor of parameers Θ = (θ, α, can be expressed as l(θ = n{log(1 θ + log[κ(α, ]} 3 log(x i + log(x i α 2 τ ( xi 2 log[1 θ(1 Φ(v i ] Birnbaum-Saunders Poisson disribuion The Birnbaum-Saunders Poisson (BSP disribuion follows by aking C(θ = e θ 1 in (4, which yields F BSP (; θ, α, = 1 0 exp{θ(1 Φ(v} 1. exp(θ 1 The densiy hazard rae funcions of he BSP disribuion are f BSP (; θ, α, = θf BS( exp[θ(1 Φ(v] exp(θ 1 h BSP (; θ, α, = θf BS( exp{θ(1 Φ(v} exp{θ(1 Φ(v} 1. 7

8 Densiy θ = 0.2 θ = 0.4 θ = 0.6 θ = 0.8 Densiy θ = 0.2 θ = 0.4 θ = 0.6 θ = (a α = 0.2 = 5.0 (b α = 2.0 = 5.0 Hazard funcion θ = 0.2 θ = 0.4 θ = 0.6 θ = 0.8 Hazard funcion θ = 0.2 θ = 0.4 θ = 0.6 θ = (c α = 0.2 = 5.0 (d α = 2.0 = 5.0 Figure 1: Plos of he BSG densiy hazard rae funcions for some parameer values. 8

9 Remark 2. The limi of he BSP hazard rae funcion as 0 is 0. In Figure 2, we plo he densiy hazard rae funcions of he BSP disribuion for seleced parameer values. This disribuion can have an upside-down bahub or an increasing failure rae funcion depending on he values of is parameers. Densiy θ = 0.5 θ = 1 θ = 3 θ = 5 Densiy θ = 0.5 θ = 1 θ = 3 θ = (a α = 0.2 = 1.0 (b α = 2.0 = 1.0 Hazard funcion θ = 0.5 θ = 1 θ = 3 θ = 5 Hazard funcion θ = 0.5 θ = 1 θ = 3 θ = (c α = 0.2 = 1.0 (d α = 2.0 = 1.0 Figure 2: Plos of he BSP densiy hazard rae funcions for some parameer values. The expressions for he densiy funcion he momens of he order saisics X 1:m,..., X m:m from a rom sample of he BSP disribuion are f i:m (x = f BS(x; α, i1 ( { i 1 exp [θφ(υ] 1 (1 k B(i, m i + 1 k e θ 1 E(X s i:m = s m j=mi+1 k=0 (1 jm+i1 (exp(θ 1 j ( j 1 m i } m+ki ( m x s1 {exp [θφ(υ] 1} j dx. j 0 9

10 Le x 1,, x n be a rom sample of size n from X BSP(θ, α,. The log-likelihood funcion for he vecor of parameers Θ = (θ, α, can be expressed as l(θ = n{log(θ log[exp(θ 1] + log[κ(α, ]} α 2 τ ( xi + θ [1 Φ(v i ] Birnbaum-Saunders logarihmic disribuion log(x i + log(x i + The cdf of he Birnbaum-Saunders logarihmic (BSL disribuion is defined by (4 wih C(θ = log(1 θ, 0 < θ < 1, which corresponds o he logarihmic disribuion. We obain F BSL (x; θ, α, = 1 The associaed densiy hazard rae funcions are log [1 θ(1 Φ(υ] log(1 θ θf BS(x f BSL (x; θ, α, = log(1 θ {1 θ[1 Φ(υ]} θf BS (x h BSL (x; θ, α, = log {1 θ[1 Φ(υ]} {1 θ[1 Φ(υ]} Remark 3. The limi of he BSL hazard rae funcion as x 0 is 0. In Figure 3, we plo he densiy hazard rae funcions of he BSL disribuion for seleced parameer values. We can verify ha his disribuion can have an upside-down bahub or an increasing failure rae funcion depending on he values of is parameers. The expressions for he densiy funcion he momens of he order saisics X 1:m,..., X m:m from a rom sample of he BSL disribuion are f i:m (x = f BS(x; α, i1 ( { i 1 log [1 θφ(υ] (1 k B(i, m i + 1 k log(1 θ k=0 } m+ki E(X s i:m = s m j=mi+1 (1 im1 [log(1 θ] j ( j 1 m i ( m j 0 [ x {log s1 1 1 θφ(υ ]} j dx. Le x 1,, x n be a rom sample of size n from X BSL(θ, α,. The log-likelihood funcion for he vecor of parameers Θ = (θ, α, can be expressed as l(θ = n {log(θ log [ log(1 θ] + log[κ(α, ]} Applicaions 1 2α 2 τ ( xi log {1 θ[1 Φ(υ i ]}. log(x i + log(x i + In his secion, we compare he fis of some special models of he BSPS class by means of wo real daa ses o show he poenialiy of he new class. In order o esimae he parameers of hese special models, we adop he maximum likelihood mehod (as discussed in Secion 4 all he compuaions were done using he subrouine NLMixed of he SAS sofware. Obviously, due o he genesis of he BS disribuion, he faigue processes are by excellence ideally modeled by his disribuion. Thus, he use of he BS BSPS disribuions for fiing hese wo daa ses is well jusified. 10

11 Densiy θ = 0.2 θ = 0.4 θ = 0.6 θ = 0.8 Densiy θ = 0.2 θ = 0.4 θ = 0.6 θ = (a α = 0.2 = 1.0 (b α = 2.0 = 1.0 Hazard θ = 0.2 θ = 0.4 θ = 0.6 θ = 0.8 Hazard θ = 0.2 θ = 0.4 θ = 0.6 θ = (c α = 0.2 = 1.0 (d α = 2.0 = 1.0 Figure 3: Plos of he BSL densiy hazard rae funcions for some parameer values. 11

12 Firs, we consider a daa se from Murhy e al. (2004 consising of he failure imes of 20 mechanical componens. The daa are: 0.067, 0.068, 0.076, 0.081, 0.084, 0.085, 0.085, 0.086, 0.089, 0.098, 0.098, 0.114, 0.114, 0.115, 0.121, 0.125, 0.131, 0.149, 0.160, The MLEs of he parameers (wih corresponding sard errors in parenheses, he value of 2l( Θ, he Kolmogorov-Smirnov saisic, Akaike informaion crierion (AIC Bayesian informaion crierion (BIC for he BSG, BSP, BSL BS models are lised in Table 2. Since he values of he AIC BIC are smaller for he BSG disribuion compared wih hose values of he oher models, he new disribuion seems o be a very compeiive model for hese daa. Table 2: Parameer esimaes, K-S saisics, AIC BIC for failure imes daa Disribuion α θ K S 2l( Θ AIC BIC BSG ( a ( ( BSP ( a ( ( BSL ( a ( ( BS ( a ( a Denoes he sard deviaion of he MLE s of α, θ. Plos of he pdf cdf of he BSG, BSP, BSL BS fied models o hese daa are displayed in Figure 4. They indicae ha he BSG disribuion is superior o he oher disribuions in erms of model fiing. Addiionally, i is eviden ha he BS disribuion presens he wors fi o he curren daa hen he proposed models ouperform his disribuion. Densiy BSG BSP BSL BS cdf BSG BSP BSL BS (a (b Figure 4: Esimaed (a pdf (b cdf for he BSG, BSP, BSL BS models for failure imes daa. We can also perform formal goodness-of-fi ess in order o verify which disribuion fis beer o hese daa. Chen Balakrishnan (1995 proposed a general approximae goodness-of-fi es for he hypohesis H 0 : X 1,..., X n wih X i following F (x; θ, i.e., under H 0, X 1,..., X n is a rom sample from a coninuous disribuion wih cumulaive disribuion F (x; θ, where he form of F is known bu he p-vecor θ is unknown. The mehod is based on he Cramér-von Mises (C M Anderson-Darling (A D saisics, in general, he smaller he values of hese saisics, he beer he fi. 12

13 Table 3 gives he values of he C M A D saisics (p-values beween parenheses for he failure imes daa. Thus, according o hese formal ess, he BSG model fis he curren daa beer han he oher models. These resuls illusrae he poenialiy of he BSG disribuion he imporance of he addiional parameer. Table 3: Goodness-of-fi ess Model Saisics C M A D BSG ( a ( a BSP ( ( BSL ( ( BS ( ( a Denoes he p-value of he es. As a second applicaion, we shall analyze a daa se from McCool (1974 on he faigue lives (in hours of 10 bearings of a cerain kind. The daa are: 152.7, 172.0, 172.5, 173.3, 193.0, 204.7, 216.5, 234.9, 262.6, Table 4 gives he MLEs (sard errors beween parenheses of he parameers, he values of 2l( Θ of he Kolmogorov-Smirnov, AIC BIC saisics for he fied BSG, BSP BS models o he second daa se. Plos of he fied BSG, BSP BS densiies are given in Figure 5. Table 4: Parameer esimaes, K-S saisics, AIC BIC for faigue life daa Disribuion α θ K S 2l( Θ AIC BIC BSG ( a ( ( BSP ( a ( ( BS ( a ( a Denoes he sard deviaion of he MLE s of α, θ. Densiy BSG BSP BS cdf BSG BSP BS (a (b Figure 5: Esimaed (a pdf (b cdf for he BSG, BSP BS models for faigue lives daa. 13

14 The figures in Table 4 indicae ha he BSG model yields a beer fi o hese daa han he oher models. Table 5 gives he values of he C M A D saisics (p-values beween parenheses for he faigue life daa. So, hese values indicae ha he null hypohesis is srongly no rejeced for he BSG disribuion, whereas he null hypohesis is weakly no rejeced for he oher models. Thus, according o hese goodness-of-fi ess, he BSG model fis he curren daa beer han he oher models. Table 5: Goodness-of-fi ess Model Saisics C M A D BSG ( a ( a BSP ( ( BS ( ( a Denoes he p-value of he es. 7. Concluding remarks The wo-parameer Birnbaum-Saunders disribuion is exended o a new class of faigue life disribuions by compounding he Birnbaum-Saunders power series disribuions, hus creaing he Birnbaum-Saunders power series disribuion. As expeced, he new class includes as special model he BS disribuion. The BSPS failure rae funcion can be increasing upside-down bahub shaped. Some mahemaical properies of he new class are provided such as he ordinary momens, quanile funcion, order saisics heir momens. The esimaion of he model parameers is approached by he mehod of maximum likelihood he observed informaion marix is derived. We fi some BSPS disribuions o wo real daa ses o show he poenialiy of he proposed class. Fuure research could be adressed o sudy he complemenary Birnbaum-Saunders power series disribuion. This class of disribuions is a suiable model in a complemenary risk problem based in he presence of laen risks which arise in several areas such as public healh, acuarial science, biomedical sudies, demography indusrial reliabiliy. Acknowledgemen The financial suppor from he Brazilian governmenal insiuions CNPq CAPES is graefully acknowledged. Appendix A. Elemens of he observed informaion marix The observed informaion marix J(Θ for he parameers θ, α is given by U J(Θ = 2 l(θ θθ U θα U θ Θ Θ = U αα U α, U 14

15 whose elemens are U θθ = n θ 2 n [ C (θ C(θ U θα = 1 α 2 θ α 2 n n U θ = 1 2α ( C ] 2 (θ C(θ [ C Φ(υ i 2 ] 2 (θφ(υ i C, (θφ(υ i φ(υ i ρ( i / C α (θφ(υ C (θφ(υ + θ α 2 φ(υ i Φ(υ i τ( i / C α (θφ(υ C (θφ(υ, C (θφ(υ C (θφ(υ θ 2α n ( i τ 2θ α 3 U αα = n α 2 ( α 2 3 α 4 ( 2 θ n α 2 φ(υ i ρ U α = 1 α 3 + θ α 2 n ( i i ( i C α (θφ(υ C (θφ(υ, φ(υ i τ( i / C α (θφ(υ C (θφ(υ C θ (θφ(υ C (θφ(υ + θ 2α C α (θφ(υ C (θφ(υ + θ α 2 C α (θφ(υ C (θφ(υ C (θφ(υ 2, θ 2α 3 ( i υ i φ(υ i + i Φ(υ i C α (θφ(υ C (θφ(υ C (θφ(υ, C αα (θφ(υ C (θφ(υ C α (θφ(υ C (θφ(υ C (θφ(υ U = + n α 2 2 θ2 2α 2 θ (2α 2 + θ 2α φ(υ i 2 τ ( i + i ( i υ i φ(υ i τ φ(υ i τ ( i 1 2α 2 2 ( i 2 2 τ 1 i 2 i ( i + ( i C (θφ(υ C (θφ(υ, ( i i i θ C (θφ(υ 2α 2 C (θφ(υ C (θφ(υ C (θφ(υ 2 C (θφ(υ C (θφ(υ + θ 4α 2 ( i C (θφ(υ φ(υ i ρ C (θφ(υ where C k ( is he firs derivaive of C ( wih respec o k C jk ( corresponds o second derivaive wih respec o j k. References References Adamidis, K. Loukas, S. (1998. A lifeime disribuion wih decreasing failure rae. Saisics Probabiliy Leers, 39, Baraka, H.M. Abdelkader, Y.H. (2004. Compuing he momens of order saisics from nonidenical rom variables. Sa. Mehods Appl., 13, Basu, A., Klein, J. (1982. Some recen developmen in compeing risks heory, vol 1. Proceedings on Survival Analysis. Birnbaum Z.W. Saunders S.C. (1969. A new family of life disribuions. Journal of Applied Probabiliy, 6, Cancho, V.G., Louzada, F., Barriga. G.D.C. (2012. The Geomeric Birnbaum-Saunders regression model wih cure rae. Journal of Saisical Planning Inference, 142, Chahki, M. Ganjali, M. (2009. On some lifeime disribuions wih decreasing failure rae. Compuaional Saisics Daa Analysis, 53, Chen, G. Balakrishnan, N. (1995. A general purpose approximae goodness-of-fi es. Journal of Qualiy Technology, 27, Cordeiro, G.M. Lemone, A.J. (2011. The -Birnbaum-Saunders disribuion: An improved disribuion for faigue life modeling. Compuaional Saisics Daa Analysis, 55,

16 Greenwood, J.A., Lwehr, J.M., Maalas, N.C. Wallis, J.R. (1979. Probabiliy weighed momens: definiion relaion o parameers of severaldisribuions expressable in inverse form. Waer Resources Research, 15, Kundu, D., Kannan, N., Balakrishnan, N. (2008. On he hazard funcion of Birnbaum-Saunders disribuion associaed inference. Compuaional Saisics Daa Analysis, 52, Kus, C. (2007. A new lifeime disribuion. Compuaional Saisics Daa Analysis, 51, Mahmoudi, E., Jafari, A.A. (2012. Generalized exponenial-power series disribuions. Compuaional Saisics Daa Analysis, 56, Marshall, A.N, Olkin, I. (1997. A new mehod for adding a parameer o a family of disribuions wih applicaions o he exponenial Weibull families. Biomerica, 84, McCool, J.I., Inferenial echniques for Weibull populaions. Aerospace Research Laboraories Repor ARL TR Wrigh Paerson Air Force Base, Dayon, OH. Morais, A.L., Barreo-Souza, W. (2011. A Compound Class of Weibull Power Series, Disribuions. Compuaional Saisics Daa Analysis, 55, Murhy, D.N.P., Xie, M., Jiang, R. (2004. Weibull Models, (Wiley Series in Probabiliy Saisics. John Wiley Sons, Inc., New Jersey. Ng, H.K.T., Kundu, D., Balakrishnan, N. (2006. Poin inerval esimaion for he wo-parameer Birnbaum-Saunders disribuion based on Type-II censored samples. Compuaional Saisics Daa Analysis, 50, Tahmasbi, R., Rezaei, S. (2008. A wo-parameer lifeime disribuion wih decreasing failure rae. Compuaional Saisics Daa Analysis, 52, Saunders, S.C. (1974. A family of rom variables closed under reciprocaion. Journal of he American Saisical Associaion, 69, Silva, R.B., Bourguignon, M., Dias, C.R.B., Cordeiro, G.M. (2013. The compound family of exended Weibull power series disribuions. Compuaional Saisics Daa Analysis 58,