On Size-Biased Double Weighted Exponential Distribution (SDWED)

This paper introduces a new distribution based on the exponential distribution, known as Size-biased Double Weighted Exponential Distribution (SDWED). Some characteristics of the new distribution are obtained. Plots for the cumulative distribution function, pdf and hazard function, tables with values of skewness and kurtosis are provided. As a motivation, the statistical application of the results to a problem of ball bearing data has been provided. It is observed that the new distribution is skewed to the right and bears most of the properties of skewed distribution. It is found that our newly proposed distribution fits better than size-biased Rayleigh and Maxwell distributions and many other distributions. Since many researchers have studied the procedure of the weighted distributions in the estates of forest, biomedicine and biostatistics etc., we hope in numerous fields of theoretical and applied sciences, the findings of this paper will be useful for the practitioners.

KEYWORDS

1. Introduction

Weighted distributions are suitable in the situation of unequal probability sampling, such as actuarial sciences, ecology, biomedicine biostatistics and survival data analysis. These distributions are applicable, when observations are recorded without any experiment, repetition and random process. The notion of weighted distributions has been used as a device for the collection of suitable model for observed data, during last 25 years. The idea is most applicable when sampling frame is not available and random sampling is not possible. Firstly the idea of weighted distributions was introduced by Fisher [1] . Cox [2] firstly provided the idea of length-biased sampling and after that Rao [3] established a unifying method that can be used for several sampling situations and can be displayed by means of the weighted distributions. Cox [4] estimated the mean of the original distribution built on length-biased data. Zelen [5] presented the concept of weighted distribution in studying cell kinetics and early discovery of disease. Warren [6] applied these distributions in forest product research. Patil and Rao [7] surveyed the idea and applications of weighted and size-biased sampling distributions. Patil and Rao [8] also discussed weighted binomial distribution to model the human families and estimation of the wildlife family size. Gupta and Keating [9] described the relationship between reliability measures of original and size-biased distribution. Arnold and Nagaraja [10] gave the idea of bivariate weighted distribution whereas Jain and Nanda [11] extended this idea and discussed multivariate aspect of weighted distribution.

Let be the pdf of the random variable x and be the unknown parameter.

The weighted distribution is defined as;

(1)

where w(x) is a weight function. When, then these distributions are termed as size-biased distribution of order m. When it is called size-biased of order 1 or say length biased distribution, whereas for it is called the area-biased distribution (Ord and Patil [12] , Patil [13] and Mahfound [14] ).

In forest product research, equilibrium and length biased distributions have been used as moment distributions. Kochar and Gupta [15] discussed the moment distributional properties in assessment with the actual distributions and derived the bound on the moments of moment distributions.

Oluyede [16] described inequalities for the reliability measures of size-biased and the original distributions. Navarro et al. [17] discussed characterization of the original and the size-biased distribution using reliability measures. Gove [18] offered the uses of size-biased distributions in forest science and ecology. Sunoj and Maya [19] established relationships among weighted and original distributions in the situation of repairable system and also characterized the sized-biased and the original distribution. Shen et al. [20] used semi-parametric transformations to model the length biased data. Hussain and Ahmad [21] presented misclassification in the size-biased modified power series distributions and its applications.

Mir and Ahmad [22] derived generalized forms of size-biased discrete distributions and discussed the practical applications in the field of Medical, Zoology and Accidental studies. Mir [23] derived size-biased Geeta distribution and size-biased consul distribution respectively, different properties are discussed and contrasts with original distributions are also done. Das and Roy [24] established size-biased form of generalized Rayleigh distribution and apply the consequences to the environmental data. They also applied the concept of size-biased sampling in the field of environmental studies

Dara [25] derived reliability measures for size-biased forms of several moment distributions as the special cases of moment distributions. Iqbal and Ahmad [26] found compound scale mixtures of limiting distribution of generalized log Pearson type VII distribution with different continuous and moment distributions. Hasnain [27] introduced a new family of distributions named as exponentiated moment exponential (EME) distribution and developed its properties. Iqbal et al. [28] found a more general class for EME distribution and built up different properties including characterization through conditional moments.

Zahida and Munir [29] worked on Weighted Weibull Distributions (WWD), Double Weibull Distributions (DWD), Weighted Double Weibull Distributions (WDWD), Double Weighted Exponential Distributions (DWED) (both in size-biased and area biased). Some basic theoretical properties of all these distributions including cumulative density function, central moments, skewness, kurtosis and moments are studied. Shannon entropy, Renyi entropy, moment generating function and information generating function of all these distributions are derived. Reliability measures including survival function, failure rates, reverse hazard rate function and Mills ratios of these distributions are also obtained. Parameters are evaluated by using method of maximum likelihood estimation along with derivation of practical examples.

The exponential distribution has a fundamental role in describing a large class of phenomena, particularly in the area of reliability theory. This distribution is commonly used to model waiting times between occurrences of rare events, lifetimes of electrical or mechanical devices. It is also used to get approximate solutions to difficult distribution problems.

2. Methodology

2.1. Size-Biased Double Weighted Exponential Distribution (SDWED)

The size-biased double weighted exponential distribution is given by:

(2)

where f(x) is the first weight and

Here and, , is the pdf of exponential distribution.

Thus the pdf of SDWED is

(3)

where is shape parameter and c is scale parameter.

Graphs of Probability Density Function

Figure 1 and Figure 2 show the probability density function of SDWED.

2.2. Distribution Function of SDWED

Distribution function of a density function is defined as:

Figure 1. The probability density function of SDWED for the indicated values of c and λ.

Figure 2. The probability density function of SDWED for the indicated values of c and λ.

After some simplification we have: (4)Figure 3 shows the commutative distribution function of SDWED.2.3. Survival FunctionThe survival function of SDWED is defined as (5)Figure 4 shows the survival function of SDWED.2.4. Hazard Rate Function of SDWEDThe hazard rate function is defined as:At and, the hazard rate function will be:

Figure 3. Cumulative Distribution Function of SDWED for the indicated values of c and λ.

Figure 4. The Survival Function of SDWED for the indicated values of c and λ.

(6)

Figure 5 shows the Hazard Rate Function of SDWED.

2.5. Reverse Hazard Rate Function

The reverse Hazard rate function of SDWED is given by

At and, the reverse hazard rate function will be: (7)Figure 6 shows the Reverse Hazard Rate Function of SDWED.

Figure 5. The Hazard Function of SDWED for the indicated values of c and λ.

Figure 6. Reverse Hazard Function of SDWED for the indicated values of c and λ.

2.6. Mills Ratio

The Mills Ratio is given by:

At and, the Mills Ratio will be: (8)Figure 7 shows the mills ratio of SDWED.2.7. Moment Generating Function of SDWEDThe moment generating function of SDWED is: (9)Using Equation (3)Applying the transformations and after simplifying:(10)2.8. Information Generating Function of SDWEDThe information generating function is defined as:

Figure 7. Mills ratio of SDWED for the indicated values of c and λ.

Using Equation (3)

Putting

and after a long simplification, the information generating function will be:

(11)

2.9. Limit and Mode of SDWED

Note that the limit of the density function given in Equation (3) is as follows:

(12)

(13)

since

and (14)

2.10. Mode of SDWED

Taking log of Equation (3) on both sides:

(15)

Differentiating Equation (15) with respect to x, we obtain:

The mode of the SDWED is obtained by solving the nonlinear equation with respect to x:

(16)

The mode of SDWED is given in Table 1.

2.11. Mean of SDWED

(17)

Table 1. Mode of SDWED.

Table 2. Mean, variance and standard deviation of SDWED.

2.12. Variance of SDWED

(18)

Table 2 shows the Mean, Variance and Standard Deviation with some values of the parameters and c.

2.13. Moments of SDWED

The rth moment of SDWED is given by

for r = 1, 2, 3, 4, the first four moments about the mean are

2.14. Moment Ratios

Table 3 shows the coefficients of skewness and kurtosis.

3. Maximum Likelihood Estimation

The maximum likelihood estimation of SDWED distribution may be defined as:

Here the independent observations are, then the likelihood function of the DWED is:

(19)

(20)

and

(21)

Equating these equations to zero, then we get:

(22)

Table 3. Coefficients of skewness and kurtosis of SDWED.

(23)

which can be solved simultaneously for and.

The asymptotic variance-covariance matrix is the inverse of

The inverse of the asymptotic covariance matrix is with

(24)

(25)

(26)

4. Numerical Example

The Ball Bearing Data Records

See for data set published in Lawless [30]

In Table 4, the approximations of the parameters are specified. For goodness-of-fit statistics Anderson-Darling and Cramer-von Mises tests have been used, SDWED model proposals the best fitting (see Figure 8 and Figure 9):

5. Conclusion

In this paper, Size-biased Double Weighted Exponential Distribution (SDWED) has

Table 4. Parameters’ estimates and goodness-of-fit statistics.

Figure 8. Size-Biased Double Weighted Exponential (Dotted Line), weighted Weibull (Solid Line), Rayleigh (Dashes Line) and Maxwell (dotted dashed Line) on the Histogram.

Figure 9. SWWD (dotted Line), SDWED (dott dashed), Maxwell (Solid Line), Rayleigh (Dashed Line), Density Estimates?cdf Estimates and Empirical cdf.

been introduced. The pdf of the SDWED has been studied as well as different reliability measures such as survival function, failure rate function or hazard function. The moments, mode, the coeff. of skewness and the coeff. of kurtosis of SDWED have been derived. For estimating the parameters of SDWED, MLE method has been used. The SDWED has been fitted to Ball Bearing data set. SDWED suggested a good fit of the data as comparing to other distributions.

Conflicts of Interest

The authors declare no conflicts of interest.

Cite this paper

Perveen, Z. , Ahmed, Z. and Ahmad, M. (2016) On Size-Biased Double Weighted Exponential Distribution (SDWED). Open Journal of Statistics, 6, 917-930. doi: 10.4236/ojs.2016.65076.