Abstract

The Weibull distribution has prominent applications in the engineering sector. However, due to its monotonic behavior of the hazard function, the Weibull model does not provide the best fit for data in many cases. This paper introduces a new family of distributions to obtain new flexible distributions. The proposed family is called a novel generalized- family. Based on this approach, an updated version of the Weibull distribution is introduced. The updated version of the Weibull distribution is called a novel generalized Weibull distribution. The proposed distribution is able to capture four different patterns of the hazard function. Some mathematical properties of the proposed method are obtained. Furthermore, the maximum likelihood estimators of the proposed family are also obtained. Moreover, a simulation study is conducted for evaluating these estimators. For illustrating the proposed model, two data sets from the engineering sector are analyzed. Based on some well-known analytical measures, it is shown that the novel generalized Weibull distribution is the best competing distribution for analyzing the engineering data sets.

1. Introduction

The two-parameter Weibull distribution is a famous statistical distribution defined on . It can be used as a suitable alternative choice for predicting, modeling, and analyzing lifetime data sets in healthcare, engineering, and other related areas. For data modeling in the engineering sector, numerous modifications of the Weibull distribution have been introduced and implemented, refer to Pedrosa et al. [1]; Bala and Napiah [2]; Huo et al. [3]; Liao et al. [4]; Shu et al. [5]; Li et al. [6]; and Bilal et al. [7].

On one side, the Weibull distribution generalizes the exponential and Rayleigh distributions whereas, on the other side, it can provide the characteristics of the other distributions. For some review studies based on the applications of the Weibull distribution, we refer to Nadarajah et al. [8]; Almalki and Nadarajah [9]; and Wais [10].

Let the random variable has the Weibull distribution with scale parameter and shape parameter , if its HF (hazard function) has the following form:

From the expression of in (1), we can see that the Weibull distribution offers data modeling with (i) decreasing HF shape, if , (ii) increasing HF shape, if , or (iii) constant HF shape, if .

In many cases, the data sets in the engineering sector follow the unimodal, bathtub, or modified unimodal HF shapes. In such cases, the Weibull distribution is not a suitable choice to implement for modeling the engineering data sets. In order to provide a close fit to such type of engineering data sets, numerous updated versions of the Weibull distribution have been produced and implemented. For example, (i) Ahmad et al. [11] studied a new Weibull–Weibull distribution for analyzing the strength of the alumina material; (ii) Zhao et al. [12] introduced TIHT–Weibull distribution for modeling the failure time of coating machines; (iii) Wang et al. [13] implemented the NG–Weibull distribution for analyzing the failure time of the electronic machine. For more recent modifications of the Weibull distribution, we refer to Klakattawi [14]; Wang et al. [15]; Tianshuai et al. [16]; Al-Sobhi [17]; Gonzalez et al. [18]; Emam and Tashkandy [19]; Prataviera [20]; and Rehman et al. [21].

To overcome the above deficiency of the Weibull distribution, we introduce a new updated form of the Weibull distribution. The new updated form of the Weibull distribution is introduced by proposing a new statistical family of distributions. The new family is called a novel generalized- (NGen-) family of distributions. The NGen- family can be used to update/increase the fitting power of the classical/traditional or other modified and generalized distributions.

Definition: a random variable V has the NGen- distributions, if its DF (distribution function) is given bywhere , and represents the baseline DF. In order to show that is a valid DF, a complete proof is provided in Proposition 1 and Proposition 2.

Proposition 1. Let’s consider the DF in Eq. (2), first, we have to show thatand

Proof. Since is a DF, we haveHence, from (5), we haveAlso, we haveAs we stated above that is a DF, we haveTherefore, from (8), we have

Proposition 2. The DF presented in Eq. (2) is a differentiable as well as a right continuous function.

Proof. The proof of Proposition 2 is very straightforward, and hence, omitted. Based on the mathematical derivations in Proposition 1 and Proposition 2, we can see that is a valid DF.
For , in link to , the PDF (probability density function) is given bywhere .
Corresponding to DF and PDF , the SF , HF , and cumulative HF (CHF) are given byandrespectively.
By using the DF in (2) as a base model, we can introduce an updated version of any base model. In this paper, we use as a DF of the Weibull model, to introduce a new version of the Weibull distribution. The updated version of the Weibull model is called a novel generalized-Weibull (NGen-Weibull) distribution. Some basic expressions of the NGen-Weibull distribution along with different plots of PDF and HF are obtained in Section 2.

2. A NGen-Weibull Distribution

Consider the DF of the Weibull distribution given by

with PDF given bywhere .

By incorporating (15) in (1), we obtain the DF of the NGen-Weibull distribution. Let V has the NGen-Weibull distribution, if its DF is given by

For , the PDF of the NGen-Weibull distribution is given by

Visual display of different plots of and is provided in Figure 1. From the plots of , we can see that the NGen-Weibull has four different patterns of PDF such as (i) reverse-J shaped (blue curve), (ii) right-skewed (red curve), (iii) left-skewed (green curve), and (vi) symmetrical (black curve).

The plots of in Figure 1 are obtained for (i) (blue curve), (ii) (red curve), (iii) (green curve), and (iv) (black curve).

Furthermore, for , the SF, HF, and CHF of the NGen-Weibull distribution are presented as follows:and

respectively.

A graphical illustration of and is presented in Figure 2. From the plots of , we can see that the NGen-Weibull distribution has the ability to model data with four different shapes of failure function. The HF shapes of NGen-Weibull distribution include (i) increasing shaped (blue curve), (ii) unimodal (green curve), (iii) modified unimodal (red curve), and (vi) decreasing (black curve).

The plots of in Figure 2 are obtained for (i) (blue curve), (ii) (red curve), (iii) (green curve), and (iv) (black curve).

The proposed model has certain advantages over the other classical and modified distributions, for example.(i)The proposed NGen-Weibull distribution is a simple extension of the Weibull model by adding one additional parameter. Most of the updated versions of the Weibull distribution have more than one additional parameter.(ii)The proposed NGen-Weibull distribution has four different shapes of the HF. There are only a few distributions that can model real-life data sets with unimodal and decreasing-increasing-decreasing shapes of HF.(iii)The proposed NGen-Weibull distribution obeys the HT properties. This fact indicates that the NGen-Weibull distribution can be a suitable candidate distribution for modeling the financial data sets. Due to the HT characteristics, the NGen-Weibull distribution can be used quite effectively for dealing with the extreme value data sets.(iv)The NGen-Weibull distribution provides the fit best to real-life engineering data sets as compared to other competing models, see Section 5.

3. Some Distributional Properties

This section offers the derivation of some distributional properties of the NGen- distributions. These properties include identifiability property, heavy-tailed property, quantile function, and moment.

3.1. The Identiability Property

In this subsection, we derive the identifiability property of the NGen- distributions. The parameter is called identifiable, if

Suppose and be the two parameters with DFs and , respectively. From the mathematical definition of the identifiability property, we have

Using (2) in (22), we get

3.2. The HT Property

A distribution with DF is called a HT model, if its SF satisfieswhere .

The HT distributions obey an important property called the regular variational property. A probability distribution is called regularly varying if it satisfieswhere the term is called the index of regular variation.

Here, we mathematical prove the regular variational property of the NGen- distributions. Based on the findings of Seneta [22], in terms of SF , we have the following:

Theorem 1. If a distribution with SFis regular varying distribution, thenis a regular varying distribution.

Proof. Suppose is finite but nonzero . Then, incorporating (2), we getSince (26) is nonzero, . Therefore, is the SF of the regular varying distribution.

3.3. The Quantile Function

The QF (quantile function) of the NGen- distributions with DF is derived aswhere and is the solution of . The expression in (27) can be implemented to generate random numbers from any special model of the NGen- family.

3.4. The Moment

Some basic distributional properties of any probability distribution can be derived and studied with the help of moments. This subsection offers the derivation of the moment of NGen- distributions. Let have the NGen- distributions with PDF , then, its moment can be obtained as

Using (12) in (28), we getwhere is the PDF of the exponentiated version of the NGen- distributions with the exponentiated parameter 2 whereas is the PDF of the exponentiated version of the NGen- distributions with the exponentiated parameter 3.

4. Estimation and Simulation

This part of the paper offers the derivation of the maximum likelihood estimators (MLEs) of the NGen- distributions with parameters and . Let be a set of sample taken from . The corresponding LF (likelihood function) is

Corresponding to (30), the log LF is

Corresponding to , the partial derivatives areand

Equating and to zero, and solving, we obtain the MLEs of .

After obtaining the expressions of the MLEs, now we provide a simulation study to assess performances of and . To carry out the evaluation of and through a simulation study, we generate RNs (random numbers) from the PDF using the inverse DF method.

The simulation results are obtained for (i) and (ii) . The simulation results are obtained using -function with -- algorithm and library. To check the performances of and , two statistical evaluation criteria such asand

were chosen. These two evaluation criteria were also calculated for .

In link to , the simulation results of the NGen-Weibull distribution are presented in Table 1 whereas the simulation results for are provided in Table 2.

From the simulation results obtained in Tables 1 and 2, it is obvious that as the size of increases.(i)The estimated values of and tend to be stable(ii)The MSEs of and decrease(iii)The biases of and tend to zero

5. Data Analyses

This section offers the illustration of the NGen-Weibull distribution using two data sets taken from the engineering sector. We apply the NGen-Weibull distribution to these data sets and compare its fitting results with the Weibull distribution, and its other versions include (i) Marshall–Olkin–Weibull (MO-Weibull), (ii) exponentiated Weibull (Exp-Weibull), and (iii) Kumaraswamy–Weibull (Kum-Weibull) distributions.

The DFs of the MO-Weibull, Exp-Weibull, and Kum-Weibull distributions are given, respectively, byand

Furthermore, we consider some analytical measures to compare the fitting results of the NGen-Weibull and other competing distributions. These measures include (i) AD (Anderson Darling) test, (ii) CM (Cramér-von Mises) test, and (iii) KS (Kolmogorov-Smirnov) test with -value. The numerical values of these measures are computed asand

respectively.

5.1. Data 1

This subsection offers the first illustration of the NGen-Weibull distribution by analyzing a data set representing the breaking stress of carbon fibres (in Gba). Previously, several authors have also studied this data set. For example, (i) Barreto-Souza et al. [23] considered this data set using the Frechet, beta Frechet (B-Frechet), and exponentiated Frechet (Exp-Frechet) distributions, (ii) Ogunde et al. [24] analyzed this data set using the Nadarajah Haghighi Gompertz (NH-Gompertz) distribution, and (iii) Eghwerido et al. [25] studied this data set using the inverse odd Weibull (IO-Weibull) and a different variant of the Weibull distribution.

Some key measures of Data 1 are given by minimum =  0.390, quartile = 1.840, median = 2.700, mean = 2.621, quartile = 3.220, maximum = 5.560, skewness = 0.3681541, kurtosis = 3.104939, variance = 1.027964, and range = 5.17. Corresponding to Data 1, some useful plots are obtained in Figure 3.

Corresponding to Data 1, the values of the MLEs , , and of the fitted models are obtained in Table 3. The profiles of the Log-LF of the MLEs of the NGen-Weibull model are plotted in Figure 4. Furthermore, the values of the analytical measures of the fitted models are presented in Table 4. Based on the numerical illustration presented in Table 4, we can see that the NGen-Weibull distribution is the best competing model.

Moreover, in support of the best fitting of the NGen-Weibull distribution using Data 1, a visual illustration of the NGen-Weibull distribution is presented in Figure 5. The plots in Figure 5 reveal that the NGen-Weibull distribution closely fits the fitted PDF, CDF, SF, HF, CFH, probability-probability (PP), and quantile-quantile (QQ) functions.

5.2. Data 2

This subsection offers the second illustration of the NGen-Weibull distribution by considering a data set representing the tensile strength, measured in GPa, of 69 carbon fibres. Previously, numerous authors have also analyzed this data set. For example, (i) Aldahlan [26] considered this data set using the Burr- exponentiated exponential (B EExp) distribution, (ii) Eyob et al. [27] studied this data set using the weighted quasi Akash (WQA) distribution, and (iii) Shukla and Shanker [28] analyzed this data set using the truncated Akash (TA) distribution.

The basic measures of Data 2 are given by minimum =  1.312, quartile = 2.150, median = 2.513, mean = 2.477, quartile = 2.816, maximum = 3.585, skewness = -0.1541874, kurtosis = 2.95123, variance = 0.2378321, and range = 2.273. Corresponding to Data 2, some basic plots are presented in Figure 6.

Corresponding to data 2, the numerical values of the MLEs , , and of the competing distributions are presented in Table 5. The profiles of the Log-LF of , and of the NGen-Weibull distribution are plotted in Figure 7. The values of the analytical measures of the competing distributions are obtained in Table 6. Based on the numerical findings in Table 6, we can see that the NGen-Weibull distribution has the smallest values of the analytical measures. This fact reveals that the NGen-Weibull model is the best model among the competing distributions.

After showing the best fitting of the NGen-Weibull distribution using Data 2 (see, Table 6), a graphical illustration is also provided in Figure 8. The plots of the fitted PDF, CDF, SF, HF, CFH, PP, and QQ functions are presented in Figure 8. These plots confirm the best fitting of the NGen-Weibull distribution.

6. Concluding Remarks

In this article, a new statistical approach for generating new probability distributions was introduced. The new statistical approach was named as NGen- distributions. Based on the proposed NGen- distributions method, an updated form of the Weibull distribution was introduced. The updated extension of the Weibull model was named as NGen-Weibull distribution. Certain distributional properties of the NGen- distributions were derived. The parameters of the NGen- distributions were estimated via implementing the maximum likelihood approach. A simulation study was also conducted to evaluate the estimators of the NGen- distributions. In the end, the NGen- distributions approach was illustrated by analyzing two data sets considered from the engineering sector. To both the engineering data sets, the NGen-Weibull distribution was applied in comparison with other competing models. Using certain statistical test (analytical measures), it is shown that the NGen-Weibull distribution was the best model.

In the future, we are motivated to implement the NGen-Weibull distribution for modeling real-life data sets in other fields such as the healthcare sector, hydrology, agriculture, metrology, geology, finance, banking, civil engineering, and education. We are also motivated to introduce the bivariate extension of the NGen-Weibull distribution for modeling the bivariate financial data sets such as import and export, among others.

Data Availability

The data sets are given in the paper.

Conflicts of Interest

The authors declare that they have no conflicts of interest.