Abstract

In this study, a new class that generates optimal univariate models called a new exponentiated-G class of distributions is developed. Numerous complementary statistical properties are derived and discussed in detail for the newly exponentiated power function (EPF) distribution. All possible shapes of the probability density and hazard rate functions are sketched for selected values of parameters. Six accredited estimation methods are discussed, and their performance is assessed and compared by a simulation study. The applicability of the new class is evaluated by analyzing the automotive engineering sector data.

1. Introduction

Modeling complicated problems is an enigma for applied researchers and practitioners. They seem to be worried about dealing with a variety of lifetime datasets that particularly follows physical and natural sciences. For this, they are searching for simple and efficient models. Consequently, a power function (PF) distribution is explored. It is a simple lifetime model as exponential and Pareto distributions. The PF distribution is a special case of the beta distribution. For more details about the classical work of the PF distribution, see Dallas [1] who developed a relationship between the Pareto and PF variables through an inverse transformation. Furthermore, for a deep understanding of the characterization of the PF distribution, see some credible work of Meniconi and Barry [2], Saran and Pandey [3], Chang [4], Tavangar [5], and Ahsanullah et al. [6].

In the most recent times, attention towards the generalization of probability distributions has grown phenomenally high. For more insight, see the trustworthy work of Cordeiro and Brito [7], Zaka and Akhter [8], Al Mutairi et al. [9], Tahir et al. [10], Shahzad et al. [11], Ahsan-ul-Haq et al. [12], Okorie et al. [13], Abdul-Moniem [14], Hassan et al. [15], Zaka et al. [16], Arshad et al. [17], Arshad et al. [18, 19], Al-Mutairi [20], Alzaatreh et al. [21], Gleaton and Lynch [22], Bourguignon et al. [23], Afify et al. [24], Tahir et al. [25], Aldahlan et al. [26], Aslam et al. [27], Balogun et al. [28], Afify et al. [29], Mansour et al. [30], Mahdavi and Kundu [31], Nassar et al. [32], Ijaz et al. [33], Klakattawi and Aljuhani [34], Afify et al. [35], Alsubie et al. [36], Ahmad et al. [37], and Nofal et al. [38].

To the best of our knowledge, the new class has not been discussed before, and it is the first time to explore the scenario particularly observed in the automobile sector. We develop a new class of distributions called the new exponentiated-G (NE-G) family and study one of its special submodels using the PF distributions as a baseline model. The studied model is called the exponentiated power function (EPF) distribution. The present study has some motivations as follows: (a) to develop new optimal models; (b) to advance the characteristics of the baseline models; (c) the density and hazard rate functions possess unimodal and bathtub-shaped curves, respectively; (d) to model the real-time scenario in the automobile sector.

This paper is outlined in the following sections. The development of the new class is presented in Section 2. A comprehensive discussion on mathematical and reliability measures is completed in Sections 3 and 4, respectively. Miscellaneous measures are discussed in Section 5. Six accredited estimation methods are discussed in Section 6. Simulation results are presented in Section 7. A lifetime application of the EPF distribution is discussed in Section 8, and finally, some conclusions are reported in Section 9.

2. The New Exponentiated-G Class

Tahir and Cordeiro [39] (Remarks 2 (ii)) developed the exponentiated generalized negative binomial G class which is defined by the CDF:

Chesneau et al. [40] reparameterized the parameters of (1) and provided the following CDF:

In this section, we provide a new generator of distributions that is very simple and capable of generating optimal alternative models. The new generator is called the new exponentiated-G class, and it is specified by the CDF:where is a CDF of any baseline (parent) distribution. It is based on a parametric vector depending on (r 1), is a shape parameter, and is a normalizing constant. , and . It is noted that the new class reduces to the baseline model with .

The probability density function (PDF) (f(x)) and the hazard rate function (HRF) (h(x)) of the new class reduce towhere

The quantile function (QF) (Q(x)) of the new class takes the form.

2.1. EPF Distribution

In this section, we discuss some useful characteristics of the EPF distribution. The PF distribution is specified by the following CDF and PDF:and

To this end, we define the analytical expressions of CDF and PDF of the new EPF distribution with two shape parameters and . The CDF of the EPF distribution has the form (for )

Its PDF reduces towhere is the normalizing constant and is a possible maximum assured life of a component.

The EPF distribution brings an additional shape parameter () that modulates the skewness and kurtosis tail weights of the baseline model. We note that a new parameter may offer a better fit to the unimodal, increasing, U-shaped, and bathtub-shaped failure rate data. The EPF distribution reduces to the baseline model (power function) for  = 1.

2.2. Asymptotic Properties of PDF and CDF

Asymptotes of the CDF and PDF at x are given by

Asymptotes of the CDF and PDF at x are given by

The derived expressions explore a dynamic effect of on the asymptotes of and

2.3. Shapes of PDF

Here, we discuss different shapes of the PDF of the EPF distribution. Figure 1 presents some different curves of the PDF for various choices of EPF parameters. We note that these curves can be decreasing, decreasing-increasing, and can be upside-down bathtub for  = 2.

Theorem 1. Let follow the EPF distribution with two shape parameters ; then, the r-th ordinary moment () of has the form

Proof. The following expression follows using (6) asThe prior expression can be rewritten by simplifying the expression as follows:After some algebra, the r-th ordinary moment of X reduces towhere
Table 1 presents some numerical results for the first-four ordinary moments (), variance =  , skewness =  , and kurtosis =  for different values of the EPF parameters.
Table 1 shows flexible and versatile behavior for moments, variance and alongside the skewness and kurtosis. The results indicate that the EPF distribution can be discussed for leptokurtic and skewed datasets.

Corollary 1. The first and second ordinary moments and the inverse moment () can be obtained by substituting and , in (7), respectively. The analytical expressions of mean, variance, and inverse moment are given, respectively, byand

Corollary 2. The factorial generating function is obtained directly followed by and it can be written for X as follows:

Theorem 2. If X ∼ EPF , then the moment generating function (MGF) () of is given by

Proof. The MGF is defined asAlso,
Hence, the MGF of X is obtained as

Theorem 3. .
If X ∼ EPF , then the characteristic function () of X is given by

Proof. The characteristic function is defined asHence, the characteristic function of X is obtained asVitality function is defined asIt is obtained for X asThe conditional moments are defined asIt is obtained for X as

3.1. Incomplete Moments and Associated Measures

Theorem 4. If X ∼ EPF , then the r-th lower incomplete moments of are

Proof. The r-th incomplete moments are defined asIt is obtained for X asHence, of X reduces to

Corollary 3. The first incomplete moment is obtained by substituting r = 1 in equation (33) asThe residual life function is defined asHence, the residual life function (RLF) and its associated CDF of X are given byandrespectively.
Furthermore, the reversed RLF is defined as . Hence, reversed RLF and its associated CDF of X take the forms:The mean RLF is defined as . It is obtained for X asThe mean inactivity time (MIT) is defined as It is obtained for X asThe strong mean inactivity time (SMIT) of a device is defined as . It is obtained for X asThe mean past lifetime (MPL) of a device is defined as MPL =  . It is obtained for X asFurthermore, the Lorenz , Bonferroni , and Zenga inequality curves have a significant role not only in the study of economics, the distribution of income, poverty, or wealth, but also they have a vital role in the fields of insurance, demography, medicine, and reliability engineering.

Theorem 5. If EPF, then the Lorenz inequality curve of is

Proof. Lorenz inequality curve is defined asIt is obtained for , using equations (17) and (34), as

Theorem 6. If EPF, then the Bonferroni inequality curve of is as follows:

Proof. The Bonferroni inequality curve is defined asIt is obtained for using and the CDF of the EPF model, as

Theorem 7. If EPF, then the Zenga inequality curve of reduces to

Proof. Zenga inequality curve is defined asIt is obtained for as

4. Reliability Function and Associated Measures

Probability distributions consider as a backbone for reliability engineering to analyze and predict the lifetime of a component/device. In this section, numerous notable reliability measures are discussed.

4.1. Survival Function

The survival function of takes the form:

4.2. Hazard Rate Function

The HRF (in demography), failure rate function (in engineering), and sometimes it is called the force of mortality (in economics). The HRF of is

4.3. Mean Time between Failures

Mean time between failures (MTBF) is defined as Hence, it is obtained for as

4.4. Cumulative HRF

The cumulative HRF is defined as Hence, the cumulative HRF of has the form:

Figure 2 presents the different curves of the EPF of HRF for various choices of its parameters. We note that it possesses increasing U-shaped and bathtub shape curves for  = 3.

4.5. Reverse HRF

The reverse HRF is defined as It is obtained for as

4.6. Odds Ratio

The odds ratio is defined as Hence, the odds ratio of is given by

4.7. Mills Ratio

The mill’s ratio is defined as Hence, the mill’s ratio of is

5. Miscellaneous Measures

5.1. Quantile Function

Theorem 8. If EPF, then the QF of is given by

Proof. The QF is defined byThe p-th QF of is obtained, by inverting the CDF (7), as

Corollary 4. The 1st quartile (Q1), median (Q2), and 3rd quartile (Q3) of are obtained by substituting p = 0.25, 0.5, and 0.75 in (61), respectively. The analytical expressions areandrespectively.

5.2. Skewness and Kurtosis

Bowley’s [41] and Moors’s [42] coefficients of skewness and kurtosis can be calculated by the following two equations:

Quartiles and octiles based on these descriptive measures provide more robust estimates than the traditional skewness and kurtosis measures. We note that these measures are almost less reactive to outliers and work more effectively for the distributions, deficient in moments. Figure 3 illustrates the skewness and kurtosis curves for the EPF distribution. We note that skewness and kurtosis are expressed as a function of . Figure 3 illustrates a positive to negative trend of skewness, and an increasing trend in kurtosis can be observed with the increase of .

5.3. Entropy Measures

Kurtosis and entropy measures have the same role in comparing the shapes and tail weights of various density functions. The entropy of a random variable is defined as a measure of uncertainty.

In this section, we have developed numerous well-known entropy measures including Rényi [43], Havrda and Charvat [44], and Mathai and Haubold [45]. For more details, see some promising work of Basit et al. [46], Dey et al. [47], and Ijaz et al. [48].

Theorem 9. If EPF, then the Rényi entropy of is

Proof. The Rényi entropy for is defined by
Using Equation (10), we can writeThen, integration givesHence, the Rényi entropy reduces toThe expression developed in (68) is quite helpful in the further computation of entropy measures of Havrda and Charvat, and Mathai and Haubold. The final expressions of Havrda and Charvat, and Mathai and Haubold entropy measures are presented in Table 2.
Table 3 presents the results of Rényi, Havrda and Charvat, and Mathai and Haubold entropy measures for some choices of model parameters for (  = 3), Set-I (), Set -II (), Set –III (), and Set-IV ().
A wide range of positive and negative values of entropy measures makes the EPF distribution more flexible and versatile.

5.4. Distribution of Order Statistics

Let be a random sample of size and their corresponding order statistics (OS) from the EPF distribution. The PDF of the i-th OS is i = 1, 2, 3, …, n.

The i-th OS density is obtained by incorporating Equations (6) and (8) in the last equation.

The minimum and maximum OS densities are obtained, respectively, by substituting in (70).

The i-th OS CDF is defined by

The i-th OS CDF of the EPF distribution reduces to

The median i-th OS PDF is

The Xm+1 median OS PDF of the EPF distribution has the form:

The i-th and j-th OS joint distribution is defined by

For the i-th and j-th OS joint distribution of the EPF model is as follows:

5.5. Bivariate and Multivariate Extensions

In this section, we develop the bivariate and multivariate extensions for the EPF distribution by following the Morgenstern family and the Clayton family

The CDF of the bivariate EPF distribution followed by the Clayton family for the random vector is

where , and .

The CDF of the bivariate EPF distribution followed by the Morgenstern family for the random vector is defined as

Let EPF (), and EPF (). Then, we set and

The CDF of the bivariate EPF distribution followed by the Clayton family for the random vector is

A simple n-dimensional extension of the last version for EPF distribution has the form:

6. Statistical Inference

In this section, we discuss six estimation techniques for the EPF parameters as follows: maximum likelihood estimators (MLEs), maximum product of spacing estimators (MPSEs), percentile estimators (PCEs), Cramér von–Mise distance estimators (CVMEs), Anderson–Darling estimators (ADEs), and right-tail Anderson–Darling estimators (RTADEs).

6.1. Maximum Likelihood Estimators

Let be a random sample of size from the EPF model; then, the likelihood function of is given by

The log  =  takes the form:

Let . The partial derivatives for the parameters and are and

The maximum likelihood estimates () of the EPF parameters can be obtained by maximizing (82) or by solving the above nonlinear equations simultaneously. These nonlinear equations although do not provide an analytical solution for the MLEs and the optimum values of and . Consequently, the Newton–Raphson type algorithm is an appropriate choice to obtain the MLEs.

6.2. Maximum Product of Spacing Estimators

The MPSEs are alternatives to the MLEs, and they are introduced by Cheng and Amin [49, 50]. Let be a uniform spacing of a random sample taken from the EPF distribution is defined by where denotes the uniform spacing, is j-th order statistics, and , and , and the MPSEs of the EPF parameters are obtained by maximizing

These estimators can also be obtained by solvingwhere

6.3. Percentile Estimators

The percentile method was introduced by Kao [51]. This method allows estimating the unknown parameters if the distribution function has a closed-form expression. Suppose be an unbiased estimator of . The PCEs of the EPF parameters are obtained by minimizingwith respect to and respectively.

6.4. Cramér von–Mise Estimators

Cramér [52] and Von Mises [53] introduced a relatively less-biased minimum distance estimator called the CVMEs. It can be obtained by making a difference between the estimates of the CDF and empirical CDF. The CVMEs of the EPF parameters are obtained by minimizingwith respect to and . Furthermore, the CVMEs follow by solving the nonlinear equations aswhere for t = 1, 2 is defined by (87).

6.5. Anderson–Darling and Right-Tail Anderson–Darling Estimators

Another type of minimum distance estimators is the ADEs. The ADEs of the EPF parameters are obtained by minimizingwith respect to and , respectively. The ADEs are also obtained by solving the following nonlinear equation aswhere (for t = 1, 2) is defined by (87).

The RTADEs of the EPF parameters can be determined by minimizingwith respect to and .

7. Simulation Experiment

In this section, we perform a simulation study to assess the behavior of different estimators in estimating the EPF parameters. We generate N = 1,000 replicates using (61) for several sample sizes n = 25, 50, and 100 with different combinations of the parameters. We calculate the average values of the estimates (AEs):and the mean square errors (MSEs):where .

The selection of the best estimation method will be made having a minimum value of MSEs. The R software (DEoptim package) is adopted to obtain the simulation results. The results of AEs and MSEs (in parenthesis) for the MLEs, MPSEs, PCEs, CVMEs, ADEs, and RTADEs are presented in Tables 47. It is noted that the AEs tend to their true parameter values, and the MSEs decrease with the increase in the sample size. This evidence is enough to favor that the estimators are unbiased asymptotically. All estimation methods perform efficiently for different combinations.

8. Application in Automobile Engineering

In this section, we analyze automobile engineering data. The data represent the time to failure (103 h) of turbocharger of one type of engine discussed by Xu et al. [54]. The observations are as follows: 1.6, 2.0, 2.6, 3.0, 3.5, 3.9, 4.5, 4.6, 4.8, 5.0, 5.1, 5.3, 5.4, 5.6, 5.8, 6.0, 6.0, 6.1, 6.3, 6.5, 6.5, 6.7, 7.0, 7.1, 7.3, 7.3, 7.3, 7.7, 7.7, 7.8, 7.9, 8.0, 8.1, 8.3, 8.4, 8.4, 8.5, 8.7, 8.8, 9.0. This dataset is analyzed by Afify et al. [55] and Nassar et al. [56].

The EPF distribution is compared with some well-known competitors, namely, the Weibull power function (W-PF) and the zero-truncated Poisson power function (ZTP-PF). Their CDFs are presented in Table 8. The criterion -log-likelihood (-LL), Akaike information criterion (AIC), along with the goodness-of-fit statistics such as Kolmogorov –Smirnov (KS) with its p-value, are adopted. Some descriptive statistics are presented in Table 9. Table 10 presents the estimates and standard errors (SEs) alongside the goodness-of-fit statistics as well. Based on the results in Table 10, we conclude that the EPF distribution provide better fit among all well-established competitors.

Furthermore, the empirical fitted density (i), estimated CDF (ii), probability-probability (PP) (iii), Kaplan–Meier survival (iv), along with the TTT plot (v), and box plot (vi) are illustrated in Figure 4, respectively. All the estimates and numerical results are calculated using the statistical software R, package adequacy model developed by Rafael et al. [57].

9. Conclusions

In this paper, we develop a new class that generates optimal univariate continuous models called the new exponentiated-G class. A special member of the proposed family called the exponentiated power function (EPF) distribution is studied in detail. Numerous statistical and reliability characteristics are discussed. Furthermore, the EPF distribution has flexible shapes for its density and hazard functions. For the estimation of EPF parameters, we followed the six accredited techniques named, MLEs, MPSEs, PCEs, CVMEs, ADEs, and RTADEs. A simulation experiment is performed to compare the performance of different estimation techniques. Our results show that the estimation techniques perform very well. The applicability of the EPF distribution is addressed using real-life data form the engineering field. The results show that the EPF distribution provides better fit as compared to other well-known competitors.

For some possible future studies, the EPF distribution can be adopted to analyze entropy measures following the works of Siddiqui et al. [59] and Rashid et al. [60].

Data Availability

This work is mainly a methodological development and has been applied on secondary data; but if required, data will be provided.

Conflicts of Interest

The authors declare that they have no conflicts of interest.