Inference for the Rayleigh Distribution Based on Progressive Type-II Fuzzy Censored Data

Size: px
Start display at page:

Download "Inference for the Rayleigh Distribution Based on Progressive Type-II Fuzzy Censored Data"

Transcription

1 Journal of Modern Appled Statstcal Methods Volue 13 Issue 1 Artcle Inference for the Raylegh Dstrbuton Based on Progressve Type-II Fuzzy Censored Data Abbas Pak Shahd Charan Unversty, Ahvaz, Iran, a-pak@scu.ac.r Ghola Al Parha Shahd Charan Unversty, Ahvaz, Iran, Parha_g@scu.ac.r Mansour Saraj Shahd Charan Unversty, Ahvaz, Iran, seraj.a@scu.ac.r Follow ths and addtonal works at: Part of the Appled Statstcs Coons, Socal and Behavoral Scences Coons, and the Statstcal Theory Coons Recoended Ctaton Pak, Abbas; Parha, Ghola Al; and Saraj, Mansour (014) "Inference for the Raylegh Dstrbuton Based on Progressve Type-II Fuzzy Censored Data," Journal of Modern Appled Statstcal Methods: Vol. 13 : Iss. 1, Artcle 19. DOI: 10.37/jas/ Avalable at: Ths Regular Artcle s brought to you for free and open access by the Open Access Journals at DgtalCoons@WayneState. It has been accepted for ncluson n Journal of Modern Appled Statstcal Methods by an authorzed edtor of DgtalCoons@WayneState.

2 Journal of Modern Appled Statstcal Methods May 014, Vol. 13, No. 1, Copyrght 014 JMASM, Inc. ISSN Inference for the Raylegh Dstrbuton Based on Progressve Type-II Fuzzy Censored Data Abbas Pak Shahd Charan Unversty Ahvaz, Iran Ghola Al Parha Shahd Charan Unversty Ahvaz, Iran Mansour Saraj Shahd Charan Unversty Ahvaz, Iran Classcal statstcal analyss of the Raylegh dstrbuton deals wth precse nforaton. However, n real world stuatons, eperental perforance results cannot always be recorded or easured precsely, but each observable event ay only be dentfed wth a fuzzy subset of the saple space. Therefore, the conventonal procedures used for estatng the Raylegh dstrbuton paraeter wll need to be adapted to the new stuaton. Ths artcle dscusses dfferent estaton ethods for the paraeters of the Raylegh dstrbuton on the bass of a progressvely type-ii censorng schee when the avalable observatons are descrbed by eans of fuzzy nforaton. They nclude the au lkelhood estaton, hghest posteror densty estaton and ethod of oents. The estaton procedures are dscussed n detal and copared va Monte Carlo sulatons n ters of ther average bases and ean squared errors. Fnally, one real data set s analyzed for llustratve purposes. Keywords: Progressve type-ii censorng, fuzzy nforaton, au lkelhood prncple, hghest posteror densty estaton Introducton The Raylegh dstrbuton was orgnally ntroduced by Lord Raylegh (1880) n the feld of acoustcs; snce ts ntroducton, any researchers have used the dstrbuton n dfferent felds of scence and technology. The Raylegh dstrbuton s frequently used to odel wave heghts n oceanography, n councaton engneerng and t also has a wde applcaton n lfete data analyss, especally n relablty theory and survval analyss. An portant characterstc of the Raylegh dstrbuton s that ts hazard rate s a lnearly ncreasng functon of te at constant rate, whch akes t a sutable odel for Abbas Pak s n the Departent of Statstcs. Eal at: a-pak@scu.ac.r. Ghola Al Parha s n the Departent of Matheatcs. Eal at: parha_g@scu.ac.r. Mansour Saraj s n the Departent of Matheatcs. Eal at: seraj.a@scu.ac.r. 87

3 INFERENCE FOR THE RAYLEIGH DISTRIBUTION the lfete of coponents/tes that age rapdly wth te. Thus, as te ncreases, the relablty functon of the Raylegh dstrbuton decreases at a uch hgher rate than the eponental relablty functon does. The probablty densty functon (pdf) and the cuulatve dstrbuton functon (cdf) of a Raylegh rando varable X can be wrtten as: f ( ) e ; 0, 0, (1) and F( ) 1 e ; 0, () respectvely. Inferences for the Raylegh dstrbuton have been dscussed by several authors. Dyer and Whsenand (1973) deonstrated the portance of ths dstrbuton n councaton engneerng. Bhattacharya and Tyag (1990) entoned that n soe clncal studes dealng wth cancer patents, the survval pattern follows the Raylegh dstrbuton. Chung (1995) obtaned the best nvarant estator and the Bayes estator of the paraeter of Raylegh dstrbuton under entropy loss. Fernandez (010) addressed the probles of estatng the paraeter, hazard rate and relablty functon of the Raylegh dstrbuton on the bass of saple quantles. Dey and Mat (01) derved Bayes estator of the Raylegh paraeter and ts assocated rsk based on etended Jeffrey s pror. In any lfe testng and relablty eperents, a saple of n tes s tested, and the eperent s ternated when all of the fal. Ths procedure ay take a long te when the lfete dstrbuton of tes has a thck tal. Moreover, f the tes are epensve, such as edcal equpent, t s costly to gather nforaton fro the whole saple. There are any stuatons where eperental unts are lost or reoved fro the test before coplete falure. For eaple, ndvduals n a clncal tral ay drop out of the study, the study ay have to be ternated early for lack of funds or the test unts ay accdentally break. In other scenaros, the eperent ay have to be ternated n order to free up testng facltes for other purposes. In vew of above, censorng s used n lfe testng to save te and cost of testng unts. The reoval of unts n a test ay be unntentonal or pre-planned. Data obtaned fro such eperents are called censored saple. There are any types of censorng schees used n lfete analyss. The two ost coon 88

4 PAK ET AL censorng schees are tered type-i and type-ii censorng schees. In the conventonal type-i censorng schee, the eperent contnues up to a prespecfed te T; the conventonal type-ii censorng schee requres the eperent to contnue untl a pre-specfed nuber of falures occur. These schees, however, do not allow reoval of unts before the ternaton of the eperent; thus, a ore general knd of censorng schee called progressve type-ii censorng s consdered, whch s as follows: Suppose that n unts are placed on a lfe test and the eperenter decdes beforehand a quantty, the nuber of unts to be faled. Now at the te of the frst falure, R1 of the reanng n 1 survvng unts are randoly reoved fro the eperent. Contnung on, at the te of the second falure, R of the reanng n R1 unts are randoly reoved fro the eperent. Fnally, at the te of the th falure, all the reanng n R1 R 1(=R) survvng unts are reoved fro the eperent. The work on progressve censorng has becoe popular n lfe-testng and relablty studes. K and Han (009) studed the proble of estatng the scale paraeter of the Raylegh dstrbuton under general progressve censorng. Krshna and Kuar (011) dscussed relablty estaton for the Lndley dstrbuton wth progressve type-ii censored data. Lee et al. (011) obtaned a Bayes estator under the Raylegh dstrbuton wth a progressve type-ii rght censored saple. Raqab and Mad (011) addressed nference for the generalzed Raylegh dstrbuton based on progressvely censored data. Az et al. (01) consdered the Bayesan estaton of the paraeter and relablty functon of Raylegh dstrbuton based on a progressvely type-ii censored saple. Rastog and Trpath (01) studed paraeter estaton of the Burr type XII dstrbuton on the bass of a progressvely type-ii censored saple. Pradhan and Kundu (009) consdered the statstcal nference of the unknown paraeters of the generalzed eponental dstrbuton n presence of progressve censorng. A recent account on progressve censorng schees can be obtaned n the onograph by Balakrshnan and Aggarwala (000) or n the ecellent revew artcle by Balakrshnan (007). The above referenced studes for estatng paraeters of dfferent lfete dstrbutons under progressve type-ii censorng are lted to precse data. However, n real world stuatons, eperents do not provde eact nforaton. For eaple, the reacton te of a person to a certan stulus n a psychologcal eperence cannot be eactly deterned, but the psychologst s able to deterne t by eans of the followng precse nforaton, such as: The te of reacton s approately 5 to 35 seconds. To deal wth the lack of precson of the data, t s necessary to ncorporate the fuzzy concept to statstcal technques. Recently, 89

5 INFERENCE FOR THE RAYLEIGH DISTRIBUTION Pak et al. (013) proposed a new ethod to deterne the au lkelhood estate of the scale paraeter of a Raylegh dstrbuton under doubly type-ii censored saple fro fuzzy data. Further, n a lfe testng eperent, soe test unts ay need to be reoved at dfferent stages n the study for varous reasons. Ths would lead to progressve censorng. The purpose of ths artcle s to develop the nferental procedures for the Raylegh dstrbuton under a progressve type-ii censorng schee when the avalable observatons are reported by eans of fuzzy nforaton. The au lkelhood estate (MLE) of the paraeter s obtaned by usng EM algorth and the hghest posteror densty (HPD) estate of the unknown paraeter s coputed. The estaton va ethod of oents s dscussed, a Monte Carlo sulaton study s presented, and a coparson of all estaton procedures developed and one real data set s analyzed for llustratve purposes. Frst, the fundaental notaton and basc defntons of fuzzy set theory used heren s revewed. Consder an eperent characterzed by a probablty space S, F, P where, F s a easurable space and P belongs to a specfed faly of probablty easures { P, } on, F. Assue that the observer cannot dstngush or transt wth eactness the outcoe n the perforance of S, but that rather the avalable observaton ay be descrbed n ters of fuzzy nforaton, whch s defned as: Defnton 1 A fuzzy event on, characterzed by a Borel easurable ebershp functon ( ) fro to [0,1], where ( ) represents the grade of ebershp of to, s called fuzzy nforaton assocated wth the eperent. The set consstng of all observable events fro the eperent fuzzy nforaton syste assocated wth t, whch s defned as: deternes a Defnton A fuzzy nforaton syste (f..s.) assocated wth the eperent s a fuzzy partton { 1,..., K} of,.e., a set of K fuzzy events on satsfyng the orthogonalty condton (see Tanaka et al., 1979): 90

6 PAK ET AL K ( ) 1, k k1 where k denotes the ebershp functon of k. Accordng to Zadeh (1968), gven the eperent, FP,,, and a f..s. assocated wth t, each probablty easure P on, F nduces a probablty easure on Defnton 3 defned as: The probablty dstrbuton on nduced by P s the appng P fro to [0,1] such that ( ) ( ) dp( ),. (3) In partcular, the condtonal densty of a contnuous rando varable Y wth p.d.f. g( y ) gven the fuzzy event can be defned as ( y) g( y) g( y ) ( u) g( u) du. (4) In order to odel precse lfetes, a generalzaton of real nubers s necessary. These lfetes can be represented by fuzzy nubers. A fuzzy nuber s a subset, denoted by, of the set of real nubers (denoted by ) and s characterzed by the so called ebershp functon (.). Fuzzy nubers satsfy the constrants (Dubos and Prade, 1980): 1. : 0,1 s Borel-easurable;. : 1; and

7 INFERENCE FOR THE RAYLEIGH DISTRIBUTION 3. The so-called cuts B :, 0 1, defned as are all closed nterval,.e., B ( ) [ a, b ], (0,1]. Aong the varous types of fuzzy nubers, the trangular and trapezodal fuzzy nubers are ost convenent and useful n descrbng fuzzy lfete data. The trangular fuzzy nuber can be defned as ( a, b, c) and ts ebershp functon s defned by the followng epresson: a a b, b a c ( ) b c, c b 0 otherwse. The trapezodal fuzzy nuber can be defned as ( a, b, c, d) wth the ebershp functon: a a b, b a 1 b c ( ) d c d, d c 0 otherwse. Mau lkelhood estaton Suppose that n dentcal unts are put on a lfe testng eperent and that the lfete dstrbuton of each unt s gven by (1). Pror to the eperent, a nuber < n s deterned and the censorng schee ( R1,..., R ) wth R 0 and R n s specfed. Let 1 1 (,..., ) denote the correspondng progressvely type-ii censored saple. The lkelhood functon for the paraeter becoes proportonal to 9

8 PAK ET AL e (1 R) 1 L( ; ) (5) Now consder the proble where s not observed precsely and only partal nforaton about s avalable n the for of fuzzy nubers ( a, c, b ), 1,...,, wth the correspondng ebershp functons ( ),..., ( ). Let c(1) c()... c( ) denote the ordered values of the 1 1 eans of these fuzzy nubers. The lfete of R survvng unts, whch are reoved fro the test after the th falure, can be encoded as fuzzy nubers z,..., 1 z R wth the ebershp functons 0 z c() z ( z), j 1,..., R. j 1 z c() The fuzzy data w ( 1,...,, z1,..., z) where z s a 1 R vector wth z ( z,..., z ) for 1,...,, s thus the set of observed lfetes. The 1 R correspondng observed data log-lkelhood functon can be obtaned by usng the epresson (3) as follows: R z O w z j 1 1 j1 L ( ; ) log e ( ) d log ze ( z) dz () 1 1 log log e ( ) d R c. (6) The au lkelhood estate of the paraeter can be obtaned by azng the log-lkelhood L ( w ; ). Equatng the partal dervatve of the O log-lkelhood (6) wth respect to to zero, the resultng equaton s: 3 e ( ) LO ( w; ) d Rc () 0. 1 e ( ) d 1 Because there s no closed for of the soluton to equaton (7), an teratve nuercal search can be used to obtan the MLE. The Epectaton Mazaton (EM) algorth s a broadly applcable approach to the teratve coputaton of au lkelhood estates and useful n a varety of ncoplete-data (7) 93

9 INFERENCE FOR THE RAYLEIGH DISTRIBUTION probles. Because the observed fuzzy data w can be seen as an ncoplete specfcaton of a coplete data, the EM algorth s applcable to obtan the au lkelhood estate of the paraeter. In the followng, the fuzzy EM algorth (see Denoeu, 011) s used to deterne the MLE of. Frst, denote the lfetes of the faled and censored unts by X ( X1,..., X ) and Z ( Z1,..., Z ), respectvely, where Z s a 1 R vector wth Z ( Z 1,..., ZR ), for 1,...,. The cobnaton of W ( X,Z ) fors the coplete lfetes and the correspondng log-lkelhood functon s denoted by L (, ) c W, then, gnorng the addtve constant, R L ( W c, ) nlog zj. (8) 1 1 j1 For the E-step, t s necessary to copute the pseudo log-lkelhood functon. It can be obtaned fro (8) as follows: n E X E Z z R log ( ) ( j j ) 1 1 j1 (9) By usng (4), the condtonal epectatons coputed as: E ( X ) and E Z z can be ( j j ) 3 e ( ) d E ( X ), 1,...,, e ( ) d 1 E ( Zj zj ) c(), 1,...,, j 1,..., R. Net, the M-step nvolves the azaton of the pseudo functon (9). Therefore, f at the h th ( h) ( h 1) stage, the estate of s, then can be obtaned by azng R * L ( W c, ) nlog E ( h) ( X ) E ( h) ( Z zj ) (10) 1 1 j1 94

10 PAK ET AL wth respect to. Fro L * c( W, ) 0, ˆ ( h1) 1 [ E ( X ) R ( c 1/ )] ( h) n ( h) () (11) The teraton process contnues untl convergence,.e., untl ( h 1) ( h) L ( w; ) - L ( w ; ) for soe pre-fed 0. O O HPD estaton Consder the hghest posteror densty (HPD) estaton of the Raylegh paraeter based on observed fuzzy saple w. As a conjugate pror for, take the Gaa ( a, b) densty wth pdf gven by b ( a) a a1 b ( ) e, 0, (1) where a 0 and b 0. Based on ths pror, the posteror densty functon of gven the data can be wrtten as follows: ( ) ( brc( ) ) a1 1 ( ) 1 w e e d (13) The ethod of HPD estaton then estates as the ode of the posteror densty w ( ); therefore, the HPD estate of can be obtaned by solvng the equaton log ( w) 0 (14) where 95

11 INFERENCE FOR THE RAYLEIGH DISTRIBUTION 3 e Rc() 1 1 e log ( w) a1 ( ) d b. (15) ( ) d However, the soluton cannot be obtaned eplctly. Theore 1 dscusses the estence and unqueness of the HPD estate of. Theore 1 Let g( ) denote the functon on the rght-hand sde of the epresson n (15). Then the root of the equaton g( ) 0 ests and s unque. Proof. Fro (15), t s seen that a1 g( ), (0, ), and consequently l g( ). Also, note that 0 a1 l g( ) l 0, (0, ). 0 0 Thus, the equaton g( ) 0 has at least one root n (0, ). To prove that the root s unque, consder the frst dervatve of g, g ( ) gven by a1 g( ) log e ( ) d 1 (16) Because the ntegrand of the second ter n (16) s a log-concave functon of, and g( ) 0. It follows that g s a strctly decreasng functon w.r.t. and hence the equaton g( ) 0 has eactly one soluton. The HPD estate of ust be derved nuercally. In the followng, the Newton-Raphson algorth to deterne the HPD estate s descrbed. The Newton-Raphson algorth s a drect approach for estatng the relevant paraeters n a lkelhood functon. In ths algorth, the soluton of the ( ) lkelhood equaton s obtaned through an teratve procedure. Let ˆ h be the paraeter value fro the h th step. Then, at the ( h 1) th step of teraton process, the updated paraeter s obtaned as 96

12 PAK ET AL ˆ ˆ ( h1) ( h) log ( w) log ( w) h h (17) where the notaton evaluated at A h, for any partal dervatve A, eans the partal dervatve ˆ( h ). The second-order dervatve of log w ( ) wth respect to the paraeter, requred for proceedng wth the Newton-Raphson ethod, s obtaned as: [ ]. (18) 5 3 n log ( w) a1 e ( ) d e ( ) d 1 e ( ) d e ( ) d The teraton process then contnues untl convergence,.e., untl for soe pre-fed 0. ˆ ˆ, ( h1) ( h) Method of oents Let X be a rando varable whch has the Raylegh dstrbuton wth pdf gven by (1). It s known that the k th oent of the Raylegh odel wth paraeter s k (19) k ( k ) (1 ). EX Equatng the frst saple oent to the correspondng populaton oent, the followng equaton can be used to fnd the estate of oent ethod: n j1 R E ( X ) E( Z zj ). (0) Because the closed for of the soluton to (0) could not be obtaned, an teratve nuercal process to obtan the paraeter estate s descrbed as: 97

13 INFERENCE FOR THE RAYLEIGH DISTRIBUTION (0) 1. Let the ntal estate of be, wth h 0.. In the ( h 1) th teraton, frst copute ( h) e ( ) d E1 E ( h) ( X ) ( h), 1,...,, e ( ) d and ( h) z z e z () z dz j ( h) j ( h) z ze z () z dz j E E ( Z z ), 1,...,, j 1,..., R. The new estate of, for eaple ( h 1), can be obtaned fro: ( h1) n 4 ( E1 R E ) Checkng convergence, f the convergence occurs then the current ( h 1) s the estate of by the ethod of oents; otherwse, set hh 1 and go to Step. Fgure 1. Fuzzy nforaton syste used to encode the sulated data 98

14 PAK ET AL Nuercal Study Eperental results llustrate how the dfferent ethods behave for varyng saple szes. Frst, for fed 1 and dfferent choces of n, and censorng schee ( R1,..., R ), progressvely censored saples fro the Raylegh dstrbuton were generated, usng the ethod proposed by Balakrshnan and Sandhu (1995), as follows: 1. Generate Z fro U(0,1) for 1,...,.. For gven values of the progressve censorng schee ( R1,..., R ), set V 1/ a Z, a R, 1,...,. j1 j 3. Set U 1 V 1V... V, 1,...,. 4. Thus, X F 1 ( U ), 1,...,, s the desred progressve type-ii censored saple fro the Raylegh dstrbuton. Each realzaton of was then fuzzfed usng the f..s. shown n Fgure 1, and the ML, HPD and oent estates (MME) of for the fuzzy saple were coputed usng the ethods provded n the precedng sectons. For coputng the HPD estate of the unknown paraeter, assue that has Gaa ( a, b ) pror. To ake the coparson eanngful, t s assued that the prors are nonnforatve, and they are ab 0. Note that n ths case the prors are nonproper also. Press (001) suggested usng very sall non-negatve values of the hyperparaeters n ths case, and t wll ake the prors proper. Ths study uses ab The results are not sgnfcantly dfferent than the correspondng results obtaned usng non-proper prors, and are not reported due to space. The average values and ean squared errors of the estates, coputed over 1,000 replcaton, are presented n Tables

15 INFERENCE FOR THE RAYLEIGH DISTRIBUTION Table 1. Average value (AV) and ean squared error (MSE) of the estates of λ for dfferent censorng schees. (n = 0) Censorng schee MLE HPD MME AV MSE AV MSE AV MSE (0,,0,1) (1,0, 0) (0,1,0,,0) (0,,0,10) (10,0, 0) (0,10,0,,0) (0,,0,5) (5,0, 0) (0,5,0,,0) Table. Average value (AV) and ean squared error (MSE) of the estates of λ for dfferent censorng schees. (n = 30) Censorng schee MLE HPD MME AV MSE AV MSE AV MSE (0,,0,1) (1,0, 0) (0,1,0,,0) (0,,0,10) (10,0, 0) (0,10,0,,0) (0,,0,5) (5,0, 0) (0,5,0,,0) Table 3. Average value (AV) and ean squared error (MSE) of the estates of λ for dfferent censorng schees. (n = 50) Censorng schee MLE HPD MME AV MSE AV MSE AV MSE (0,,0,1) (1,0, 0) 1, (0,1,0,,0) (0,,0,10) (10,0, 0) (0,10,0,,0) (0,,0,5) (5,0, 0) (0,5,0,,0)

16 PAK ET AL Several ponts are clear fro the eperent: Even for sall saple szes, the perforances of the estates are satsfactory n ters of AVs and MSEs. For all the ethods, t s observed that for fed n as ncreases, the MSEs of the estates decrease. Aong the three estaton procedures developed n the paper, the HPD procedure gves the ost precse paraeter estates as shown by MSEs n Tables 1-3. Applcaton eaple To deonstrate the applcaton of proposed ethods to real data, consder the data collected durng the eperent reported by Pak et al. (013). In ths eperent, a saple of 5 ball bearngs s placed on a lfe test. A ball bearng ay work perfectly over a certan perod but be breakng for soe te and fnally be unusable at a certan te. So, the observed falure tes of the ball bearngs are reported by fuzzy nubers ( a,, b ), where a 0.05 and b 0.03 wth the ebershp functons ( a) a, a ( ) b b, b 1,...,5. Progressvely censored saples of sze 10 were consdered fro these fuzzy data usng three dfferent saplng schees, naely: Schee 1: R1... R 1 0 and R 15. Schee : R1 15 and R... R 0. Schee 3: R1... R 1 1 and R 6. the estate of the paraeter was then coputed usng the ML, HPD and oent ethods. For coputng the HPD estates, t was assued that has Gaa ( a, b) pror, ncludng the non-nforatve gaa pror,.e. ab 0 and nforatve gaa pror,.e. ab. All the results are suarzed n Table

17 INFERENCE FOR THE RAYLEIGH DISTRIBUTION Table 4. ML, HPD and oent estates of the paraeter for Eaple Schee MLE MME HPD (a = b = 0) HPD (a = b = ) Conclusons Soe work has been done n the past on the estaton of the paraeter of Raylegh dstrbuton based on coplete and censored saples, but tradtonally t s assued that the data avalable are perfored n eact nubers. In real world stuatons, however, soe collected lfete data ght be precse and are represented n the for of fuzzy nubers. Therefore, sutable statstcal ethodology s needed to handle these data. Ths artcle proposed dfferent procedures for estatng the paraeter of Raylegh dstrbuton under progressve type-ii censorng when the avalable observatons are descrbed by eans of fuzzy nforaton. They are au lkelhood estaton (MLE), hghest posteror densty (HPD) estaton and ethod of oents (MME). A sulaton study was conducted to assess the perforance of these procedures. Based on the results of the sulaton study, t ay be observed that, the perforance of the HPD estates s generally best followed by the MLE and MME. Thus, t would see reasonable to recoend the use of the HPD procedure for estatng the unknown paraeter fro the Raylegh dstrbuton. References Az, R., Yaghae, F., & Az, D. (01). Coparson of Bayesan estaton ethods for Raylegh progressve censored data under the dfferent asyetrc loss functon. Internatonal Journal of Appled Matheatcal Research, 1(4): Balakrshnan, N. (007). Progressve censorng ethodology: an apprasal. Test, 16(): Balakrshnan, N. & Aggarwala, R. (000). Progressve censorng, theory, ethods and applcatons. Boston, MA: Brkhauser. 30

18 PAK ET AL Balakrshnan, N. & Sandhu, R. A. (1995). A sple algorth for generatng progressvely Type-II censored saples. The Aercan Statstcan, 49(): Bhattacharya, S. K., & Tyag, R. K. (1990). Bayesan survval analyss based on the Raylegh odel. Trabajos de Estadstca, 5: Chung, Y. (1995). Estaton of scale paraeter fro Raylegh dstrbuton under entropy loss. Journal of Appled Matheatcs and Coputng, (1): Denoeu, T. (011). Mau lkelhood estaton fro fuzzy data usng the EM algorth. Fuzzy Sets and Systes, 183(1): Dey, S., & Mat, S. S. (01). Bayesan estaton of the paraeter of Raylegh dstrbuton under the etended Jeffrey's pror. Electronc Journal of Appled Statstcal Analyss, 5(1), Dubos, D., & Prade, H. (1980). Fuzzy Sets and Systes: Theory and Applcatons. Acadec Press, New York. Dyer, D. D. & Whsenand, C. W. (1973). Best lnear unbased estator of the paraeter of the Raylegh dstrbuton. IEEE Transactons on Relablty R- : 7-34, Fernandez, A. J. (010). Bayesan estaton and predcton based on Raylegh saple quantles. Qualty & Quantty, 44: K, C., and Han, K. (009). Estaton of the scale paraeter of the Raylegh dstrbuton under general progressve censorng. Journal of the Korean Statstcal Socety, 38: Lee, W., Wu, J., Hong, M., Ln, L., & Chan, R. (011). Assessng the lfete perforance nde of Raylegh products based on the Bayesan estaton under progressve type II rght censored saples. Journal of Coputatonal and Appled Matheatcs, 35: Krshna, H. & Kuar, K. (011). Relablty estaton n Lndley dstrbuton wth progressvely type-ii rght censored saple. Matheatcs and Coputers n Sulaton, 8: Pak, A., Parha, G. A., & Saraj, M. (013). On estaton of Raylegh scale paraeter under doubly type-ii censorng fro precse data. Journal of Data Scence, 11: Pradhan, B., & Kundu, D. (009). On progressvely censored generalzed eponental dstrbuton. Test, 18:

19 INFERENCE FOR THE RAYLEIGH DISTRIBUTION Raqab, M. Z., & Mad, T. M. (011). Inference for the generalzed Rayleg dstrbuton based on progressvely censored data. Journal of Statstcal Plannng and Inference, 141: Rastog, M. K., & Trpath, Y. M. (01). Estatng the paraeters of Burr dstrbuton under progressve type II censorng. Statstcak Methodology, 9: Raylegh, J. W. S. (1880). On the resultant of a large nuber of vbraton of the soepth and of arbtrary phase. Phlosophcal Magazne, 5th seres, 10: Press, S. J. (001). The subjectvty of scentsts and the Bayesan approach. New York: Wley. Tanaka, H. Okuda, T. & Asa, K. (1979). Fuzzy nforaton and decson n statstcal odel. In M. M. Gupta et al., Eds. Advances n Fuzzy Sets Theory and Applcatons, pp Asterda: North-Holland Publshng Co. Zadeh, L. A. (1968). Probablty easures of fuzzy events. Journal of Matheatcal Analyss and Applcatons 10:

Reliability estimation in Pareto-I distribution based on progressively type II censored sample with binomial removals

Reliability estimation in Pareto-I distribution based on progressively type II censored sample with binomial removals Journal of Scentfc esearch Developent (): 08-3 05 Avalable onlne at wwwjsradorg ISSN 5-7569 05 JSAD elablty estaton n Pareto-I dstrbuton based on progressvely type II censored saple wth bnoal reovals Ilhan

More information

System in Weibull Distribution

System in Weibull Distribution Internatonal Matheatcal Foru 4 9 no. 9 94-95 Relablty Equvalence Factors of a Seres-Parallel Syste n Webull Dstrbuton M. A. El-Dacese Matheatcs Departent Faculty of Scence Tanta Unversty Tanta Egypt eldacese@yahoo.co

More information

PARAMETER ESTIMATION IN WEIBULL DISTRIBUTION ON PROGRESSIVELY TYPE- II CENSORED SAMPLE WITH BETA-BINOMIAL REMOVALS

PARAMETER ESTIMATION IN WEIBULL DISTRIBUTION ON PROGRESSIVELY TYPE- II CENSORED SAMPLE WITH BETA-BINOMIAL REMOVALS Econoy & Busness ISSN 1314-7242, Volue 10, 2016 PARAMETER ESTIMATION IN WEIBULL DISTRIBUTION ON PROGRESSIVELY TYPE- II CENSORED SAMPLE WITH BETA-BINOMIAL REMOVALS Ilhan Usta, Hanef Gezer Departent of Statstcs,

More information

BAYESIAN AND NON BAYESIAN ESTIMATION OF ERLANG DISTRIBUTION UNDER PROGRESSIVE CENSORING

BAYESIAN AND NON BAYESIAN ESTIMATION OF ERLANG DISTRIBUTION UNDER PROGRESSIVE CENSORING www.arpapress.co/volues/volissue3/ijrras 3_8.pdf BAYESIAN AND NON BAYESIAN ESTIMATION OF ERLANG DISTRIBUTION UNDER PROGRESSIVE CENSORING R.A. Bakoban Departent of Statstcs, Scences Faculty for Grls, Kng

More information

Bayesian estimation using MCMC approach based on progressive first-failure censoring from generalized Pareto distribution

Bayesian estimation using MCMC approach based on progressive first-failure censoring from generalized Pareto distribution Aercan Journal of Theoretcal and Appled Statstcs 03; (5): 8-4 Publshed onlne August 30 03 (http://www.scencepublshnggroup.co/j/ajtas) do: 0.648/j.ajtas.03005.3 Bayesan estaton usng MCMC approach based

More information

Estimation of Reliability in Multicomponent Stress-Strength Based on Generalized Rayleigh Distribution

Estimation of Reliability in Multicomponent Stress-Strength Based on Generalized Rayleigh Distribution Journal of Modern Appled Statstcal Methods Volue 13 Issue 1 Artcle 4 5-1-014 Estaton of Relablty n Multcoponent Stress-Strength Based on Generalzed Raylegh Dstrbuton Gadde Srnvasa Rao Unversty of Dodoa,

More information

Applied Mathematics Letters

Applied Mathematics Letters Appled Matheatcs Letters 2 (2) 46 5 Contents lsts avalable at ScenceDrect Appled Matheatcs Letters journal hoepage: wwwelseverco/locate/al Calculaton of coeffcents of a cardnal B-splne Gradr V Mlovanovć

More information

Statistical inference for generalized Pareto distribution based on progressive Type-II censored data with random removals

Statistical inference for generalized Pareto distribution based on progressive Type-II censored data with random removals Internatonal Journal of Scentfc World, 2 1) 2014) 1-9 c Scence Publshng Corporaton www.scencepubco.com/ndex.php/ijsw do: 10.14419/jsw.v21.1780 Research Paper Statstcal nference for generalzed Pareto dstrbuton

More information

Estimation in Step-stress Partially Accelerated Life Test for Exponentiated Pareto Distribution under Progressive Censoring with Random Removal

Estimation in Step-stress Partially Accelerated Life Test for Exponentiated Pareto Distribution under Progressive Censoring with Random Removal Journal of Advances n Matheatcs and Coputer Scence 5(): -6, 07; Artcle no.jamcs.3469 Prevously known as Brtsh Journal of Matheatcs & Coputer Scence ISSN: 3-085 Estaton n Step-stress Partally Accelerated

More information

BAYESIAN CURVE FITTING USING PIECEWISE POLYNOMIALS. Dariusz Biskup

BAYESIAN CURVE FITTING USING PIECEWISE POLYNOMIALS. Dariusz Biskup BAYESIAN CURVE FITTING USING PIECEWISE POLYNOMIALS Darusz Bskup 1. Introducton The paper presents a nonparaetrc procedure for estaton of an unknown functon f n the regresson odel y = f x + ε = N. (1) (

More information

XII.3 The EM (Expectation-Maximization) Algorithm

XII.3 The EM (Expectation-Maximization) Algorithm XII.3 The EM (Expectaton-Maxzaton) Algorth Toshnor Munaata 3/7/06 The EM algorth s a technque to deal wth varous types of ncoplete data or hdden varables. It can be appled to a wde range of learnng probles

More information

Excess Error, Approximation Error, and Estimation Error

Excess Error, Approximation Error, and Estimation Error E0 370 Statstcal Learnng Theory Lecture 10 Sep 15, 011 Excess Error, Approxaton Error, and Estaton Error Lecturer: Shvan Agarwal Scrbe: Shvan Agarwal 1 Introducton So far, we have consdered the fnte saple

More information

Several generation methods of multinomial distributed random number Tian Lei 1, a,linxihe 1,b,Zhigang Zhang 1,c

Several generation methods of multinomial distributed random number Tian Lei 1, a,linxihe 1,b,Zhigang Zhang 1,c Internatonal Conference on Appled Scence and Engneerng Innovaton (ASEI 205) Several generaton ethods of ultnoal dstrbuted rando nuber Tan Le, a,lnhe,b,zhgang Zhang,c School of Matheatcs and Physcs, USTB,

More information

Statistical analysis of Accelerated life testing under Weibull distribution based on fuzzy theory

Statistical analysis of Accelerated life testing under Weibull distribution based on fuzzy theory Statstcal analyss of Accelerated lfe testng under Webull dstrbuton based on fuzzy theory Han Xu, Scence & Technology on Relablty & Envronental Engneerng Laboratory, School of Relablty and Syste Engneerng,

More information

PROBABILITY AND STATISTICS Vol. III - Analysis of Variance and Analysis of Covariance - V. Nollau ANALYSIS OF VARIANCE AND ANALYSIS OF COVARIANCE

PROBABILITY AND STATISTICS Vol. III - Analysis of Variance and Analysis of Covariance - V. Nollau ANALYSIS OF VARIANCE AND ANALYSIS OF COVARIANCE ANALYSIS OF VARIANCE AND ANALYSIS OF COVARIANCE V. Nollau Insttute of Matheatcal Stochastcs, Techncal Unversty of Dresden, Gerany Keywords: Analyss of varance, least squares ethod, odels wth fxed effects,

More information

Bayesian Estimation and Prediction of Generalized Pareto Distribution Based on Type II Censored Samples

Bayesian Estimation and Prediction of Generalized Pareto Distribution Based on Type II Censored Samples ISSN 1684-8403 Journal of Statstcs Volue, 015. pp. 139-165 Bayesan Estaton and Predcton of Generalzed Pareto Dstrbuton Based on Type II Censored Saples Abstract Navd Feroze 1, Muhaad Asla and Azhar Salee

More information

CHAPT II : Prob-stats, estimation

CHAPT II : Prob-stats, estimation CHAPT II : Prob-stats, estaton Randoness, probablty Probablty densty functons and cuulatve densty functons. Jont, argnal and condtonal dstrbutons. The Bayes forula. Saplng and statstcs Descrptve and nferental

More information

1 Review From Last Time

1 Review From Last Time COS 5: Foundatons of Machne Learnng Rob Schapre Lecture #8 Scrbe: Monrul I Sharf Aprl 0, 2003 Revew Fro Last Te Last te, we were talkng about how to odel dstrbutons, and we had ths setup: Gven - exaples

More information

Determination of the Confidence Level of PSD Estimation with Given D.O.F. Based on WELCH Algorithm

Determination of the Confidence Level of PSD Estimation with Given D.O.F. Based on WELCH Algorithm Internatonal Conference on Inforaton Technology and Manageent Innovaton (ICITMI 05) Deternaton of the Confdence Level of PSD Estaton wth Gven D.O.F. Based on WELCH Algorth Xue-wang Zhu, *, S-jan Zhang

More information

Designing Fuzzy Time Series Model Using Generalized Wang s Method and Its application to Forecasting Interest Rate of Bank Indonesia Certificate

Designing Fuzzy Time Series Model Using Generalized Wang s Method and Its application to Forecasting Interest Rate of Bank Indonesia Certificate The Frst Internatonal Senar on Scence and Technology, Islac Unversty of Indonesa, 4-5 January 009. Desgnng Fuzzy Te Seres odel Usng Generalzed Wang s ethod and Its applcaton to Forecastng Interest Rate

More information

COS 511: Theoretical Machine Learning

COS 511: Theoretical Machine Learning COS 5: Theoretcal Machne Learnng Lecturer: Rob Schapre Lecture #0 Scrbe: José Sões Ferrera March 06, 203 In the last lecture the concept of Radeacher coplexty was ntroduced, wth the goal of showng that

More information

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4) I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes

More information

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family IOSR Journal of Mathematcs IOSR-JM) ISSN: 2278-5728. Volume 3, Issue 3 Sep-Oct. 202), PP 44-48 www.osrjournals.org Usng T.O.M to Estmate Parameter of dstrbutons that have not Sngle Exponental Famly Jubran

More information

Least Squares Fitting of Data

Least Squares Fitting of Data Least Squares Fttng of Data Davd Eberly Geoetrc Tools, LLC http://www.geoetrctools.co/ Copyrght c 1998-2014. All Rghts Reserved. Created: July 15, 1999 Last Modfed: February 9, 2008 Contents 1 Lnear Fttng

More information

PGM Learning Tasks and Metrics

PGM Learning Tasks and Metrics Probablstc Graphcal odels Learnng Overvew PG Learnng Tasks and etrcs Learnng doan epert True dstrbuton P* aybe correspondng to a PG * dataset of nstances D{d],...d]} sapled fro P* elctaton Network Learnng

More information

RELIABILITY ASSESSMENT

RELIABILITY ASSESSMENT CHAPTER Rsk Analyss n Engneerng and Economcs RELIABILITY ASSESSMENT A. J. Clark School of Engneerng Department of Cvl and Envronmental Engneerng 4a CHAPMAN HALL/CRC Rsk Analyss for Engneerng Department

More information

ON THE NUMBER OF PRIMITIVE PYTHAGOREAN QUINTUPLES

ON THE NUMBER OF PRIMITIVE PYTHAGOREAN QUINTUPLES Journal of Algebra, Nuber Theory: Advances and Applcatons Volue 3, Nuber, 05, Pages 3-8 ON THE NUMBER OF PRIMITIVE PYTHAGOREAN QUINTUPLES Feldstrasse 45 CH-8004, Zürch Swtzerland e-al: whurlann@bluewn.ch

More information

On Pfaff s solution of the Pfaff problem

On Pfaff s solution of the Pfaff problem Zur Pfaff scen Lösung des Pfaff scen Probles Mat. Ann. 7 (880) 53-530. On Pfaff s soluton of te Pfaff proble By A. MAYER n Lepzg Translated by D. H. Delpenc Te way tat Pfaff adopted for te ntegraton of

More information

Denote the function derivatives f(x) in given points. x a b. Using relationships (1.2), polynomials (1.1) are written in the form

Denote the function derivatives f(x) in given points. x a b. Using relationships (1.2), polynomials (1.1) are written in the form SET OF METHODS FO SOUTION THE AUHY POBEM FO STIFF SYSTEMS OF ODINAY DIFFEENTIA EUATIONS AF atypov and YuV Nulchev Insttute of Theoretcal and Appled Mechancs SB AS 639 Novosbrs ussa Introducton A constructon

More information

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method Appled Mathematcal Scences, Vol. 7, 0, no. 47, 07-0 HIARI Ltd, www.m-hkar.com Comparson of the Populaton Varance Estmators of -Parameter Exponental Dstrbuton Based on Multple Crtera Decson Makng Method

More information

International Journal of Mathematical Archive-9(3), 2018, Available online through ISSN

International Journal of Mathematical Archive-9(3), 2018, Available online through   ISSN Internatonal Journal of Matheatcal Archve-9(3), 208, 20-24 Avalable onlne through www.ja.nfo ISSN 2229 5046 CONSTRUCTION OF BALANCED INCOMPLETE BLOCK DESIGNS T. SHEKAR GOUD, JAGAN MOHAN RAO M AND N.CH.

More information

Chapter One Mixture of Ideal Gases

Chapter One Mixture of Ideal Gases herodynacs II AA Chapter One Mxture of Ideal Gases. Coposton of a Gas Mxture: Mass and Mole Fractons o deterne the propertes of a xture, we need to now the coposton of the xture as well as the propertes

More information

Bayesian predictive Configural Frequency Analysis

Bayesian predictive Configural Frequency Analysis Psychologcal Test and Assessment Modelng, Volume 54, 2012 (3), 285-292 Bayesan predctve Confgural Frequency Analyss Eduardo Gutérrez-Peña 1 Abstract Confgural Frequency Analyss s a method for cell-wse

More information

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

More information

LECTURE :FACTOR ANALYSIS

LECTURE :FACTOR ANALYSIS LCUR :FACOR ANALYSIS Rta Osadchy Based on Lecture Notes by A. Ng Motvaton Dstrbuton coes fro MoG Have suffcent aount of data: >>n denson Use M to ft Mture of Gaussans nu. of tranng ponts If

More information

On an Extension of Stochastic Approximation EM Algorithm for Incomplete Data Problems. Vahid Tadayon 1

On an Extension of Stochastic Approximation EM Algorithm for Incomplete Data Problems. Vahid Tadayon 1 On an Extenson of Stochastc Approxmaton EM Algorthm for Incomplete Data Problems Vahd Tadayon Abstract: The Stochastc Approxmaton EM (SAEM algorthm, a varant stochastc approxmaton of EM, s a versatle tool

More information

Integral Transforms and Dual Integral Equations to Solve Heat Equation with Mixed Conditions

Integral Transforms and Dual Integral Equations to Solve Heat Equation with Mixed Conditions Int J Open Probles Copt Math, Vol 7, No 4, Deceber 214 ISSN 1998-6262; Copyrght ICSS Publcaton, 214 www-csrsorg Integral Transfors and Dual Integral Equatons to Solve Heat Equaton wth Mxed Condtons Naser

More information

Collaborative Filtering Recommendation Algorithm

Collaborative Filtering Recommendation Algorithm Vol.141 (GST 2016), pp.199-203 http://dx.do.org/10.14257/astl.2016.141.43 Collaboratve Flterng Recoendaton Algorth Dong Lang Qongta Teachers College, Haou 570100, Chna, 18689851015@163.co Abstract. Ths

More information

What is LP? LP is an optimization technique that allocates limited resources among competing activities in the best possible manner.

What is LP? LP is an optimization technique that allocates limited resources among competing activities in the best possible manner. (C) 998 Gerald B Sheblé, all rghts reserved Lnear Prograng Introducton Contents I. What s LP? II. LP Theor III. The Splex Method IV. Refneents to the Splex Method What s LP? LP s an optzaton technque that

More information

Chapter 12 Lyes KADEM [Thermodynamics II] 2007

Chapter 12 Lyes KADEM [Thermodynamics II] 2007 Chapter 2 Lyes KDEM [Therodynacs II] 2007 Gas Mxtures In ths chapter we wll develop ethods for deternng therodynac propertes of a xture n order to apply the frst law to systes nvolvng xtures. Ths wll be

More information

Centroid Uncertainty Bounds for Interval Type-2 Fuzzy Sets: Forward and Inverse Problems

Centroid Uncertainty Bounds for Interval Type-2 Fuzzy Sets: Forward and Inverse Problems Centrod Uncertanty Bounds for Interval Type-2 Fuzzy Sets: Forward and Inverse Probles Jerry M. Mendel and Hongwe Wu Sgnal and Iage Processng Insttute Departent of Electrcal Engneerng Unversty of Southern

More information

Change Point Estimation for Pareto Type-II Model

Change Point Estimation for Pareto Type-II Model Journal of Modern Appled Statstcal Methods Volume 3 Issue Artcle 5--04 Change Pont Estmaton for Pareto Type-II Model Gyan Prakash S. N. Medcal College, Agra, U. P., Inda, ggyanj@yahoo.com Follow ths and

More information

Multipoint Analysis for Sibling Pairs. Biostatistics 666 Lecture 18

Multipoint Analysis for Sibling Pairs. Biostatistics 666 Lecture 18 Multpont Analyss for Sblng ars Bostatstcs 666 Lecture 8 revously Lnkage analyss wth pars of ndvduals Non-paraetrc BS Methods Maxu Lkelhood BD Based Method ossble Trangle Constrant AS Methods Covered So

More information

Slobodan Lakić. Communicated by R. Van Keer

Slobodan Lakić. Communicated by R. Van Keer Serdca Math. J. 21 (1995), 335-344 AN ITERATIVE METHOD FOR THE MATRIX PRINCIPAL n-th ROOT Slobodan Lakć Councated by R. Van Keer In ths paper we gve an teratve ethod to copute the prncpal n-th root and

More information

Gadjah Mada University, Indonesia. Yogyakarta State University, Indonesia Karangmalang Yogyakarta 55281

Gadjah Mada University, Indonesia. Yogyakarta State University, Indonesia Karangmalang Yogyakarta 55281 Reducng Fuzzy Relatons of Fuzzy Te Seres odel Usng QR Factorzaton ethod and Its Applcaton to Forecastng Interest Rate of Bank Indonesa Certfcate Agus aan Abad Subanar Wdodo 3 Sasubar Saleh 4 Ph.D Student

More information

Fermi-Dirac statistics

Fermi-Dirac statistics UCC/Physcs/MK/EM/October 8, 205 Fer-Drac statstcs Fer-Drac dstrbuton Matter partcles that are eleentary ostly have a type of angular oentu called spn. hese partcles are known to have a agnetc oent whch

More information

1 Definition of Rademacher Complexity

1 Definition of Rademacher Complexity COS 511: Theoretcal Machne Learnng Lecturer: Rob Schapre Lecture #9 Scrbe: Josh Chen March 5, 2013 We ve spent the past few classes provng bounds on the generalzaton error of PAClearnng algorths for the

More information

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

2E Pattern Recognition Solutions to Introduction to Pattern Recognition, Chapter 2: Bayesian pattern classification

2E Pattern Recognition Solutions to Introduction to Pattern Recognition, Chapter 2: Bayesian pattern classification E395 - Pattern Recognton Solutons to Introducton to Pattern Recognton, Chapter : Bayesan pattern classfcaton Preface Ths document s a soluton manual for selected exercses from Introducton to Pattern Recognton

More information

AN ANALYSIS OF A FRACTAL KINETICS CURVE OF SAVAGEAU

AN ANALYSIS OF A FRACTAL KINETICS CURVE OF SAVAGEAU AN ANALYI OF A FRACTAL KINETIC CURE OF AAGEAU by John Maloney and Jack Hedel Departent of Matheatcs Unversty of Nebraska at Oaha Oaha, Nebraska 688 Eal addresses: aloney@unoaha.edu, jhedel@unoaha.edu Runnng

More information

Maximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models

Maximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models ECO 452 -- OE 4: Probt and Logt Models ECO 452 -- OE 4 Mamum Lkelhood Estmaton of Bnary Dependent Varables Models: Probt and Logt hs note demonstrates how to formulate bnary dependent varables models for

More information

Quantum Particle Motion in Physical Space

Quantum Particle Motion in Physical Space Adv. Studes Theor. Phys., Vol. 8, 014, no. 1, 7-34 HIKARI Ltd, www.-hkar.co http://dx.do.org/10.1988/astp.014.311136 Quantu Partcle Moton n Physcal Space A. Yu. Saarn Dept. of Physcs, Saara State Techncal

More information

Small-Sample Equating With Prior Information

Small-Sample Equating With Prior Information Research Report Sall-Saple Equatng Wth Pror Inforaton Sauel A Lvngston Charles Lews June 009 ETS RR-09-5 Lstenng Learnng Leadng Sall-Saple Equatng Wth Pror Inforaton Sauel A Lvngston and Charles Lews ETS,

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

Solutions for Homework #9

Solutions for Homework #9 Solutons for Hoewor #9 PROBEM. (P. 3 on page 379 n the note) Consder a sprng ounted rgd bar of total ass and length, to whch an addtonal ass s luped at the rghtost end. he syste has no dapng. Fnd the natural

More information

Maximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models

Maximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models ECO 452 -- OE 4: Probt and Logt Models ECO 452 -- OE 4 Maxmum Lkelhood Estmaton of Bnary Dependent Varables Models: Probt and Logt hs note demonstrates how to formulate bnary dependent varables models

More information

Interval Estimation of Stress-Strength Reliability for a General Exponential Form Distribution with Different Unknown Parameters

Interval Estimation of Stress-Strength Reliability for a General Exponential Form Distribution with Different Unknown Parameters Internatonal Journal of Statstcs and Probablty; Vol. 6, No. 6; November 17 ISSN 197-73 E-ISSN 197-74 Publshed by Canadan Center of Scence and Educaton Interval Estmaton of Stress-Strength Relablty for

More information

On the number of regions in an m-dimensional space cut by n hyperplanes

On the number of regions in an m-dimensional space cut by n hyperplanes 6 On the nuber of regons n an -densonal space cut by n hyperplanes Chungwu Ho and Seth Zeran Abstract In ths note we provde a unfor approach for the nuber of bounded regons cut by n hyperplanes n general

More information

Least Squares Fitting of Data

Least Squares Fitting of Data Least Squares Fttng of Data Davd Eberly Geoetrc Tools, LLC http://www.geoetrctools.co/ Copyrght c 1998-2015. All Rghts Reserved. Created: July 15, 1999 Last Modfed: January 5, 2015 Contents 1 Lnear Fttng

More information

International Journal of Engineering Research and Modern Education (IJERME) Impact Factor: 7.018, ISSN (Online): (

International Journal of Engineering Research and Modern Education (IJERME) Impact Factor: 7.018, ISSN (Online): ( CONSTRUCTION AND SELECTION OF CHAIN SAMPLING PLAN WITH ZERO INFLATED POISSON DISTRIBUTION A. Palansamy* & M. Latha** * Research Scholar, Department of Statstcs, Government Arts College, Udumalpet, Tamlnadu

More information

Department of Economics, Niigata Sangyo University, Niigata, Japan

Department of Economics, Niigata Sangyo University, Niigata, Japan Appled Matheatcs, 0, 5, 777-78 Publshed Onlne March 0 n ScRes. http://www.scrp.org/journal/a http://d.do.org/0.6/a.0.507 On Relatons between the General Recurrence Forula o the Etenson o Murase-Newton

More information

Outline. Prior Information and Subjective Probability. Subjective Probability. The Histogram Approach. Subjective Determination of the Prior Density

Outline. Prior Information and Subjective Probability. Subjective Probability. The Histogram Approach. Subjective Determination of the Prior Density Outlne Pror Inforaton and Subjectve Probablty u89603 1 Subjectve Probablty Subjectve Deternaton of the Pror Densty Nonnforatve Prors Maxu Entropy Prors Usng the Margnal Dstrbuton to Deterne the Pror Herarchcal

More information

4 Analysis of Variance (ANOVA) 5 ANOVA. 5.1 Introduction. 5.2 Fixed Effects ANOVA

4 Analysis of Variance (ANOVA) 5 ANOVA. 5.1 Introduction. 5.2 Fixed Effects ANOVA 4 Analyss of Varance (ANOVA) 5 ANOVA 51 Introducton ANOVA ANOVA s a way to estmate and test the means of multple populatons We wll start wth one-way ANOVA If the populatons ncluded n the study are selected

More information

Non-Mixture Cure Model for Interval Censored Data: Simulation Study ABSTRACT

Non-Mixture Cure Model for Interval Censored Data: Simulation Study ABSTRACT Malaysan Journal of Mathematcal Scences 8(S): 37-44 (2014) Specal Issue: Internatonal Conference on Mathematcal Scences and Statstcs 2013 (ICMSS2013) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES Journal

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

1. Statement of the problem

1. Statement of the problem Volue 14, 010 15 ON THE ITERATIVE SOUTION OF A SYSTEM OF DISCRETE TIMOSHENKO EQUATIONS Peradze J. and Tsklaur Z. I. Javakhshvl Tbls State Uversty,, Uversty St., Tbls 0186, Georga Georgan Techcal Uversty,

More information

CIS526: Machine Learning Lecture 3 (Sept 16, 2003) Linear Regression. Preparation help: Xiaoying Huang. x 1 θ 1 output... θ M x M

CIS526: Machine Learning Lecture 3 (Sept 16, 2003) Linear Regression. Preparation help: Xiaoying Huang. x 1 θ 1 output... θ M x M CIS56: achne Learnng Lecture 3 (Sept 6, 003) Preparaton help: Xaoyng Huang Lnear Regresson Lnear regresson can be represented by a functonal form: f(; θ) = θ 0 0 +θ + + θ = θ = 0 ote: 0 s a dummy attrbute

More information

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests Smulated of the Cramér-von Mses Goodness-of-Ft Tests Steele, M., Chaselng, J. and 3 Hurst, C. School of Mathematcal and Physcal Scences, James Cook Unversty, Australan School of Envronmental Studes, Grffth

More information

Introducing Entropy Distributions

Introducing Entropy Distributions Graubner, Schdt & Proske: Proceedngs of the 6 th Internatonal Probablstc Workshop, Darstadt 8 Introducng Entropy Dstrbutons Noel van Erp & Peter van Gelder Structural Hydraulc Engneerng and Probablstc

More information

Parametric fractional imputation for missing data analysis. Jae Kwang Kim Survey Working Group Seminar March 29, 2010

Parametric fractional imputation for missing data analysis. Jae Kwang Kim Survey Working Group Seminar March 29, 2010 Parametrc fractonal mputaton for mssng data analyss Jae Kwang Km Survey Workng Group Semnar March 29, 2010 1 Outlne Introducton Proposed method Fractonal mputaton Approxmaton Varance estmaton Multple mputaton

More information

Chapter 20 Duration Analysis

Chapter 20 Duration Analysis Chapter 20 Duraton Analyss Duraton: tme elapsed untl a certan event occurs (weeks unemployed, months spent on welfare). Survval analyss: duraton of nterest s survval tme of a subject, begn n an ntal state

More information

Finite Fields and Their Applications

Finite Fields and Their Applications Fnte Felds and Ther Applcatons 5 009 796 807 Contents lsts avalable at ScenceDrect Fnte Felds and Ther Applcatons www.elsever.co/locate/ffa Typcal prtve polynoals over nteger resdue rngs Tan Tan a, Wen-Feng

More information

The Parity of the Number of Irreducible Factors for Some Pentanomials

The Parity of the Number of Irreducible Factors for Some Pentanomials The Party of the Nuber of Irreducble Factors for Soe Pentanoals Wolfra Koepf 1, Ryul K 1 Departent of Matheatcs Unversty of Kassel, Kassel, F. R. Gerany Faculty of Matheatcs and Mechancs K Il Sung Unversty,

More information

Interactive Bi-Level Multi-Objective Integer. Non-linear Programming Problem

Interactive Bi-Level Multi-Objective Integer. Non-linear Programming Problem Appled Mathematcal Scences Vol 5 0 no 65 3 33 Interactve B-Level Mult-Objectve Integer Non-lnear Programmng Problem O E Emam Department of Informaton Systems aculty of Computer Scence and nformaton Helwan

More information

Approximate Technique for Solving Class of Fractional Variational Problems

Approximate Technique for Solving Class of Fractional Variational Problems Appled Matheatcs, 5, 6, 837-846 Publshed Onlne May 5 n ScRes. http://www.scrp.org/journal/a http://dx.do.org/.436/a.5.6578 Approxate Technque for Solvng Class of Fractonal Varatonal Probles Ead M. Soloua,,

More information

Double Acceptance Sampling Plan for Time Truncated Life Tests Based on Transmuted Generalized Inverse Weibull Distribution

Double Acceptance Sampling Plan for Time Truncated Life Tests Based on Transmuted Generalized Inverse Weibull Distribution J. Stat. Appl. Pro. 6, No. 1, 1-6 2017 1 Journal of Statstcs Applcatons & Probablty An Internatonal Journal http://dx.do.org/10.18576/jsap/060101 Double Acceptance Samplng Plan for Tme Truncated Lfe Tests

More information

Study of the possibility of eliminating the Gibbs paradox within the framework of classical thermodynamics *

Study of the possibility of eliminating the Gibbs paradox within the framework of classical thermodynamics * tudy of the possblty of elnatng the Gbbs paradox wthn the fraework of classcal therodynacs * V. Ihnatovych Departent of Phlosophy, Natonal echncal Unversty of Ukrane Kyv Polytechnc Insttute, Kyv, Ukrane

More information

Numerical Solution of Nonlinear Singular Ordinary Differential Equations Arising in Biology Via Operational Matrix of Shifted Legendre Polynomials

Numerical Solution of Nonlinear Singular Ordinary Differential Equations Arising in Biology Via Operational Matrix of Shifted Legendre Polynomials Aercan Journal of Coputatonal and Appled Matheatcs: ; (): -9 DOI:.93/.aca..4 Nuercal Soluton of Nonlnear Sngular Ordnary Dfferental Equatons Arsng n Bology Va Operatonal Matr of Shfted Legendre Polynoals

More information

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng

More information

Markov Chain Monte-Carlo (MCMC)

Markov Chain Monte-Carlo (MCMC) Markov Chan Monte-Carlo (MCMC) What for s t and what does t look lke? A. Favorov, 2003-2017 favorov@sens.org favorov@gal.co Monte Carlo ethod: a fgure square The value s unknown. Let s saple a rando value

More information

Durban Watson for Testing the Lack-of-Fit of Polynomial Regression Models without Replications

Durban Watson for Testing the Lack-of-Fit of Polynomial Regression Models without Replications Durban Watson for Testng the Lack-of-Ft of Polynomal Regresson Models wthout Replcatons Ruba A. Alyaf, Maha A. Omar, Abdullah A. Al-Shha ralyaf@ksu.edu.sa, maomar@ksu.edu.sa, aalshha@ksu.edu.sa Department

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look

More information

Composite Hypotheses testing

Composite Hypotheses testing Composte ypotheses testng In many hypothess testng problems there are many possble dstrbutons that can occur under each of the hypotheses. The output of the source s a set of parameters (ponts n a parameter

More information

A Hybrid Variational Iteration Method for Blasius Equation

A Hybrid Variational Iteration Method for Blasius Equation Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method

More information

Handling Overload (G. Buttazzo, Hard Real-Time Systems, Ch. 9) Causes for Overload

Handling Overload (G. Buttazzo, Hard Real-Time Systems, Ch. 9) Causes for Overload PS-663: Real-Te Systes Handlng Overloads Handlng Overload (G Buttazzo, Hard Real-Te Systes, h 9) auses for Overload Bad syste desgn eg poor estaton of worst-case executon tes Sultaneous arrval of unexpected

More information

1.3 Hence, calculate a formula for the force required to break the bond (i.e. the maximum value of F)

1.3 Hence, calculate a formula for the force required to break the bond (i.e. the maximum value of F) EN40: Dynacs and Vbratons Hoework 4: Work, Energy and Lnear Moentu Due Frday March 6 th School of Engneerng Brown Unversty 1. The Rydberg potental s a sple odel of atoc nteractons. It specfes the potental

More information

A new Approach for Solving Linear Ordinary Differential Equations

A new Approach for Solving Linear Ordinary Differential Equations , ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

Chapter 11: Simple Linear Regression and Correlation

Chapter 11: Simple Linear Regression and Correlation Chapter 11: Smple Lnear Regresson and Correlaton 11-1 Emprcal Models 11-2 Smple Lnear Regresson 11-3 Propertes of the Least Squares Estmators 11-4 Hypothess Test n Smple Lnear Regresson 11-4.1 Use of t-tests

More information

Outline. Review Numerical Approach. Schedule for April and May. Review Simple Methods. Review Notation and Order

Outline. Review Numerical Approach. Schedule for April and May. Review Simple Methods. Review Notation and Order Sstes of Ordnar Dfferental Equatons Aprl, Solvng Sstes of Ordnar Dfferental Equatons Larr Caretto Mecancal Engneerng 9 Nuercal Analss of Engneerng Sstes Aprl, Outlne Revew bascs of nuercal solutons of

More information

Solving Fuzzy Linear Programming Problem With Fuzzy Relational Equation Constraint

Solving Fuzzy Linear Programming Problem With Fuzzy Relational Equation Constraint Intern. J. Fuzz Maeatcal Archve Vol., 0, -0 ISSN: 0 (P, 0 0 (onlne Publshed on 0 Septeber 0 www.researchasc.org Internatonal Journal of Solvng Fuzz Lnear Prograng Proble W Fuzz Relatonal Equaton Constrant

More information

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0 MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector

More information

Our focus will be on linear systems. A system is linear if it obeys the principle of superposition and homogenity, i.e.

Our focus will be on linear systems. A system is linear if it obeys the principle of superposition and homogenity, i.e. SSTEM MODELLIN In order to solve a control syste proble, the descrptons of the syste and ts coponents ust be put nto a for sutable for analyss and evaluaton. The followng ethods can be used to odel physcal

More information

MIMA Group. Chapter 2 Bayesian Decision Theory. School of Computer Science and Technology, Shandong University. Xin-Shun SDU

MIMA Group. Chapter 2 Bayesian Decision Theory. School of Computer Science and Technology, Shandong University. Xin-Shun SDU Group M D L M Chapter Bayesan Decson heory Xn-Shun Xu @ SDU School of Computer Scence and echnology, Shandong Unversty Bayesan Decson heory Bayesan decson theory s a statstcal approach to data mnng/pattern

More information

Numerical Solution of Ordinary Differential Equations

Numerical Solution of Ordinary Differential Equations Numercal Methods (CENG 00) CHAPTER-VI Numercal Soluton of Ordnar Dfferental Equatons 6 Introducton Dfferental equatons are equatons composed of an unknown functon and ts dervatves The followng are examples

More information

Research on Network Security Prediction Method Based on Kalman Filtering Fusion Decision Entropy Theory

Research on Network Security Prediction Method Based on Kalman Filtering Fusion Decision Entropy Theory Vol.39 (SecTech 206) pp.32-326 http://dx.do.org/0.4257/astl.206.39.66 Research on Networ Securty Predcton Method Based on Kalan Flterng Fuson Decson Entropy Theory Lang Huang Xnhao Chen and Xnsheng La

More information

Fuzzy Boundaries of Sample Selection Model

Fuzzy Boundaries of Sample Selection Model Proceedngs of the 9th WSES Internatonal Conference on ppled Mathematcs, Istanbul, Turkey, May 7-9, 006 (pp309-34) Fuzzy Boundares of Sample Selecton Model L. MUHMD SFIIH, NTON BDULBSH KMIL, M. T. BU OSMN

More information

CS 2750 Machine Learning. Lecture 5. Density estimation. CS 2750 Machine Learning. Announcements

CS 2750 Machine Learning. Lecture 5. Density estimation. CS 2750 Machine Learning. Announcements CS 750 Machne Learnng Lecture 5 Densty estmaton Mlos Hauskrecht mlos@cs.ptt.edu 539 Sennott Square CS 750 Machne Learnng Announcements Homework Due on Wednesday before the class Reports: hand n before

More information

BOOTSTRAP METHOD FOR TESTING OF EQUALITY OF SEVERAL MEANS. M. Krishna Reddy, B. Naveen Kumar and Y. Ramu

BOOTSTRAP METHOD FOR TESTING OF EQUALITY OF SEVERAL MEANS. M. Krishna Reddy, B. Naveen Kumar and Y. Ramu BOOTSTRAP METHOD FOR TESTING OF EQUALITY OF SEVERAL MEANS M. Krshna Reddy, B. Naveen Kumar and Y. Ramu Department of Statstcs, Osmana Unversty, Hyderabad -500 007, Inda. nanbyrozu@gmal.com, ramu0@gmal.com

More information