Objective Priors for Estimation of Extended Exponential Geometric Distribution
|
|
- Shona Boyd
- 6 years ago
- Views:
Transcription
1 Journal of Modern Appled Statstcal Methods Volume 13 Issue Artcle Objectve Prors for Estmaton of Extended Exponental Geometrc Dstrbuton Pedro L. Ramos Unversdade de São Paulo, São Paulo, Brazl, pedrolramos@hotmal.com Fernando A. Moala Unversdade Estadual Paulsta, São Paulo, Brazl, femoala@fct.unesp.br Jorge A. Achcar Unversdade de São Paulo, São Paulo, Brazl, achcar@fmrp.usp.br Follow ths and addtonal works at: Part of the Appled Statstcs Commons, Socal and Behavoral Scences Commons, and the Statstcal Theory Commons Recommended Ctaton Ramos, Pedro L.; Moala, Fernando A.; and Achcar, Jorge A. (014) "Objectve Prors for Estmaton of Extended Exponental Geometrc Dstrbuton," Journal of Modern Appled Statstcal Methods: Vol. 13 : Iss., Artcle 1. DOI: 10.37/jmasm/ Avalable at: Ths Regular Artcle s brought to you for free and open access by the Open Access Journals at DgtalCommons@WayneState. It has been accepted for ncluson n Journal of Modern Appled Statstcal Methods by an authorzed edtor of DgtalCommons@WayneState.
2 Journal of Modern Appled Statstcal Methods November 014, Vol. 13, No. 1, Copyrght 014 JMASM, Inc. ISSN Objectve Prors for Estmaton of Extended Exponental Geometrc Dstrbuton Pedro L. Ramos Unversdade de São Paulo São Paulo, Brazl Fernando A. Moala Unversdade Estadual Paulsta São Paulo, Brazl Jorge A. Achcar Unversdade de São Paulo São Paulo, Brazl A Bayesan analyss was developed wth dfferent nonnformatve pror dstrbutons such as Jeffreys, Maxmal Data Informaton, and Reference. The am was to nvestgate the effects of each pror dstrbuton on the posteror estmates of the parameters of the extended exponental geometrc dstrbuton, based on smulated data and a real applcaton. Keywords: Extended exponental geometrc dstrbuton, Jeffreys, MDIP, Reference, Bayesan, nonnformatve, pror Introducton Adamds & Loukas (005) ntroduced an extenson of the exponental geometrc dstrbuton (Adamds & Loukas, 1998), namng t as an extended exponental geometrc (EEG) dstrbuton, to analyze lfetme data. Ths dstrbuton provdes ncreasng or decreasng hazard functons, dependng on the values of ts parameters. In ths way, EEG gves a great flexblty of ft for the data. If T s a random varable denotng the lfetme of a component wth an extended exponental geometrc (EEG) dstrbuton, then the probablty densty s gven by: f t, e t t 11 e, (1) wth t > 0 and parameters γ > 0 and λ > 0. Let us denote ths dstrbuton as EEG( γ, λ ). Dr. Ramos s n the Faculty of Medcne of Rberão Preto. Emal hm at: pedrolramos@hotmal.com. Dr. Moala s a Professor n the Faculty of Scence and Technology. Emal hm at: femoala@fct.unesp.br. Dr. Achcar s a Professor n the Insttute of Mathematcs and Computer Scence. Emal hm at: achcar@fmrp.usp.br. 6
3 RAMOS ET AL. The survval and hazard functons of EEG( γ, λ ) dstrbuton, for a fxed tme t, s gven by t e S t;, and ht;, t 1 1 e 1 1 e t, () respectvely. The mean and varance of the EEG dstrbuton are gven, respectvely, by ET 1,1,1 and var 1,,1 1,1,1 T (3) where Ψ( z, s, a ) s known as Lerch transcendental functon (Erdely et al., 1953), gven by s1 au k 1 u e z z, s, a du ; z 1; a, s 0 s. u s 1 ze 0 k0 a k Adamds et al. (005) and Ktdamrongsuk (010) gave addtonal propertes of the EEG dstrbuton. Fgures 1 and present dfferent forms for the densty, survval and hazard functons for the EEG dstrbuton consderng dfferent values of γ and λ. The motvaton here s to present a Bayesan analyss when there s lttle pror knowledge avalable or that reflects manly the nformaton from the sample.. In ths stuaton, t s mportant to use nonnformatve prors, however, t can be dffcult to choose a pror dstrbuton that represent one of ths stuatons. Thus, the man am of ths paper s to choose a nonnformatve pror dstrbuton s for the parameters parameters λ and γ of the EEG dstrbuton and to study the effects of these dfferent prors n the resultng posteror dstrbutons, especally n stuatons of small sample szes, a common stuaton n applcatons. 7
4 OBJECTIVE PRIORS FOR ESTIMATION OF GEOMETRIC DISTRIBUTION Fgure 1. Probablty densty functons for the EEG dstrbuton wth values for the scale and shape parameters gven respectvely, by λ = 0.5, 1.0, 1.5 and γ = 0.5, 1.5,.0,.5, 3.0, 3.5. Fgure. Survval functons and hazard functons for the EEG dstrbuton wth values for the scale and shape parameters gven respectvely, by λ = 0.5, 1.0 and γ = 0.5, 1, 3, 5. Commonly used nonnformatve pror dstrbutons are derved, such as unform (Bayes, 1763; Laplace, 1774), Jeffreys (1967), reference prors (Bernardo, 1979; Berger & Bernardo, 199), and the uncommon MDIP pror (Zellner, 1977, 1984). A smulaton study s conducted comparng ther performance n terms of ther summares and coverage rates of credble ntervals. 8
5 RAMOS ET AL. Numercal ntegraton based on stochastc smulaton methods as the Markov Chan Monte Carlo (MCMC) wll be used to smulate samples of the margnal posteror dstrbuton of nterest. In partcular, we wll be usng the Metropols- Hastngs algorthm to obtan the posteror summares of nterest (see Gelfand & Smth, 1990 or Chb & Greenberg, 1995). Methodology Maxmum Lkelhood Estmaton Let X 1,, X n be a random sample from EEG( γ, λ ) dstrbuton wth densty (1). The lkelhood functon for the parameters γ and λ s gven by n Lx;, exp x 1 1 e 1 1 n n x, (4) where γ > 0 and λ > 0. The logarthm of the lkelhood functon (4) s gven by n n x. (5) l x;, nlog x log 1 1 e 1 1 By settng l(x; γ, λ ) / γ = 0 and l(x; γ, λ ) / λ = 0 and after some algebrac manpulatons, we obtan the lkelhood equatons n n x n xe x 1 0 and x 1 1 e n x n e 0 x 1 1 e (6) whose solutons provde the maxmum lkelhood estmators of the parameters γ and λ. Note that the solutons of the lkelhood equatons (6) cannot be obtaned analytcally and hence numercal approaches need to be used. Adamds et al. (005) propose to use the EM algorthm (Dempster et al., 1977) to solve the nonlnear equatons (6) and fnd the MLE of γ and λ. The EM teratons are gven by 9
6 OBJECTIVE PRIORS FOR ESTIMATION OF GEOMETRIC DISTRIBUTION ˆ j1 j ˆ ˆ j j j n 1 1 exp x n ˆ exp ˆ x j j n x 1 1 ˆ exp ˆ 1 x j n. j j 1 ˆ ˆ ˆ x exp 1 (7) Ktdamrongsuk (010) shows n detal the computatons of the expected Fsher nformaton matrx Ι(γ,λ) of the EEG dstrbuton, gven by I,, n 1 log I11, 3 1 n 1 log n (8) wth where I 11, (9) n 31 1 L 1, n 1 L 1 1 log, p Lr p 1, 0 p 1 r s a polygarthmca functon (Erdely et al., 1953). The maxmum lkelhood estmates for γ and λ are based for small sample problems. In the case of large samples they become unbased and asymptotcally effcent. Such estmates are asymptotcally Normal dstrbuted wth jont dstrbuton gven by 30
7 RAMOS ET AL. Bayesan Analyss ˆ 1 ˆ, ~ N,, I,, when n. (10) In ths secton we consder the Bayesan estmaton of the unknown parameters λ and γ. Frst, a pror dstrbuton whch expresses lttle nformaton on the parameters γ and λ can be obtaned from unform denstes, whch do not favor any partcular value of λ and γ. In ths case, the jont pror dstrbuton for λ and γ s gven by, constant. (11) U Another wdely-used method to specfy pror nformaton s through the product of ndependent gamma dstrbutons for each parameter λ and γ, snce γ > 0 and λ > 0, that s, γ ~ Gamma(α 1, β 1 ) and λ ~ Gamma(α, β ), where Gamma(a,b) denotes a gamma dstrbuton wth mean a/b and varance α/b ; and α 1, α, β 1 and β are known hyperparameters. Thus, the jont pror dstrbuton for λ and γ s gven by 1 1 1, exp. (1) G 1 Assume α 1 = α = β 1 = β = 0.01, that s, a non-nformatve pror gven by (1). An another well-known exstng non-nformatve pror, whch represents a stuaton wth lttle a pror nformaton on the parameters was ntroduced by Jeffreys (1967), also known as the Jeffreys rule. The Jeffreys pror has been wdely used due to the nvarance property for one to one transformatons of the parameters. The Jeffreys pror s defned as, det I,. (13) J where Ι(γ, λ) s the Fsher nformaton matrx defned n (8) and (9). From the equaton (13), we get 31
8 OBJECTIVE PRIORS FOR ESTIMATION OF GEOMETRIC DISTRIBUTION J, 3 L log 1 log 3 J, 1 ; 0 1, log 3 log L (14) 1 ; 1. It s nterestng to observe t was found n (14) ndependent prors for the parameters λ and γ, but ths jont pror has a dependence structure. Zellner (1977, 1984) proposed a non-nformatve pror based on the Shannon's entropy (1948). The dea s to maxmze the nformaton from the data n relaton to the pror nformaton on the parameters. Ths non-nformatve pror dstrbuton known as "Maxmal Data Informaton Pror" s obtaned from the soluton of the equaton Z, exp log f x, f x, dx. (15) 0 For the EEG dstrbuton gven n (1) the resultng non-nformatve pror s gven by log Z, exp 1,1,1, 1 (16) where Ψ(z,s,a) s defned n (5). The proposed Zellner pror dstrbuton (15) has lmted nvarance propertes, where nvarance s only verfed under lnear transformatons of the vector (γ, λ) and not for all dfferentable one by one transformatons. Bernardo (1979) and Berger & Bernardo (199) use the Kullback-Lebler dstance between the posteror dstrbuton p (θ x) the pror dstrbuton π (θ) to maxmze the nformaton from the data n relaton to the known pror nformaton for the parameters to fnd a nonnformatve pror. Addtonal nformaton about the reference pror can be found n Bernardo (005). 3
9 RAMOS ET AL. An mportant feature n ths approach s the dfferent treatment for nterest and nusance parameters when θ s a vector of parameters. In the presence of nusance parameters, a typcal case n ths paper, one must establsh an ordered parameterzaton wth the parameter of nterest sngled out and then follow the procedure below. The algorthm of Berger and Bernardo (199) to derve the reference pror can be descrbed n four steps, as follows. We wll present here the two-parameters case n detals. Let θ = (θ 1, θ ) be the two parameters vector; θ 1 wll be consdered the parameter of nterest and θ s the nusance parameter. The algorthm used to obtan the reference pror s gven by Step 1: gven by Fnd the condtonal reference pror π (θ θ 1 ), assumng that θ 1 s I,, (17) 1 1 where I (θ 1, θ ) s the term of order (,) of the nformaton Fsher matrx. Step : Normalze π (θ θ 1 ). If π (θ θ 1 ) s mproper, choose a sequence of sets Ω 1 Ω Ω, where π (θ θ 1 ) s proper. Fnd p m c m 1 1 (18) d m 1 c 1. (19) 1 m 1 1 m Step 3: Fnd the reference pror for θ 1. The result s gven by the soluton of 1 det I 1, m1 exp pm 1 log d. I 1, m (0) 33
10 OBJECTIVE PRIORS FOR ESTIMATION OF GEOMETRIC DISTRIBUTION Step 4: parameter Fnd the pror dstrbuton for (θ 1, θ ), when θ s the nusance c m 1 m 1 1, lm 1, m * * cm1 m1 (1) * where 1 s any fxed pont wthn the postve densty for all π m. For the EEG dstrbuton gven n (1), the reference pror when λ s the parameter of nterest s gven by, 1 L log L () The reference pror when γ s the parameter of nterest s gven by, for 0 1 and, for 1 3 L 1 1 log 1 log L 3log 3 log (3) Fnally, derve the pror dstrbutons for the parameters the resultng jont posteror dstrbutons for γ and λ s proportonal to the product of the lkelhood functon (4) and the pror dstrbutons π (γ, λ) gven n (11), (1), (14), (16), () and (3), that s n n n x,, exp 1 1 p x x e. (4)
11 RAMOS ET AL. By usng any pror dstrbuton proposed s not possble to derve the margnal posteror dstrbutons n an analytcal form for the parameters γ and λ. Thus, to obtan the posteror nformaton on the parameters of nterest as the pont estmator and Bayes credblty ntervals, we use MCMC algorthms to smulate samples of the values of γ and λ from the jont posteror dstrbutons. Results Two applcatons of the theoretcal results dscussed n the prevous sectons are presented. The frst nvolves a comparson of the estmaton methods based on smulated data; the second shows an applcaton of the EEG dstrbuton to real data. Analyss va numercal smulaton In ths example, some smulatons are performed va the Monte Carlo method. The goal s to study the effect of dfferent non-nformatve pror dstrbutons on the posteror summares and also to compare these results wth the obtaned results usng classcal nference analyss. Posteror summares of nterest are evaluated usng Monte Carlo Markov Chan (MCMC) methods. The nfluence of sample sze on the accuracy of the obtaned estmators s also examned. The followng procedure was adopted: 1. Determne the values of γ and λ.. Specfy the sample sze n. 3. Generate values of a dstrbuton EEG(γ,λ) wth sze n. 4. Usng the data obtaned n Step 3, calculate the estmates for the parameters γ and λ usng MCMC n the Bayesan approach and MLE n the classcal approach. 5. Repeat the steps 3 and 4 N tmes. Consder two set of the true values for the parameter (γ, λ) gven by (γ, λ) = (0.5, ) and (γ, λ) = (, 4) representng decreasng and ncreasng hazard functons, respectvely. The smulated data are generate from EEG dstrbuton wth the parameter values above for dfferent sample szes, as n = 10, 5 and 50. Tables 1 and show the posteror mean and medan, respectvely, by consderng the non-nformatve prors proposed n ths paper for the parameters. The maxmum lkelhood estmates (MLE) are also avalable. 35
12 OBJECTIVE PRIORS FOR ESTIMATION OF GEOMETRIC DISTRIBUTION Table 1. Posteror medans and MLE for λ = and γ = ½ for 1000 samples of szes 10, 5 and 50. λ = Jeffreys MDIP Ref. λ Ref. γ Unform Gamma MLE n=10.98(1.6).80(0.84).96(0.87).88(0.81) 3.61(1.54) 1.15(1.06) 3.44(.61) n=5.6(1.10) 3.00(1.05).93(0.88).9(0.83).95(1.00).01(1.06).54(1.6) n=50.17(0.76).5(0.70).58(0.67).71(0.67).56(0.7).07(0.77).6(0.80) γ = ½ Jeffreys MDIP Ref. λ Ref. γ Unform Gamma MLE n= (0.43) 1.31(0.77) 1.1(0.17) 1.13(0.14) 1.50(0.77) 0.35(0.49) 1.51(.0) n=5 0.87(0.39) 0.9(0.4) 0.95(0.3) 0.99(0.1) 0.97(0.46) 0.56(0.41) 0.81(0.68) n= (0.3) 0.71(0.7) 0.81(0.7) 0.89(0.6) 0.76(0.3) 0.56(0.9) 0.64(0.36) Table. Posteror medans and MLE for λ = 4 and γ = for 1000 samples of szes 10, 5 and 50. λ = 4 Jeffreys MDIP Ref. λ Ref. γ Unform Gamma MLE n= (1.1) 4.6(1.50) 5.31(1.55) 4.73(0.99) 5.08(1.5) 5.08(1.5) 4.86(.03) n=5 4.56(0.81) 4.0(0.77) 3.8(0.71) 3.90(0.73) 4.77(0.87) 3.50(0.84) 4.41(1.4) n= (0.57) 3.98(0.61) 3.51(0.41) 3.56(0.40) 3.40(0.64) 3.67(0.65) 4.(0.94) γ = Jeffreys MDIP Ref. λ Ref. γ Unform Gamma MLE n= (0.50) 1.81(0.40) 3.5(0.4) 3.33(0.40) 3.6(0.4) 3.6(0.4) 3.65(3.3) n=5.70(0.49) 1.76(0.34) 1.6(0.5) 1.74(0.7).94(0.53) 1.5(0.56).87(.06) n= (0.45) 1.90(0.41) 1.53(0.9) 1.60(0.3) 1.95(0.45) 1.70(0.5).43(1.16) From Tables 1 and, t s observed that when the hazard functon s decreasng (0 < γ < 1) the pror dstrbuton gven by product of ndependent gamma dstrbutons gves the best estmaton for the parameters whle for the ncreasng hazard functon (γ > 1) the MDIP pror dstrbuton provdes the best one for all sample szes consdered. A crteron for comparson of the pror dstrbutons conssts on checkng the frequentst coverage probabltes of the posteror ntervals. We therefore compare the frequency at whch the true values of γ and λ are ncluded n ther 95% posteror ntervals. Ths frequency should be close to 95% for large numbers of repeated experments. 36
13 RAMOS ET AL. Table 3. Coverage probabltes for λ = and γ = ½ for 1000 samples of szes 10, 5 and 50. λ = Jeffreys MDIP Ref. λ Ref. γ Unform Gamma MLE n= % 99.10% 95.70% 96.70% 90.60% 96.50% 95.0% n=5 95.0% 91.60% 89.70% 90.00% 93.40% 97.0% 95.00% n= % 94.30% 95.00% 91.30% 9.0% 96.60% 94.70% γ = ½ Jeffreys MDIP Ref. λ Ref. γ Unform Gamma MLE n= % 98.00% 95.70% 94.80% 90.60% 97.10% 94.60% n= % 97.80% 97.90% 97.70% 95.00% 97.80% 95.50% n= % 97.90% 96.00% 95.50% 93.70% 97.70% 95.40% Table 4. Coverage probabltes for λ = 4 and γ = for 1000 samples of szes 10, 5 and 50. λ = 4 Jeffreys MDIP Ref. λ Ref. γ Unform Gamma MLE n= % 96.60% 88.10% 98.00% 91.0% 91.0% 96.60% n= % 98.0% 96.80% 97.00% 94.00% 95.60% 94.30% n= % 98.40% 96.00% 97.60% 98.50% 96.60% 95.00% γ = Jeffreys MDIP Ref. λ Ref. γ Unform Gamma MLE n= % 99.50% 93.80% 97.00% 96.60% 96.60% 9.50% n= % 99.10% 99.50% 99.90% 98.50% 97.10% 9.00% n= % 98.90% 97.00% 99.50% 98.80% 96.50% 93.90% An example wth lterature data Now consder a lfetme dataset related to an electrcal nsulator subjected to constant stress and stran, ntroduced by Lawless (198). The dataset does not have censored values and represents the lfetme (n mnutes) to falure: 0.96, 4.15, 0.19, 0.78, 8.01, 31.75, 7.35, 6.50, 8.7, 33.91, 3.5, 16.03, 4.85,.78, 4.67, 1.31, 1.06, and Assume that the EEG dstrbuton s approprated to analyze ths dataset, and then t wll be compared wth other lfetme dstrbutons such as Webull, Gamma, and Lognormal. As shown, the effcency of the dfferent non-nformatve pror dstrbutons changes wth the shapes of the hazard functons, therefore, to get good nferences on parameters of nterest t s necessary to have some pror nformaton on how the hazard functon behaves for the Lawless data set. In ths way, Barlow & Campo (1975) proposed a smple graphcal technque that has been wdely used to verfy the behavor of the rsk functon called TTT plot (total tme for testng). The graph s constructed wth the plot of the consecutve quanttes [r/n, G (r/n)], where G (r/n) s a functon gven by 37
14 OBJECTIVE PRIORS FOR ESTIMATION OF GEOMETRIC DISTRIBUTION r 1 t n r t r G r / n, t where r = 1,..., n and t (), = 1,, n are the order statstcs n the sample. Usng the TTT curve n an emprcal scale, one can determne the shape of the hazard functon for the lfetme data. A dagonal lne ndcates that the data have a constant hazard functon; f the curve s convex the rsk s decreasng; f t s concave, there s an ndcaton that the rsk s ncreasng; f frst s convex and after ths s concave then there s an ndcaton that there s a bathtub shape for the hazard functon; f t s frst concave and after ths convex, there s an ndcaton of nverse form of the bath for the hazard functon. The Fgure 3 shows how to verfy the behavor of the hazard functon. n 1 Fgure 3. TTT plots for dfferent dstrbutons ndcatng the shape of the hazard functon. Some TTT transformatons can be studed to solve other problems. Nar et al. (008) show some of these transformatons appled n survval analyss. Fgure 4 shows the TTT plot for the Lawless data set. 38
15 RAMOS ET AL. Fgure 4. TTT plot for the dataset lfetme related to an electrcal nsulator subjected to constant stress and stran (Lawless data). It s observed n Fgure 4 that the TTT plot s convex; then t can be concluded that the rsk s a decreasng functon. When the hazard functon s decreasng, t was observed from the results of secton 4.1, that non-nformatve prors obtaned through the product of ndependent gamma dstrbutons s the best pror wth lttle pror nformaton about the parameters of nterest. The jont posteror dstrbuton of λ and γ (4) s obtaned by replacng π(γ, λ) by (1). It s necessary to use numercal methods to extract nformaton from the margnal posteror dstrbutons λ and γ. MCMC methods are used to smulate samples for the jont posteror dstrbuton; that s, also for the margnal posteror dstrbutons of nterest. It was generated 110,000 teratons wth a burn-n of 10,000 values and jumps of sze 10; so we get chans of the margnal posteror dstrbutons for λ and γ of sze 10,000 obtaned usng MCMC methods. To verfy the convergence of the chans, we have used Geweke (199) dagnoss, whch ndcated the convergence of the two chans. The convergence and autocorrelatons s also observed n the trace-plots of the smulated seres gven n Fgure 5. To verfy the performance of other lfetme dstrbutons we also consder as non-nformatve pror, the product of ndependent gamma dstrbutons γ ~ Gamma (0.01, 0.01), λ ~ Gamma (0.01, 0.01) assumng the followng lfetme dstrbutons: EGE, Webull, Gamma and Lognormal dstrbuton. The results are compled n Tables 5 and 6, respectvely. 39
16 OBJECTIVE PRIORS FOR ESTIMATION OF GEOMETRIC DISTRIBUTION Fgure 5. Trace-plots and autocorrelaton graphs for the generated values of λ and γ. Table 5. Posteror estmates (means) for γ and λ consderng dfferent probablty dstrbutons for the Lawless data set. Parameters EGE Webull Gamma Log-normal γ (0.0161) 0.769(0.1356) 0.675(0.1870) (0.3969) λ (0.1649) 7.850(4.50) (0.0183).6410(1.000) Table 6. Obtaned results for the DIC, BIC and AIC crtera for the dfferent probablty dstrbutons for the Lawless data set. Crtera EGE Webull Gamma Log-normal DIC BIC AIC Based on any of the crtera used by the table t can be concluded that EGE was the best ft to the offered data. Concluson The use of extended exponental geometrc (EEG) dstrbutons showed a good flexblty of ft for lfetme data applcatons and could be an alternatve 40
17 RAMOS ET AL. dstrbuton to other usual dstrbutons n Survval analyss. The great number of exstng non-nformatve pror dstrbutons can cause dffcultes n the choce of an adequate pror wth lttle nformaton a pror, manly when these pror dstrbutons do not produce smlar posteror summares. In ths way, the development of a general theory for the constructon of non-nformatve pror dstrbutons s an mportant topc to be nvestgated by researchers n the Bayesan nference. The results showed the effects of dfferent non-nformatve pror dstrbutons related to the changes n the rsk functon usng extended exponental geometrc (EEG) dstrbutons. Therefore, we recommend the product of gamma dstrbutons Gamma (0.01, 0.01) when the hazard functon s decreasng and the nonnformatve MDIP pror dstrbuton when the hazard functon s ncreasng. Wth these choces of pror we surely get better nferences for the parameters. References Achcar, J. A, Moala, F. A, & Boleta, J. (March, 010). Generalzed Exponental Dstrbuton: A Bayesan approach usng MCMC methods. Poster presented at the 10th Bayesan Statstcs Brazlan Meetng, Ro de Janero, Brazl. Adamds, K, Dmtrakopoulou, T, & Loukas, S. (005). On an Extenson of the Exponental Geometrc Dstrbuton. Statstcs and Probablty Letters, 73, Adamds, K, & Loukas, S. (1998). A Lfetme Dstrbuton wth Decreasng Falure Rate. Statstcs and Probablty Letters, 39, Barlow, R. E., & Campo, R. A. (1975). Total Tme on Test processes and applcatons to falure data analyss. In R. E. Barlow, J. B. Fussel, & N. D. Sngpurwalla (Eds.), Relablty and Fault Tree Analyss: Theoretcal and Appled Aspects of System Relablty and Safety Assessment ( ). Phladelpha, PA: Socety for Industral and Appled Mathematcs. Bayes, T. R. (1958). Essay towards solvng a problem n the doctrne of changes. Reprnted n Bometrka, 45, , Berger, J. O, & Bernardo, J. M. (199). On the Development of the Reference Pror Method. Fourth Valenca Internatonal Meetng on Bayesan Statstcs, Span. 41
18 OBJECTIVE PRIORS FOR ESTIMATION OF GEOMETRIC DISTRIBUTION Bernardo, J. M. (1979). Reference Posteror Dstrbutons for Bayesan Inference. Journal Royal Statstcal Socety, 41(), Bernardo, J. M. (005). Reference analyss. In D.K. Dey and C. R. Rao, (Eds.), Handbook of Statstcs, 5 (17-90). Amsterdam: Elsever. Box, G. E. P., & Tao, G. C. (1973). Bayesan nference n statstcal analyss. Readng, MA: Addson-Wesley. Chb, S., & Greenberg, E. (1995). Understandng The Metropols-Hastng Algorthm. The Amercan Statstcan, 49(4), Degroot, M. H., & Schervsh, M. J. (00). Probablty and Statstcs. Readng, MA: Addson-Wesley. Dempster, A. P., Lard, N. M., & Rubn, D. B. (1977). Maxmum Lkelhood from Incomplete Data va the EM Algorthm (wth dscusson). Journal of the Royal Statstcal Socety: Seres B (Methodologcal), 39(1), Erdely, A., Maguns, W., Oberhettnger, F., & Trcom, F.G. (1953). Hgher Transcendental Functons. New York: McGraw-Hll. Gelfand, A. E., & Smth, F. M. (1990). Samplng-based approaches to calculatng margnal denstes. Journal of the Amercan Statstcal Assocaton, 85, Geweke, J. (199). Evaluatng the accuracy of samplng-based approaches to calculatng posteror moments. In J. M. Bernado, J. O. Berger, A. P. Dawd, & A. F. M. Smth (Eds.), Bayesan Statstcs 4. Oxford, UK: Clarendon Press. Gupta, R. D., & Kundu, D. (1999). Generalzed exponental dstrbutons. Australan and New Zealand Journal of Statstcs, 41, Gupta, R. D., & Kundu, D. (001). Generalzed exponental dstrbutons: dfferent methods of estmaton. Journal of Statstcal Computaton and Smulaton, 69, Gupta, R. D., & Kundu, D. (007). Generalzed exponental dstrbuton: exstng results and some recent developments. Journal of Statstcal Plannng and Inference, 137(11), do: /j.jsp Gupta, R. D., & Kundu, D. (008). Generalzed exponental dstrbuton; Bayesan Inference. Computatonal Statstcs and Data Analyss, 5(4), Hamada, M., Wlson, A. G., Reese, C. S., & Martz, H. F. (008). Bayesan Relablty. New York: Sprnger. Jeffreys, H. (1967). Theory of probablty (3rd ed.). London: Oxford Unversty Press. 4
19 RAMOS ET AL. Ktdamrongsuk, P. (010). Dscrmnatng Between the Extended Exponental Geometrc Dstrbuton and the Gamma Dstrbuton. (Unpublshed Doctoral Dssertaton). Assumpton Unversty of Thaland. Laplace, P. (1774). Mémore sur la probablté des causes par les é venemens. Mem. Acad. R. Sc. Presentés par Dvers Savans, v. 6, p (translated n Statstcal Scence, 1, ). Lawless, J. F. (198). Statstcal models and methods for lfetme data. New York: John Wley and Sons. Marshall, A.W., & Olkn, I. (1997). A New Method for Addng a Parameter to a Famly of Dstrbutons wth Applcaton to the Exponental and Webull Famles. Bometrka, 84, Moala, F. A. (010). Bayesan analyss for the Webull parameters by usng nonnformatve pror dstrbutons. Advances and Applcatons n Statstcs, 14(), Nar, N. U., Sankaran, P. G., & Kumar, B. V., (008). Total tme on test transforms of order n and ts mplcatons n relablty analyss. Journal of Appled Probablty, 45(4), Raqab, M. Z. (00). Inferences for generalzed exponental dstrbuton based on record statstcs. Journal of Statstcal Plannng and Inference, 104, Raqab, M. Z., & Ahsanullah, M. (001). Estmaton of the locaton and scale parameters of generalzed exponental dstrbuton based on order statstcs. Journal of Statstcal Computaton and Smulaton, 69, Sarhan, A. M. (007). Analyss of ncomplete, censored data n competng rsks models wth generalzed exponental dstrbutons. IEEE Transactons on Relablty, 56(1), Shannon, C. E. (1948). A Mathematcal theory of communcaton. Bell System Techncal Journal, 7, Zellner, A. (1977). Maxmal Data Informaton Pror Dstrbutons. In A. Aykac & C. Brumat (Eds.), New Methods n the applcatons of Bayesan Methods. Amsterdam: North-Holland Publshng Co. Zellner, A. (1984). Maxmal Data Informaton Pror Dstrbutons: Basc Issues n Econometrcs. Chcago, IL: Unversty of Chcago Press. Zheng, G. (00). Fsher nformaton matrx n type-ii censored data from exponentated exponental famly. Bometrcal Journal, 44,
International Journal of Industrial Engineering Computations
Internatonal Journal of Industral Engneerng Computatons 4 (03) 47 46 Contents lsts avalable at GrowngScence Internatonal Journal of Industral Engneerng Computatons homepage: www.growngscence.com/ec Modelng
More informationOn 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 informationStatistical 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 informationComparison 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 informationInterval 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 informationCOMPARISON 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 informationEcon107 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 informationParametric 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 informationComputing MLE Bias Empirically
Computng MLE Bas Emprcally Kar Wa Lm Australan atonal Unversty January 3, 27 Abstract Ths note studes the bas arses from the MLE estmate of the rate parameter and the mean parameter of an exponental dstrbuton.
More informationChange 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 informationStat260: Bayesian Modeling and Inference Lecture Date: February 22, Reference Priors
Stat60: Bayesan Modelng and Inference Lecture Date: February, 00 Reference Prors Lecturer: Mchael I. Jordan Scrbe: Steven Troxler and Wayne Lee In ths lecture, we assume that θ R; n hgher-dmensons, reference
More informationA 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 informationParameters Estimation of the Modified Weibull Distribution Based on Type I Censored Samples
Appled Mathematcal Scences, Vol. 5, 011, no. 59, 899-917 Parameters Estmaton of the Modfed Webull Dstrbuton Based on Type I Censored Samples Soufane Gasm École Supereure des Scences et Technques de Tuns
More informationANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)
Econ 413 Exam 13 H ANSWERS Settet er nndelt 9 deloppgaver, A,B,C, som alle anbefales å telle lkt for å gøre det ltt lettere å stå. Svar er gtt . Unfortunately, there s a prntng error n the hnt of
More informationJoint Statistical Meetings - Biopharmaceutical Section
Iteratve Ch-Square Test for Equvalence of Multple Treatment Groups Te-Hua Ng*, U.S. Food and Drug Admnstraton 1401 Rockvlle Pke, #200S, HFM-217, Rockvlle, MD 20852-1448 Key Words: Equvalence Testng; Actve
More informationComputation of Higher Order Moments from Two Multinomial Overdispersion Likelihood Models
Computaton of Hgher Order Moments from Two Multnomal Overdsperson Lkelhood Models BY J. T. NEWCOMER, N. K. NEERCHAL Department of Mathematcs and Statstcs, Unversty of Maryland, Baltmore County, Baltmore,
More informationGlobal Sensitivity. Tuesday 20 th February, 2018
Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values
More informationMarkov Chain Monte Carlo (MCMC), Gibbs Sampling, Metropolis Algorithms, and Simulated Annealing Bioinformatics Course Supplement
Markov Chan Monte Carlo MCMC, Gbbs Samplng, Metropols Algorthms, and Smulated Annealng 2001 Bonformatcs Course Supplement SNU Bontellgence Lab http://bsnuackr/ Outlne! Markov Chan Monte Carlo MCMC! Metropols-Hastngs
More informationSimulated 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 informationEstimation of the Mean of Truncated Exponential Distribution
Journal of Mathematcs and Statstcs 4 (4): 84-88, 008 ISSN 549-644 008 Scence Publcatons Estmaton of the Mean of Truncated Exponental Dstrbuton Fars Muslm Al-Athar Department of Mathematcs, Faculty of Scence,
More informationNon-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 information2E 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 information3.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 informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationUsing 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 informationDERIVATION OF THE PROBABILITY PLOT CORRELATION COEFFICIENT TEST STATISTICS FOR THE GENERALIZED LOGISTIC DISTRIBUTION
Internatonal Worshop ADVANCES IN STATISTICAL HYDROLOGY May 3-5, Taormna, Italy DERIVATION OF THE PROBABILITY PLOT CORRELATION COEFFICIENT TEST STATISTICS FOR THE GENERALIZED LOGISTIC DISTRIBUTION by Sooyoung
More informationRELIABILITY 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 informationLOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin
Proceedngs of the 007 Wnter Smulaton Conference S G Henderson, B Bller, M-H Hseh, J Shortle, J D Tew, and R R Barton, eds LOW BIAS INTEGRATED PATH ESTIMATORS James M Calvn Department of Computer Scence
More informationComposite 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 informationDouble 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 informationA 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 information4 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 informationChapter 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 informationPower law and dimension of the maximum value for belief distribution with the max Deng entropy
Power law and dmenson of the maxmum value for belef dstrbuton wth the max Deng entropy Bngy Kang a, a College of Informaton Engneerng, Northwest A&F Unversty, Yanglng, Shaanx, 712100, Chna. Abstract Deng
More informationConjugacy and the Exponential Family
CS281B/Stat241B: Advanced Topcs n Learnng & Decson Makng Conjugacy and the Exponental Famly Lecturer: Mchael I. Jordan Scrbes: Bran Mlch 1 Conjugacy In the prevous lecture, we saw conjugate prors for the
More informationLecture Notes on Linear Regression
Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume
More informationStat 543 Exam 2 Spring 2016
Stat 543 Exam 2 Sprng 206 I have nether gven nor receved unauthorzed assstance on ths exam. Name Sgned Date Name Prnted Ths Exam conssts of questons. Do at least 0 of the parts of the man exam. I wll score
More informationParametric fractional imputation for missing data analysis
Secton on Survey Research Methods JSM 2008 Parametrc fractonal mputaton for mssng data analyss Jae Kwang Km Wayne Fuller Abstract Under a parametrc model for mssng data, the EM algorthm s a popular tool
More informationMATH 829: Introduction to Data Mining and Analysis The EM algorithm (part 2)
1/16 MATH 829: Introducton to Data Mnng and Analyss The EM algorthm (part 2) Domnque Gullot Departments of Mathematcal Scences Unversty of Delaware Aprl 20, 2016 Recall 2/16 We are gven ndependent observatons
More informationMarkov Chain Monte Carlo Lecture 6
where (x 1,..., x N ) X N, N s called the populaton sze, f(x) f (x) for at least one {1, 2,..., N}, and those dfferent from f(x) are called the tral dstrbutons n terms of mportance samplng. Dfferent ways
More informationANOMALIES OF THE MAGNITUDE OF THE BIAS OF THE MAXIMUM LIKELIHOOD ESTIMATOR OF THE REGRESSION SLOPE
P a g e ANOMALIES OF THE MAGNITUDE OF THE BIAS OF THE MAXIMUM LIKELIHOOD ESTIMATOR OF THE REGRESSION SLOPE Darmud O Drscoll ¹, Donald E. Ramrez ² ¹ Head of Department of Mathematcs and Computer Studes
More informationModule 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 informationStat 543 Exam 2 Spring 2016
Stat 543 Exam 2 Sprng 2016 I have nether gven nor receved unauthorzed assstance on ths exam. Name Sgned Date Name Prnted Ths Exam conssts of 11 questons. Do at least 10 of the 11 parts of the man exam.
More informationHere is the rationale: If X and y have a strong positive relationship to one another, then ( x x) will tend to be positive when ( y y)
Secton 1.5 Correlaton In the prevous sectons, we looked at regresson and the value r was a measurement of how much of the varaton n y can be attrbuted to the lnear relatonshp between y and x. In ths secton,
More informationSmall Area Interval Estimation
.. Small Area Interval Estmaton Partha Lahr Jont Program n Survey Methodology Unversty of Maryland, College Park (Based on jont work wth Masayo Yoshmor, Former JPSM Vstng PhD Student and Research Fellow
More informationTest for Intraclass Correlation Coefficient under Unequal Family Sizes
Journal of Modern Appled Statstcal Methods Volume Issue Artcle 9 --03 Test for Intraclass Correlaton Coeffcent under Unequal Famly Szes Madhusudan Bhandary Columbus State Unversty, Columbus, GA, bhandary_madhusudan@colstate.edu
More informationBayesian Planning of Hit-Miss Inspection Tests
Bayesan Plannng of Ht-Mss Inspecton Tests Yew-Meng Koh a and Wllam Q Meeker a a Center for Nondestructve Evaluaton, Department of Statstcs, Iowa State Unversty, Ames, Iowa 5000 Abstract Although some useful
More informationA Note on Bound for Jensen-Shannon Divergence by Jeffreys
OPEN ACCESS Conference Proceedngs Paper Entropy www.scforum.net/conference/ecea- A Note on Bound for Jensen-Shannon Dvergence by Jeffreys Takuya Yamano, * Department of Mathematcs and Physcs, Faculty of
More informationBOOTSTRAP 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 informationAssignment 5. Simulation for Logistics. Monti, N.E. Yunita, T.
Assgnment 5 Smulaton for Logstcs Mont, N.E. Yunta, T. November 26, 2007 1. Smulaton Desgn The frst objectve of ths assgnment s to derve a 90% two-sded Confdence Interval (CI) for the average watng tme
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationDETERMINATION OF UNCERTAINTY ASSOCIATED WITH QUANTIZATION ERRORS USING THE BAYESIAN APPROACH
Proceedngs, XVII IMEKO World Congress, June 7, 3, Dubrovn, Croata Proceedngs, XVII IMEKO World Congress, June 7, 3, Dubrovn, Croata TC XVII IMEKO World Congress Metrology n the 3rd Mllennum June 7, 3,
More informationOn the correction of the h-index for career length
1 On the correcton of the h-ndex for career length by L. Egghe Unverstet Hasselt (UHasselt), Campus Depenbeek, Agoralaan, B-3590 Depenbeek, Belgum 1 and Unverstet Antwerpen (UA), IBW, Stadscampus, Venusstraat
More informationThe EM Algorithm (Dempster, Laird, Rubin 1977) The missing data or incomplete data setting: ODL(φ;Y ) = [Y;φ] = [Y X,φ][X φ] = X
The EM Algorthm (Dempster, Lard, Rubn 1977 The mssng data or ncomplete data settng: An Observed Data Lkelhood (ODL that s a mxture or ntegral of Complete Data Lkelhoods (CDL. (1a ODL(;Y = [Y;] = [Y,][
More informationRockefeller College University at Albany
Rockefeller College Unverst at Alban PAD 705 Handout: Maxmum Lkelhood Estmaton Orgnal b Davd A. Wse John F. Kenned School of Government, Harvard Unverst Modfcatons b R. Karl Rethemeer Up to ths pont n
More informationDegradation Data Analysis Using Wiener Process and MCMC Approach
Engneerng Letters 5:3 EL_5_3_0 Degradaton Data Analyss Usng Wener Process and MCMC Approach Chunpng L Hubng Hao Abstract Tradtonal relablty assessment methods are based on lfetme data. However the lfetme
More informationLINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity
LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 30 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 2 Remedes for multcollnearty Varous technques have
More informationBayesian 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 informationIntroduction to Generalized Linear Models
INTRODUCTION TO STATISTICAL MODELLING TRINITY 00 Introducton to Generalzed Lnear Models I. Motvaton In ths lecture we extend the deas of lnear regresson to the more general dea of a generalzed lnear model
More informationStatistics for Managers Using Microsoft Excel/SPSS Chapter 13 The Simple Linear Regression Model and Correlation
Statstcs for Managers Usng Mcrosoft Excel/SPSS Chapter 13 The Smple Lnear Regresson Model and Correlaton 1999 Prentce-Hall, Inc. Chap. 13-1 Chapter Topcs Types of Regresson Models Determnng the Smple Lnear
More informationA PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS
HCMC Unversty of Pedagogy Thong Nguyen Huu et al. A PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS Thong Nguyen Huu and Hao Tran Van Department of mathematcs-nformaton,
More informationTHE ROYAL STATISTICAL SOCIETY 2006 EXAMINATIONS SOLUTIONS HIGHER CERTIFICATE
THE ROYAL STATISTICAL SOCIETY 6 EXAMINATIONS SOLUTIONS HIGHER CERTIFICATE PAPER I STATISTICAL THEORY The Socety provdes these solutons to assst canddates preparng for the eamnatons n future years and for
More informationChapter 3. Two-Variable Regression Model: The Problem of Estimation
Chapter 3. Two-Varable Regresson Model: The Problem of Estmaton Ordnary Least Squares Method (OLS) Recall that, PRF: Y = β 1 + β X + u Thus, snce PRF s not drectly observable, t s estmated by SRF; that
More informationAn (almost) unbiased estimator for the S-Gini index
An (almost unbased estmator for the S-Gn ndex Thomas Demuynck February 25, 2009 Abstract Ths note provdes an unbased estmator for the absolute S-Gn and an almost unbased estmator for the relatve S-Gn for
More informationA Note on Test of Homogeneity Against Umbrella Scale Alternative Based on U-Statistics
J Stat Appl Pro No 3 93- () 93 NSP Journal of Statstcs Applcatons & Probablty --- An Internatonal Journal @ NSP Natural Scences Publshng Cor A Note on Test of Homogenety Aganst Umbrella Scale Alternatve
More informationHidden Markov Models & The Multivariate Gaussian (10/26/04)
CS281A/Stat241A: Statstcal Learnng Theory Hdden Markov Models & The Multvarate Gaussan (10/26/04) Lecturer: Mchael I. Jordan Scrbes: Jonathan W. Hu 1 Hdden Markov Models As a bref revew, hdden Markov models
More informationDr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur
Analyss of Varance and Desgn of Experment-I MODULE VIII LECTURE - 34 ANALYSIS OF VARIANCE IN RANDOM-EFFECTS MODEL AND MIXED-EFFECTS EFFECTS MODEL Dr Shalabh Department of Mathematcs and Statstcs Indan
More informationOn the Influential Points in the Functional Circular Relationship Models
On the Influental Ponts n the Functonal Crcular Relatonshp Models Department of Mathematcs, Faculty of Scence Al-Azhar Unversty-Gaza, Gaza, Palestne alzad33@yahoo.com Abstract If the nterest s to calbrate
More informationPOWER AND SIZE OF NORMAL DISTRIBUTION AND ITS APPLICATIONS
Jurnal Ilmah Matematka dan Penddkan Matematka (JMP) Vol. 9 No., Desember 07, hal. -6 ISSN (Cetak) : 085-456; ISSN (Onlne) : 550-04; https://jmpunsoed.com/ POWER AND SIZE OF NORMAL DISTRIBION AND ITS APPLICATIONS
More informationPopulation Design in Nonlinear Mixed Effects Multiple Response Models: extension of PFIM and evaluation by simulation with NONMEM and MONOLIX
Populaton Desgn n Nonlnear Mxed Effects Multple Response Models: extenson of PFIM and evaluaton by smulaton wth NONMEM and MONOLIX May 4th 007 Carolne Bazzol, Sylve Retout, France Mentré Inserm U738 Unversty
More informationOn Distribution Function Estimation Using Double Ranked Set Samples With Application
Journal of Modern Appled Statstcal Methods Volume Issue Artcle -- On Dstrbuton Functon Estmaton Usng Double Raned Set Samples Wth Applcaton Wald A. Abu-Dayyeh Yarmou Unversty, Irbd Jordan, abudayyehw@yahoo.com
More informationFirst Year Examination Department of Statistics, University of Florida
Frst Year Examnaton Department of Statstcs, Unversty of Florda May 7, 010, 8:00 am - 1:00 noon Instructons: 1. You have four hours to answer questons n ths examnaton.. You must show your work to receve
More informationSTAT 3008 Applied Regression Analysis
STAT 3008 Appled Regresson Analyss Tutoral : Smple Lnear Regresson LAI Chun He Department of Statstcs, The Chnese Unversty of Hong Kong 1 Model Assumpton To quantfy the relatonshp between two factors,
More informationChapter 5 Multilevel Models
Chapter 5 Multlevel Models 5.1 Cross-sectonal multlevel models 5.1.1 Two-level models 5.1.2 Multple level models 5.1.3 Multple level modelng n other felds 5.2 Longtudnal multlevel models 5.2.1 Two-level
More information1. Inference on Regression Parameters a. Finding Mean, s.d and covariance amongst estimates. 2. Confidence Intervals and Working Hotelling Bands
Content. Inference on Regresson Parameters a. Fndng Mean, s.d and covarance amongst estmates.. Confdence Intervals and Workng Hotellng Bands 3. Cochran s Theorem 4. General Lnear Testng 5. Measures of
More informationEstimation: Part 2. Chapter GREG estimation
Chapter 9 Estmaton: Part 2 9. GREG estmaton In Chapter 8, we have seen that the regresson estmator s an effcent estmator when there s a lnear relatonshp between y and x. In ths chapter, we generalzed the
More informationISSN X Robust bayesian inference of generalized Pareto distribution
Afrka Statstka Vol. 112), 2016, pages 1061 1074. DOI: http://dx.do.org/10.16929/as/2016.1061.92 Afrka Statstka ISSN 2316-090X Robust bayesan nference of generalzed Pareto dstrbuton Fatha Mokran 1, Hocne
More informationA quantum-statistical-mechanical extension of Gaussian mixture model
A quantum-statstcal-mechancal extenson of Gaussan mxture model Kazuyuk Tanaka, and Koj Tsuda 2 Graduate School of Informaton Scences, Tohoku Unversty, 6-3-09 Aramak-aza-aoba, Aoba-ku, Senda 980-8579, Japan
More informationLINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity
LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 31 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 6. Rdge regresson The OLSE s the best lnear unbased
More informationEconomics 130. Lecture 4 Simple Linear Regression Continued
Economcs 130 Lecture 4 Contnued Readngs for Week 4 Text, Chapter and 3. We contnue wth addressng our second ssue + add n how we evaluate these relatonshps: Where do we get data to do ths analyss? How do
More informationLinear Approximation with Regularization and Moving Least Squares
Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...
More informationTurbulence classification of load data by the frequency and severity of wind gusts. Oscar Moñux, DEWI GmbH Kevin Bleibler, DEWI GmbH
Turbulence classfcaton of load data by the frequency and severty of wnd gusts Introducton Oscar Moñux, DEWI GmbH Kevn Blebler, DEWI GmbH Durng the wnd turbne developng process, one of the most mportant
More informationUNIVERSITY OF TORONTO Faculty of Arts and Science. December 2005 Examinations STA437H1F/STA1005HF. Duration - 3 hours
UNIVERSITY OF TORONTO Faculty of Arts and Scence December 005 Examnatons STA47HF/STA005HF Duraton - hours AIDS ALLOWED: (to be suppled by the student) Non-programmable calculator One handwrtten 8.5'' x
More informationarxiv:cs.cv/ Jun 2000
Correlaton over Decomposed Sgnals: A Non-Lnear Approach to Fast and Effectve Sequences Comparson Lucano da Fontoura Costa arxv:cs.cv/0006040 28 Jun 2000 Cybernetc Vson Research Group IFSC Unversty of São
More informationCredit Card Pricing and Impact of Adverse Selection
Credt Card Prcng and Impact of Adverse Selecton Bo Huang and Lyn C. Thomas Unversty of Southampton Contents Background Aucton model of credt card solctaton - Errors n probablty of beng Good - Errors n
More informationSOME NEW CHARACTERIZATION RESULTS ON EXPONENTIAL AND RELATED DISTRIBUTIONS. Communicated by Ahmad Reza Soltani. 1. Introduction
Bulletn of the Iranan Mathematcal Socety Vol. 36 No. 1 (21), pp 257-272. SOME NEW CHARACTERIZATION RESULTS ON EXPONENTIAL AND RELATED DISTRIBUTIONS M. TAVANGAR AND M. ASADI Communcated by Ahmad Reza Soltan
More informationObjective Bayesian Analysis for the Lomax Distribution arxiv: v1 [stat.me] 26 Feb 2016
Objectve Bayesan Analyss for the Lomax Dstrbuton arxv:62.845v [stat.me] 26 Feb 26 Paulo Ferrera, Jhon Gonzales, Vera Tomazella, Rcardo Ehlers, Francsco Louzada, Evelny Slva February 26 Abstract In ths
More informationQuantum and Classical Information Theory with Disentropy
Quantum and Classcal Informaton Theory wth Dsentropy R V Ramos rubensramos@ufcbr Lab of Quantum Informaton Technology, Department of Telenformatc Engneerng Federal Unversty of Ceara - DETI/UFC, CP 6007
More informationChapter 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 information4.3 Poisson Regression
of teratvely reweghted least squares regressons (the IRLS algorthm). We do wthout gvng further detals, but nstead focus on the practcal applcaton. > glm(survval~log(weght)+age, famly="bnomal", data=baby)
More informationTesting for outliers in nonlinear longitudinal data models based on M-estimation
ISS 1746-7659, England, UK Journal of Informaton and Computng Scence Vol 1, o, 017, pp107-11 estng for outlers n nonlnear longtudnal data models based on M-estmaton Huhu Sun 1 1 School of Mathematcs and
More informationAn adaptive SMC scheme for ABC. Bayesian Computation (ABC)
An adaptve SMC scheme for Approxmate Bayesan Computaton (ABC) (ont work wth Prof. Mke West) Department of Statstcal Scence - Duke Unversty Aprl/2011 Approxmate Bayesan Computaton (ABC) Problems n whch
More informationLECTURE 9 CANONICAL CORRELATION ANALYSIS
LECURE 9 CANONICAL CORRELAION ANALYSIS Introducton he concept of canoncal correlaton arses when we want to quantfy the assocatons between two sets of varables. For example, suppose that the frst set of
More informationx i1 =1 for all i (the constant ).
Chapter 5 The Multple Regresson Model Consder an economc model where the dependent varable s a functon of K explanatory varables. The economc model has the form: y = f ( x,x,..., ) xk Approxmate ths by
More informationHidden Markov Models
CM229S: Machne Learnng for Bonformatcs Lecture 12-05/05/2016 Hdden Markov Models Lecturer: Srram Sankararaman Scrbe: Akshay Dattatray Shnde Edted by: TBD 1 Introducton For a drected graph G we can wrte
More informationsince [1-( 0+ 1x1i+ 2x2 i)] [ 0+ 1x1i+ assumed to be a reasonable approximation
Econ 388 R. Butler 204 revsons Lecture 4 Dummy Dependent Varables I. Lnear Probablty Model: the Regresson model wth a dummy varables as the dependent varable assumpton, mplcaton regular multple regresson
More informationWhy Bayesian? 3. Bayes and Normal Models. State of nature: class. Decision rule. Rev. Thomas Bayes ( ) Bayes Theorem (yes, the famous one)
Why Bayesan? 3. Bayes and Normal Models Alex M. Martnez alex@ece.osu.edu Handouts Handoutsfor forece ECE874 874Sp Sp007 If all our research (n PR was to dsappear and you could only save one theory, whch
More informationLecture 16 Statistical Analysis in Biomaterials Research (Part II)
3.051J/0.340J 1 Lecture 16 Statstcal Analyss n Bomaterals Research (Part II) C. F Dstrbuton Allows comparson of varablty of behavor between populatons usng test of hypothess: σ x = σ x amed for Brtsh statstcan
More informationAndreas C. Drichoutis Agriculural University of Athens. Abstract
Heteroskedastcty, the sngle crossng property and ordered response models Andreas C. Drchouts Agrculural Unversty of Athens Panagots Lazards Agrculural Unversty of Athens Rodolfo M. Nayga, Jr. Texas AMUnversty
More informationChapter 13: Multiple Regression
Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to
More information