The length-biased weighted exponentiated inverted Weibull distribution

Abstract: Length-Biased distributions are a special case of the more general form known as weighted distributions. We can exploit the conceptuality of Length-Biased distribution in the development of appropriate models for lifetime data. Its method is adjusting the original probability density function from real data and the expectation of those data. This modification can lead to correct conclusions of the models. Therefore, we introduced the Length-Biased version of the weighted Exponentiated inverted Weibull distribution in this paper. Various properties and the expressions for moments, coefficient of skewness, coefficient of kurtosis, moment generating function, hazard rate function, etc. are derived. The maximum likelihood estimates of the parameters of the proposed distribution are determined. The study results suggest that this distribution is an efficacious model in life time data analysis and other related fields.


PUBLIC INTEREST STATEMENT
The theory of weighted distributions provides a unifying approach for model specification and data interpretation problems. The current work is on length biased form of weighted exponential inverted distribution which is a better model for modeling of uncensored data. Weibull (1951), the Swedish Scientist interposed the Weibull distribution. It is the most widely used distribution for the analyzation of life time data. This distribution stipulates the immense encroachment of reliability and quality control such as ball bearing auto mobile components concrete bridges demography actuarial study and electrical insulation. It is also utilized in biological and medical applications. The weighted distributions speculation gives a compositional formulation for the model confinement and data illustration problems. Fisher (1934) and Rao (1965) first familiarized the idea of weighted distributions. In Rao's paper, he recognized several situations that can be modeled by weighted distributions. These situations refer to instances where the recoded observations cannot be considered as a random sample from the original distributions. This may occur due to non observability of some events or damage caused to the original observation resulting in a reduced value, or adoption of a sampling procedure which gives unequal chances to the units in the original.

Introduction
The weighted version of the bivariate logarithmic series distribution is presented by Gupta and Tripathi (1990) and a size biased sampling and related invariant weighted distributions are studied by Patil and Ord (1976). Zelen and Feinleib (1969) introduced the application examples for weighted distributions to represent what he broadly perceived as Length-Biased sampling. The Weibull distribution is a member of the family of extreme value distributions. These distributions are the limit distributions of the smallest or the greatest value, respectively, in a sample with sample size n → ∞ and is often used in the field of life data analysis due to its flexibility. The Weibull distribution includes the exponential and the Rayleigh distributions as sub models. The use of the distribution in reliability and quality control work was advocated by many authors following Kao (1959) and Berrettoni (1964). The usefulness and applications of parametric distributions including Weibull, Rayleigh are seen in various areas including reliability, renewal theory, and branching processes which can be seen in papers by several authors including Patil and Rao (1978), Gupta and Kirmani (1990), Gupta and Keating (1985), Oluyede (1999). Examples for weighted distribution see El-Shaarawi, Walter, and Patil (2002), Oluyede and George (2000). For more important results of weighted distribution you can see also Ghitany and Al-Mutairi (2008). Kersey (2010) introduced the weighted inverse Weibull distribution and beta-inverse Weibull distribution. Del Castillo and Pérez-Casany (1998) introduced new exponential families that come from the concept of weighted distribution, that include and generalize the Poisson distribution. Weighted distributions in general and Length-Biased distributions in particular are very useful and convenient for the analysis of lifetime data and the Length-Biased distribution representation was primitively familiar by Cox (1962). A table for some basic distributions and their Length-Biased forms is given such as Lognormal by Ratnaparkhi and Naik-Nimbalkar (2012), Nanuwong and Bodhisuwan (2014) presented the Length-Biased beta Pareto distribution. Khattree (1989) proposed a useful result by giving a relationship between the original random variable X and its Length-Biased form Y when X is either Inverse Gaussian or Gamma distribution. Das and Roy (2011a) developed the Length-Biased form of the Weighted Generalized Rayleigh distribution (WGRD) known as Length-Biased WGRD and also developed the Length-Biased form of weighted Weibull distribution by Das and Roy (2011b). Al-Khadim and Hussein (2014) proposed the length biased form of weighted Exponential and Rayleigh distribution. Seenoi, Supapa, and Bodhisuwan (2014) and Al-Khadim and Hussein (2014) developed the Length-Biased form of Exponentiated inverted Weibull distribution.
Considering the importance of weighted and Length-Biased distribution, we present the Length-Biased weighted Exponentiated Inverted Weibull Distribution (LBWEIWD) and the sub models which are the special cases of our proposed distribution. Various useful mathematical properties of the LBWEIWD are derived in the next sections. We also present a numerical example of the proposed distribution considering the real life data-set. This paper is organized as follows. Section 2 defines some basic materials and in Section 3, we provide the derivation of PDF and CDF of LBWEIWD. Some particular cases are obtained in Sections 4 and 5 defines the shape and mode of the LBWEIWD. Section 6 discusses the different statistical properties of our proposed model. Estimation of the unknown parameters of the proposed model is carried out in Section 7. The real data-set has been analyzed in Sections 8 and 9 gives some brief conclusion.

Materials and methods
Weighted distribution concept can be traced from the study of Fisher and Rao. If we have a non negative random variable X, with its PDF f(X) then the PDF of the weighted random variable X w is obtained as where w(x) is a non negative weight function and [w( When the weight function depends on the lengths of units of interest (i.e. w(x) = x) the resulting distribution is called Length-Biased. In this case the PDF of a Length-Biased random variable x is defined as Let T be a random variable with Weibull distribution its PDF takes the form where > 0 and β > 0 are the scale and shape parameters, respectively. The distribution includes especially the exponential and Rayleigh distributions as special cases where β = 1 and β = 2, respectively. The Weibull distribution has a distribution function of the form The survival function and hazard rate are, respectively, The hazard rate has several forms e.g. decreasing, constant, increasing, w shape, etc. The Weibull hazard rate is monotone increasing if β > 1 monotone decreasing if β < 1 and constant when β = 1.
Different PDF and CDF plots of the Weibull distribution are shown in Figures 1(a-b) and 2(a-b), respectively. Now consider the Exponentiated Inverted Weibull distribution and it was proposed by Flaih, Elsalloukh, Mendi, and Milanova (2012). It is the generalization of the Inverted Weibull distribution and it has two shape parameters. The probability density function and cumulative distribution function (CDF) of EIWD are given by (Figures 3 and 4): The corresponding survival function, hazard rate function and reversed function are, respectively, as: The kth moment of the EIW distribution is given as the following: Some special cases of EIWD: EIWD has the Inverted Weibull distribution and the Exponentiated Inverted Exponential distribution as its special cases.
Hence the EIW distribution is a generality of the EIED as well as IW distribution. And some basic properties distribution is estimated through the maximum likelihood estimation method (MLE).

Derivation of the Length-Biased Weighted Exponentiated Inverted Weibull Distribution
In this section, we drive the LBWEIWD and also its CDF. We also drive the shape of PDF and CDF of the LBWEIWD at various choices of parametric values.
Recall the probability density function of the EIWD And the weight function used is as follows: Substitute (3.1) and (3.2) in (2.1) then we get The density function (3.3) is known as Weighted Exponentiated Inverted Weibull Distribution.
Again substituting w(x) = x in (2.1) and using the density function of (3.3) we obtain the density function of the LBWEIWD of the form The CDF of (3.4) is defined as The plots of the LBWEIWD at various choices of parameters are shown in Figure 5(a-c) for comments. Figure 5 interprets the geometric behavior of PDF of the LBWEIWD. Figure 5(a-c) depict that as the parameter c changes for representative values of the parameters = 1, = 2 for c = 1, 2, 3, as the parameter θ changes for representative values of the parameters = 2, c = 2 for = 1, 2, 3 and as the parameter β changes for representative values of the parameters = 2, c = 3 for = 1, 2, 3, respectively. We can also notice that as the values of c, and β in Figure 5(a-c), respectively, increase, the "height" of the probability density function of LBWEIWD becomes lower and the PDF is more skewed right or it can be said that the LBWEIWD is positively skewed. The increasing parameters c, and β control the shape and skewness of the density. Now some plots of CDF of the LBWEIWD with specific parameter values are shown in Figure 6. Figure 6 shows the different CDF plots of the LBWEIWD for various combinations of parameters. We observe that the CDF of the proposed distribution increase slowly as the values of parameters increases. All the curves given in Figure 6(a-c) approaches to 1 when x ⟶ ∞.

Some particular cases of LBWEIWD
This section presents some sub-models that deduced from Equation (

The shape and mode of LBWEIWD
The shape and mode of LBWEIWD is confabulated in this section. The shape of the density function in (3.4) can be categorized by studding this function defined over the positive real line [0, ∞) and the behavior of its derivative as follows:

Limit and derivative of the function
The limit of the density function in (3.4) is as follows:

Mode
For these limits, we conclude the PDF of the LBWEIWD when w(x) = x −c and w(x) = x has one mode say x 0 as follows: (1) Taking logarithm of g(x) in (3.4) we get (2) Derivation for the x gives Its second derivative with respect to x is: (3) The mode of LBWEIWD is obtained by solving the non linear Equation (5.4) with respect to x and equating the first derivative to 0 leads to Now for different values of parameters , , and c the mode (x 0 ) of LBWEIWD is calculated and presented in Table 1.

Statistical properties of LBWEIWD
In this section, some basic properties involving moments, mean, variance, coefficient of Variation, coefficient of Skewness, and coefficient of Kurtosis, survival function, hazard rate, reversed hazard rate, and moment generating function are presented.

Moments of LBWEIWD:
If X follows the LBWEIWD with parameters β, θ, and c then kth moment of X, say ′ k , is given as From the kth moment of the LBWEIWD, the mean, variance and coefficient of Variation, Skewness (CS), and Kurtosis (CK), respectively, are as follows: Now for different values of parameters, Table 2 provided the values of the mean, variance, standard deviation (STD), coefficient of Variation (CV), coefficient of Skewness (CS), coefficient of Kurtosis (CK).

Moment Generating Function:
The moment generating function (mgf) of the LBWEIWD is given by

Hazard Rate Function:
The hazard rate function of the LBWEIWD is given by The hazard rate of the LBWEIWD with specific parameters values is displayed in Figure 8.

Reversed Hazard Rate Function:
The reversed hazard rate function of the LBWEIWD is given by The plot of reversed hazard rate with specific parameter values is made in Figure 9.

Parameter estimation
In this section estimation of the three parameters , and c of Length-Biased Weighted Exponentiated Inverted Weibull Distribution are obtained by the method of moments and maximum likelihood method:

Method of moments (MOM)
This procedure follows by equating the population moments to the sample moments then we have the system of three equations; Solving the system (i), (ii), and (ii), the moment's estimates of , and c can be determined. These approximations are mostly used as initial values for the MLE procedure when no closed form exists for the MLE and normal equations need to be resolved iteratively. (6.10)

Maximum-likelihood estimators (MLE)
Let suppose that sample was drawn from (3.4) and we drive the non linear equations for finding the maximum likelihood estimators of the parameters. It is assumed that X follows LBWEIWD (β, θ, c) and (β, θ, c) T denotes the parameter vector. The estimation of the parameters for the LBWEIWD is done via the MLE; we will provide the MLE procedure as follows: The likelihood function L( , , c; x) of the LBWEIWD , , and c , and then find the value of , , and c which maximize L( , , c; x) is given by The log likelihood function of n observations of X can be written as Differentiating (7.2) with respect to , and c, respectively, as follows: Solving the Equation (7.3) log L(X 1 , X 2 … X n , , , c) = 0 yields the maximum likelihood Solving the Equation (7.4) log L(X 1 , X 2 … X n , , , c) = 0 yields the maximum likelihood Solving the Equation (7.5) log L(X 1 , X 2 … X n , , , c) c = 0 yields the maximum likelihood To find out the likelihood estimators of , and c, we have to resolve the given system of non linear Equations (7.6)-(7.8) with respect to , and c but the non linear Equations (7.6)-(7.8) do not seem to have a closed form solution and must be solved iteratively to gain the estimate of the parameters or the system of the three equations can be solved consecutively. We take the MLE ̂,̂and ĉ, respectively, by resolving iteratively to zero. These equations are not in closed form and we can apply a numerical technique to attain the expectation of them, such as Newton-Raphson method to obtain the solution.
Example As in Section 8, we consider uncensored data on distance between cracks in a pipe dataset given in Table 4.The data can be modeled by LBWEIWD and also we estimate the parameters , and c using Newton-Raphson method for simultaneously three equations, by taking the initial estimate ̂ = 0.5, ̂ = 2, and ĉ = 1.5 then the estimate of values are: ̂ = 1.4256, ̂ = 100.7943, ĉ = 1.4857.

Data analysis
This section presents the comparison analysis of the LBWEIW, LBEIW, and EIW distributions applying to the real data-sets. We consider uncensored data on distance between cracks in a pipe dataset. The data is given in Table 4. We can also perform formal goodness-of-fit test in order to verify which distribution fits to the data. We apply the Kolmogorov-Smirnov test (KS test) for the goodnessof-fit purpose. Table 5 lists the MLE estimates of the parameters β, θ, and c values of the test statistics which is KS test. The p-value of KS for the LBWEIWD is 0.9388. The results in Table 5 showed that the LBWEIWD fits the data as well as LBIW, Weibull, LBEIW, and EIW distribution. In fact, based on the value of the KS-statistic, we observe that the LBWEIWD provides the best fit for this data among all the models considered. So it is evident that the LBWEIWD is the best distribution and is a strong competitor to other distributions commonly used in literature for fitting lifetime data as given below.

Conclusions
In this article, the Length-Biased version of weighted Exponentiated inverted Weibull distribution is introduced. A new three parameter distribution is proposed. At first the PDF of the weighted Exponentiated inverted Weibull distribution has been obtained considering weight as w(x) = x −c then from PDF of WEIWD, LBWEIWD have been obtained considering weight as w(x) = x. Next to characterize the distribution of random variable x of LBWEIWD, three functions have been introduced namely survival function or mass function and the failure rate function or hazard function. The moment's coefficient of variation coefficient of skewness coefficient or kurtosis of LBWEIWD has been derived, for estimation of parameters of LBWEIWD the maximum likelihood estimators have been used.