International Journal of Physical Sciences

Transcription

1 Vol. 9(4), pp , 28 February, 2014 DOI: /IJPS ISSN Copyright 2014 Author(s) retain the copyright of this article International Journal of Physical Sciences Full Length Research Paper Lehmann Type II weighted Weibull distribution N. I. Badmus 1 *, T. A. Bamiduro 2 and S. G. Ogunobi 3 1 Department of Statistics, Abraham Adesanya Polytechnic, Ijebu-Igbo, Nigeria. 2 Department of Mathematical Sciences, Redeemer's University, Redemption Camp, Ogun State, Nigeria. 3 Department of Science Laboratory Tehnology, Abraham Adesanya Polytechnic, Ijebu-Igbo, Nigeria. Accepted 5 February, 2014 We upgrade on the method used on weighted Weibull model proposed by Azzalini (1985) using the logit of Beta function by Jones (2004) to have Lehmann Type II weighted Weibull model with a view to obtaining a distribution that is more better than both weighted Weibull and Weibull distribution in terms of estimate of their characteristics and their parameter. The weighted Weibull distribution is proposed by slightly modifying the method of Azzalini (1985) on weighted distribution with additional shape and scale parameters. Some basic properties of the newly proposed distribution including moments and moment generating function, survival rate function, hazard rate function, asymptotic behaviours, and the estimation of parameters have been studied. Key words: Asymptotic, hazard rate, Lehmann Type II, moments, weighted-weibull. INTRODUCTION The Weibull distribution is a well known common distribution and has been a powerful probability distribution in reliability analysis, while weighted distributions are used to adjust the probabilities of the events as observed and recorded. The Weibull distribution can also be used as an alternative to Gamma and Log-normal distribution in reliability engineering and life testing. Numerous authors/researchers have been worked on Weibull and eponential distribution in literature. In a nut shell, Shahbaz et al. (2010) applied Azzalini s method with the Weibull distribution that produced a new class of weighted Weibull distribution in which they obtained the probability density and cumulative density function as given below:, > 0 (1) The graph of the probability density function of Equation (1) with is given Figure 1. In the literature, researchers/authors on Lehmann Type II are very few. Figure 1 is the pdf graph of the weighted Weibull by Shabarz (2010). The aim of this article is to introduce this distribution and study on its statistical properties. This article presents the proposed distribution Lehmann Type II weighted Weibull distribution. Thereafter, moments and moment generating function is studied. This is followed by a critical discussion of the estimation of parameters showing the empirical distribution of data and the study was concluded. (2) *Corresponding author. idowuolasunkanmi8169@yahoo.com Author(s) agree that this article remain permanently open access under the terms of the Creative Commons Attribution License 4.0 International License

2 weipdf Int. J. Phys. Sci. PDF of Weighted Weibull a=4,ß=3 Figure 1. The probability density function of weighted Weibull distribution with function from Equation (1) is given in Equation (2). and the cumulative density MATERIALS AND METHODS The proposed Lehmann Type Ii weighted Weibull distribution The logit of beta function (the link function of the beta generalized distribution) is introduced by Jones (2004), since then etensive work has been done using the logit of beta distribution in literature. For instance, Gupta and Kundu (1990) proposed a generalized eponential distribution that provides an alternative to eponential and Weibull distribution. Lee et al. (2007), Eugene and Famoye (2002) and Nadarajah and Kotz (2005) also used the logit of beta distribution and they provided an etension on eponential distribution. The logit of beta distribution used by Famoye and Olugbenga (2005) to introduce the Beta-Weibull distribution alongside its major properties. The transformation is called the Lehmann Type I distribution. There is a dual transformation defined by referred to as the Lehmann Type II distribution introduced by Cordeiro et al. (2011) among others. Figures 2, 3 and 4 shows pdf, cdf and hazard rate graphs at several values of parameters. Figure 2 clearly shows that it is rightly skewed, Figure 3 rise upward to right and Figure 4 shows the level at which the hazard rate by rightly upward breaking points. Now, let be a random variable form of the distribution with parameters and defined (1) and (2) using the logit of beta by Jones (2004), we then have (3) Both probability density and cumulative density function (pdf and cdf) of Lehmann Type II weighted Weibull distribution is obtained by substituting pdf and cdf of weighted Weibull distribution in Equations (1) and (2) into Equation (3) and set we then obtain (4) Where,. Equation (5) is the pdf of Lehmann Type II weighted Weibull distribution. setting (5) substituting into Equation (4), we get (6) (6)

3 LWWpdf LWWpdf LWWpdf Badmus et al. 73 The PDF of Lehmann Type II Weighted Weibull a=1, b=5,a=0.5,ß=2 The PDF of Lehmann Type II Weighted Weibull a=1,b=3,a=3,ß=4.5 The PDF of Lehmann Type II Weighted Weibull a=1,b=3,a=4,ß=3 Figure 2. The pdf of LWW distribution at various values of the parameters and it is clear that indeed it is rightly skewed. Figure 2 shows the pdf of LWW distribution at various values of the p clear that

4 LWWcdf LWWcdf LWWcdf Int. J. Phys. Sci. The CDF of Lehmann Type II Weighted Weibull b=3.5 The CDF of Lehmann Type II Weighted Weibull a=1,b=3,a=4,ß=2 The CDF of Lehmann Type II Weighted Weibull a=1,b=3,a=1,ß=1 Figure 3. The CDF of LWW with a=1, b=3, give the same pattern and shape.. However, at given various values of the parameters, the graphs

5 LWWhaz haz 0e+00 2e+05 4e+05 6e+05 8e+05 Badmus et al. 75 Hazard Rate of LWW a=1, b=0.5,a=1,ß=2 The Hazard Rate of LWW Distribution (b=1.8) Figure 4. The hazard rate function of Lehmann Type II weighted Weibull distribution with several values of the parameter. Therefore, Equation (4) can be re-written as From logit of beta link function, we show that Equation (7) equal to one (7) Figure 2 show the graph of the Lehmann Type II weighted Weibull distribution at different values of the parameters. Cumulative density function (CDF) We have defined the probability density function (pdf) of in Equation (6), then Equation (6) can be given as Then, recall that (8) and using the differential equation method on have Therefore, we and cdf is given as (9)

6 76 Int. J. Phys. Sci. The graph of cumulative distribution of Lehmann Type II Weighted Weibull is shown Figure 3. The survival rate function of the distribution which is investigated as follows. Now from Epression (4), we have The survival rate function of the is defined by taking the limit (10) The hazard rate function The hazard rate function of a random variable with the pdf and cdf is given by This shows that the asymptotic behaviours when, and it is sure that the distribution has at least a mode. = 0 However, with and respectively defined in Equations (4) and (9), the hazard rate function can be epressed as: or (11) Where W is the epression in Equation (5) To show that we get the following where, = MOMENTS AND MOMENT GENERATING FUNCTION Shittu and Adepoju (2013) described that when a random variable following a generalized beta generated distribution that is, is a constant and where is the inverse of CDF of the weighted Weibull distribution, since distribution is a special form when a=c=1.we then derive the moment generating function (mgf) of the proposed distribution and the general rth moment of a beta generated distribution is defined by (12) Also, using the taylor series epansion around the point to obtain (13) For simplification, we take the limit of the following: When Cordeiro et al. (2011) described an alternative series epansion in their paper for in terms of where k follows the parent distribution then for z = 0, 1, = 0, Equation (11) above tends to and equal to zero (Figure 4). Asymptotic behaviours Following the steps in Hazard Rate Function above, by taking the They further discussed another mgf of for generated beta distribution Where, Therefore, (14) (15)

7 Badmus et al. 77 Substituting both pdf and cdf (F() and f()) of the weighted Weibull distribution into Equation (15), we obtain We can also get the standard deviation (SD) as (23) (16) Others like coefficient of variation, skewness and kurtosis can easily be obtained. Equation (17) becomes the mgf of Lehmann Type II weighted Weibull distribution and by setting b=1 and i = 0, the same epression (17) is reduced to become the parent distribution. To obtain the rth moment of, the weighted Weibull distribution by Shahbaz et al. (2010) is given by (17) ESTIMATION OF PARAMETER The maimum likelihood estimate (MLEs) of the parameter of distribution following Cordeiro et al (2011) given the log-likelihood function for, where by setting. Shittu and Adepoju (2013) in their study discussed that to be a vector of parameter and given by (24) Equation (16) can be re-written as (18) Recall that, a = c = 1 (24) becomes (25) Moreover, the rth moment of the proposed distribution can also be written from the above equation (18) as given below: where, and (19) Again, setting b = 1 in Equation (19) leads to the rth moment of the weighted Weibull distribution by Shahbaz et al. (2010) and is given by While the rth moment of the proposed distribution is For determining the MLE of Equation (26) with respect to ( (26), we took the partial derivative of ) as follows: (27) (20) +(b- From Equation (20), one can easily obtain the mean about the origin e.g when r = 1and the second moment when r = 2, etc. The first moment of is obtain 1) (28) (21) +(b-1) The second moment can also be obtained as follows: (22) (29) where, The above equations can be solved using Newton Ralphson method to obtain the the MLE of ( ) respectively. Taking second derivatives of Equations 27, 28 and 29 with respect to the parameters above we can derive the interval estimate and hypothesis tests on the model parameter and inverse of fisher s information matri needed.

8 78 Int. J. Phys. Sci. RESULTS AND DISCUSSION Application to life time data set Here we used a data set studied by Shahbaz et al (2010) on life components in years to compare the Lehmann Type II weighted Weibull and weighted Weibull distribution. The data contains Life (Grouped data say (0 1.0, ,,> 5.0) and frequency all together is R software is used to determine the maimum likelihood estimates and the log-likelihood for the Lehmann Type II weighted Weibull distribution are: while the maimum likelihood estimates and the loglikelihood for the weighted Weibull distribution are: and, where denote log-likelihood of both Lehmann Type II weighted Weibull distribution and weighted Weibull distribution. The asymptotic covariance matri of the maimum likelihood estimates for the Lehmann Type II weighted Weibull distribution, which is generated from the inverse of fisher s information matri and is given by REFERENCES Azzalini A (1985). A class of distributions which include the normal ones, Scand. J. Stat. 12: Shahbaz S, Shahbaz MQ, Butt NZ (2010). A class of weighted Weibull distribution in Pakistan J. Statistics Operation Res. VI(1): Cordeiro GM, Aleandra C, Ortega MM, Edwin Sarabia JM (2011). Generalized Beta Generated distributions. ICMA Centre. Discussion Papers in Finance DP pp Gupta RD, Kundu D (1990). The role of weighted distribution in stochastic modeling. Communications in statistics, 19: Eugene NC, Famoye F (2002). Beta-Normal distribution and its applications. Communications in Statistics Theory and Methods 31: Famoye F, Lee C, Olugbenga O (2005). The beta-weibull distribution. J. Statistical Theory Applications, 4(2): Jones MC (2004). Families of distributions arising from distributions of order statistics test 13:1-43. Lee C, Famoye F, Olumolade O (2007). Beta-Weibull distribution: Some properties and applications to censored data. J. modern Appl. statistical methods, 6: Nadarajah S, Kotz S (2005). The beta eponential distribution. Reliability Eng. Syst. Safety, 91: Shittu OI, Adepoju AK (2013). On the Beta-Nakagami Distribution. Progress in Appl. Maths. 5(1): Following the result above, this shows that the new proposed distribution can take care even more skew data than the weighted Weibull distribution because of the additional shape parameter b, and with the hope that the model will be wider application to many areas of research e.g economics, finance, environmental, biomedical among others. Conclusion In a nut shell, we studied the statistical properties of the proposed distribution e.g moments, moment generating function, estimation of parameters with R software are presented in this paper. We also upgraded with an additional parameter to the eisting two parameters in the weighted Weibull distribution and with the method of maimum likelihood the parameters of the model were estimated in order to give way for the derivation of fisher information matri. The data used indicates that Lehmann Type II weighted Weibull distribution has a better representation of data and more fleible than weighted Weibull distribution. Conflict of Interests The author(s) have not declared any conflict of interests.