A continuous probability model allowing occurrence of zero values: Polynomial-exponential
|
|
- Priscilla Price
- 5 years ago
- Views:
Transcription
1 A cotiuous probability model allowig occurrece of zero values: Polyomial-epoetial Christophe Cheseau, Hassa Bakouch, Pedro Ramos To cite this versio: Christophe Cheseau, Hassa Bakouch, Pedro Ramos. A cotiuous probability model allowig occurrece of zero values: Polyomial-epoetial. 7. <hal-57656> HAL Id: hal Submitted o 3 Aug 7 HAL is a multi-discipliary ope access archive for the deposit ad dissemiatio of scietific research documets, whether they are published or ot. The documets may come from teachig ad research istitutios i Frace or abroad, or from public or private research ceters. L archive ouverte pluridiscipliaire HAL, est destiée au dépôt et à la diffusio de documets scietifiques de iveau recherche, publiés ou o, émaat des établissemets d eseigemet et de recherche fraçais ou étragers, des laboratoires publics ou privés.
2 A cotiuous probability model allowig occurrece of zero values: Polyomial-epoetial Christophe Cheseau, Hassa S. Bakouch, Pedro L. Ramos 3 Uiversité de Cae Basse-Normadie, LMNO, Départemet de Mathématiques, UFR de Scieces, F-43, Cae. Departmet of Mathematics, Faculty of Sciece, Tata Uiversity, Tata, Egypt 3 Istitute of Mathematical Sciece ad Computig, Uiversity of São Paulo, São Carlos, Brazil ARTICLE HISTORY Compiled August, 7 ABSTRACT This paper deals with a ew two-parameter lifetime distributio with icreasig, decreasig ad costat hazard rate. This distributio allows the occurrece of zero values ad ivolves the epoetial, liear epoetial ad other combiatios of Weibull distributios as submodels. May statistical properties of the distributio are derived. Maimum likelihood estimatio of the parameters is ivestigated with a simulatio study for performace of the estimators. Two real data sets are aalyzed for illustrative purposes ad it is oted that the distributio is a highly alterative to the gamma, Weibull, Logormal ad epoetiated epoetial distributios. KEYWORDS Hazard rate; Two-parameter distributios; Reliability ad statistical measures; Maimum Likelihood Estimatio; Data applicatios.. Itroductio with motivatios I aalysis of the lifetime data, mootoe hazard rates are commo. Such data ca be modelled usig the log-ormal, Weibull ad gamma distributios. The Weibull distributio is more popular tha log-ormal ad gamma because the survival ad hazard rate fuctios of the last two distributios have ot a closed form ad hece umerical itegratios are required. Gupta ad Kudu [6] itroduced the epoetiated epoetial EE) distributio as a etesio to the epoetial distributio ad also as a alterative to the gamma distributio. Further developmets o the epoetiated epoetial distributio ca be see i [7]. I may practical applicatios, cotiuous probability models that allow occurrece of zero values have vast importace, for eample i forecast models whe we observe the mothly raifall precipitatio, it is commo i dry periods the o occurrece of precipitatio, therefore the occurrece of zero values ca be observed i differet measures, such as the average, maimum ad miimum. I survival aalysis, we may Correspodig author: Pedro Luiz Ramos, pedrolramos@usp.br
3 observe data with istataeous failure due costructio problem. Aother eample of zero occurrece is the hydrologic data i arid ad semiarid regios, like aual peak flow discharges. Moreover, zero occurrece ca be met i may areas, e.g. maufacturig defects, medical cosultatios, hydrology, ecology ad ecoometrics. The above four distributios log-ormal, Weibull, gamma ad epoetiated epoetial) do ot provide the characteristic of zero occurrece. Therefore, we itroduce a ew cotiuous probability model that allows occurrece of zero values ad it shall be called Polyomial-epoetial PE) distributio. This model represets a powerful alterative to the metioed distributios above ad also a etesio to epoetial, liear epoetial ad other combiatios of Weibull distributios as it shall be see later. The model is specified i terms of the cumulative distributio fuctio cdf) as: F ) e λ α,, + )\{}, where α >, λ > ad F ) if < ) ad Figure presets plots of the cdf for differet values of α ad λ. Note that, there is a discotiuity whe, to avoid that we ca cosidered the cotiuous etesio for give by F ) { e λ α,, + )\{}, e λα,. ) Moreover, a sigificat accout of mathematical properties for the ew distributio are provided ad the hazard rate fuctio has costat, icreasig or decreasig shape, which makig the PE distributio a alterative to the metioed distributios above. Aother attractive feature of the PE distributio is that it has closed form epressios for its cdf ad hazard rate fuctio, which is ot the case, for the log-ormal ad gamma distributios. Further, the distributio has several particular sub-models. For α, the PE distributio gives the epoetial ad whe α, it reduces to the oe parameter liear epoetial distributio. Also, whe α is a iteger, we ca epress F ) as F ) ep λ α ), k >, ad hece the survival fuctio k S) F ) of the PE distributio is the product of the survival fuctios of Weibull distributios with parameters λ, ), λ, ),..., ad λ, α) with respect to value of k, which meas the distributio has a epoetial geeral polyomial. Also, we ca ote F ) give by ) has the form: F ) GH)), where G) is the cdf associated to the epoetial distributio with parameter λ > ad H) α with H) α) is a positive icreasig fuctio with H) ad lim + H) +. The iferetial procedure for the parameters of PE distributio is preseted usig the maimum likelihood estimatio. It is show that oe of the estimators ca be obtaied i closed-form ad this allows us to obtai the estimates solvig a simple oe o-liear equatio. The performace of the MLEs are compared usig etesive umerical simulatios. The paper is orgaized as follows. I Sectio we itroduce the PE distributio. Sectio 3 is devoted to some of its mathematical properties. Estimatios of the parameters via the maimum likelihood method are ivestigated i Sectio 4. A simulatio
4 aalysis is give i Sectio 5. I Sectio 6 we apply our proposed model i two real data sets. Fially i Sectio 7 we coclude the paper. F) α.5 α. α.5 α. α.8 α.6 F) α. α.5 α. α.99 α.95 α Figure. λ. Cumulative desity fuctio shapes for PE distributio cosiderig differet values of α ad fied. Polyomial-epoetial distributio The associated probability desity fuctio pdf) of the cdf give by the equatio ) is f) λ αα+ α + ) α + ) e λ α,, + )\{}, ) with the cotiuous etesio for : f) λαα+) e λα. Figure presets some plot of the PE distributio for differet values of α ad λ, ad showig various shapes of the desity fuctio with left skewess. f) α3., λ.5 α., λ.5 α., λ.5 α.8, λ.5 α.4, λ.5 α., λ.5 f) α5, λ. α5, λ.5 α5, λ. α5, λ.5 α5, λ. α5, λ Figure. Desity fuctio shapes for PE distributio cosiderig differet values of α ad λ. As we kow, may distributios such as log-ormal, Weibull, gamma ad epoetiated epoetial, to list a few, do ot allow occurrece of zero values. I this regard, 3
5 the followig remark shows that the PE distributio ca be used as a model with occurrece of zero values. Let us observe that f) λ > for all α > ad λ >. Further, it follows from equatio 3) that f) λα α e λα as. Therefore, the upper tail behavior of the pdf is a product of a polyomial power ad a epoetial polyomial power decay, both of them deped oly o α. Obviously, larger values of α lead to faster decay of the upper tail, which iterprets α as a shape parameter. The hazard rate fuctio hrf) is give by h) λ αα+ α + ) α + ),, + )\{}, with the cotiuous etesio for give by h) λαα+). The study of the behavior of the hazard fuctio is ot a easy task. Glaser s [5] lemma is difficult to be implemeted sice ηt) d dt logft)) does ot has a simple form. However, from graphical aalysis we observed that the hazard fuctio presets a decreasig hazard rate for α > ad λ >, icreasig hazard rate for α < ad λ > ad costat rate for α, ad for this purpose some plots of the hazard fuctio with various values for the parameters α ad λ are preseted i Figure 3. h) 3 4 α3. α. α. α.8 α.4 α. h) α. α.5 α. α.99 α.95 α Figure 3. Hazard fuctio shapes for PF distributio for λ ad cosiderig differet values of α. Moreover, ote that, h) λ as, ad h) λα α as. Hece, we coclude that, the lower tail of the hazard rate fuctio is a costat, while its upper tail is a polyomial which allows for icreasig, decreasig ad costat hazard rate shapes. Remark. Usig the idicator fuctio: A ) if A ad elsewhere, we have followig aalytic epressios for F ), f) ad h): F ) e λ α ),+ )\{}) e λα) {}),+ ) ), 4
6 f) λ αα+ α + ) α + ) e λ ) λαα + ) {}) e λα,+ ) ), ),+ )\{}) α 3) ad h) λ αα+ α + ) α + ) ),+ )\{}) λαα + ) These aalytic epressios will be useful i the et. ) {}),+ ) ). 3. Mathematical properties 3.. Some useful epasios The result below presets a polyomial epasio of the cdf F ) give by ). Propositio 3.. We have the followig epasio, for, + )\{}, where F ) k l j ] A k,l,j [ ) k k+αl+j,) ) + αl j,+ ) ), A k,l,j k l ) k j ) k! )l+j λ k. Proof of Propositio 3.. First of all, let us ow ivestigate a epasio for e λ α by distiguishig the case [, ) ad the case >. If [, ), ote that α α ) give e λ α k ) λ α k k! k l j k l ) k j. The epoetial ad biomial series k k! λ)k k α ) k ) k ) k! )l+j λ) k k+αl+j. If >, ote that α α ). It follows from the epoetial ad biomial series that e λ α k ) λ α k k! k l j k l ) k j 5 k k! λk α ) k ) k! )l+j λ k αl j. ) k
7 Hece e λ α Therefore F ) e λ α k l j ) A k,l,j ) k k+αl+j,) ) + αl j,+ ) ). k l j This complete the proof of Propositio 3.. ) A k,l,j ) k k+αl+j,) ) + αl j,+ ) ). The result below presets a epasio of the pdf f) give by ) via polyomial ad the epoetial fuctio e λ, which will be importat to esure the permutatio of sum ad itegral i several probabilistic quatities. Propositio 3.. We have the followig epasio, for, + )\{}, f) where k l j ] B k,l,j [ ) k R k,l,j ),) ) + S k,l,j ),+ ) ) e λ, B k,l,j k l ) k j ) k! )l+j λ k+, 4) ad R k,l,j ) α k+α )l+j+α+ α + ) k+α )l+j+α + k+α )l+j S k,l,j ) α k+α )l j+α α + ) k+α )l j +α + k+α )l j. Proof of Propositio 3.. Let us observe that λ α α ) λ + λ λ α + λ. Therefore, we ca epress the pdf f) as f) λα α+ α + ) α + )e λ Let us ow ivestigate a epasio for ad >. ) e λ α α ) e λ. by distiguishig [, ) 6
8 If [, ), it follows from the epoetial ad biomial series that α ) e λ k! ) k k ) k λ α k! λ)k k α ) k ) k+ k l j k l ) k j ) k! )l+j λ) k k+α )l+j. Hece If >, epoetial ad biomial series give α ) e λ ) k k! ) λ α k! λk k α ) k k k k l j k l ) k j ) k+ ) k! )l+j λ k k+α )l j. α ) e λ ) ) k k l j k! )l+j λ k k l j [ ] ) k k+α )l+j,) ) + k+α )l j,+ ) ). Owig to this equality, we obtai the desired epasio: f) k l j ] B k,l,j [ ) k k+α )l+j,) ) + k+α )l j,+ ) ) α α+ α + ) α + ) e λ ) B k,l,j ) k R k,l,j ),) ) + S k,l,j ),+ ) ) e λ. k l j This eds the proof of Propositio Momets ad momet geeratig fuctio Here ad after, we cosider a radom variable X followig the PEα, λ) distributio with α > ad λ >. We defie the upper icomplete gamma fuctio as Γs, ) + t s e t dt ad the lower icomplete gamma fuctio as γs, ) ts e t dt, s >,. 7
9 Let r. Usig the otatios ad the result of Propositio 3., the r-momets of X is give by EX r ) r f)d k l j B k,l,j ) k U k,l,j,r + V k,l,j,r ), where U k,l,j,r r R k,l,j )e λ d ad V k,l,j,r + r S k,l,j )e λ d. We have U k,l,j,r α + r+k+α )l+j+α+ e λ d α + ) r+k+α )l+j e λ d γr + k + α )l + j + α +, λ) α λ r+k+α )l+j+α+ α + ) γr + k + α )l + j +, λ) + λ r+k+α )l+j+. O the other had, we have V k,l,j,r α r+k+α )l j+α e λ d α + ) r+k+α )l j e λ d Γr + k + α )l j + α, λ) α λ r+k+α )l j+α α + ) Γr + k + α )l j, λ) + λ r+k+α )l j. r+k+α )l+j+α e λ d + The momet geeratig fuctio of X is give by, for t < λ, Mt) Ee tx ) e t f)d k l j γr + k + α )l + j + α +, λ) λ r+k+α )l+j+α+ r+k+α )l j +α e λ d Γr + k + α )l j + α, λ) λ r+k+α )l j +α B k,l,j ) k U k,l,j t) + V k,l,j t) ), where Uk,l,j t) et R k,l,j )e λ d ad Vk,l,j t) + e t S k,l,j )e λ d. We have Uk,l,j t) α k+α )l+j+α+ e λ t) d α + ) + k+α )l+j e λ t) d γk + α )l + j + α +, λ t) α λ t) k+α )l+j+α+ α + ) γk + α )l + j +, λ t) + λ t) k+α )l+j+. k+α )l+j+α e λ t) d γk + α )l + j + α +, λ t) λ t) k+α )l+j+α+ 8
10 O the other had, we have V + k,l,j t) α k+α )l j+α e λ t) d α + ) + + k+α )l j e λ t) d Γk + α )l j + α, λ t) α λ t) k+α )l j+α α + ) Γk + α )l j, λ t) + λ t) k+α )l j. + k+α )l j +α e t λ) d Γk + α )l j + α, λ t) λ t) k+α )l j +α 3.3. O other meas ad momets The followig result proposes a epasio of the primitive t r f)d, with t >. It will be useful i the et. Propositio 3.3. For ay r ad t >, we have t r f)d k l j where B k,l,j is defied by 4), ad ) B k,l,j ) k Uk,l,j,r t) + V k,l,j,r t) [,+ )t), 5) Uk,l,j,r + k + α )l + j + α +, λ mit, )) t) αγr λ r+k+α )l+j+α+ γr + k + α )l + j + α +, λ mit, )) α + ) + λ r+k+α )l+j+α+ γr + k + α )l + j +, λ mit, )) λ r+k+α )l+j+ Vk,l,j,r + k + α )l j + α, λt) t) αγr + λ r+k+α )l j+α Γr + k + α )l j, λt) λ r+k+α )l j. α + ) Γr + k + α )l j + α, λt) λ r+k+α )l j +α The proof of Propositio 3.3 follows from Propositio 3. with Uk,l,j,r t) mit,) r R k,l,j )e λ d ad Vk,l,j,r t) t r S k,l,j )e λ d. The epressios of these itegrals i terms of upper icomplete gamma fuctio ad the lower icomplete gamma fuctio is obtaied proceedig as Subsectio 3.. Several crucial coditioal momets use the itegral t r f)d for various values of r. The most useful of them are preseted below. For ay t >, The r-th coditioal momets of X is give by, EX r X > t) F t) + t r f)d EX r ) F t) t ) r f)d. 9
11 The r-th reversed momets of X is give by Let µ EX). EX r X t) F t) The mea deviatios of X about µ is give by t δ E X µ ) µf µ) r f)d. µ f)d The mea deviatios of X about the media M is give by M η E X M ) µ f)d. Residual life parameters ca be also determied usig EX r ) ad t r f)d for several values of r. I particular, The mea residual life is defied as Kt) EX t X > t) St) ad the variace residual life is give by V t) V arx t X > t) EX ) St) The mea reversed residual life is defied as EX) t Lt) Et X X t) t F t) ad the variace reversed residual life is give by t ) f)d t ) f)d t tkt) [Kt)]. t f)d W t) V art X X t) tlt) [Lt)] t + F t) t f)d Stress-stregth reliability Let X be a radom variable followig the PEα, λ ) distributio with pdf deoted by f X ) ad Y be a radom variable followig the PEα, λ ) distributio with cdf deoted by F Y ), with α >, λ > ad λ >. The the stress-stregth reliability is defied by R P X > Y ). Sice the itegral o, + ) of the pdf ) with
12 parameters α, λ + λ ) deoted by f ) is equal to oe, we have R P X > Y ) + λ λ + λ + f X )F Y )d α α+ α + ) α + λ ) e λ+λ) α d + f )d λ λ + λ. This result is of iterest i a parametric estimatio cotet; oly λ ad λ eed to be estimated to have a estimatio of R the maimum likelihood estimators for λ ad λ yield the maimum likelihood estimator for R by the plug-i method) Order statistics distributios We ow itroduce order statistics ad preset some of their properties i our mathematical framework geeral results ca be foud, for istace, i [4]). Let X, X,..., X be i.i.d. radom variables followig the PEα, λ) distributio with α > ad λ >. Let us cosider its order statistics as X :, X :,..., X :, i.e., for ay i {,..., }, X i: {X,..., X } with X :... X : so X : X ) ifx,..., X ) ad X : X ) supx,..., X )). Let us ow preset some importat distributios related to X :, X :,..., X :. Some importat of them ivolvig our distributio are preseted below. The geeral epressio of the cdf of X i: is give by F Xi: )! i )! i)! i ) i ) k k i + k [F )]i+k, R. k Hece, for ay, + )/{}, we have F Xi: ) For the case, we have F Xi: )! i )! i)! i ) i ) k [ e λ k i + k k! i )! i)! i ) i ) k k i + k k The geeral epressio of the pdf of X i: is give by f Xi: ) Thus, for ay, + )/{}, we have f Xi: ) α ] i+k. [ e λα] i+k.! i )! i)! [F )]i [ F )] i f), R.! [ i )! i)! λ e λ α ] i α α+ α + ) α + ) e λ i+) α.
13 For the case, we have! [ f Xi: ) e λα] i λαα + ) e λα i+). i )! i)! As i Remark, oe ca epress f Xi: ) i a oe form usig,+ )\{} ) ad {} ). For i < j ad i < j, the geeral epressio of the joit pdf of X i:, X j: ) is give by f Xi:,X j:) i, j )! i )! j)! j i ) [F i)] i [F j ) F i )] j i [ F j )] j f i )f j ). For the case i, j ), + ) /{, )}, we have f Xi:,X j:) i, j )! i )! j)! j i ) λ αα+ i [ e λi α i i ] i [ e λi α i i α + ) α i + )αα+ j α + ) α j + ) i ) j ) e λi α i α e λj j j i λ j+)j ] j i α j j. The epressio of f Xi:,X j:) i, j ) for i, j ), + ) /{, )} ca be set i a similar maer, usig the values of F ) ad f). I the followig propositio, we provide the asymptotic distributios of the etreme values X : ad X :, ad show that they are epoetial ad Gumbel distributios, respectively, which adapt the stadards of the asymptotic distributio of etremes. Propositio 3.4. Let X ) be a sequece of i.i.d. radom variables followig the P Eα, λ) distributio, the X ) coverges i distributio to a radom variable X havig the epoetial distributio ) of parameter λ. ) X) α log) λ coverges i distributio to a radom variable X havig the Gumbel distributio of parameters ad λ. Proof of Propositio 3.4. Let us prove the two poits i tur. Sice X,..., X are i.i.d., usig stadard mathematical argumets, for, + )/{}, the cdf of X ) is give by F X) ) F )) e λ ) α. So lim + F X) ) e λ F X ). This eds the proof of the first poit. Agai, sice X,..., X are i.i.d., usig stadard mathematical argumets, for
14 ) { log) λ, + / F X α log) ) λ ) F log) λ + log) }, the cdf of X α ) log) λ λ ) )) α e is give by λ+ log) λ ) α + log) λ log) + α λ ). Therefore, whe +, several equivaleces give F X α log) ) λ ) e log e λ ) e e λ. Hece lim + F X α log) ) λ ) e e λ F X ). The secod poit is proved Stochastic orderig The orderig mechaism i life time distributios ca be illustrate by the cocept of stochastic orderig. See, for istace, []. This subsectio presets the basic of this cocept, with a result usig the proposed distributio. A radom variable X is said to be stochastically smaller tha a radom variable Y i the stochastic order X st Y ) if the associated cdfs satisfy: F X ) F Y ) for all. hazard rate order X hr Y ) if the associated hrfs satisfy: h X ) h Y ) for all. likelihood ratio order X lr Y ) if the ratio of the associated pdfs give by fx) f Y ) decreases i. Importat equivaleces eist; whe the supports of X ad Y have a commo fiite left ed-poit, the we have: X lr Y X hr Y X st Y. Propositio 3.5. Let X be a radom variable followig the PEα, λ ) distributio with pdf deoted by f X ) ad Y be a radom variable followig the PEα, λ ) distributio with pdf deoted by f Y ), with α >, λ > ad λ >. If λ > λ, the we have X lr Y. Proof of Propositio 3.5. For ay, + )/{}, we have f X ) f Y ) λ e λ λ) α λ ad fx) f Y ) λ λ e λ λ)α. As metioed i Itroductio, the fuctio H) α is icreasig fuctio of. Ideed, we have H ) +αα+ +α) α ) ad a study of fuctio shows that + α α+ + α) α for ay α >. Hece, if λ > λ, fx) f Y ) decreases i ad X lr Y. Propositio 3.5 is proved. Remark. Oe ca prove that λ > λ implies that X hr Y,which follows immediately from the defiitio ): for ay, + )/{}, we have h X ) 3
15 λ α α+ α+) α + α λ α+ α+) α + ) h ) Y ), ad for, we have h X ) λ αα+) λαα+) h Y ) Record values distributios Let us ow focus our attetio o record values ad preset some of their properties i our mathematical cotet geeral results ca be foud, for istace, i []). Let X, X,..., X be i.i.d. radom variables followig the PEα, λ) distributio with α > ad λ >. We defie a sequece of record times U) as follows: U), U) mi{j; j > U ), X j > X U ) } for. We defie the i-th upper record value by R i X Ui), with R X. The geeral epressio of the cdf of R i is give by i [ log F ))] k F Ri ) F )), R. k! k Hece, for ay, + )/{}, we have For the case, we have F Ri ) e λ α i k i F Ri ) e λα [ ] λ α k. k! k The geeral epressio of the pdf of R i is give by λα) k. k! f Ri ) [ log F ))]i f), R. i )! Hece, for ay, + )/{}, we have f Ri ) Note that, for, we have i )! λi i α ) i α α+ α + ) α + ) ) i+ e λ α. f Ri ) α + i )! λα)i e λα. The geeral epressio of the joit pdf of R,..., R ) is give by, for,..., ) R with <... <, f R,...,R ),..., ) f ) h k ). k 4
16 For the case,..., ), + ) /{,..., )}, we have + αα+ α + ) α + f R,...,R ),..., ) λ ) e λ α k α α+ k α + ) α k + k ). The epressio of f R,...,R ),..., ) for,..., ), + ) /{,..., )} ca be set i a similar maer, usig the values of f) ad h). For i < j ad i < j, the geeral epressio of the joit pdf of R i, R j ) is give by f Ri,R j) i, j ) [ log F i))] i i )! For the case i, j ), + ) /{, )}, we have f Ri,R j) i, j ) [ log F i))] i i )! [ α i )! λj i i i αα+ i [ )] j i log F i) F j) h i )f j ). j i )! [ log )] j i F i) F j) j i )! ] i j i )! h i )f j ) [ j α j j i α i ] j i i α + ) α i + )αα+ j α + ) α j + ) α i ) j ) e λj j j. The epressio of f Ri,R j) i, j ) for i, j ), + ) /{, )} ca be set i a similar maer, usig the values of F ), f) ad h). 4. Maimum likelihood estimatio Let X, X,..., X be a radom sample with commo distributio the PEα, λ) distributio with α > ad λ >. Let θ α, λ) be the parameter vector ad,,..., be the observed values. The the likelihood fuctio associated to,..., is give by Lθ) i λ αα+ i α + ) α i + i ) e λi α i i ),+ )\{} i) λαα + ) For the set of simplicity, let us set u i,+ )\{} i ), v i {} i ), ad v i. The log-likelihood fuctio ca be epressed as i lθ) logλ) λ u i + i α α+ i α + ) α i u i log + ) i ) i α i u i i i + v i logλ) + logα) + logα + ) log)). i i 5 e λα ) {} i). i u i +
17 The oliear log-likelihood equatios lθ) θ are give by lθ) α u i α i )[ + α log i )] log i ) i α α+ i i α + ) α i + λ + α + ) v i, α + i α i u i i i log i) i 6) ad lθ) λ λ i u i i α i i. 7) Note that after some algebraic maipulatios i 7) we have that λ i u i α i. 8) i i Replacig 8) i 6) the maimum likelihood estimates of α ad λ are determied by solvig the oe liear equatio 6). Sice it does ot admit ay eplicit solutio, umerical procedures ca be used. Uder mild coditios the maimum likelihood estimators are asymptotically ormal, with a asymptotic variace-covariace matri depedig o the Fisher iformatio matri. Crucial quatities to determie the etries of this matri are the secod partial derivatives of the log-likelihood fuctio give by l θ) α u i α log i ) i )[ + α log i )] log i ))α α+ i α + ) α i + ) i α α+ i i α + ) α i + ) u i α α i [ i )[ + α log i )] log i )] i α α+ i i α + ) α i + ) α ) i λ u i i i log i)) α + α + ) v i, i i ad l θ) λ λ, lθ) λ α α i u i i i log i). i 5. Simulatio aalysis I this sectio a simulatio study is preseted to compare the efficiecy of the maimum likelihood method. The compariso is performed by computig the Bias ad the 6
18 mea square errors MSE) give by Biasα i ) N Biasλ i ) N N ˆα i α, i N ˆλ i λ, i MSEα i ) N MSEλ i ) N N ˆα i α), j N ˆλ i λ), j where N is the umber of estimates obtaied through the MLE. The 95% coverage probability of the asymptotic cofidece itervals are also evaluated. Here we epect that the most efficiet estimatio method returs both Bias ad MSE closer to zero. Additioally, for a large umber of eperimets, usig a 95% cofidece level, the frequecies of itervals that covered the true values of α ad λ should be closer to 95%. The programs ca be obtaied, upo request. The values of the PE were geerated cosiderig the followig algorithm: ) Geerate U i Uiform, ), i,..., ; ) Fid i from the solutio of F i ) u i, i,..., ; The simulatio study is performed uder the assumptio.5, ) ad 4, ), N, ad, 35,..., 46). The chose values allow us to obtai data with both icreasig α < ) ad decreasig α > ) hazard rate. It is importat to poit out that, the results of this simulatio study were similar for differet choices of α ad λ. Figures 4 ad 5 preset the Bias, the MSE ad the coverage probability with a 95% cofidece level of the estimates obtaied through the MLE for differet samples of size. Bias α) MSE α) CP α) Bias λ) MSE λ)....3 CP λ) MLE Figure 4. MLE. Bias, MSEs related from the estimates of α 3 ad λ for N simulated samples uder the From the obtaied results, we ca coclude that as there is a icrease of both Bias ad MSE ted to zero, i.e., the estimator are asymptotic efficiecy. Moreover, the coverage probability of the cofidece levels ted to the omial value assumed.95. Therefore, the MLE showed to be a good estimator for the parameters of the PE distributio. 7
19 Bias α) MSE α) CP α) Bias λ)..4.8 MSE λ)..4.8 CP λ) MLE Figure 5. MLE. Bias, MSEs related from the estimates of α.9 ad λ. for N simulated samples uder the 6. Applicatio to real data I this sectio, we illustrate the fleibility of our proposed distributio by cosiderig two real data sets. The results obtaied from the PE distributio are compared with oes of the Weibull, Gamma, Logormal ad the EE distributios, ad oparametric survival fuctio Here, differet discrimiatio criterio are cosidered based o log likelihood fuctio. Let k be the umber of parameters to be fitted ad ˆθ the MLEs of θ, the discrimiatio criterio methods are respectively: Akaike iformatio criterio AIC lˆθ; ) + k; Corrected Akaike iformatio criterio AICC AIC + k k + ))/ k ); Haa-Qui iformatio criterio HQIC lˆθ; ) + k log log)); Cosistet Akaike iformatio criterio CAIC lˆθ; ) + k log) + ). The best model is the oe which provides the miimum values of these criteria. The Kolmogorov-Smirov KS) test is also cosidered aimig to check the goodess of the fit for the models. This procedure is widely kow ad based o the KS statistic D sup F ) F ; θ), where sup is the supremum of the set of distaces, F ) is the empirical distributio fuctio ad F ; θ) is the cdf of the fitted distributio. Uder a sigificace level of 5% if the data comes from F ; θ) ull hypothesis), the hypothesis is rejected if the P-value is smaller tha.5. The et subsectios give a descriptio of the used data ad their aalysis uder the metioed distributios above. 6.. Air coditioig system data The data have bee preseted by Proscha [9] ad further aalyzed by Adamidis ad Loukas []. Table cosists of the umber of successive failures of the air coditioig system of each member of a fleet of 3 Boeig 7 jet airplaes. Table displays the MLEs, stadard-error ad 95% cofidece itervals for α ad λ. Table 3 presets the results of AIC, AICc, HQIC, CAIC criteria, for the compared 8
20 Table. Data set related to the umber of successive failures of the air coditioig system of each member of a fleet of 3 Boeig 7 jet airplaes distributios. Table. MLE, Stadard-error ad 95% cofidece itervals for α ad λ. θ MLE S. error CI 95% θ) α ;.3) λ ;.8) Table 3. Results of AIC, AICc, HQIC, CAIC criteria ad the p-value for the KS test for the compared distributios cosiderig the umber of successive failures of the air coditioig system of each member of a fleet of 3 Boeig 7 jet airplaes. Test PE Weibull Gamma Logormal EE AIC AICc CAIC HQIC P-value I Figure 6, we have the TTT-plot, the survival fuctio adjusted by the compared distributios ad the o-parametric survival fuctio. Comparig the empirical survival fuctio with the adjusted models we observe a goodess of the fit for the PE distributio, which is cofirmed from differet discrimiatio criterio methods as the PE distributio has the miimum value for all statistics ad the largest for the P-value. Cosequetly, we coclude that the data related to the umber of successive failures of the air coditioig system of each member of a fleet of 3 Boeig 7 jet airplaes ca be described by the PE distributio. 9
21 St) Empirical PE Weibull Gamma Logormal EE ht) Failure Failure Figure 6. Survival fuctio adjusted by the compared distributios ad a o-parametric method cosiderig the data sets related to the umber of successive failures of the air coditioig system of each member of a fleet of 3 Boeig 7 jet airplaes. 6.. Mothly raifall data I this subsectio, we cosidered the data set firstly preseted i Bakouch et al. [3]. The data set is related to the total mothly raifall durig September at São Carlos located i southeaster Brazil. Such city has a active idustrial profile ad high agricultural importace where the study of the behavior of dry ad wet periods has proved to be strategic ad ecoomically sigificat its developmet. Table 4 presets the data related to the total mothly raifall mm) durig September at São Carlos. Table 4. The data set related to the total mothly mm) raifall durig September at São Carlos Nadarajah ad Haghighi [8] observed that maimum likelihood estimate of the shape parameter is o-uique for the Gamma, Weibull ad Geeralized epoetial distributios if data set cosists of zeros ad therefore oe of these three distributios ca fit this kid of data set. O the other had the PE distributio is defied as, which allow us to use the origial values i the presece of zero. Table 5 displays the MLE, stadard-error ad 95% cofidece itervals for α ad λ. Table 6 presets the results of the P-value for the KS test for the compared distributios. Table 5. MLE, Stadard-error ad 95% cofidece itervals for α ad λ. θ MLE S. error CI 95% θ) α ;.69) λ.44..;.9)
22 Table 6. Results of KS test for the compared distributios cosiderig the data set related to the total mothly raifall durig September at São Carlos. Test PE Weibull Gamma Logormal EE P-value I the Figure 7, the survival fuctio adjusted by the compared distributios ad the Kapla-Meier estimator. St) Empirical PE Weibull Gamma Logormal EE ht) Failure Failure Figure 7. Survival fuctio adjusted by the compared distributios ad a o-parametric method cosiderig the data set related to the total mothly raifall durig September at São Carlos. The adjusted models whe compared to the empirical survival show a goodess of the fit for the PE distributio. Additioally, this result is corroborated by the P-value of the KS test. Therefore, our proposed distributio ca be used to describe the data related to the total mothly raifall durig September at São Carlos.
23 7. Cocludig remarks I this paper, we itroduced a ew two-parameter distributio called polyomial epoetial distributio, which geeralizes the ordiary epoetial, liear epoetial ad other combiatios of Weibull distributios, because survival fuctio of the PE distributio represets the product of the survival fuctios of Weibull distributios with parameters λ, ), λ, ),..., ad λ, α). The ew distributio could be a alterative model for lifetime data, specially for the presece of istataeous failures iliers), sice stadard distributios such as Gamma, Weibull, Logormal ad epoetiated epoetial may ot be suitable. We provided a mathematical treatmet of the ew distributio. The estimatio of parameters was discussed by the maimum likelihood approach. Simulatio studies were performed to assess the performace of the maimum likelihood estimators. We fitted the proposed distributio to two real data sets ad compared its fit to those of commoly kow lifetime distributios, establishig that the ew model ca be a good competitor for the latter. We hope that the proposed distributio may be used i wide applicatios as well as lifetime modelig. Future studies ca be ivestigated by usig other baselie fuctios G), see Itroductio sectio. Refereces [] K. Adamidis ad S. Loukas, A lifetime distributio with decreasig failure rate, Statistics & Probability Letters ), pp [] M. Ahsaullah, Record statistics, Nova Sciece Publishers, 995. [3] H. Bakouch, S. Dey, P. Ramos, ad F. Louzada, Biomial-epoetial distributio: Differet estimatio methods with weather applicatios, TEMA to appear) 7). [4] H.A. David ad H. Nagaraja, Wiley series i probability ad statistics, Order Statistics, Third Editio 3), pp [5] R.E. Glaser, Bathtub ad related failure rate characterizatios, Joural of the America Statistical Associatio 75 98), pp [6] R. Gupta ad D. Kudu, Epoetiated epoetial family: A alterative to gamma ad weibull distributios, Biometrical Joural 43 ), pp [7] R. Gupta ad D. Kudu, Geeralized epoetial distributio: Eistig results ad some recet developmet, J. Stat. Pla. If. 37 7), pp [8] S. Nadarajah ad F. Haghighi, A etesio of the epoetial distributio, Statistics 45 ), pp [9] F. Proscha, Theoretical eplaatio of observed decreasig failure rate, Techometrics 5 963), pp [] M. Shaked, J.G. Shathikumar, ad Y. Tog, Stochastic orders ad their applicatios, SIAM Review ), pp
Confidence interval for the two-parameter exponentiated Gumbel distribution based on record values
Iteratioal Joural of Applied Operatioal Research Vol. 4 No. 1 pp. 61-68 Witer 2014 Joural homepage: www.ijorlu.ir Cofidece iterval for the two-parameter expoetiated Gumbel distributio based o record values
More informationOn the behavior at infinity of an integrable function
O the behavior at ifiity of a itegrable fuctio Emmauel Lesige To cite this versio: Emmauel Lesige. O the behavior at ifiity of a itegrable fuctio. The America Mathematical Mothly, 200, 7 (2), pp.75-8.
More informationPOWER AKASH DISTRIBUTION AND ITS APPLICATION
POWER AKASH DISTRIBUTION AND ITS APPLICATION Rama SHANKER PhD, Uiversity Professor, Departmet of Statistics, College of Sciece, Eritrea Istitute of Techology, Asmara, Eritrea E-mail: shakerrama009@gmail.com
More informationThe Sampling Distribution of the Maximum. Likelihood Estimators for the Parameters of. Beta-Binomial Distribution
Iteratioal Mathematical Forum, Vol. 8, 2013, o. 26, 1263-1277 HIKARI Ltd, www.m-hikari.com http://d.doi.org/10.12988/imf.2013.3475 The Samplig Distributio of the Maimum Likelihood Estimators for the Parameters
More informationResampling Methods. X (1/2), i.e., Pr (X i m) = 1/2. We order the data: X (1) X (2) X (n). Define the sample median: ( n.
Jauary 1, 2019 Resamplig Methods Motivatio We have so may estimators with the property θ θ d N 0, σ 2 We ca also write θ a N θ, σ 2 /, where a meas approximately distributed as Oce we have a cosistet estimator
More informationGoodness-Of-Fit For The Generalized Exponential Distribution. Abstract
Goodess-Of-Fit For The Geeralized Expoetial Distributio By Amal S. Hassa stitute of Statistical Studies & Research Cairo Uiversity Abstract Recetly a ew distributio called geeralized expoetial or expoetiated
More informationImprovement of Generic Attacks on the Rank Syndrome Decoding Problem
Improvemet of Geeric Attacks o the Rak Sydrome Decodig Problem Nicolas Arago, Philippe Gaborit, Adrie Hauteville, Jea-Pierre Tillich To cite this versio: Nicolas Arago, Philippe Gaborit, Adrie Hauteville,
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.436J/15.085J Fall 2008 Lecture 19 11/17/2008 LAWS OF LARGE NUMBERS II THE STRONG LAW OF LARGE NUMBERS
MASSACHUSTTS INSTITUT OF TCHNOLOGY 6.436J/5.085J Fall 2008 Lecture 9 /7/2008 LAWS OF LARG NUMBRS II Cotets. The strog law of large umbers 2. The Cheroff boud TH STRONG LAW OF LARG NUMBRS While the weak
More informationMATH 320: Probability and Statistics 9. Estimation and Testing of Parameters. Readings: Pruim, Chapter 4
MATH 30: Probability ad Statistics 9. Estimatio ad Testig of Parameters Estimatio ad Testig of Parameters We have bee dealig situatios i which we have full kowledge of the distributio of a radom variable.
More informationComparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series
Applied Mathematical Scieces, Vol. 7, 03, o. 6, 3-337 HIKARI Ltd, www.m-hikari.com http://d.doi.org/0.988/ams.03.3430 Compariso Study of Series Approimatio ad Covergece betwee Chebyshev ad Legedre Series
More informationA New Lifetime Distribution For Series System: Model, Properties and Application
Joural of Moder Applied Statistical Methods Volume 7 Issue Article 3 08 A New Lifetime Distributio For Series System: Model, Properties ad Applicatio Adil Rashid Uiversity of Kashmir, Sriagar, Idia, adilstat@gmail.com
More informationMaximum likelihood estimation from record-breaking data for the generalized Pareto distribution
METRON - Iteratioal Joural of Statistics 004, vol. LXII,. 3, pp. 377-389 NAGI S. ABD-EL-HAKIM KHALAF S. SULTAN Maximum likelihood estimatio from record-breakig data for the geeralized Pareto distributio
More informationDouble Stage Shrinkage Estimator of Two Parameters. Generalized Exponential Distribution
Iteratioal Mathematical Forum, Vol., 3, o. 3, 3-53 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/.9/imf.3.335 Double Stage Shrikage Estimator of Two Parameters Geeralized Expoetial Distributio Alaa M.
More informationRandom Variables, Sampling and Estimation
Chapter 1 Radom Variables, Samplig ad Estimatio 1.1 Itroductio This chapter will cover the most importat basic statistical theory you eed i order to uderstad the ecoometric material that will be comig
More informationCEE 522 Autumn Uncertainty Concepts for Geotechnical Engineering
CEE 5 Autum 005 Ucertaity Cocepts for Geotechical Egieerig Basic Termiology Set A set is a collectio of (mutually exclusive) objects or evets. The sample space is the (collectively exhaustive) collectio
More informationConvergence of random variables. (telegram style notes) P.J.C. Spreij
Covergece of radom variables (telegram style otes).j.c. Spreij this versio: September 6, 2005 Itroductio As we kow, radom variables are by defiitio measurable fuctios o some uderlyig measurable space
More informationANOTHER WEIGHTED WEIBULL DISTRIBUTION FROM AZZALINI S FAMILY
ANOTHER WEIGHTED WEIBULL DISTRIBUTION FROM AZZALINI S FAMILY Sulema Nasiru, MSc. Departmet of Statistics, Faculty of Mathematical Scieces, Uiversity for Developmet Studies, Navrogo, Upper East Regio, Ghaa,
More informationLecture 2: Monte Carlo Simulation
STAT/Q SCI 43: Itroductio to Resamplig ethods Sprig 27 Istructor: Ye-Chi Che Lecture 2: ote Carlo Simulatio 2 ote Carlo Itegratio Assume we wat to evaluate the followig itegratio: e x3 dx What ca we do?
More informationGUIDELINES ON REPRESENTATIVE SAMPLING
DRUGS WORKING GROUP VALIDATION OF THE GUIDELINES ON REPRESENTATIVE SAMPLING DOCUMENT TYPE : REF. CODE: ISSUE NO: ISSUE DATE: VALIDATION REPORT DWG-SGL-001 002 08 DECEMBER 2012 Ref code: DWG-SGL-001 Issue
More informationThis section is optional.
4 Momet Geeratig Fuctios* This sectio is optioal. The momet geeratig fuctio g : R R of a radom variable X is defied as g(t) = E[e tx ]. Propositio 1. We have g () (0) = E[X ] for = 1, 2,... Proof. Therefore
More informationA statistical method to determine sample size to estimate characteristic value of soil parameters
A statistical method to determie sample size to estimate characteristic value of soil parameters Y. Hojo, B. Setiawa 2 ad M. Suzuki 3 Abstract Sample size is a importat factor to be cosidered i determiig
More informationLecture 19: Convergence
Lecture 19: Covergece Asymptotic approach I statistical aalysis or iferece, a key to the success of fidig a good procedure is beig able to fid some momets ad/or distributios of various statistics. I may
More informationEstimation for Complete Data
Estimatio for Complete Data complete data: there is o loss of iformatio durig study. complete idividual complete data= grouped data A complete idividual data is the oe i which the complete iformatio of
More informationGini Index and Polynomial Pen s Parade
Gii Idex ad Polyomial Pe s Parade Jules Sadefo Kamdem To cite this versio: Jules Sadefo Kamdem. Gii Idex ad Polyomial Pe s Parade. 2011. HAL Id: hal-00582625 https://hal.archives-ouvertes.fr/hal-00582625
More informationMOMENT-METHOD ESTIMATION BASED ON CENSORED SAMPLE
Vol. 8 o. Joural of Systems Sciece ad Complexity Apr., 5 MOMET-METHOD ESTIMATIO BASED O CESORED SAMPLE I Zhogxi Departmet of Mathematics, East Chia Uiversity of Sciece ad Techology, Shaghai 37, Chia. Email:
More informationTURBULENT FUNCTIONS AND SOLVING THE NAVIER-STOKES EQUATION BY FOURIER SERIES
TURBULENT FUNCTIONS AND SOLVING THE NAVIER-STOKES EQUATION BY FOURIER SERIES M Sghiar To cite this versio: M Sghiar. TURBULENT FUNCTIONS AND SOLVING THE NAVIER-STOKES EQUATION BY FOURIER SERIES. Iteratioal
More informationProperties and Hypothesis Testing
Chapter 3 Properties ad Hypothesis Testig 3.1 Types of data The regressio techiques developed i previous chapters ca be applied to three differet kids of data. 1. Cross-sectioal data. 2. Time series data.
More informationAn Introduction to Randomized Algorithms
A Itroductio to Radomized Algorithms The focus of this lecture is to study a radomized algorithm for quick sort, aalyze it usig probabilistic recurrece relatios, ad also provide more geeral tools for aalysis
More informationTopic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2016 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
More informationOptimization Results for a Generalized Coupon Collector Problem
Optimizatio Results for a Geeralized Coupo Collector Problem Emmauelle Aceaume, Ya Busel, E Schulte-Geers, B Sericola To cite this versio: Emmauelle Aceaume, Ya Busel, E Schulte-Geers, B Sericola. Optimizatio
More information( θ. sup θ Θ f X (x θ) = L. sup Pr (Λ (X) < c) = α. x : Λ (x) = sup θ H 0. sup θ Θ f X (x θ) = ) < c. NH : θ 1 = θ 2 against AH : θ 1 θ 2
82 CHAPTER 4. MAXIMUM IKEIHOOD ESTIMATION Defiitio: et X be a radom sample with joit p.m/d.f. f X x θ. The geeralised likelihood ratio test g.l.r.t. of the NH : θ H 0 agaist the alterative AH : θ H 1,
More informationLecture 27: Optimal Estimators and Functional Delta Method
Stat210B: Theoretical Statistics Lecture Date: April 19, 2007 Lecture 27: Optimal Estimators ad Fuctioal Delta Method Lecturer: Michael I. Jorda Scribe: Guilherme V. Rocha 1 Achievig Optimal Estimators
More informationA Generalized Gamma-Weibull Distribution: Model, Properties and Applications
Marquette Uiversity e-publicatios@marquette Mathematics, Statistics ad Computer Sciece Faculty Research ad Publicatios Mathematics, Statistics ad Computer Sciece, Departmet of --06 A Geeralized Gamma-Weibull
More informationA proposed discrete distribution for the statistical modeling of
It. Statistical Ist.: Proc. 58th World Statistical Cogress, 0, Dubli (Sessio CPS047) p.5059 A proposed discrete distributio for the statistical modelig of Likert data Kidd, Marti Cetre for Statistical
More information32 estimating the cumulative distribution function
32 estimatig the cumulative distributio fuctio 4.6 types of cofidece itervals/bads Let F be a class of distributio fuctios F ad let θ be some quatity of iterest, such as the mea of F or the whole fuctio
More informationThe standard deviation of the mean
Physics 6C Fall 20 The stadard deviatio of the mea These otes provide some clarificatio o the distictio betwee the stadard deviatio ad the stadard deviatio of the mea.. The sample mea ad variace Cosider
More informationInternational Journal of Mathematical Archive-5(7), 2014, Available online through ISSN
Iteratioal Joural of Mathematical Archive-5(7), 214, 11-117 Available olie through www.ijma.ifo ISSN 2229 546 USING SQUARED-LOG ERROR LOSS FUNCTION TO ESTIMATE THE SHAPE PARAMETER AND THE RELIABILITY FUNCTION
More informationProduct measures, Tonelli s and Fubini s theorems For use in MAT3400/4400, autumn 2014 Nadia S. Larsen. Version of 13 October 2014.
Product measures, Toelli s ad Fubii s theorems For use i MAT3400/4400, autum 2014 Nadia S. Larse Versio of 13 October 2014. 1. Costructio of the product measure The purpose of these otes is to preset the
More informationTesting the number of parameters with multidimensional MLP
Testig the umber of parameters with multidimesioal MLP Joseph Rykiewicz To cite this versio: Joseph Rykiewicz. Testig the umber of parameters with multidimesioal MLP. ASMDA 2005, 2005, Brest, Frace. pp.561-568,
More informationChapter 3. Strong convergence. 3.1 Definition of almost sure convergence
Chapter 3 Strog covergece As poited out i the Chapter 2, there are multiple ways to defie the otio of covergece of a sequece of radom variables. That chapter defied covergece i probability, covergece i
More informationEECS564 Estimation, Filtering, and Detection Hwk 2 Solns. Winter p θ (z) = (2θz + 1 θ), 0 z 1
EECS564 Estimatio, Filterig, ad Detectio Hwk 2 Sols. Witer 25 4. Let Z be a sigle observatio havig desity fuctio where. p (z) = (2z + ), z (a) Assumig that is a oradom parameter, fid ad plot the maximum
More informationThe Goldbach conjectures
The Goldbach cojectures Jamel Ghaouchi To cite this versio: Jamel Ghaouchi. The Goldbach cojectures. 2015. HAL Id: hal-01243303 https://hal.archives-ouvertes.fr/hal-01243303 Submitted o
More informationGoodness-Of-Fit For The Generalized Exponential Distribution. Abstract
Goodess-Of-Fit For The Geeralized Expoetial Distributio By Amal S. Hassa stitute of Statistical Studies & Research Cairo Uiversity Abstract Recetly a ew distributio called geeralized expoetial or expoetiated
More information62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +
62. Power series Defiitio 16. (Power series) Give a sequece {c }, the series c x = c 0 + c 1 x + c 2 x 2 + c 3 x 3 + is called a power series i the variable x. The umbers c are called the coefficiets of
More informationBayesian and E- Bayesian Method of Estimation of Parameter of Rayleigh Distribution- A Bayesian Approach under Linex Loss Function
Iteratioal Joural of Statistics ad Systems ISSN 973-2675 Volume 12, Number 4 (217), pp. 791-796 Research Idia Publicatios http://www.ripublicatio.com Bayesia ad E- Bayesia Method of Estimatio of Parameter
More information6.3 Testing Series With Positive Terms
6.3. TESTING SERIES WITH POSITIVE TERMS 307 6.3 Testig Series With Positive Terms 6.3. Review of what is kow up to ow I theory, testig a series a i for covergece amouts to fidig the i= sequece of partial
More informationLet us give one more example of MLE. Example 3. The uniform distribution U[0, θ] on the interval [0, θ] has p.d.f.
Lecture 5 Let us give oe more example of MLE. Example 3. The uiform distributio U[0, ] o the iterval [0, ] has p.d.f. { 1 f(x =, 0 x, 0, otherwise The likelihood fuctio ϕ( = f(x i = 1 I(X 1,..., X [0,
More informationECONOMETRIC THEORY. MODULE XIII Lecture - 34 Asymptotic Theory and Stochastic Regressors
ECONOMETRIC THEORY MODULE XIII Lecture - 34 Asymptotic Theory ad Stochastic Regressors Dr. Shalabh Departmet of Mathematics ad Statistics Idia Istitute of Techology Kapur Asymptotic theory The asymptotic
More informationDepartment of Mathematics
Departmet of Mathematics Ma 3/103 KC Border Itroductio to Probability ad Statistics Witer 2017 Lecture 19: Estimatio II Relevat textbook passages: Larse Marx [1]: Sectios 5.2 5.7 19.1 The method of momets
More informationLecture 6 Simple alternatives and the Neyman-Pearson lemma
STATS 00: Itroductio to Statistical Iferece Autum 06 Lecture 6 Simple alteratives ad the Neyma-Pearso lemma Last lecture, we discussed a umber of ways to costruct test statistics for testig a simple ull
More informationJournal of Multivariate Analysis. Superefficient estimation of the marginals by exploiting knowledge on the copula
Joural of Multivariate Aalysis 102 (2011) 1315 1319 Cotets lists available at ScieceDirect Joural of Multivariate Aalysis joural homepage: www.elsevier.com/locate/jmva Superefficiet estimatio of the margials
More informationBootstrap Intervals of the Parameters of Lognormal Distribution Using Power Rule Model and Accelerated Life Tests
Joural of Moder Applied Statistical Methods Volume 5 Issue Article --5 Bootstrap Itervals of the Parameters of Logormal Distributio Usig Power Rule Model ad Accelerated Life Tests Mohammed Al-Ha Ebrahem
More informationThe (P-A-L) Generalized Exponential Distribution: Properties and Estimation
Iteratioal Mathematical Forum, Vol. 12, 2017, o. 1, 27-37 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/imf.2017.610140 The (P-A-L) Geeralized Expoetial Distributio: Properties ad Estimatio M.R.
More informationG. R. Pasha Department of Statistics Bahauddin Zakariya University Multan, Pakistan
Deviatio of the Variaces of Classical Estimators ad Negative Iteger Momet Estimator from Miimum Variace Boud with Referece to Maxwell Distributio G. R. Pasha Departmet of Statistics Bahauddi Zakariya Uiversity
More informationFACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING. Lectures
FACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING Lectures MODULE 5 STATISTICS II. Mea ad stadard error of sample data. Biomial distributio. Normal distributio 4. Samplig 5. Cofidece itervals
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 3 9/11/2013. Large deviations Theory. Cramér s Theorem
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/5.070J Fall 203 Lecture 3 9//203 Large deviatios Theory. Cramér s Theorem Cotet.. Cramér s Theorem. 2. Rate fuctio ad properties. 3. Chage of measure techique.
More informationA Simple Proof of the Shallow Packing Lemma
A Simple Proof of the Shallow Packig Lemma Nabil Mustafa To cite this versio: Nabil Mustafa. A Simple Proof of the Shallow Packig Lemma. Discrete ad Computatioal Geometry, Spriger Verlag, 06, 55 (3), pp.739-743.
More informationComparison of Minimum Initial Capital with Investment and Non-investment Discrete Time Surplus Processes
The 22 d Aual Meetig i Mathematics (AMM 207) Departmet of Mathematics, Faculty of Sciece Chiag Mai Uiversity, Chiag Mai, Thailad Compariso of Miimum Iitial Capital with Ivestmet ad -ivestmet Discrete Time
More informationMAT1026 Calculus II Basic Convergence Tests for Series
MAT026 Calculus II Basic Covergece Tests for Series Egi MERMUT 202.03.08 Dokuz Eylül Uiversity Faculty of Sciece Departmet of Mathematics İzmir/TURKEY Cotets Mootoe Covergece Theorem 2 2 Series of Real
More informationMathematical Modeling of Optimum 3 Step Stress Accelerated Life Testing for Generalized Pareto Distribution
America Joural of Theoretical ad Applied Statistics 05; 4(: 6-69 Published olie May 8, 05 (http://www.sciecepublishiggroup.com/j/ajtas doi: 0.648/j.ajtas.05040. ISSN: 6-8999 (Prit; ISSN: 6-9006 (Olie Mathematical
More informationPrecise Rates in Complete Moment Convergence for Negatively Associated Sequences
Commuicatios of the Korea Statistical Society 29, Vol. 16, No. 5, 841 849 Precise Rates i Complete Momet Covergece for Negatively Associated Sequeces Dae-Hee Ryu 1,a a Departmet of Computer Sciece, ChugWoo
More informationA new distribution-free quantile estimator
Biometrika (1982), 69, 3, pp. 635-40 Prited i Great Britai 635 A ew distributio-free quatile estimator BY FRANK E. HARRELL Cliical Biostatistics, Duke Uiversity Medical Ceter, Durham, North Carolia, U.S.A.
More informationExpectation and Variance of a random variable
Chapter 11 Expectatio ad Variace of a radom variable The aim of this lecture is to defie ad itroduce mathematical Expectatio ad variace of a fuctio of discrete & cotiuous radom variables ad the distributio
More informationTopic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2018 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
More informationCoefficient of variation and Power Pen s parade computation
Coefficiet of variatio ad Power Pe s parade computatio Jules Sadefo Kamdem To cite this versio: Jules Sadefo Kamdem. Coefficiet of variatio ad Power Pe s parade computatio. 20. HAL Id: hal-0058658
More informationESTIMATION AND PREDICTION BASED ON K-RECORD VALUES FROM NORMAL DISTRIBUTION
STATISTICA, ao LXXIII,. 4, 013 ESTIMATION AND PREDICTION BASED ON K-RECORD VALUES FROM NORMAL DISTRIBUTION Maoj Chacko Departmet of Statistics, Uiversity of Kerala, Trivadrum- 695581, Kerala, Idia M. Shy
More informationof the matrix is =-85, so it is not positive definite. Thus, the first
BOSTON COLLEGE Departmet of Ecoomics EC771: Ecoometrics Sprig 4 Prof. Baum, Ms. Uysal Solutio Key for Problem Set 1 1. Are the followig quadratic forms positive for all values of x? (a) y = x 1 8x 1 x
More informationProbability and statistics: basic terms
Probability ad statistics: basic terms M. Veeraraghava August 203 A radom variable is a rule that assigs a umerical value to each possible outcome of a experimet. Outcomes of a experimet form the sample
More informationChapter 6 Infinite Series
Chapter 6 Ifiite Series I the previous chapter we cosidered itegrals which were improper i the sese that the iterval of itegratio was ubouded. I this chapter we are goig to discuss a topic which is somewhat
More information1.010 Uncertainty in Engineering Fall 2008
MIT OpeCourseWare http://ocw.mit.edu.00 Ucertaity i Egieerig Fall 2008 For iformatio about citig these materials or our Terms of Use, visit: http://ocw.mit.edu.terms. .00 - Brief Notes # 9 Poit ad Iterval
More informationProbability and Statistics
robability ad Statistics rof. Zheg Zheg Radom Variable A fiite sigle valued fuctio.) that maps the set of all eperimetal outcomes ito the set of real umbers R is a r.v., if the set ) is a evet F ) for
More informationMath 152. Rumbos Fall Solutions to Review Problems for Exam #2. Number of Heads Frequency
Math 152. Rumbos Fall 2009 1 Solutios to Review Problems for Exam #2 1. I the book Experimetatio ad Measuremet, by W. J. Youde ad published by the by the Natioal Sciece Teachers Associatio i 1962, the
More information7.1 Convergence of sequences of random variables
Chapter 7 Limit Theorems Throughout this sectio we will assume a probability space (, F, P), i which is defied a ifiite sequece of radom variables (X ) ad a radom variable X. The fact that for every ifiite
More informationDistribution of Random Samples & Limit theorems
STAT/MATH 395 A - PROBABILITY II UW Witer Quarter 2017 Néhémy Lim Distributio of Radom Samples & Limit theorems 1 Distributio of i.i.d. Samples Motivatig example. Assume that the goal of a study is to
More informationProblem Set 4 Due Oct, 12
EE226: Radom Processes i Systems Lecturer: Jea C. Walrad Problem Set 4 Due Oct, 12 Fall 06 GSI: Assae Gueye This problem set essetially reviews detectio theory ad hypothesis testig ad some basic otios
More informationLast time: Moments of the Poisson distribution from its generating function. Example: Using telescope to measure intensity of an object
6.3 Stochastic Estimatio ad Cotrol, Fall 004 Lecture 7 Last time: Momets of the Poisso distributio from its geeratig fuctio. Gs () e dg µ e ds dg µ ( s) µ ( s) µ ( s) µ e ds dg X µ ds X s dg dg + ds ds
More informationIntroductory statistics
CM9S: Machie Learig for Bioiformatics Lecture - 03/3/06 Itroductory statistics Lecturer: Sriram Sakararama Scribe: Sriram Sakararama We will provide a overview of statistical iferece focussig o the key
More informationSince X n /n P p, we know that X n (n. Xn (n X n ) Using the asymptotic result above to obtain an approximation for fixed n, we obtain
Assigmet 9 Exercise 5.5 Let X biomial, p, where p 0, 1 is ukow. Obtai cofidece itervals for p i two differet ways: a Sice X / p d N0, p1 p], the variace of the limitig distributio depeds oly o p. Use the
More informationTopic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2018 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
More informationECO 312 Fall 2013 Chris Sims LIKELIHOOD, POSTERIORS, DIAGNOSING NON-NORMALITY
ECO 312 Fall 2013 Chris Sims LIKELIHOOD, POSTERIORS, DIAGNOSING NON-NORMALITY (1) A distributio that allows asymmetry differet probabilities for egative ad positive outliers is the asymmetric double expoetial,
More informationRecord Values from T-X Family of. Pareto-Exponential Distribution with. Properties and Simulations
Applied Mathematical Scieces, Vol. 3, 209, o., 33-44 HIKARI Ltd, www.m-hikari.com https://doi.org/0.2988/ams.209.879 Record Values from T-X Family of Pareto-Epoetial Distributio with Properties ad Simulatios
More informationCS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 5
CS434a/54a: Patter Recogitio Prof. Olga Veksler Lecture 5 Today Itroductio to parameter estimatio Two methods for parameter estimatio Maimum Likelihood Estimatio Bayesia Estimatio Itroducto Bayesia Decisio
More information17. Joint distributions of extreme order statistics Lehmann 5.1; Ferguson 15
17. Joit distributios of extreme order statistics Lehma 5.1; Ferguso 15 I Example 10., we derived the asymptotic distributio of the maximum from a radom sample from a uiform distributio. We did this usig
More informationLecture 3 The Lebesgue Integral
Lecture 3: The Lebesgue Itegral 1 of 14 Course: Theory of Probability I Term: Fall 2013 Istructor: Gorda Zitkovic Lecture 3 The Lebesgue Itegral The costructio of the itegral Uless expressly specified
More information7.1 Convergence of sequences of random variables
Chapter 7 Limit theorems Throughout this sectio we will assume a probability space (Ω, F, P), i which is defied a ifiite sequece of radom variables (X ) ad a radom variable X. The fact that for every ifiite
More informationLecture Note 8 Point Estimators and Point Estimation Methods. MIT Spring 2006 Herman Bennett
Lecture Note 8 Poit Estimators ad Poit Estimatio Methods MIT 14.30 Sprig 2006 Herma Beett Give a parameter with ukow value, the goal of poit estimatio is to use a sample to compute a umber that represets
More informationKolmogorov-Smirnov type Tests for Local Gaussianity in High-Frequency Data
Proceedigs 59th ISI World Statistics Cogress, 5-30 August 013, Hog Kog (Sessio STS046) p.09 Kolmogorov-Smirov type Tests for Local Gaussiaity i High-Frequecy Data George Tauche, Duke Uiversity Viktor Todorov,
More informationThis exam contains 19 pages (including this cover page) and 10 questions. A Formulae sheet is provided with the exam.
Probability ad Statistics FS 07 Secod Sessio Exam 09.0.08 Time Limit: 80 Miutes Name: Studet ID: This exam cotais 9 pages (icludig this cover page) ad 0 questios. A Formulae sheet is provided with the
More informationModeling and Performance Analysis with Discrete-Event Simulation
Simulatio Modelig ad Performace Aalysis with Discrete-Evet Simulatio Chapter 5 Statistical Models i Simulatio Cotets Basic Probability Theory Cocepts Useful Statistical Models Discrete Distributios Cotiuous
More information1 of 7 7/16/2009 6:06 AM Virtual Laboratories > 6. Radom Samples > 1 2 3 4 5 6 7 6. Order Statistics Defiitios Suppose agai that we have a basic radom experimet, ad that X is a real-valued radom variable
More informationEE 4TM4: Digital Communications II Probability Theory
1 EE 4TM4: Digital Commuicatios II Probability Theory I. RANDOM VARIABLES A radom variable is a real-valued fuctio defied o the sample space. Example: Suppose that our experimet cosists of tossig two fair
More informationA RANK STATISTIC FOR NON-PARAMETRIC K-SAMPLE AND CHANGE POINT PROBLEMS
J. Japa Statist. Soc. Vol. 41 No. 1 2011 67 73 A RANK STATISTIC FOR NON-PARAMETRIC K-SAMPLE AND CHANGE POINT PROBLEMS Yoichi Nishiyama* We cosider k-sample ad chage poit problems for idepedet data i a
More informationDirection: This test is worth 150 points. You are required to complete this test within 55 minutes.
Term Test 3 (Part A) November 1, 004 Name Math 6 Studet Number Directio: This test is worth 10 poits. You are required to complete this test withi miutes. I order to receive full credit, aswer each problem
More informationStatistical Inference (Chapter 10) Statistical inference = learn about a population based on the information provided by a sample.
Statistical Iferece (Chapter 10) Statistical iferece = lear about a populatio based o the iformatio provided by a sample. Populatio: The set of all values of a radom variable X of iterest. Characterized
More informationRank tests and regression rank scores tests in measurement error models
Rak tests ad regressio rak scores tests i measuremet error models J. Jurečková ad A.K.Md.E. Saleh Charles Uiversity i Prague ad Carleto Uiversity i Ottawa Abstract The rak ad regressio rak score tests
More informationGoodness-of-Fit Tests and Categorical Data Analysis (Devore Chapter Fourteen)
Goodess-of-Fit Tests ad Categorical Data Aalysis (Devore Chapter Fourtee) MATH-252-01: Probability ad Statistics II Sprig 2019 Cotets 1 Chi-Squared Tests with Kow Probabilities 1 1.1 Chi-Squared Testig................
More informationEcon 325/327 Notes on Sample Mean, Sample Proportion, Central Limit Theorem, Chi-square Distribution, Student s t distribution 1.
Eco 325/327 Notes o Sample Mea, Sample Proportio, Cetral Limit Theorem, Chi-square Distributio, Studet s t distributio 1 Sample Mea By Hiro Kasahara We cosider a radom sample from a populatio. Defiitio
More informationThe Choquet Integral with Respect to Fuzzy-Valued Set Functions
The Choquet Itegral with Respect to Fuzzy-Valued Set Fuctios Weiwei Zhag Abstract The Choquet itegral with respect to real-valued oadditive set fuctios, such as siged efficiecy measures, has bee used i
More informationIntroduction to Extreme Value Theory Laurens de Haan, ISM Japan, Erasmus University Rotterdam, NL University of Lisbon, PT
Itroductio to Extreme Value Theory Laures de Haa, ISM Japa, 202 Itroductio to Extreme Value Theory Laures de Haa Erasmus Uiversity Rotterdam, NL Uiversity of Lisbo, PT Itroductio to Extreme Value Theory
More informationμ are complex parameters. Other
A New Numerical Itegrator for the Solutio of Iitial Value Problems i Ordiary Differetial Equatios. J. Suday * ad M.R. Odekule Departmet of Mathematical Scieces, Adamawa State Uiversity, Mubi, Nigeria.
More information