Exponentiated Kumaraswamy-Dagum distribution with applications to income and lifetime data
|
|
- Ethelbert Small
- 5 years ago
- Views:
Transcription
1 Huang and Oluyede Journal of Statstcal Dstrbutons and Applcatons :8 RESEARCH Open Access Exponentated Kumaraswamy-Dagum dstrbuton wth applcatons to ncome and lfetme data Shujao Huang and Broderck O Oluyede * *Correspondence: boluyede@georgasouthern.edu Department of Mathematcal Scences Georga Southern Unversty Statesboro GA USA Abstract A new famly of dstrbutons called exponentated Kumaraswamy-Dagum (EKD) dstrbuton s proposed and studed. Ths famly ncludes several well known sub-models such as Dagum (D) Burr III (BIII) Fsk or Log-logstc (F or LLog) and new sub-models namely Kumaraswamy-Dagum (KD) Kumaraswamy-Burr III (KBIII) Kumaraswamy-Fsk or Kumaraswamy-Log-logstc (KF or KLLog) exponentated Kumaraswamy-Burr III (EKBIII) and exponentated Kumaraswamy-Fsk or exponentated Kumaraswamy-Log-logstc (EKF or EKLLog) dstrbutons. Statstcal propertes ncludng seres representaton of the probablty densty functon hazard and reverse hazard functons moments mean and medan devatons relablty Bonferron and Lorenz curves as well as entropy measures for ths class of dstrbutons and the sub-models are presented. Maxmum lkelhood estmates of the model parameters are obtaned. Smulaton studes are conducted. Examples and applcatons as well as comparsons of the EKD and ts sub-dstrbutons wth other dstrbutons are gven. Mathematcs Subject Classfcaton (2000): 62E10; 62F30 Keywords: Dagum dstrbuton; Exponentated Kumaraswamy-Dagum dstrbuton; Maxmum lkelhood estmaton 1 Introducton Camlo Dagum proposed the dstrbuton whch s referred to as Dagum dstrbuton n Ths proposal enable the development of statstcal dstrbutons used to ft emprcal ncome and wealth data that could accommodate heavy tals n ncome and wealth dstrbutons. Dagum s proposed dstrbuton has both Type-I and Type-II specfcaton where Type-I s the three parameter specfcaton and Type-II deals wth four parameter specfcaton. Ths dstrbuton s a specal case of generalzed beta dstrbuton of the second knd (GB2) McDonald (1984) McDonald and Xu (1995) when the parameter q = 1 where the probablty densty functon (pdf) of the GB2 dstrbuton s gven by: ay ap 1 f GB2 (y; a b p q) = b ap B(p q)[1 + ( y ) a fory > 0. b ] p+q See Kleber and Kotz (2003) for detals. Note that a > 0 p > 0 q > 0 are the shape parameters and b s the scale parameter and B(p q) = Ɣ(p)Ɣ(q) Ɣ(p+q) s the beta functon. Kleber (2008) traced the geness of Dagum dstrbuton and summarzed several 2014 Huang and Oluyede; lcensee Sprnger. Ths s an Open Access artcle dstrbuted under the terms of the Creatve Commons Attrbuton Lcense ( whch permts unrestrcted use dstrbuton and reproducton n any medum provded the orgnal work s properly credted.
2 Huang and Oluyede Journal of Statstcal Dstrbutons and Applcatons :8 Page 2 of 20 statstcal propertes of ths dstrbuton. Domma et al. (2011) obtaned the maxmum lkelhood estmates of the parameters of Dagum dstrbuton for censored data. Domma and Condno (2013) presented the beta-dagum dstrbuton. Cordero et al. (2013) proposed the beta exponentated Webull dstrbuton. Cordero et al. (2010) ntroduced and studed some mathematcal propertes of the Kumaraswamy Webull dstrbuton. Oluyede and Rajasoorya (2013) developed the Mc-Dagum dstrbuton and presented ts statstcal propertes. See references theren for addtonal results. The pdf and cumulatve dstrbuton functon (cdf) of Dagum dstrbuton are gven by: g D (x; λ β δ) = βλδx δ 1 ( ) β 1 (1) and G D (x; λ β δ) = ( ) β (2) for x > 0 where λ s a scale parameter δ and β are shape parameters. Dagum (1977) refers to hs model as the generalzed logstc-burr dstrbuton. The k th raw or non central moments are gven by E (X k) ( = βλ k δ B β + k δ 1 k ) δ for k <δandλ β>0 where B( ) s the beta functon. The q th percentle s x q = λ 1 δ ( q 1 β 1) 1 δ. In ths paper we present generalzatons of the Dagum dstrbuton va Kumaraswamy dstrbuton and ts exponentated verson. Ths leads to the exponentated Kumaraswamy Dagum dstrbuton. The motvaton for the development of ths dstrbuton s the modelng of sze dstrbuton of personal ncome and lfetme data wth a dverse model that takes nto consderaton not only shape and scale but also skewness kurtoss and tal varaton. Also the EKD dstrbuton and ts sub-models has desrable features of exhbtng a non-monotone falure rate thereby accommodatng dfferent shapes for the hazard rate functon and should be an attractve choce for survval and relablty data analyss. Ths paper s organzed as follows. In secton 3 we present the exponentated Kumaraswamy-Dagum dstrbuton and ts sub models as well as seres expanson hazard and reverse hazard functons. Moments moment generatng functon Lorenz and Bonferron curves mean and medan devatons and relablty are obtaned n secton 4. Secton 5 contans results on the dstrbuton of the order statstcs and Reny entropy. Estmaton of model parameters va the method of maxmum lkelhood s presented n secton 6. In secton 7 varous smulatons are conducted for dfferent sample szes. Secton 8 contans examples and applcatons of the EKD dstrbuton and ts sub-models followed by concludng remarks. 2 Methods results and dscussons Methods results and dscussons for the class of EKD dstrbutons are presented n sectons 3 to 8. These sectons nclude the sub-models seres expanson of the pdf closed form expressons for the hazard and reverse hazard functons moments moment generatng functon Bonferron and Lorenz curves relablty mean and medan devatons
3 Huang and Oluyede Journal of Statstcal Dstrbutons and Applcatons :8 Page 3 of 20 dstrbuton of order statstcs and entropy as well as estmaton of model parameters and applcatons. 3 The exponentated Kumaraswamy-Dagum dstrbuton In ths secton we present the proposed dstrbuton and ts sub-models. Seres expanson hazard and reverse hazard functons are also studed n ths secton. 3.1 Kumaraswamy-Dagum dstrbuton Kumaraswamy (1980) ntroduced a two-parameter dstrbuton on (0 1). Itscdfsgven by G(x) = ( x ψ) φ x (0 1) for ψ>0andφ>0. For an arbtrary cdf F(x) wth pdf f (x) = df(x) dx the famly of Kumaraswamy-G dstrbutons wth cdf G k (x) s gven by G K (x) = ( F ψ (x) ) φ for ψ>0andφ>0. By lettng F(x) = G D (x) we obtan the Kumaraswamy-Dagum (KD) dstrbuton wth cdf ( ) G KD (x) = G ψ φ D (x). 3.2 The EKD dstrbuton In general the EKD dstrbuton s G EKD (x) = [F KD (x)] θ wheref KD (x) s a baselne (Kum-Dagum) cdf θ > 0 wth the correspondng pdf gven by g EKD (x) = θ[f KD (x)] θ 1 f KD (x). For large values of x andforθ>1(< 1) the multplcatve factor θ[f KD (x)] θ 1 > 1(< 1) respectvely. The reverse statement holds for smaller values of x. Consequently ths mples that the ordnary moments of g EKD (x) are larger (smaller) than those of f KD (x) when θ>1(< 1). Replacng the dependent parameter βψ by α the cdf and pdf of the EKD dstrbuton are gven by G EKD (x; α λ δ φ θ) = { ) α ] } φ θ (3) and g EKD (x; α λ δ φ θ) = αλδφθx δ 1 ( ) α 1 [ ( ) α ] φ 1 { ) α ] } φ θ 1 (4) for α λ δ φ θ>0 and x > 0 respectvely. The quantle functon of the EKD dstrbuton s n closed form [ ( ) G 1 (q) = x EKD q = λ 1 δ q θ 1 1 ] 1 α φ 1 1 δ. (5) Plots of the pdf for some combnatons of values of the model parameters are gven n Fgure 1. The plots ndcate that the EKD pdf can be decreasng or rght skewed. The EKD dstrbuton has a postve asymmetry.
4 Huang and Oluyede Journal of Statstcal Dstrbutons and Applcatons :8 Page 4 of 20 Fgure 1 Graph of pdfs. 3.3 Some sub-models Sub-models of EKD dstrbuton for selected values of the parameters are presented n ths secton. 1. When θ = 1 we obtan Kumaraswamy-Dagum dstrbuton wth cdf: [ G(x; α λ δ φ) = ( ) α ] φ for α λ δ φ>0 and x > When φ = θ = 1 we obtan Dagum dstrbuton wth cdf: G(x; α λ δ) = ( ) α for α λ δ>0 and x > When λ = 1 we obtan exponentated Kumaraswamy-Burr III dstrbuton wth cdf: G(x; α δ φ θ) = { 1 + x δ) α ] } φ θ for α δ φ θ>0and x > When λ = θ = 1 we obtan Kumaraswamy-Burr III dstrbuton wth cdf: [ G(x; α δ φ) = ( 1 + x δ) α ] φ for α δ φ>0 and x > When λ = φ = θ = 1 we obtan Burr III dstrbuton wth cdf: G(x; α δ) = ( 1 + x δ) α for α δ>0and x > When α = 1 we obtan exponentated Kumaraswamy-Fsk or Kumaraswamy-Log-logstc dstrbuton wth cdf: { G(x; λ δ φ θ) = ) 1 ] } φ θ for λ δ φ θ>0 and x > 0.
5 Huang and Oluyede Journal of Statstcal Dstrbutons and Applcatons :8 Page 5 of When α = θ = 1 we obtan Kumaraswamy-Fsk or Kumaraswamy-Log-logstc dstrbuton wth cdf: [ G(x; λ δ φ) = ( ) 1 ] φ for λ δ φ>0 and x > When α = φ = θ = 1 we obtan Fsk or Log-logstc dstrbuton wth cdf: G(x; λ δ) = ( ) 1 for λ δ>0 and x > Seres expanson of EKD dstrbuton We apply the seres expanson ( z) b 1 = j=0 ( 1) j Ɣ(b) Ɣ(b j)j! zj (6) for b > 0and z < 1 to obtan the seres expanson of the EKD dstrbuton. By usng equaton (6) g EKD (x) = =0 j=0 ω( j)x δ 1 ( ) α(j+1) 1 (7) where ω( j) = αλδφθ ( 1)+j Ɣ(θ)Ɣ(φ+φ) Ɣ(θ )Ɣ(φ+φ j)!j!. Note that n the Dagum(α δ λ) dstrbuton α and δ are shape parameters and λ s a scale parameter. In the Exponentated-Kumaraswamy(ψ φ θ) dstrbuton ψ s a skewness parameter φ s a tal varaton parameter and the parameter θ characterzes the skewness kurtoss and tal of the dstrbuton. Consequently for the EKD(α λ δ φ θ) dstrbuton α s shape and skewness parameter δ s shape parameter λ s a scale parameter φ s a tal varaton parameter and the parameter θ characterzes the skewness kurtoss and tal of the dstrbuton. 3.5 Hazard and reverse hazard functon The hazard functon of the EKD dstrbuton s h EKD (x) = g (x) EKD G EKD (x) = αλδφθx δ 1 ( ) α 1 [ ( ) α ] φ 1 [ {1 ( ) α ] φ } θ 1 ( { ) α ] } ) φ θ 1. (8) Plots of the hazard functon are presented n Fgure 2. The plots show varous shapes ncludng monotoncally decreasng unmodal and bathtub followed by upsde down bathtub shapes wth fve combnatons of the values of the parameters. Ths attractve flexblty makes the EKD hazard rate functon useful and sutable for non-monotone emprcal hazard behavors whch are more lkely to be encountered or observed n real lfe stuatons. Unfortunately the analytcal analyss of the shape of both the densty (except for zero modal when αδ 1 and unmodal f αδ > 1 both for φ = θ = 1) and
6 Huang and Oluyede Journal of Statstcal Dstrbutons and Applcatons :8 Page 6 of 20 Fgure 2 Graphs of hazard functons. hazard rate functon seems to be very complcated. We could not determne any specfc rules for the shapes of the hazard rate functon. The reverse hazard functon of the EKD dstrbuton s τ EKD (x) = g EKD (x) G EKD (x) = αλδφθx δ 1 ( ) α 1 [ ( ) α ] φ 1 { ) α ] } φ 1. (9) 4 Moments moment generatng functon Bonferron and Lorenz curves mean and medan devatons and relablty In ths secton we present the moments moment generatng functon Bonferron and Lorenz curves mean and medan devatons as well as the relablty of the EKD dstrbuton. The moments of the sub-models can be readly obtaned from the general results. 4.1 Moments and moment generatng functon Let t = ( ) 1 n equaton (7) then the s th raw moment of the EKD dstrbuton s gven by E(X s ) = = = 0 =0 j=0 =0 j=0 x s g EKD (x)dx ω( j)λ δ s 1 1 ( δ B α(j + 1) + s δ 1 s ) δ ( ω( j s)b α(j + 1) + s δ 1 s ) δ (10) where ω( j s) = αλ δ s φθ ( 1)+j Ɣ(θ)Ɣ(φ+φ) Ɣ(θ )Ɣ(φ+φ j)!j! ands <δ.
7 Huang and Oluyede Journal of Statstcal Dstrbutons and Applcatons :8 Page 7 of 20 The moment generatng functon of the EKD dstrbuton s gven by ( M(t) = ω( j r) tr r! B α(j + 1) + r δ 1 r ) δ for r <δ. r=0 =0 j=0 4.2 Bonferron and Lorenz curves Bonferron and Lorenz curves are wdely used tool for analyzng and vsualzng ncome nequalty. Lorenz curve L(p) can be regarded as the proporton of total ncome volume accumulated bythoseunts wth ncome lower than or equal to the volume a and Bonferron curve B(p) s the scaled condtonal mean curve that s rato of group mean ncome of the populaton. Plots of Bonferron and Lorenz curves are gven n Fgure 3. Let I(a) = a 0 x g EKD (x)dx and μ = E(X) then Bonferron and Lorenz curves are gven by B(p) = I(q) I(q) and L(p) = pμ μ respectvely for 0 p 1 and q = G 1 EKD (p). The mean of the EKD dstrbuton s obtaned from equaton (10) wth s = 1 and the quantle functon s gven n equaton (5). Consequently ( I(a) = ω( j1)b t(a) α(j + 1) + 1 δ 1 1 ) (11) δ =0 j=0 for δ > 1 where t(a) = (1 + λa δ ) 1 andb G(x) (c d) = G(x) 0 t c 1 ( t) d 1 dt for G(x) < 1 s ncomplete Beta functon. 4.3 Mean and medan devatons If X has the EKD dstrbuton we can derve the mean devaton about the mean μ = E(X) and the medan devaton about the medan M from δ 1 = x μ g EKD (x)dx and δ 2 = x M g EKD (x)dx 0 respectvely. The mean μ s obtaned from equaton (10) wth s = 1 and the medan M s gvenbyequaton(5)whenq = Fgure 3 Graphs of Bonferron and Lorenz curves.
8 Huang and Oluyede Journal of Statstcal Dstrbutons and Applcatons :8 Page 8 of 20 The measure δ 1 and δ 2 can be calculated by the followng relatonshps: δ 1 = 2μ G EKD (μ) 2μ + 2T(μ) and δ 2 = 2T(M) μ where T(a) = a x g EKD (x)dx follows from equaton (11) that s T(a) = ω( j1) B α(j + 1) + 1 δ 1 1 δ =0 j=0 ) B t(a) ( α(j + 1) + 1 δ 1 1 δ 4.4 Relablty The relablty R = P(X 1 > X 2 ) when X 1 and X 2 have ndependent EKD(α 1 λ 1 δ 1 φ 1 θ 1 ) and EKD(α 2 λ 2 δ 2 φ 2 θ 2 ) dstrbutons s gven by R = = 0 g 1 (x)g 2 (x)dx ζ( j k l) =0 j=0 k=0 l=0 0 )]. x δ 1 1 ( 1 + λ 1 x δ 1 ) α 1 (j+1) 1 ( 1 + λ2 x δ 2 ) α 2 l dx ( 1) +j+k+l Ɣ(θ 1 )Ɣ(φ 1 +φ 1 )Ɣ(θ 2 +1)Ɣ(φ 2 k+1) Ɣ(θ )Ɣ(φ 1 +φ j)ɣ(θ 2 +1 k)ɣ(φ 2 k+1 l)!j!k!l!. where ζ( j k l) = α 1 λ 1 δ 1 φ 1 θ 1 If λ = λ 1 = λ 2 and δ = δ 1 = δ 2 then relablty can be reduced to ζ( j k l) R = λδ[α 1 (j + 1) + α 2 l]. j=1 k=1 l=1 5 Order statstcs and entropy In ths secton the dstrbuton of the k th order statstc and Reny entropy (Reny 1960) for the EKD dstrbuton are presented. The entropy of a random varable s a measure of varaton of the uncertanty. 5.1 Order statstcs The pdf of the k th order statstcs from a pdf f (x) s f (x) f k:n (x) = B(k n k + 1) Fk 1 (x)[ F(x)] n k ( ) n = k f (x)f k 1 (x)[ F(x)] n k. (12) k Usng equaton (6) the pdf of the k th order statstc from EKD dstrbuton s gven by g k:n (x) = K( j p k) x δ 1 ( ) α αp 1 =0 j=0 p=0 where K( j p k) = 5.2 Entropy ( 1)+j+p Ɣ(n k+1)ɣ(θk+θ)ɣ(φj+φ) Ɣ(n k+1 )Ɣ(θk+θ j)ɣ(φj+φ p)!j!p! k( n k) αλδφθ. Reny entropy of a dstrbuton wth pdf f (x) s defned as { } I R (τ) = ( τ) 1 log f τ (x)dx τ>0 τ = 1. R Usng equaton (6) Reny entropy of EKD dstrbuton s gven by I R (τ) = ( τ) 1 log ( 1) +j Ɣ(θτ τ + 1)Ɣ(φτ τ + φ + 1) Ɣ(θτ τ + )Ɣ(φτ τ + φ + j)! j! =0 j=0 ( α τ λ τ δ + 1 δ δ τ 1 φ τ θ τ B ατ + αj + τ τ + τ 1 δ δ ) ]. for ατ + αj + 1 τ δ obtaned. > 0andτ + τ 1 δ > 0. Reny entropy for the sub-models can be readly
9 Huang and Oluyede Journal of Statstcal Dstrbutons and Applcatons :8 Page 9 of 20 6 Estmaton of model parameters In ths secton we present estmates of the parameters of the EKD dstrbuton va method of maxmum lkelhood estmaton. The elements of the score functon are presented. There are no closed form solutons to the nonlnear equatons obtaned by settng the elements of the score functon to zero. Thus the estmates of the model parameters must be obtaned va numercal methods. 6.1 Maxmum lkelhood estmaton Let x = (x 1 x n ) T be a random sample of the EKD dstrbuton wth unknown parameter vector = (α λ δ φ θ) T. The log-lkelhood functon for s n l( ) = n (ln α + ln λ + ln δ + ln φ + ln θ) (δ + 1) ln x (α + 1) + (θ 1) n ( ln { n ln ) + (φ 1) n ln ) ] } α φ [ ( ) α ]. (13) The partal dervatves of l( ) wth respect to the parameters are ( ) α ( ) l α = n n ( ) n α ln ln + (φ 1) ( ) α ) ] α φ 1 ( ) α ( ) n ln (θ 1)φ ) ] α φ l λ = n n λ (α + 1) (θ 1)φα n x δ + (φ 1)α ( ) α 1 n x δ ( ) α ) ] α φ 1 ( ) α 1 x δ ) ] α φ l δ = n n n δ x δ ln x ln x + (α + 1)λ ( ) α 1 n x δ ln x (φ 1)αλ ( ) α n + (θ 1)φαλ l φ = n n φ + and ln 1+λx δ ) α ] (θ 1) { l θ = n n θ + ln ) α ] φ 1 ( n ) α 1 x δ ln x ) ] α φ ) ] α φ ln ) ] } α φ ) α ] ) ] α φ
10 Huang and Oluyede Journal of Statstcal Dstrbutons and Applcatons :8 Page 10 of 20 respectvely. The MLE of the parameters α λ δ φ andθ say ˆα ˆλ ˆδ ˆφ and ˆθ mustbe obtaned by numercal methods. 6.2 Asymptotc confdence ntervals In ths secton we present the asymptotc confdence ntervals for the parameters of the EKD dstrbuton. The expectatons n the Fsher Informaton Matrx (FIM) can be obtaned numercally. Let ˆ = ( ˆα ˆλ ˆδ ˆφ ˆθ) be the maxmum lkelhood estmate of = (α λ δ φ θ). Under the usual regularty condtons and that the parameters are n the nteror of the parameter space but not on the boundary we have: d n( ˆ ) N 5 (0 I 1 ( )) wherei( ) s the expected Fsher nformaton matrx. The asymptotc behavor s stll vald f I( ) s replaced by the observed nformaton matrx evaluated at ˆ thatsj( ˆ ). The multvarate normal dstrbuton N 5 (0 J( ˆ ) 1) where the mean vector 0 = ( ) T can be used to construct confdence ntervals and confdence regons for the ndvdual model parameters and for the survval and hazard rate functons. The approxmate 100( η)% two-sded confdence ntervals for α λ δ φ and θ are gven by: α ± Z η Iαα 1 ( ) 2 λ ± Z η I 1 2 λλ ( ) δ ± Z η I 1 2 δδ ( ) φ ± Z η I 1 2 φφ ( ) θ ± Z η I 1 2 θθ ( ) respectvely where Z η s the upper η th 2 2 percentle of a standard normal dstrbuton. We can use the lkelhood rato (LR) test to compare the ft of the EKD dstrbuton wth ts sub-models for a gven data set. For example to test θ = 1 the LR statstc s ( )) ( ( ))] ω = 2 ln L ˆα ˆλ ˆδ ˆφ ˆθ ln L α λ δ φ 1 where ˆα ˆλ ˆδ ˆφ and ˆθ are the unrestrcted estmates and α λ δ and φ are the restrcted estmates.thelrtestrejectsthenullhypothessfω>χ 2whereχ2 denote the upper d d 100d% pont of the χ 2 dstrbuton wth 1 degrees of freedom. 7 Smulaton study In ths secton we examne the performance of the EKD dstrbuton by conductng varous smulatons for dfferent szes (n= ) va the subroutne NLP n SAS. We smulate 2000 samples for the true parameters values I : α = 2 λ = 1 δ = 3 φ = 2 θ = 2andII : α = 1 λ = 1 δ = 1 φ = 1 θ = 1. Table 1 lsts the means MLEs of the fve model parameters along wth the respectve root mean squared errors (RMSE). From the results we can verfy that as the sample sze n ncreases the mean estmates of the parameters tend to be closer to the true parameter values snce RMSEs decay toward zero. 8 Applcaton: EKD and sub-dstrbutons In ths secton applcatons based on real data as well as comparson of the EKD dstrbuton wth ts sub-models are presented. We provde examples to llustrate the flexblty of the EKD dstrbuton n contrast to other models ncludng the exponentated
11 Huang and Oluyede Journal of Statstcal Dstrbutons and Applcatons :8 Page 11 of 20 Table 1 Monte Carlo smulaton results: mean estmates and RMSEs I II n Parameter Mean RMSE Mean RMSE 200 α λ δ φ θ α λ δ φ θ α λ δ φ θ α λ δ φ θ Kumaraswamy-Webull (EKW) and beta-kumaraswamy-webull (BKW) dstrbutons for data modelng. The pdfs of EKW and BKW dstrbutons are [ f EKW (x) = θabcλ c x c 1 e (λx)c e (λx)c] a 1 { [ e (λx)c] a} b 1 [ { [ e (λx)c] ] a} b θ 1 and f BKW (x) = 1 B(a b) αβcλc x c 1 e (λx)c [ e (λx)c] α 1 { [ e (λx)c] α} βb 1 [ { [ e (λx)c] ] α} β a 1 respectvely. The frst data set conssts of the number of successve falures for the ar condtonng system of each member n a fleet of 13 Boeng 720 jet arplanes (Proschan 1963). The data s presented n Table 2. The second data set conssts of the salares of 818 professonal baseball players for the year 2009 (USA TODAY). The thrd data set represents the poverty rate of 533 dstrcts wth more than students n 2009 (Dgest of Educaton Statstcs d11/tables/dt11_096.asp ). These data sets are modeled by the EKD dstrbuton and compared wth the correspondng sub-models the Kumaraswamy-Dagum and Dagum
12 Huang and Oluyede Journal of Statstcal Dstrbutons and Applcatons :8 Page 12 of 20 Table 2 Ar condtonng system data dstrbutons and as well as EKW BKW dstrbutons. Table 3 gves a descrptve summary of each sample. The ar condtonng system sample has far more varablty and the baseball player salary sample has the lowest varablty. The maxmum lkelhood estmates (MLEs) of the parameters are computed by maxmzng the objectve functon va the subroutne NLMIXED n SAS. The estmated values of the parameters (standard error n parenthess) -2 Log-lkelhood statstc Akake Informaton Crteron AIC = 2p 2ln(L) Bayesan Informaton Crteron BIC = p ln(n) 2ln(L) and Consstent Akake Informaton Crteron AICC = AIC + 2 p(p+1) n p 1 where L = L( ˆ ) s the value of the lkelhood functon evaluated at the parameter estmates n s the number of observatons and p s the number of estmated parameters for the EKD dstrbuton and ts sub-dstrbutons are tabulated. See Table 4 Table 5 and Table 6. Plots of the ftted EKD KD D and the hstogram of the data are gven n Fgure 4. The probablty plots (Chambers et al. 1983) conssts of plots of the observed probabltes aganst the probabltes predcted by the ftted model are also presented n Fgure 5. For the EKD dstrbuton we plotted for example { G(x (j) ) = 1 + ˆλx ˆδ (j) ) ˆα ] ˆφ } ˆθ aganst j n+0.25 j = n wherex (j) are the ordered values of the observed data. A measure of closeness of the plot to the dagonal lne gven by the sum of squares SS = n j=1 [ G(x (j) ) ( j n )] 2 Table 3 Descrptve statstcs Data Mean Medan Mode SD Varance Skewness Kurtoss Mn. Max. I II III
13 Table 4 Estmaton of models for ar condtonng system data Estmates Statstcs Dstrbuton α λ δ φ θ -2 log lkelhood AIC AICC BIC SS EKD (1.2347) (0.4174) (0.0459) (5.8028) (0.0089) KD (2.1177) (3.0727) (0.1253) ( ) - D (0.1749) ( ) (0.0663) - - a b c λ θ EKW (0.8783) (0.0183) (0.1448) (0.0224) (0.1136) a b α β c λ BKW (1.2507) (0.0875) (1.6573) (1.9807) (0.0343) (0.0388) Huang and Oluyede Journal of Statstcal Dstrbutons and Applcatons :8 Page 13 of 20
14 Table 5 Estmaton of models for baseball player salary data Estmates Statstcs Dstrbuton α λ δ φ θ -2 log lkelhood AIC AICC BIC SS EKD ( ) ( ) (0.0557) (0.0044) (0.0327) KD ( ) ( ) (0.037) (0.0036) - D ( ) (0.0058) (0.0301) - - a b c λ θ EKW (2.0692) (0.0266) (0.0756) (4.9198) (0.2098) a b α β c λ BKW (0.6879) (0.0039) (0.2069) (0.4549) (0.006) (2.4559) Huang and Oluyede Journal of Statstcal Dstrbutons and Applcatons :8 Page 14 of 20
15 Table 6 Estmaton of models for poverty rate data Estmates Statstcs Dstrbuton α λ δ φ θ -2 log lkelhood AIC AICC BIC SS EKD ( ) (0.32) (0.0714) ( ) ( ) KD ( ) (0.0963) (0.0555) ( ) - D (0.2034) (105.94) (0.0784) - - a b c λ θ EKW (0.0944) (1.8026) (2.2276) (0.0199) ( ) a b α β c λ BKW (0.0069) (0.0431) (0.0584) (0.0017) (0.2564) (0.0075) Huang and Oluyede Journal of Statstcal Dstrbutons and Applcatons :8 Page 15 of 20
16 Huang and Oluyede Journal of Statstcal Dstrbutons and Applcatons :8 Page 16 of 20 Fgure 4 Ftted PDF for data sets. was calculated for each plot. The plot wth the smallest SS corresponds to the model wth ponts that are closer to the dagonal lne. Plots of the emprcal and estmated survval functons for the models are also presented n Fgure 6. For the ar condtonng system data ntal values α = 1 λ = 2 δ = 0.6 φ = 3 θ = 1 are used n SAS code for EKD model. The LR statstcs for the test of the hypothess H 0 : KD aganst H a : EKD and H 0 : D aganst H a : EKD are 1.9 (p-value = 0.17) and 13.4 (p-value = ). Consequently KD dstrbuton s the best dstrbuton based on the LR statstc. The KD dstrbuton gves smaller SS value than Dagum dstrbuton and slghtly bgger than EKD. For the non nested models the values of AIC and AICC for KD and EKW models are very close however the BIC value for KD dstrbuton s slghtly smaller than the correspondng value for the EKW dstrbuton. We conclude that KD model compares favorably wth the EKW dstrbuton and thus provdes a good ft for the ar condtonng system data. For the baseball player salary data set ntal values for EKD model n SAS code are α = 70 λ = 0.01 δ = φ = 0.1 θ = 1. The EKD dstrbuton s a better ft than KD and Dagum dstrbutons for ths data as well as the other dstrbutons. The values of the statstcs AIC AICC and BIC for KD dstrbuton are smaller compared to the non nested dstrbutons. The LR statstcs for the test of the hypotheses H 0 : KD aganst H a : EKD and H 0 : D aganst H a : EKD are 93.1 (p-value< ) and 361.5
17 Huang and Oluyede Journal of Statstcal Dstrbutons and Applcatons :8 Page 17 of 20 Fgure 5 Observed probablty vs predcted probablty for data sets. (p-value < ). Consequently we reject the null hypothess n favor of the EKD dstrbuton and conclude that the EKD dstrbuton s sgnfcantly better than the KD and Dagum dstrbutons based on the LR statstc. The value of AIC AICC and BIC statstcs are lower for the EKD dstrbuton when compared to those for the EKW and BKW dstrbutons. For poverty rate data ntal values for EKD model are α = 73 λ = 0.1 δ = 0.15 φ = 60 θ = The LR statstc for the test of the hypotheses H 0 : KD aganst H a : EKD and H 0 : D aganst H a : EKD are 8.2 (p-value = ) and 81.1 (p-value < ) respectvely. The values of AIC AICC and BIC statstcs shows EKD dstrbutons s a better model and the SS value of EKD model s comparatvely smaller than the correspondng values for the KD and D dstrbutons. Consequently we conclude that EKD dstrbuton sthebestftforthepovertyratedata. 9 Conclusons We have proposed and presented results on a new class of dstrbutons called the EKD dstrbuton. Ths class of dstrbutons have applcatons n ncome and lfetme data analyss. Propertes of ths class of dstrbutons ncludng the seres expanson of pdfs cdfs moments hazard functon reverse hazard functon ncome nequalty measures such as Lorenz and Bonferron curves are derved. Reny entropy order statstcs relablty mean
18 Huang and Oluyede Journal of Statstcal Dstrbutons and Applcatons :8 Page 18 of 20 Fgure 6 Emprcal survval functon for data sets. and medan devatons are presented. Estmaton of the parameters of the models and applcatons are also gven. Future work nclude MCMC methods wth censored data and regresson problems wth concomtant nformaton. Appendx RCodes nstall.packages( stats4 ) nstall.packages( bbmle ) lbrary(stats4) lbrary(bbmle) lbrary(stats) # Defne Functon #defneekdpdf g=functon(alphalambdadeltaphthetax){ y=alpha*lambda*delta*ph*theta*((x)ˆ(-delta-1))*((1+lambda*(xˆ(-delta)))ˆ(-alpha-1))*((1-((1+lambda*(xˆ(-delta)))ˆ(- alpha)))ˆ(ph-1))*((1-((1-((1+lambda*(xˆ(-delta)))ˆ(-alpha)))ˆ(ph)))ˆ(theta-1)) return(y) } #defneekdcdf
19 Huang and Oluyede Journal of Statstcal Dstrbutons and Applcatons :8 Page 19 of 20 G=functon(alphalambdadeltaphthetax){ y=(1-((1-((1+lambda*(xˆ(-delta)))ˆ(-alpha)))ˆ(ph)))ˆ(theta) return(y) } #defneekdhazard h=functon(alphalambdadeltaphthetax){ y=g(alphalambdadeltaphthetax)/(1-g(alphalambdadeltaphthetax)) return(y) } #defneekdquantle quantle=functon(alphalambdadeltaphthetaq){ ((lambda)ˆ(1/delta))*((((1-((1-((q)ˆ(1/theta)))ˆ(1/ph)))ˆ(-1/alpha))-1)ˆ(-1/delta)) } # defne EKD moments. #note:k<delta moments=functon(alphalambdadeltaphthetak){ f=functon(alphalambdadeltaphthetakx){(xˆk)*(g(alphalambdadeltaphthetax))} y=ntegrate(flower=0upper=infsubdvsons=100000alpha=alphalambda=lambdadelta=deltaph=phtheta=thetak=k) return(y) } #defneekdi(a) Ia=functon(alphalambdadeltaphthetaa){ n=length(a) y=0 for( n 1:n){ y[]=ntegrate(functon(alphalambdadeltaphthetax){x*g(alphalambdadeltaphthetax)}lower=0upper=a[] subdvsons=100000alpha=alphalambda=lambdadelta=deltaph=phtheta=theta)$value } return(y) } # defne EKD bonferron # note: p s between (01) bonferron=functon(alphalambdadeltaphthetap){ q=quantle(alphalambdadeltaphthetap) mu=moments(alphalambdadeltaphtheta1)$value y=(ia(alphalambdadeltaphthetaq))/(p*mu) return(y) } #defneekdlorenz # note: p s between (01) lorenz=functon(alphalambdadeltaphthetap){ q=quantle(alphalambdadeltaphthetap) mu=moments(alphalambdadeltaphtheta1)$value y=(ia(alphalambdadeltaphthetaq))/(mu) return(y) }
20 Huang and Oluyede Journal of Statstcal Dstrbutons and Applcatons :8 Page 20 of 20 Competng nterests The authors declare that they have no competng nterests. Authors contrbutons SH provded the R codes conducted the smulatons and the applcatons. BOO proposed developed the statstcal propertes of the new famly of dstrbutons and drafted the manuscrpt. All authors read and approved the fnal manuscrpt. Authors nformaton Shujao Huang s a graduate student at Georga Southern Unversty and Broderck O. Oluyede s Professor of Mathematcs and Statstcs at Georga Southern Unversty. Receved: 22 November 2013 Accepted: 11 March 2014 Publshed: 16 June 2014 References Chambers J Cleveland W Klener B Tukey P: Graphcal Methods for Data Analyss. Chapman and Hall London (1983) Cordero GM Ortega EMM Nadarajah S: The Kumaraswamy Webull dstrbuton wth applcaton to falure data. J. Frankln Inst (2010) Cordero GM Gomes AE de-slva CQ Ortega EMM: The Beta Exponentated Webull Dstrbuton. J. Stat. Comput. Smulat. 83(1) (2013) Dagum CA: New model of personal ncome dstrbuton: specfcaton and estmaton. Econome Applque e (1977) Domma F Condno F: The Beta-Dagum dstrbuton: defnton and propertes. Communcatons n Statstcs-Theory and Methods. 44(22) (2013) Domma F Gordano S Zenga M: Maxmum lkelhood estmaton n Dagum dstrbuton wth censored samples. J. Appl. Stat. 38(21) (2011) Kleber C: A Gude to the Dagum Dstrbutons. In: Duangkamon C (ed.) Modelng Income Dstrbutons and Lorenz Curve Seres: Economc Studes n Inequalty Socal Excluson and Well-Beng 5. Sprnger New York (2008) Kleber C Kotz S: Statstcal sze dstrbutons n economcs and actuaral scences. Wley New York (2003) Kumaraswamy P: Generalzed probablty densty functon for double-bounded random process. J. Hydrol (1980) McDonald B: Some generalzed functons for the sze dstrbuton of ncome. Econometrca. 52(3) (1984) McDonald B Xu J: A generalzaton of the beta dstrbuton wth applcaton. J. Econometrcs. 69(2) (1995) Oluyede BO Rajasoorya S: The Mc-Dagum Dstrbuton and Its Statstcal Propertes wth Applcatons. Asan J. Mathematcs and Applcatons. 2013(85) (2013) Proschan F: Theoretcal explanaton of observed decreasng falure rate. Technometrcs (1963) Reny A: On measures of entropy and nformaton. Berkeley Symp. Math. Stat. Probablty. 1(1) (1960) do: / Cte ths artcle as: Huang and Oluyede: Exponentated Kumaraswamy-Dagum dstrbuton wth applcatons to ncome and lfetme data. Journal of Statstcal Dstrbutons and Applcatons :8. Submt your manuscrpt to a journal and beneft from: 7 Convenent onlne submsson 7 Rgorous peer revew 7 Immedate publcaton on acceptance 7 Open access: artcles freely avalable onlne 7 Hgh vsblty wthn the feld 7 Retanng the copyrght to your artcle Submt your next manuscrpt at 7 sprngeropen.com
Marshall-Olkin Log-Logistic Extended Weibull Distribution : Theory, Properties and Applications
Journal of Data Scence 15(217), 691-722 Marshall-Olkn Log-Logstc Extended Webull Dstrbuton : Theory, Propertes and Applcatons Lornah Lepetu 1, Broderck O. Oluyede 2, Bokanyo Makubate 3, Susan Foya 4 and
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 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 informationStatistical Properties of the Mc-Dagum and Related Distributions
Georga Southern Unversty Dgtal Commons@Georga Southern Electronc Theses and Dssertatons Graduate Studes, Jack N. Avertt College of Sprng 2013 Statstcal Propertes of the Mc-Dagum and Related Dstrbutons
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 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 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 informationMaximum 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 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 informationThe Order Relation and Trace Inequalities for. Hermitian Operators
Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence
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 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 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 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 information/ n ) are compared. The logic is: if the two
STAT C141, Sprng 2005 Lecture 13 Two sample tests One sample tests: examples of goodness of ft tests, where we are testng whether our data supports predctons. Two sample tests: called as tests of ndependence
More informationMaximum 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 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 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 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 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 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 informationLecture 6 More on Complete Randomized Block Design (RBD)
Lecture 6 More on Complete Randomzed Block Desgn (RBD) Multple test Multple test The multple comparsons or multple testng problem occurs when one consders a set of statstcal nferences smultaneously. For
More informationASYMPTOTIC PROPERTIES OF ESTIMATES FOR THE PARAMETERS IN THE LOGISTIC REGRESSION MODEL
Asymptotc Asan-Afrcan Propertes Journal of Estmates Economcs for and the Econometrcs, Parameters n Vol. the Logstc, No., Regresson 20: 65-74 Model 65 ASYMPTOTIC PROPERTIES OF ESTIMATES FOR THE PARAMETERS
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 informationPredictive Analytics : QM901.1x Prof U Dinesh Kumar, IIMB. All Rights Reserved, Indian Institute of Management Bangalore
Sesson Outlne Introducton to classfcaton problems and dscrete choce models. Introducton to Logstcs Regresson. Logstc functon and Logt functon. Maxmum Lkelhood Estmator (MLE) for estmaton of LR parameters.
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 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 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 informationWhich estimator of the dispersion parameter for the Gamma family generalized linear models is to be chosen?
STATISTICS Dalarna Unversty D-level Master s Thess 007 Whch estmator of the dsperson parameter for the Gamma famly generalzed lnear models s to be chosen? Submtted by: Juan Du Regstraton Number: 8096-T084
More informationStatistics II Final Exam 26/6/18
Statstcs II Fnal Exam 26/6/18 Academc Year 2017/18 Solutons Exam duraton: 2 h 30 mn 1. (3 ponts) A town hall s conductng a study to determne the amount of leftover food produced by the restaurants n the
More informationDagum Distribution: Properties and Different Methods of Estimation
Internatonal Journal of Statstcs and Probablty; Vol. 6, No. 2; March 27 ISSN 927-732 E-ISSN 927-74 Publshed by Canadan Center of Scence and Educaton Dagum Dstrbuton: Propertes and Dfferent Methods of Estmaton
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 information2016 Wiley. Study Session 2: Ethical and Professional Standards Application
6 Wley Study Sesson : Ethcal and Professonal Standards Applcaton LESSON : CORRECTION ANALYSIS Readng 9: Correlaton and Regresson LOS 9a: Calculate and nterpret a sample covarance and a sample correlaton
More informationDepartment of Statistics University of Toronto STA305H1S / 1004 HS Design and Analysis of Experiments Term Test - Winter Solution
Department of Statstcs Unversty of Toronto STA35HS / HS Desgn and Analyss of Experments Term Test - Wnter - Soluton February, Last Name: Frst Name: Student Number: Instructons: Tme: hours. Ads: a non-programmable
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 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 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 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 informationThe binomial transforms of the generalized (s, t )-Jacobsthal matrix sequence
Int. J. Adv. Appl. Math. and Mech. 6(3 (2019 14 20 (ISSN: 2347-2529 Journal homepage: www.jaamm.com IJAAMM Internatonal Journal of Advances n Appled Mathematcs and Mechancs The bnomal transforms of the
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 informationLinear Regression Analysis: Terminology and Notation
ECON 35* -- Secton : Basc Concepts of Regresson Analyss (Page ) Lnear Regresson Analyss: Termnology and Notaton Consder the generc verson of the smple (two-varable) lnear regresson model. It s represented
More informationOn the spectral norm of r-circulant matrices with the Pell and Pell-Lucas numbers
Türkmen and Gökbaş Journal of Inequaltes and Applcatons (06) 06:65 DOI 086/s3660-06-0997-0 R E S E A R C H Open Access On the spectral norm of r-crculant matrces wth the Pell and Pell-Lucas numbers Ramazan
More informationStatistical analysis using matlab. HY 439 Presented by: George Fortetsanakis
Statstcal analyss usng matlab HY 439 Presented by: George Fortetsanaks Roadmap Probablty dstrbutons Statstcal estmaton Fttng data to probablty dstrbutons Contnuous dstrbutons Contnuous random varable X
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 informationNew Extended Weibull Distribution
Crculaton n Computer Scence Vol No6 pp: (4-9) July 7 https:doorg6ccs-7-5- New Extended Webull Dstrbuton Zubar Ahmad Research Scholar: Department of Statstcs Quad--Aam Unversty 45 Islamabad 44 Pastan Zawar
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 informationThe Jacobsthal and Jacobsthal-Lucas Numbers via Square Roots of Matrices
Internatonal Mathematcal Forum, Vol 11, 2016, no 11, 513-520 HIKARI Ltd, wwwm-hkarcom http://dxdoorg/1012988/mf20166442 The Jacobsthal and Jacobsthal-Lucas Numbers va Square Roots of Matrces Saadet Arslan
More informationMaximum Likelihood Estimation of Lorenz Curves using Alternative Parametric Model
Metodološk zvezk, Vol. 1, No. 1, 2004, 109-118 Maxmum Lkelhood Estmaton of Lorenz Curves usng Alternatve Parametrc Model Ibrahm M. Abdalla 1 and Mohamed Y. Hassan 2 Abstract In ths paper the Lorenz curve
More informationDO NOT OPEN THE QUESTION PAPER UNTIL INSTRUCTED TO DO SO BY THE CHIEF INVIGILATOR. Introductory Econometrics 1 hour 30 minutes
25/6 Canddates Only January Examnatons 26 Student Number: Desk Number:...... DO NOT OPEN THE QUESTION PAPER UNTIL INSTRUCTED TO DO SO BY THE CHIEF INVIGILATOR Department Module Code Module Ttle Exam Duraton
More informationThe Exponentiated-Log-logistic Geometric Distribution: Dual Activation
The Exponentated-Log-logstc Geometrc Dstrbuton: Dual Actvaton Natale V. R. Mendoza a Edwn M. M. Ortega b and Gauss M. Cordero c a Departamento de Cêncas Exatas Unversdade de São Paulo Av. Pádua Das 11
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 Generalization of the Weibull Distribution with Applications
Journal of Modern Appled Statstcal Methods Volume 5 Issue 2 Artcle 47 --26 A Generalzaton of the Webull Dstrbuton wth Applcatons Maalee Almhedat Unversty of Petra, Amman, Jordan, malmhedat@uop.edu.jo Carl
More informationSTAT 405 BIOSTATISTICS (Fall 2016) Handout 15 Introduction to Logistic Regression
STAT 45 BIOSTATISTICS (Fall 26) Handout 5 Introducton to Logstc Regresson Ths handout covers materal found n Secton 3.7 of your text. You may also want to revew regresson technques n Chapter. In ths handout,
More informationThe RS Generalized Lambda Distribution Based Calibration Model
Internatonal Journal of Statstcs and Probablty; Vol. 2, No. 1; 2013 ISSN 1927-7032 E-ISSN 1927-7040 Publshed by Canadan Center of Scence and Educaton The RS Generalzed Lambda Dstrbuton Based Calbraton
More informationa. (All your answers should be in the letter!
Econ 301 Blkent Unversty Taskn Econometrcs Department of Economcs Md Term Exam I November 8, 015 Name For each hypothess testng n the exam complete the followng steps: Indcate the test statstc, ts crtcal
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 informationAsymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation
Nonl. Analyss and Dfferental Equatons, ol., 4, no., 5 - HIKARI Ltd, www.m-har.com http://dx.do.org/.988/nade.4.456 Asymptotcs of the Soluton of a Boundary alue Problem for One-Characterstc Dfferental Equaton
More informationStatistics Chapter 4
Statstcs Chapter 4 "There are three knds of les: les, damned les, and statstcs." Benjamn Dsrael, 1895 (Brtsh statesman) Gaussan Dstrbuton, 4-1 If a measurement s repeated many tmes a statstcal treatment
More informationStatistics for Economics & Business
Statstcs for Economcs & Busness Smple Lnear Regresson Learnng Objectves In ths chapter, you learn: How to use regresson analyss to predct the value of a dependent varable based on an ndependent varable
More informationOn the Transmuted Additive Weibull Distribution
AUSTRIAN JOURNAL OF STATISTICS Volume 42 213, Number 2, 117 132 On the Transmuted Addtve Webull Dstrbuton Ibrahm Elbatal 1 and Gokarna Aryal 2 1 Mathematcal Statstcs, Caro Unversty, Egypt 2 Mathematcs,
More informationCS-433: Simulation and Modeling Modeling and Probability Review
CS-433: Smulaton and Modelng Modelng and Probablty Revew Exercse 1. (Probablty of Smple Events) Exercse 1.1 The owner of a camera shop receves a shpment of fve cameras from a camera manufacturer. Unknown
More informationSee Book Chapter 11 2 nd Edition (Chapter 10 1 st Edition)
Count Data Models See Book Chapter 11 2 nd Edton (Chapter 10 1 st Edton) Count data consst of non-negatve nteger values Examples: number of drver route changes per week, the number of trp departure changes
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 informationSTATISTICS QUESTIONS. Step by Step Solutions.
STATISTICS QUESTIONS Step by Step Solutons www.mathcracker.com 9//016 Problem 1: A researcher s nterested n the effects of famly sze on delnquency for a group of offenders and examnes famles wth one to
More informationComparison of Regression Lines
STATGRAPHICS Rev. 9/13/2013 Comparson of Regresson Lnes Summary... 1 Data Input... 3 Analyss Summary... 4 Plot of Ftted Model... 6 Condtonal Sums of Squares... 6 Analyss Optons... 7 Forecasts... 8 Confdence
More informationA Comparative Study for Estimation Parameters in Panel Data Model
A Comparatve Study for Estmaton Parameters n Panel Data Model Ahmed H. Youssef and Mohamed R. Abonazel hs paper examnes the panel data models when the regresson coeffcents are fxed random and mxed and
More informationECONOMICS 351*-A Mid-Term Exam -- Fall Term 2000 Page 1 of 13 pages. QUEEN'S UNIVERSITY AT KINGSTON Department of Economics
ECOOMICS 35*-A Md-Term Exam -- Fall Term 000 Page of 3 pages QUEE'S UIVERSITY AT KIGSTO Department of Economcs ECOOMICS 35* - Secton A Introductory Econometrcs Fall Term 000 MID-TERM EAM ASWERS MG Abbott
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 informationStatistical Inference. 2.3 Summary Statistics Measures of Center and Spread. parameters ( population characteristics )
Ismor Fscher, 8//008 Stat 54 / -8.3 Summary Statstcs Measures of Center and Spread Dstrbuton of dscrete contnuous POPULATION Random Varable, numercal True center =??? True spread =???? parameters ( populaton
More informationStatistics for Business and Economics
Statstcs for Busness and Economcs Chapter 11 Smple Regresson Copyrght 010 Pearson Educaton, Inc. Publshng as Prentce Hall Ch. 11-1 11.1 Overvew of Lnear Models n An equaton can be ft to show the best lnear
More informationSTAT 511 FINAL EXAM NAME Spring 2001
STAT 5 FINAL EXAM NAME Sprng Instructons: Ths s a closed book exam. No notes or books are allowed. ou may use a calculator but you are not allowed to store notes or formulas n the calculator. Please wrte
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 informationChapter 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 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 informationStat 642, Lecture notes for 01/27/ d i = 1 t. n i t nj. n j
Stat 642, Lecture notes for 01/27/05 18 Rate Standardzaton Contnued: Note that f T n t where T s the cumulatve follow-up tme and n s the number of subjects at rsk at the mdpont or nterval, and d s the
More informationLecture 6: Introduction to Linear Regression
Lecture 6: Introducton to Lnear Regresson An Manchakul amancha@jhsph.edu 24 Aprl 27 Lnear regresson: man dea Lnear regresson can be used to study an outcome as a lnear functon of a predctor Example: 6
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 informationChapter 8 Indicator Variables
Chapter 8 Indcator Varables In general, e explanatory varables n any regresson analyss are assumed to be quanttatve n nature. For example, e varables lke temperature, dstance, age etc. are quanttatve n
More informationISQS 6348 Final Open notes, no books. Points out of 100 in parentheses. Y 1 ε 2
ISQS 6348 Fnal Open notes, no books. Ponts out of 100 n parentheses. 1. The followng path dagram s gven: ε 1 Y 1 ε F Y 1.A. (10) Wrte down the usual model and assumptons that are mpled by ths dagram. Soluton:
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 informationModeling and Simulation NETW 707
Modelng and Smulaton NETW 707 Lecture 5 Tests for Random Numbers Course Instructor: Dr.-Ing. Magge Mashaly magge.ezzat@guc.edu.eg C3.220 1 Propertes of Random Numbers Random Number Generators (RNGs) must
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 informationGoodness of fit and Wilks theorem
DRAFT 0.0 Glen Cowan 3 June, 2013 Goodness of ft and Wlks theorem Suppose we model data y wth a lkelhood L(µ) that depends on a set of N parameters µ = (µ 1,...,µ N ). Defne the statstc t µ ln L(µ) L(ˆµ),
More informationChapter 14 Simple Linear Regression
Chapter 4 Smple Lnear Regresson Chapter 4 - Smple Lnear Regresson Manageral decsons often are based on the relatonshp between two or more varables. Regresson analss can be used to develop an equaton showng
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 informationTransfer 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 informationNegative Binomial Regression
STATGRAPHICS Rev. 9/16/2013 Negatve Bnomal Regresson Summary... 1 Data Input... 3 Statstcal Model... 3 Analyss Summary... 4 Analyss Optons... 7 Plot of Ftted Model... 8 Observed Versus Predcted... 10 Predctons...
More informationDepartment of Quantitative Methods & Information Systems. Time Series and Their Components QMIS 320. Chapter 6
Department of Quanttatve Methods & Informaton Systems Tme Seres and Ther Components QMIS 30 Chapter 6 Fall 00 Dr. Mohammad Zanal These sldes were modfed from ther orgnal source for educatonal purpose only.
More informationChapter 3 Describing Data Using Numerical Measures
Chapter 3 Student Lecture Notes 3-1 Chapter 3 Descrbng Data Usng Numercal Measures Fall 2006 Fundamentals of Busness Statstcs 1 Chapter Goals To establsh the usefulness of summary measures of data. 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 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 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 informationResearch Article Green s Theorem for Sign Data
Internatonal Scholarly Research Network ISRN Appled Mathematcs Volume 2012, Artcle ID 539359, 10 pages do:10.5402/2012/539359 Research Artcle Green s Theorem for Sgn Data Lous M. Houston The Unversty of
More informationA Solution of the Harry-Dym Equation Using Lattice-Boltzmannn and a Solitary Wave Methods
Appled Mathematcal Scences, Vol. 11, 2017, no. 52, 2579-2586 HIKARI Ltd, www.m-hkar.com https://do.org/10.12988/ams.2017.79280 A Soluton of the Harry-Dym Equaton Usng Lattce-Boltzmannn and a Soltary Wave
More informationBasic Statistical Analysis and Yield Calculations
October 17, 007 Basc Statstcal Analyss and Yeld Calculatons Dr. José Ernesto Rayas Sánchez 1 Outlne Sources of desgn-performance uncertanty Desgn and development processes Desgn for manufacturablty A general
More informationLINEAR REGRESSION ANALYSIS. MODULE VIII Lecture Indicator Variables
LINEAR REGRESSION ANALYSIS MODULE VIII Lecture - 7 Indcator Varables Dr. Shalabh Department of Maematcs and Statstcs Indan Insttute of Technology Kanpur Indcator varables versus quanttatve explanatory
More informationMore metrics on cartesian products
More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of
More informationCorrelation and Regression. Correlation 9.1. Correlation. Chapter 9
Chapter 9 Correlaton and Regresson 9. Correlaton Correlaton A correlaton s a relatonshp between two varables. The data can be represented b the ordered pars (, ) where s the ndependent (or eplanator) varable,
More information