The Iverse Weghted Ldley Dstrbuto: Propertes, Estmato ad a Applcato o a Falure Tme Data Pedro L. Ramos 1 Fracsco Louzada 1 Tacaa K.O. Shmzu 1 Ale O. Luz 1 arxv:1711.09338v1 [math.st] 26 Nov 2017 1 Isttute of Mathematcal Scece ad Computg, Uversty of São Paulo, São Carlos, Brazl ARTICLE HISTORY Compled November 28, 2017 ABSTRACT I ths paper a ew dstrbuto s proposed. Ths ew model provdes more flexblty to modelg data wth upsde-dow bathtub hazard rate fucto. A sgfcat accout of mathematcal propertes of the ew dstrbuto s preseted. The maxmum lkelhood estmators for the parameters the presece of complete ad cesored data are preseted. Two correctve approaches are cosdered to derve modfed estmators that are bas-free to secod order. A umercal smulato s carred out to exame the effcecy of the bas correcto. Fally, a applcato usg a real data set s preseted order to llustrate our proposed dstrbuto. KEYWORDS Iverse weghted Ldley dstrbuto; Maxmum Lkelhood Estmato; Bas correcto; Radom cesorg. 1. Itroducto I recet years, several ew dstrbutos have bee troduced lterature for descrbg real problems. A mportat dstrbuto was preseted by Ldley [16] the cotext of fducal statstcs ad Bayes theorem. Ghtay et al. [11] argued that the Ldley dstrbuto provdes flexble mathematcal propertes ad outled that may cases ths dstrbuto outperforms the expoetal dstrbuto. Sce the, ew geeralzatos of Ldley dstrbuto have bee proposed such as the geeralzed Ldley [26], exteded Ldley [3], ad Power Ldley [9] dstrbuto. The study of weght dstrbutos provde ew compreheso of stadard dstrbutos ad cotrbutes addg more flexblty for fttg data [18]. Ghtay et al. [10] preseted a two-parameter weghted Ldley (WL dstrbuto whch has bathtub ad creasg hazard rate. The WL dstrbuto has probablty desty fucto (PDF gve by f(t φ, λ = λ φ+1 (φ + λγ(φ tφ 1 (1 + te λt, (1 for all t > 0, φ > 0 ad λ > 0 where Γ(φ = 0 e x x φ 1 dx s the gamma fucto. Mazuchel et al. [17] compared the fte sample propertes of the parameters of the WL Correspodg author: Pedro Luz Ramos, Emal: pedrolramos@usp.br
dstrbuto umercal smulatos usg four methods. Wag ad Wag [25] preseted bas-corrected MLEs ad argued that the proposed estmators are strogly recommeded over other estmators wthout bas-correcto. Al [2] cosdered a Bayesa approach ad derved several formatve ad oformatve prors uder dfferet loss fuctos. Ramos ad Louzada [20] troduced three parameters geeralzed weghted Ldley dstrbuto. I ths study, a ew two-parameter dstrbuto wth upsde-dow bathtub hazard rate s proposed, hereafter, verse weghted Ldley (IWL dstrbuto. Ths ew model ca be rewrtte as the verse of the WL dstrbuto. A sgfcat accout of mathematcal propertes for the IWD dstrbuto s preseted such as momets, survval propertes ad etropy fuctos. The maxmum lkelhood estmators of the parameters ad ts asymptotc propertes are obtaed. Further, two correctve approaches are dscussed to derve modfed MLEs that are bas-free to secod order. The frst has a aalytcal expresso derved by Cox ad Sell (12 ad the secod s based o the bootstrap resamplg method (see Efro [8] for more detals, whch ca be used to reduce bas. Smlar correctve approaches has bee cosdered by may authors for other dstrbutos, e.g., Cordero et al. [5], Lemote [15], Temour ad Nadarajah [24], Gles et al. [12], Ramos et al. [19], Schwartz et al. [22] ad Reath et al. [21]. I addto, the MLEs the presece of radomly cesored data s preseted. Approxmated bas-corrected MLEs for cesored data are also dscussed. A umercal smulato s performed to exame the effect of the bas correctos the MLEs for complete ad cesored data. The ew dstrbuto s a useful geeralzato of the verse Ldley dstrbuto [23] ad ca be represeted by a two-compoet mxture model. Mxture models play a mportat role statstcs for descrbg heterogeety (see, Aale [1]. Therefore, the IWL dstrbuto ca be used to descrbe data sets the presece of heterogeety. For stace, we ca be terested descrbg the lfetme of compoets that are composed of ew ad repared products, however, oly the falure tme s observed ad the groups are latet varables. I ths case, the proposed dstrbuto, as a mxture dstrbuto, ca express the heterogeety the data. I relablty, ths model may be used to descrbe the lfetme of compoets assocated wth a hgh falure rate after short repar tme. I studes volvg the lfetme of patets ths model ca be useful to descrbe the course of a dsease, where ther mortalty rate reaches a peak ad the decles as the tme crease,.e., problems where ther hazard fucto has upsde-dow bathtub shape. I order to llustrate our proposed methodology, we cosdered a real data set related to falure tme of devces of a arle compay. Such study s mportat order to prevet customer dssatsfacto ad customer attrto, ad cosequetly to avod customer loss. I ths cotext, the choce of the dstrbuto that fts better ths data s fudametal for the compay reduces ts costs. We showed that the verse weghted Ldley dstrbuto fts better tha other well-kow dstrbutos for ths data set. The paper s orgazed as follows. Secto 2 troduces the verse weghted Ldley dstrbuto. Secto 3 presets the propertes of the IWL dstrbuto such as momets, survval propertes ad etropy. Secto 4 dscusses the feretal procedure based o MLEs for complete ad cesored data. A bas correcto approach s also preseted for complete ad cesored data. Secto 5 descrbes two correctve approaches to reduce the bas the MLEs for complete ad cesored data. Secto 6 presets a smulato study to verfy the performace of the proposed estmators. Secto 7 llustrates the relevace of our proposed methodology a real lfetme data. Secto 8 summarzes the preset study. 2
2. Iverse Weghted Ldley dstrbuto A o-egatve radom varable T follows the IWL dstrbuto wth parameters φ > 0 ad λ > 0 f ts PDF s gve by f(t φ, λ = λ φ+1 ( (φ + λγ(φ t φ 1 1 + 1 e λt 1. (2 t Note that f φ = 1, the IWL dstrbuto reduces to the verse Ldley dstrbuto [23]. The IWL dstrbuto ca be expressed as a two-compoet mxture f(t φ, λ = pf 1 (t φ, λ + (1 pf 2 (t φ, λ, where p = λ/(λ + φ ad T j IG(φ + j 1, λ, for j = 1, 2,.e., f j (t λ, φ s Iverse Gamma dstrbuto, gve by f j (t φ, λ = λ φ+j 1 Γ(φ + j 1 t φ j e λt 1. Therefore, the IWL dstrbuto s a mxture dstrbuto ad ca express the heterogeety the data. Proposto 2.1. Let T IWL(φ, λ the X = 1/T follows a weghted Ldley dstrbuto [10]. Proof. Defe the trasformato X = g(t = 1 T the the resultg trasformato s ( f X (x = f T g 1 (x d dx g 1 (x = λ φ+1 (φ + λγ(φ xφ+1 (1 + x e λx 1 x 2 = λ φ+1 (φ + λγ(φ xφ 1 (1 + x e λx. Fgure 1 gves examples from the shapes of the desty fucto for dfferet values of φ ad λ. The cumulatve dstrbuto fucto from the IWL dstrbuto s gve by F (t φ, λ = Γ ( φ, λt 1 (λ + φ + (λt 1 φ e λt 1, (λ + φγ(φ where Γ(x, y = x wy 1 e x dw s the upper complete gamma. 3
f(t 0 2 4 6 φ=3, λ=1.5 φ=5, λ=1.5 φ=6, λ=1.5 φ=7, λ=1.5 φ=8, λ=1.5 f(t 0.0 0.5 1.0 1.5 2.0 φ=0.1, λ=0.5 φ=0.3, λ=0.5 φ=0.5, λ=0.5 φ=0.8, λ=0.5 φ=1.0, λ=0,5 0.0 0.2 0.4 0.6 0.8 1.0 1.2 0.0 0.5 1.0 1.5 2.0 t t Fgure 1. Desty fucto shapes for IWL dstrbuto ad cosderg dfferet values of φ ad λ. 3. Propertes of IWL Dstrbuto I ths secto, we provde a sgfcat accout of mathematcal propertes of the ew dstrbuto. 3.1. Momets Momets play a mportat role statstcs. They ca be used may applcatos, for stace the frst momet of the PDF s the well kow mea, whle the secod momet s used to obta the varace, skewess ad kurtoss are also obtaed from the momets. I the followg, we wll derve the momets for the IWL dstrbuto. Proposto 3.1. For the radom varable T wth IWL dstrbuto, the r-th momet s gve by µ r = E[T r ] = λ r (φ + λ r, where φ > r. (3 (λ + φ(φ 1(φ 2 (φ r Proof. Note that f W IG(φ, λ dstrbuto the the r-th momet from the radom varable W s gve by E (φ,λ [W r ] = λr Γ(φ r Γ(φ = λ r, where φ > r. (φ 1(φ 2... (φ r 4
Sce the IWL dstrbuto ca be expressed as a two-compoet mxture, we have µ r = E[T r ] = = = 0 t r f(t φ, λdt = pe (φ,λ [W r ] + (1 pe (φ+1,λ [W r ] λ Γ(φ r φ Γ(φ + 1 r + = λr (λ + φ rγ(φ r (λ + φ Γ(φ (λ + φ Γ(φ + 1 (λ + φγ(φ λ r (φ + λ r, where φ > r. (λ + φ(φ 1(φ 2 (φ r Proposto 3.2. The r-th cetral momet for the radom varable T s gve by M r = E[T µ] r = = r =0 ( ( r r =0 ( r ( µ r E[T ] r ( λ(φ + λ 1 (λ + φ(φ 1 λ (φ + λ (λ + φ(φ 1(φ 2 (φ Proof. The result follows drectly from the proposto 3.1. Proposto 3.3. A radom varable T wth IWL dstrbuto, has the mea ad varace gve by µ = λ(φ + λ 1 (λ + φ(φ 1,. (4 σ 2 = λ ((φ 2 + λ 2(φ 1 (φ + λ 1 2 (φ 2 (λ + φ(φ 2(φ 1 2. Proof. From (3 ad cosderg r = 1, t follows that µ 1 = µ. The secod result follows from (4 cosderg r = 2 ad wth some algebrac operato the proof s completed. 3.2. Survval Propertes Survval aalyss has become a popular brach of statstcs wth wde rage of applcatos. Although may fuctos related to survval aalyss ca be derved for ths model, ths secto we wll preset the most commo fuctos. The survval fucto of IWL dstrbuto represetg the probablty of a observato does ot fal utl a specfed tme t s gve by S(t φ, λ = γ ( φ, λt 1 (λ + φ (λt 1 φ e λt 1, (λ + φγ(φ 5
where γ(y, x = x 0 wy 1 e w dw s the lower complete gamma fucto. The hazard fucto of T s gve by h(t φ, λ = λ φ+1 t φ 1 ( 1 + t 1 e λt 1 γ (φ, λt 1 (λ + φ (λt 1 φ e λt 1. (5 Ths model has upsde-dow bathtub hazard rate. The followg Lemma s useful to prove such result. Lemma 3.4. Glaser [13]: Let T be a o-egatve cotuous radom varable wth twce dfferetable PDF f(t, hazard rate fucto h(t ad η(t = t log f(t. The f η(t has a upsde-dow bathtub shape, h(t has a upsde-dow bathtub shape. Theorem 3.5. The hazard fucto (5 s upsde-dow bathtub for all φ > 0 ad λ > 0. Proof. For IWL dstrbuto we have t follows that η(t = φ t + 2 t 1 (t + 1 λ t 2, η (t = φ t 2 2 t 2 + 1 (t + 1 2 + 2λ t 3. The study of the behavour of η (t s ot smple. However usg the Wolfram Alpha software, we ca check that for all φ > 0 ad λ > 0, η (t s creasg (0, ξ(t φ, λ ad decreasg (ξ(t φ, λ,,.e., η (t = 0 at ξ(t φ, λ, where ξ(t φ, λ s a very large fucto computed to the Wolfram Alpha (avalable upo request. Therefore, η(t ad cosequetly h(t has upsde-dow bathtub shape. Ths propertes make the IWL dstrbuto a useful model for relablty data. Fgure 2 gves examples of dfferet shapes for the hazard fucto. h(t 0 5 10 15 φ=3, λ=1.5 φ=5, λ=1.5 φ=6, λ=1.5 φ=7, λ=1.5 φ=8, λ=1.5 h(t 0.0 0.5 1.0 1.5 2.0 2.5 φ=0.1, λ=0.5 φ=0.3, λ=0.5 φ=0.5, λ=0.5 φ=0.8, λ=0.5 φ=1.0, λ=0,5 0 1 2 3 4 5 0 1 2 3 4 5 t t Fgure 2. Hazard fucto shapes for IWL dstrbuto ad cosderg dfferet values of φ ad λ. 6
Proposto 3.6. The mea resdual lfe fucto r(t φ, λ of the IWL dstrbuto s gve by r(t φ, λ = 1 S(t t yf(y λ, φdy t = λγ ( φ, λt 1 + λ 2 γ ( φ, λt 1 γ (φ, λt 1 (λ + φ (λt 1 φ e λt 1 t. Proof. Note that, for the Iverse Gamma dstrbuto we have that t yf j (y φ, λdy = Usg the followg relatoshp r(t φ, λ = 1 [ p S(t t λ Γ(φ + j 1 γ [ φ + j 2, λt 1], j = 1, 2. yf 1 (y λ, φdy + (1 p ad after some algebrac mapulatos, the proof s completed. t ] yf 2 (y λ, φdy t, 3.3. Etropy I formato theory, etropy has played a cetral role as a measure of the ucertaty assocated wth a radom varable. Shao s etropy s oe of the most mportat metrcs formato theory. The Shao s Etropy from IWL dstrbuto s gve by solvg the followg equato H S (φ, λ = 0 ( log λ φ+1 (φ + λγ(φ t φ 1 ( 1 + 1 e λt 1 f(t φ, λdt. (6 t Proposto 3.7. A radom varable T wth IWL dstrbuto, has the Shao s Etropy gve by H S (φ, λ = log(λ + φ + log Γ(φ + λφ+1 Ω(φ, λ (λ + φγ(φ. φ(λ + φ + 1 (λ + φ where Ω(φ, λ = 0 (x + 1 log(x + 1x( φ 1e λx dx. Proof. From the equato (6 we have ( 1 (φ + 1 λ + φ + ψ(φ H S (φ, λ = (φ+1 log λ log(λ+φ log Γ(φ λe [ t 1 ] (φ+1e[log(t ]+E [ log(1 + t 1 ]. Sce E[log(t] = log(λ 1 λ + φ ψ(φ, ad E [ t 1 ] (φ + 1 = 1 λ λ + φ 7 φ(λ + φ + 1 =. λ(λ + φ
The E [ log(1 + t 1 ] = H S (φ, λ = (φ + 1 λ φ+1 (λ + φγ(φ 0 (x + 1 log(x + 1x ( φ 1e λx dx. ( 1 λ + φ + ψ(φ φ(λ + φ + 1 log(λ + φ log Γ(φ (λ + φ λφ+1 Ω(φ, λ (λ + φγ(φ. 4. Iferece I ths secto, we preset the maxmum lkelhood estmator of the parameters φ ad λ of the IWL dstrbuto. Addtoally, MLEs cosderg radomly cesored data are also dscussed. 4.1. Maxmum Lkelhood Estmato Amog the statstcal ferece methods, the maxmum lkelhood method s wdely used due to ts better asymptotc propertes. Uder the maxmum lkelhood method, the estmators are obtaed from maxmzg the lkelhood fucto. Let T 1,..., T be a radom sample such that T IWL(φ, µ. I ths case, the lkelhood fucto from (2 s gve by L(θ; t = { λ (φ+1 (φ + λ Γ(φ t φ 1 =1 } =1 { (1 + 1t exp λ The log-lkelhood fucto l(θ; t = log L(θ; t s gve by l(θ; t = (φ + 1 log λ log(λ + φ log Γ(φ λ =1 1 t (φ + 1 } 1. t =1 log(t. (7 From the expressos φ l(θ; t = 0, λl(θ; t = 0, we get the lkelhood equatos log(λ log(t λ + φ ψ(φ = 0, =1 =1 (φ + 1 λ =1 1 t λ + φ = 0, where ψ(k = k log Γ(k = Γ (k Γ(k s the dgamma fucto. After some algebrac 8
mapulato the soluto of λ MLE s gve by ˆλ MLE = ( ˆφ 2 ( MLE (ξ(t 1 + ˆφMLE (ξ(t 1 + 4 ξ(t ˆφ2 MLE + ˆφ MLE 2ξ(t, where ξ(t = =1 (t 1 ad ˆφ MLE ca be obtaed solvg the olear system log(ˆλ MLE log(t =1 ˆλ MLE + ˆφ MLE ψ( ˆφ MLE = 0. (8 These results are a smple modfcato of the results obtaed for Ghtay et al. [11] for the WL dstrbuto. Uder mld codtos the ML estmates are asymptotcally ormal dstrbuted wth a bvarate ormal dstrbuto gve by ( ˆφ, ˆλ N 2 [(φ, λ, I 1 (φ, λ] for, where the elemets of the Fsher formato matrx I(φ, λ are gve by h 11 (φ, λ = (λ + φ 2 + ψ (φ, h 12 (φ, λ = h 21 (φ, λ = λ (λ + φ 2, h 22 (φ, λ = (φ + 1 λ 2 (λ + φ 2, ad ψ (k = 2 k log Γ(k s the trgamma fucto. A terestg property of the IWL dstrbuto s that the observed matrx formato s equal to the expected formato matrx. 4.2. Radom Cesorg I survval aalyss ad dustral lfetme testg, radom cesorg schemes have bee receved specal atteto. Suppose that the th dvdual has a lfetme T ad a cesorg tme C, moreover the radom cesorg tmes C s are depedet of T s ad ther dstrbuto does ot deped o the parameters, the the data set s (t, δ, where t = m(t, C ad δ = I(T C. Ths type of cesorg have as specal case the type I ad II cesorg mechasm. The lkelhood fucto for θ s gve by L(θ, t = f(t θ δ S(t θ 1 δ. =1 Let T 1,, T be a radom sample of IWL dstrbuto, the lkelhood fucto 9
cosderg data wth radom cesorg s gve by λ d(φ+1 ( L(λ, φ t = (λ + φ Γ(φ =1 ( (1 + t 1 t φ 1 (λ + φγ(φ, λt 1 e λt 1 δ. ( λt 1 φ e λt 1 1 δ (9 The logarthm of the lkelhood fucto (9 s gve by l(λ, φ t = (φ + 1 + + δ log(t λ =1 ( (1 δ log =1 =1 δ log(1 + t 1. =1 (λ + φγ(φ, λt 1 δ t 1 + d(φ + 1 log(λ log(φ + λ (λt 1 φ e λt 1 log (Γ(φ From l(λ, φ t/ λ = 0 ad l(λ, φ t/ φ = 0, the lkelhood equatos are gve as follows ( (1 δ γ(φ, λt 1 + (λ + φ ( λt 1 φ 1 e λt 1 φλ φ 1 t φ e λt 1 ( λt 1 φ+1 e λt 1 ( =1 (λ + φγ(φ, λt 1 ( λt 1 φ = e λt 1 d(φ + 1 + λ + φ λ =1 δ t 1, (10 =1 ( (1 δ γ(φ, λt 1 + (λ + φψ(φ, λt 1 ( λt 1 ( (λ + φγ(φ, λt 1 ( λt 1 φ e λt 1 φ log(λt 1 e λt 1 = d log(λ + λ + φ + ψ(φ + δ log(t 1, (11 where Ψ(k, x = γ(k, x/ k ca be computed umercally. Numercal methods are requred order to fd the soluto of these o-lear equatos. =1 5. Bas correcto for the maxmum lkelhood estmators I ths secto, we dscuss modfed MLEs based o two correctve approaches that are bas-free to secod order. Frstly a correctve aalytcal approach s preseted tha the bootstrap resamplg method s preseted. 10
5.1. A correctve approach Cosder the lkelhood fucto L(θ; t wth a p-dmesoal vector of parameters θ. Thus, the jot cumulats of the dervatves of l(θ; t ca be wrtte by ( 2 l(θ; t h j (θ = E θ θ j ( 2 l(θ; t h j,l (θ = E. θ θ j ( 3 l(θ; t, h jl (θ = E θ θ j θ l ad l(θ; t, for, j, l = 1,..., p. θ l Cosequetly, the dervatves of such cumulats are gve by h (l j (θ = h j(θ θ l, for, j, l = 1,..., p. The bas of θ m studed by Cox ad Sell [7] for depedet sample wthout ecessarly be detcally dstrbuted ca be wrtte by Bas(ˆθ m = p p p s m (θs jl (θ (h j,l (θ + 0.5h jl (θ + O( 2, (12 =1 j=1 k=1 where s j s the (, j-th elemet of the verse of Fsher s formato matrx of ˆθ, K = { h j }. Cordero ad Kle [6] proved that eve f the data are depedet the expresso (12 ca be re-wrtte as Bas(ˆθ m = p s m (θ =1 p p j=1 k=1 ( s jl (θ h (l j (θ 0.5h jl(θ + O( 2. (13 Let a l j = h(l j 1 2 h(l j ad defe the matrx A = [A(1 A (2... A (p ] wth A (l = {a (l j }, for, j, l = 1,..., p. Thus, the expresso for the bas of ˆθ ca be expressed as A bas corrected MLE for ˆθ s obtaed as Bas(ˆθ m = K 1 A.vec(K 1 + O( 2. (14 ˆθ CMLE = ˆθ K 1 A.vec(K 1, (15 where ˆθ s the MLE of the parameter θ, ˆK = K θ=ˆθ ad  = A θ=ˆθ. The bas of ˆθ CMLE s ubased O( 2. For the IWL dstrbuto the hgher-order dervatves ca be easly obtaed sce they do ot volve t, thus, we have h 111 (θ = h (1 2 11 (θ = (λ + φ 3 ψ (φ, h 122 (θ = h 221 (θ = h 212 (θ = h (2 12 h 222 (θ = h (2 2 2(φ + 1 22 (θ = 3 (λ + φ λ 3 2 (θ = h(1 22 (θ = (λ + φ 3 λ 2, ad 11
where ψ (k = h 211 (θ = h 112 (θ = h 121 (θ = h (1 2 12 (θ = h(2 11 (θ = (λ + φ 3, 3 k log Γ(k. The matrx K s gve by K = [ (λ+φ 2 ψ (φ λ + (λ+φ 2 λ + (λ+φ 2 (φ+1 λ 2 + (λ+φ 2 To obta the matrx A of (14, we preset the elemets of A (1 ad the elemets of A (2 are a (1 11 = h(1 11 1 2 h 111 = (λ + φ 3 ψ (φ, 2 a (1 12 = a(1 21 = h(1 12 1 2 h 112 = (λ + φ 3, a (1 22 = h(1 22 1 2 h 221 = (λ + φ 3 2λ 2, a (2 11 = h(2 11 1 2 h 112 = (λ + φ 3, a (2 12 = a(2 21 = h(2 12 1 2 h 122 = (λ + φ 3 2λ 2, a (2 22 = h(2 22 1 2 h (φ + 1 222 = (λ + φ 3 λ 3. Thus, the matrx A = [A (1 A (2 ] s expressed by A = ( 1 (λ+φ ψ (φ 3 2 1 (λ+φ 1 3 (λ+φ 1 3 (λ+φ 1 3 2λ 2 1 (λ+φ 3 1 (λ+φ 3 1 2λ 2 1 (λ+φ 3 1 2λ 2 1 (λ+φ 3 Fally, the bas-corrected maxmum lkelhood estmators are gve by ]. (φ+1 λ 3 ( ( ˆφCMLE ˆφMLE = ˆλ CMLE ˆλ ˆK 1 Â.vec( ˆK 1, (16 MLE. where ˆK = K φ= ˆφ,λ=ˆλ ad  = A φ= ˆφ,λ=ˆλ. It s mportat to pot out that, sce the hgher-order do ot volve t, they are the same of the WL dstrbuto [25]. A bas corrected approach ca be cosdered for cesored data. Although the Fsher formato matrx related to the MLEs (9 does ot preset closed-form expressos, we ca cosder the bas corrected preseted (5.1. I ths case, approxmated bascorrected maxmum lkelhood estmates (ACMLE are archved by ( ( ˆφACMLE ˆφMLE = ˆλ ACMLE ˆλ ˆK 1 Â.vec( ˆK 1, MLE where ˆK = K φ= ˆφMLE,λ=ˆλ MLE,  = A φ= ˆφMLE,λ=ˆλ MLE ad ˆφ MLE ad ˆλ MLE are the 12
solutos of (10 ad (11. However, the bas of ˆθ ACMLE s ot a ubased estmator wth O( 2. 5.2. Bootstrap resamplg method I what follows we cosder the bootstrap resamplg method proposed by Efro [8] to reduce the bas of the MLEs. Such method cossts geeratg pseudo-samples from the orgal sample to estmate the bas of the MLEs. Thus, the bas-corrected MLEs s gve by subtracto of the estmated bas wth the orgal MLEs. Let y = (y 1,..., y be a sample wth observatos radomly selected from the radom varable Y whch has the dstrbuto fucto F = F ν (y. Thus, let the parameter ν be a fucto of F gve by ν = t(f. Fally, let ˆν be a estmator of ν based o y,.e., ˆν = s(y. The pseudo-samples y = (y 1,..., y s obtaed from the orgal sample y through resamplg wth replacemet. The bootstrap replcates of ˆν s calculated, where ˆν = s(y ad the emprcal cdf (ecdf of ˆν s used to estmate Fˆν (cdf of ˆν. Let B F (ˆν, ν be the bas of the estmator ˆν = s(y gve by B F (ˆν, ν = E F [ˆν, ν] = E F [s(y] ν(f. Note that the subscrpt of the expectato F dcates that s take wth respect to F. The bootstrap estmators of the bas were obtaed by replacg F wth Fˆν, where F geerated the orgal sample. Therefore, the bootstrap bas estmate s gve by ˆB Fˆν (ˆν, ν = E Fˆν [ ˆν ] ˆν. If we have B bootstrap samples (y (1, y (2,..., y (B whch are geerated depedetly from the orgal sample y ad the respectve bootstrap estmates (ˆν (1, ˆν (2,..., ˆν (B are calculated, the t s achevable to determe the bootstrap expectatos E Fˆν [ ˆν ] approxmately by ˆν (. = 1 B B ˆν (. =1 Therefore, the bootstrap bas estmate based o B replcatos of ˆν s ˆB F (ˆν, ν = ˆν (. ˆν, whch results the bas corrected estmators obtaed through by bootstrap resamplg method that s gve by ν B = ˆν ˆB F (ˆν, ν = 2ˆν ˆν (.. I our case, we have ν B deoted by ˆθ BOOT = ( ˆφ BOOT, ˆλ BOOT. 6. Smulato Aalyss I ths secto a smulato study s preseted to compare the effcecy of the maxmum lkelhood method ad the bas correcto approaches the presece of complete ad cesored data. The proposed comparsos are performed by computg the mea 13
relatve errors (MRE ad the relatve mea square errors (RMSE gve by MRE = 1 N N j=1 ˆθ,j θ, RMSE = 1 N N (ˆθ,j θ 2, for = 1, 2, j=1 θ 2 where N s the umber of estmates obtaed through the MLE, CMLE ad the bootstrap approach. The 95% coverage probablty of the asymptotc cofdece tervals are also evaluated. Cosderg ths approach, we expected that the most effcet estmato method returs the MREs closer to oe wth smaller RMSEs. Moreover, for a large umber of expermets, usg a 95% cofdece level, the frequeces of tervals that covered the true values of θ should be closer to 95%. Followg Reath et al. [21] we used B=1,000 for the bootstrap method. The programs ca be obtaed, upo request. The radom sample of the IWL were geerated cosderg the followg algorthm: (1 Geerate U Uform(0, 1, = 1,..., ; (2 Geerate X IG(φ, λ, = 1,..., ; (3 Geerate Y IG(φ + 1, λ, = 1,..., ; (4 If U p = λ/(λ + φ, the set T = X, otherwse, set T = Y, = 1,...,. 6.1. Complete Data The smulato study s performed cosderg the values: θ = ((0.5, 2,(2, 4, N = 30, 000 ad = (20, 25,..., 130. It s mportat to pot out that, smlar results were acheved for dfferet choces of φ ad λ. The uroot procedure avalable R s cosdered to fd the soluto of the o-lear equato (8. The bas correcto s computed drectly from (16. Fgures 3 ad 4 preset the MRE, RMSE ad the coverage probablty wth a 95% cofdece level related to the MLE, CMLE ad the bootstrap uder dfferet values of. From Fgures 3 ad 4, we observed that the estmates of φ ad λ are asymptotcally ubased,.e., the MREs ted to oe whe creases ad the RMSEs decrease to zero for large. The CMLE preset superor performace tha the bootstrap approach for both parameters for ay sample szes. Takg to accout the results of the smulato studes, the maxmum lkelhood estmators combed wth the correctve bas approach dscussed Secto 5.1 should be cosdered for estmatg the parameters of the IWL dstrbuto. 6.2. Cesored Data I ths secto, we cosdered the MLES the presece of radom cesored data. The cesored data s geerated followg the same procedure preseted by Goodma et al. [14]. I our case, we preseted two scearos where we obtaed approxmately 0.3 ad 0.5 proportos of cesored data,.e., 30% ad 50% of cesorshp. The smulato study s performed cosderg θ = (2, 4, N = 2, 000 ad = (10, 15,..., 130 The maxmum lkelhood estmates were computed usg the log-lkelhood fuctos (4.2 wth the maxlk package avalable R. The soluto for the maxmum was uque for all tal values. Fgures 5 ad 6 preset the MRE, the RMSE ad the coverage probablty wth a 95% cofdece level related to the MLE, CMLE ad the bootstrap uder dfferet values of. 14
MRE (φ 0.90 1.00 1.10 1.20 MSE (φ 0.0 0.5 1.0 1.5 2.0 CP (φ 0.90 0.92 0.94 0.96 MLE CLME Bootstrap MRE (λ 0.90 1.00 1.10 1.20 MSE (λ 0 1 2 3 4 CP (λ 0.90 0.92 0.94 0.96 Fgure 3. MREs, MSEs related to the estmates of λ = 2 ad φ = 4 for N = 30, 000 smulated samples, cosderg dfferet values of. MRE (φ 0.90 1.00 1.10 1.20 MSE (φ 0.00 0.05 0.10 0.15 CP (φ 0.90 0.92 0.94 0.96 MLE CLME Bootstrap MRE (λ 0.90 1.00 1.10 1.20 MSE (λ 0 1 2 3 4 CP (λ 0.90 0.92 0.94 0.96 Fgure 4. MREs, MSEs related to the estmates of λ = 0.5 ad φ = 2 for N = 30, 000 smulated samples, cosderg dfferet values of. 15
MRE (φ 0.90 1.00 1.10 1.20 MSE (φ 0.0 0.5 1.0 1.5 2.0 CP (φ 0.90 0.92 0.94 0.96 MLE CLME Bootstrap MRE (λ 0.90 1.00 1.10 1.20 MSE (λ 0 1 2 3 4 CP (λ 0.90 0.92 0.94 0.96 Fgure 5. MREs, RMSEs related to the estmates of λ = 2 ad φ = 4 for N = 10, 000 smulated samples, cosderg dfferet values of ad 30% of cesorshp. MRE (φ 0.90 1.00 1.10 1.20 MSE (φ 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 CP (φ 0.90 0.92 0.94 0.96 MLE CLME Bootstrap MRE (λ 0.90 1.00 1.10 1.20 MSE (λ 0 1 2 3 4 CP (λ 0.90 0.92 0.94 0.96 Fgure 6. MREs, RMSEs related to the estmates of λ = 2 ad φ = 4 for N = 2, 000 smulated samples, cosderg dfferet values of ad 50% of cesorshp. 16
As show Fgures 5 ad 6 the proposed ACMLE retured more accurate estmates for both parameters whe compared wth the bootstrap approach or the MLEs. Takg to accout the results of the smulato studes, the approxmated corrected bas approach combed wth the maxmum lkelhood estmators should be cosder for estmatg the parameters of the IWL dstrbuto the presece of cesorshp. 7. Applcato I ths secto, recall the real data set brefly preseted Secto 1. The aalyze of the dstrbuto that better ft the proposed data s relevat to avod hgher costs for the compay. Table 1 presets the data related to falure tme of ( days of 194 devces a arcraft (+ dcates the presece of cesorshp. Table 1. Data set related to the falure tme of 194 devces a arcraft. 43 29 37 88 5 14 9 43+ 1 78 1 77 17 100 3 119+ 22 3 8 80 1 19 157+ 65 34 13 62+ 2 1 1 2 3 6 1 2 5 7 6 1 1 4 1 1 1 2 7 2 1 1 2 1 1 7 1 1 4 1 4 2 4 5 5 4 3 2 2 2 3 3 9 1 6 9 2 5 7 4 2 1 2 2 3 11 8 3 1 2 2 2 2 2 1 3 20+ 8 8 197 20 14 7 29 7 16 34 25 10 80 42 32 1 3 1 12 7 7 39+ 60 53 32 9 8 1 1 27 2 4 8 13 7 7 1 19 7 12 19 5 18 1 4 18 20 9 14 13 70 18 3 7 20 3 11 10 3 38+ 278 13 79 145+ 19 2 18 2 65 14 31 10 19 5 9 45 13 5 1 1 31 35 34 4 3 5 12 140+ 106 5 40 130+ 21 19 7 10 91 193 64 85+ The results obtaed from the IWL dstrbuto were compared to the Webull, Gamma, Logormal, Logstc, Iverse Webull ad Iverse Ldley dstrbuto ad the oparametrc survval curve adjusted usg the Kapla-Meer estmator. Itally, order to verfy the behavor of the emprcal hazard fucto t wll be cosdered the TTT-plot (total tme o test proposed by Barlow ad Campo [4]. The TTT-plot s acheved through the cosecutve plot of the values [r/, G(r/] where G(r/ = ( r =1 t + ( rt (r / =1 t, r = 1,...,, = 1,...,, ad t s the order statstcs. If the curve s cocave (covex, the hazard fucto s creasg (decreasg, whe t starts covex ad the cocave (cocave ad the covex the hazard fucto wll have a bathtub (verse bathtub shape. Dfferet dscrmato crtero methods based o log-lkelhood fucto evaluated at the MLEs were also cosdered. The dscrmato crtero methods are respectvely: Akake formato crtero (AIC computed through AIC = 2l( ˆθ; x + 2k, Corrected Akake formato crtero AICC = AIC +[2 k (k + 1/( k 1], Haa-Qu formato crtero HQIC = 2 l( ˆθ; x + 2 k log (log( ad the cosstet Akake formato crtero CAIC = AIC +k log( k, where k s the umber of parameters to be ftted ad ˆθ the estmates of θ. The best model s the oe whch provdes the mmum values of those crtera. Sce the data has radom cesorg mechasm, cosequetly the equatos (10 17
ad (11 were used to compute the MLEs. Table 2 dsplays the MLEs, stadard-error ad 95% cofdece tervals for φ ad λ. Table 3 presets the results of AIC, AICC, HQIC, CAIC crtera, for dfferet probablty dstrbutos. Table 2. MLE, Stadard-error ad 95% cofdece tervals for φ ad λ θ MLE SE CI 95% (θ φ 0.643 0.059 (0.527; 0.760 λ 2.825 0.296 (2.245; 3.405 Table 3. Results of AIC, AICC, HQIC, CAIC crtera for dfferet probablty dstrbutos cosderg the data set related to the falure tme of 194 of devces a arcraft. Test IW. Ldley Webull Gamma Logormal Logstc I Webull I Ldley AIC 1392.66 1452.37 1474.44 1408.44 1818.42 1392.70 1418.75 AICC 1392.73 1452.43 1474.50 1408.51 1818.48 1392.76 1416.78 HQIC 1395.31 1455.02 1477.08 1411.09 1821.06 1395.34 1418.08 CAIC 1401.20 1460.91 1482.97 1416.98 1826.95 1401.23 1427.29 Fgure 7 presets the TTT-plot, the survval fucto adjusted by dfferet dstrbutos ad the Kapla-Meer estmator ad the hazard fucto adjusted by the IWL dstrbuto. G(r/ 0.0 0.2 0.4 0.6 0.8 1.0 S(t 0.0 0.2 0.4 0.6 0.8 1.0 Kapla Meer IW Ldley Webull Gamma Logormal Logstc I. Webull I. Ldley h(t 0.00 0.05 0.10 0.15 0.0 0.2 0.4 0.6 0.8 1.0 r/ 0 50 100 150 200 250 Tme 0 50 100 150 200 250 Tme Fgure 7. TTT-plot, survval fucto adjusted by dfferet dstrbutos ad the KaplaMeer estmator ad the hazard fucto adjusted by the IWL dstrbuto cosderg data set related to the falure tme of 194 devces. Comparg the emprcal survval fucto wth the adjusted models we observed a goodess of the ft for the verse weghted Ldley dstrbuto. Ths result s also cofrmed by the dfferet dscrmato crtero methods cosdered sce IWL dstrbuto has the mmum value. Based o the TTT-plot there s a dcato that the hazard fucto has upsde-dow bathtub falure rate ths result s cofrmed by the adjusted hazard fucto. Therefore, from the proposed methodology the data related to the falure tme of 194 devces a arcraft ca be descrbed by the verse weghted Ldley dstrbuto. 18
8. Cocludg Remarks I ths paper, a ew dstrbuto called verse weghted Ldley s proposed ad ts mathematcal propertes were studed detal. The maxmum lkelhood estmators of the parameters ad ther asymptotc propertes were obtaed, we also preseted two correctve approaches to derve a modfed MLEs that are bas-free to secod order, as well as the MLEs the presece of radomly cesored data. The smulato study showed that the CMLE ad ACLME preset extremely effcet estmators for both parameters for ay sample szes. The practcal mportace of the IWL dstrbuto was reported a real applcato, whch our ew dstrbuto retured better fttg comparso wth other well-kow dstrbutos. Ackowledgemets The authors are very grateful to the revewers for ther helpful ad useful commets that mproved the mauscrpt. The authors researchers are partally supported by the Brazla Isttutos: CNPq, CAPES ad FAPESP. Refereces [1] O.O. Aale, Heterogeety survval aalyss, Statstcs medce 7 (1988, pp. 1121 1137. [2] S. Al, O the bayesa estmato of the weghted ldley dstrbuto, Joural of Statstcal Computato ad Smulato 85 (2015, pp. 855 880. [3] H.S. Bakouch, B.M. Al-Zahra, A.A. Al-Shomra, V.A. March, ad F. Louzada, A exteded ldley dstrbuto, Joural of the Korea Statstcal Socety 41 (2012, pp. 75 85. [4] R.E. Barlow ad R.A. Campo, Total tme o test processes ad applcatos to falure data aalyss., Tech. Rep., DTIC Documet, 1975. [5] G.M. Cordero, E.C. Da Rocha, J.G.C. Da Rocha, ad F. Crbar-Neto, Bas-corrected maxmum lkelhood estmato for the beta dstrbuto, Joural of Statstcal Computato ad Smulato 58 (1997, pp. 21 35. [6] G.M. Cordero ad R. Kle, Bas correcto arma models, Statstcs & Probablty Letters 19 (1994, pp. 169 176. [7] D.R. Cox ad E.J. Sell, A geeral defto of resduals, Joural of the Royal Statstcal Socety. Seres B (Methodologcal (1968, pp. 248 275. [8] B. Efro, Bootstrap methods: aother look at the jackkfe, The Aals of Statstcs 7 (1979, pp. 1 26. [9] M. Ghtay, D. Al-Mutar, N. Balakrsha, ad L. Al-Eez, Power ldley dstrbuto ad assocated ferece, Computatoal Statstcs & Data Aalyss 64 (2013, pp. 20 33. [10] M. Ghtay, F. Alqallaf, D. Al-Mutar, ad H. Husa, A two-parameter weghted ldley dstrbuto ad ts applcatos to survval data, Mathematcs ad Computers Smulato 81 (2011, pp. 1190 1201. [11] M. Ghtay, B. Ateh, ad S. Nadarajah, Ldley dstrbuto ad ts applcato, Mathematcs ad computers smulato 78 (2008, pp. 493 506. [12] D.E. Gles, H. Feg, ad R.T. Godw, O the bas of the maxmum lkelhood estmator for the two-parameter lomax dstrbuto, Commucatos Statstcs-Theory ad Methods 42 (2013, pp. 1934 1950. [13] R.E. Glaser, Bathtub ad related falure rate characterzatos, Joural of the Amerca Statstcal Assocato 75 (1980, pp. 667 672. 19
[14] M.S. Goodma, Y. L, ad R.C. Twar, Survval aalyss wth chage pot hazard fuctos (2006. [15] A.J. Lemote, Improved pot estmato for the kumaraswamy dstrbuto, Joural of Statstcal Computato ad Smulato 81 (2011, pp. 1971 1982. [16] D.V. Ldley, Fducal dstrbutos ad bayes theorem, Joural of the Royal Statstcal Socety. Seres B (Methodologcal (1958, pp. 102 107. [17] J. Mazuchel, F. Louzada, ad M. Ghtay, Comparso of estmato methods for the parameters of the weghted ldley dstrbuto, Appled Mathematcs ad Computato 220 (2013, pp. 463 471. [18] G.P. Patl ad C.R. Rao, Weghted dstrbutos ad sze-based samplg wth applcatos to wldlfe populatos ad huma famles, Bometrcs (1978, pp. 179 189. [19] P.L. Ramos, F. Louzada, ad E. Ramos, A effcet, closed-form map estmator for akagam-m fadg parameter, IEEE Commucatos Letters 20 (2016, pp. 2328 2331. [20] P. Ramos ad F. Louzada, The geeralzed weghted ldley dstrbuto: Propertes, estmato ad applcatos, Coget Mathematcs 3 (2016, p. 1256022. [21] J. Reath, J. Dog, ad M. Wag, Improved parameter estmato of the log-logstc dstrbuto wth applcatos, Computatoal Statstcs (2016, pp. 1 18. [22] J. Schwartz, R.T. Godw, ad D.E. Gles, Improved maxmum-lkelhood estmato of the shape parameter the akagam dstrbuto, Joural of Statstcal Computato ad Smulato 83 (2013, pp. 434 445. [23] V.K. Sharma, S.K. Sgh, U. Sgh, ad V. Agwal, The verse ldley dstrbuto: a stress-stregth relablty model wth applcato to head ad eck cacer data, Joural of Idustral ad Producto Egeerg 32 (2015, pp. 162 173. [24] M. Temour ad S. Nadarajah, Bas corrected mles for the webull dstrbuto based o records, Statstcal Methodology 13 (2013, pp. 12 24. [25] M. Wag ad W. Wag, Bas-corrected maxmum lkelhood estmato of the parameters of the weghted ldley dstrbuto, Commucatos Statstcs-Smulato ad Computato 46 (2017, pp. 530 545. [26] H. Zakerzadeh ad A. Dolat, geeralzed ldley dstrbuto, Joural of Mathematcal Exteso (2009. 20