A Generalized Rayleigh Distribution and Its Application
Total Page:16
File Type:pdf, Size:1020Kb
Biometrics & Biostatistics International Journal Research Article Open Access A generalized Rayleigh distribution and its application Abstract Volume 8 Issue 4 - 2019 This paper studies properties and applications of the generalized Rayleigh-truncated 1 2 1 Lishamol Tomy, Jiju Gillariose negative binomial distribution established by Jiju & Lishamol. For the concerned objective, 1Department of Statistics, Deva Matha College, India we applied the model to a real life data set and its performance is compared with that of 2Department of Statistics, St.Thomas College, India other four-parameter generalized Rayleigh distributions which are derived using different generators. Correspondence: Lishamol Tomy, Department of Statistics, Deva Matha College, Kuravilangad, Kerala-686633, India, Keywords: hazard rate function, marshall-olkin family, maximum likelihood estimation, Email rayleigh-truncated negative binomial distribution Received: May 31, 2019 | Published: July 08, 2019 Introduction distributed random variables with survival function F(x) and N be a truncated negative binomial random variable, independent of X ’s, 2 i The Rayleigh distribution was introduced by Rayleigh and with parameters 0<α <1 and β >0, such that originally proposed in the fields of acoustics and optics. It has emerged β β +− as a special case of the Weibull distribution. The Rayleigh distribution α n1 n P(N = n) = β (1−α ) ;n = 1,2,... (3) has widely used in communication theory to describe hourly median 1−α β −1 and instantaneous peak power of received radio signals. Moreover, It If U = min(X ,X ,...X ), then the survival function of U is has received a considerable attention from engineers and physicists N 12 N N for modeling wave propagation, radiation, synthetic aperture radar α β αβ +α −β −∈αβ images, and other related phenomena. There have been many forms G(x, , ) =β [F(x) F(x)] 1 x R, , > 0 (4) 1−α { } for the Rayleigh distribution to provide flexibility for modeling data. Vod 3,4 proposed a powerful extension of the Rayleigh distribution and Similarly, if α >1 and N is a truncated negative binomial 1 studied its properties. Its probability density function (pdf) is given by random variable with parameters and β >0, then α −−2λλ 2 2 +− 1 (x/ θ )2 2θ xe WN = max(X 12 ,X ,...X N) also has the same survival function f (x,θλ , ) = , x,θ > 0, λ >− 1. (1) Γ+(λ 1) given in (4). If α →1 in (4), then G(x,αβ , )→ F(x) . If β = 1 , where θ is the scale parameter, λ is the shape parameter and then this family reduces to the Marshall-Olkin family of distributions. ∞ Γ(a) = ta1−− e t dt is the complete gamma function. Its survival The pdf of survival function given in equation (4) is ∫0 (1−α ) βα β f (x) function is αβ 2 g(x, , ) = β +1 (5) Γ+(λθ 1,(x / ) ) β F(x,θλ , ) = 1−+−=γλ ( 1,(x / θ )2 ), x, θ > 0, λ > 1. (1−+αα ){ F(x) F(x)} Γ+(λ 1) Recently, several authors have used this approach to introduce new (2) distributions. Jayakumar & Sankaran10 defined a generalized uniform x 8 11 where Γ(a,x) = ta1−− e t dt is the incomplete gamma function, distribution using the approach of Nadarajah et al., Babu introduced ∫0 Weibull truncated negative binomial distribution. Further, Jayakumar hence statistical software’s can be used for various values of θ and & Sankaran12 introduced generalized exponential truncated negative λ . In literature, there are many studies based on extensions of this binomial distribution and studied its properties. Also, Jose & Sivadas13 generalized Rayleigh (GR) distribution using different generators. used the family given by equation (4) to introduce the negative Cordeiro et al.,5 derived four-parameter beta-GR distribution, Gomes binomial Marshall-Olkin Rayleigh distribution. et al.,6 proposed the four-parameter Kumaraswamy-GR distribution, and MirMostafaee et al.,7 introduced Marshall-Olkin extended GR The contents of this paper are organized as follows. Section distribution respectively. 2 deals with a generalized Rayleigh-truncated negative binomial distribution and its properties. Section 3 gives a real-life application. In modern era, the literature has suggested several ways of The concluding remarks are given in Section 4. extending well-known distributions to generate a more flexible of distributions. Recently new generator approach introduced by Generalized Rayleigh-truncated negative Nadarajah et al.,8 as a generalization of the family of Marshall-Olkin extended (MOE) distributions by Marshall & Olkin.9 This approach binomial distribution deals with the shape parameter induction in parent (or baseline) In this section, we are focused on generalized Rayleigh truncated distribution to explore tail properties and to improve goodness-of- negative binomial (GR-TNB) distribution introduced by Jiju & 1 fits. Let X12 ,X ,... be a sequence of independent and identically Lishamol. The distribution is derived by using generator approach of Submit Manuscript | http://medcraveonline.com Biom Biostat Int J. 2019;8(4):139‒143. 139 © 2019 Lishamol et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and build upon your work non-commercially. Copyright: A generalized Rayleigh distribution and its application ©2019 Lishamol et al. 140 Nadarajah et al.,8 They have examined various statistical properties function of X follows GR-TNB distribution, it can be expressed as of this distribution including estimation of parameters and have 1 showed that this distribution is more flexible comparing to other − 1 (1−− u)(1α β ) β generalizations of the rayleigh distribution. The survival function of Q(u) =θα QG β +−1 the GR-TNB distribution is given by 1−α α β α 22−β G(x,αβλθ , , , ) =β [ γλ (+ 1,(x / θ ) ) ++ αγλ ( 1,(x / θ ) )] − 1 , x,αβθ , , > 0, λ >− 1 ⋅ −α { } where u is generated from the uniform (0, 1) distribution and QG () 1 is the (standardized) gamma quantile function with shape parameter (6) λ +1 and unit scale parameter (Figure 1). corresponding pdf is Simulation study 2 (1−α ) βαβ 2 θ −2( λλ + 1) x 2 +− 1 e (x/ θ ) In this section, we carry out Monte Carlo simulation study to αβλθ αβθ λ − g(x, , , , ) = β +1 , x, , , > 0, > 1 β 22 assess the performance of the maximum likelihood estimates (MLE). Γ+(λ 1)(1 − α ){ γλ ( + 1,(x / θ ) ) + αγλ ( + 1,(x / θ ) )} The results are obtained from generating 1000 samples from the GR- TNB distribution. For each replication, a random sample of size n = (7) 50, 100 and 200 is drawn from the GR-TNB distribution. The GR- TNB random number generation was performed using the quantile and the hazard rate function (hrf) of the GR-TNB distribution function of GR-TNB distribution and the parameters are estimated by becomes using the method of MLE by using package nlm in R, we get MLEs ˆ ˆ ˆ 2 xx22β , θ and λ for fixed α = 0.2 or α = 0.02 . The evaluation of the (1−αβθ ) 2−−2λλ 2 x 2 +− 1 e (x/ θ ) [γλ (+++ 1, ) αγλ ( 1, )]−1 22performance is based on the bias and the mean squared errors (MSE) h(x,α , β , λθ , ) = θθ Γ+(λ 1)(1 − [ γλ ( + 1,(x / θ )22 ) + αγλ ( + 1,(x / θ ) )]β ) defined as follows: Average bias of the simulated N estimates of R : (8) N 1 ˆ The GR-TNB distribution enfolds some sub models such as MOE ∑(Ri − R) GR distribution, MOE half normal distribution and MOE Reyleigh N i=1 distribution. Figure 1 shows the plots for pdf and hrf for GR-TNB Average Mean square error of the simulated N estimates of R : distribution with various parameter values. As seen in Figure 1, the pdf N and the hrf of the GR-TNB distribution have several different shapes 1 ˆ 2 ∑(Ri − R) according to the values of the parameters. This shows that GR-TNB N i=1 distribution is more flexible than Rayleigh distribution. The quantile Figure 1 Graphs of pdf and hrf the GR-TNB distribution for different values of α , β , λ and θ . where R is the true value of parameters β ,θ and λ and also Data analysis N is the number of replications. The initial values of parameter are In this section, we consider a real data set on breaking stress β = 1.5 ,θ = 0.1 and λ = 0.2 . The results of our simulation study are of carbon fibres of 50 mm length (GPa) to assess the flexibility of summarized. From this table, we can see that the bias and MSE of the the GR-TNB distribution over some well-known generalizations of MLEs converge to zero when the sample size is increased (Table 1). Rayleigh distribution. The data have been previously used by Nichols & Padgettcit.14 The data are as follows: Citation: Lishamol T, Jiju G. A generalized Rayleigh distribution and its application. Biom Biostat Int J. 2019;8(4):139‒143. DOI: 10.15406/bbij.2019.08.00282 Copyright: A generalized Rayleigh distribution and its application ©2019 Lishamol et al. 141 0.39, 0.85, 1.08, 1.25, 1.47, 1.57, 1.61, 1.61, 1.69, 1.80, 1.84 ,1.87, We compare the results of GR-TNB distribution with following 1.89, 2.03, 2.03, 2.05, 2.12, 2.35, 2.41, 2.43, 2.48, 2.50, 2.53, 2.55, four-parameter generalizations of Rayleigh distribution which are 2.55, 2.56, 2.59, 2.67, 2.73, 2.74, 2.79, 2.81, 2.82, 2.85, 2.87, 2.88, generalized by using different generators: 2.93, 2.95, 2.96, 2.97, 3.09, 3.11, 3.11, 3.15, 3.15, 3.19, 3.22, 3.22, 1.