Next Article in Journal
YOLO-ABD: A Multi-Scale Detection Model for Pedestrian Anomaly Behavior Detection
Next Article in Special Issue
Reheating Constraints and the H0 Tension in Quintessential Inflation
Previous Article in Journal
Energy Efficiency for Faster-than-Nyquist Data Transmission Using Processing Algorithms with Decision Feedback
Previous Article in Special Issue
Multidimensional Representation of Semantic Relations between Physical Theories, Fundamental Constants and Units of Measurement with Formal Concept Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Generalization of the Inverse Generalized Weibull Distribution with Different Methods of Estimation and Applications in Medicine and Engineering

Department of Statistics, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia
*
Author to whom correspondence should be addressed.
Symmetry 2024, 16(8), 1002; https://doi.org/10.3390/sym16081002
Submission received: 19 June 2024 / Revised: 30 July 2024 / Accepted: 1 August 2024 / Published: 7 August 2024

Abstract

:
Limitations inherent to existing statistical distributions in capturing the complexities of real-world data often necessitate the development of novel models. This paper introduces the new exponential generalized inverse generalized Weibull (NEGIGW) distribution. The NEGIGW distribution boasts significant flexibility with symmetrical and asymmetrical shapes, allowing its hazard rate function to be adapted to many failure patterns observed in various fields such as medicine, biology, and engineering. Some statistical properties of the NEGIGW distribution, such as moments, quantile function, and Renyi entropy, are studied. Three methods are used for parameter estimation, including maximum likelihood, maximum product of spacing, and percentile methods. The performance of the estimation methods is evaluated via Monte Carlo simulations. The NEGIGW distribution excels in its ability to fit real-world data accurately. Five medical and engineering datasets are applied to demonstrate the superior fit of NEGIGW distribution compared to competing models. This compelling evidence suggests that the NEGIGW distribution is promising for lifetime data analysis and reliability assessments across different disciplines.

1. Introduction

Probability distributions are essential for data modeling in many disciplines, including economics, engineering, biology, business, and the medical sciences. Therefore, several extended distributions have been proposed to enhance the functionality and adaptability of the density and hazard rate functions to model data diversity. The techniques for extending distributions include compounding, adding parameters, composing, and transforming. For instance, the beta-generated method by [1], the Kumaraswamy-generated method by [2], and the transformed-transformer approach by [3] among others.
A new lifetime family named the exponential-X (NLTE-X) was developed by [4]. This family is based on the T-X generator with T E x p ( 1 ) and W ( F ( x ) ) = log 1 F ( x ) e θ F ( x ) . The cumulative distribution function (CDF) and probability density function (PDF) of NLTE-X are given by:
G ( x ; θ , ξ ) = 1 1 F ( x ; Θ ) e θ F ( x ; Θ ) ; θ > 0 , x > 0 ,
g ( x ; θ , ξ ) = f ( x ) { 1 + θ F ¯ ( x ; Θ ) } e θ F ( x ; Θ ) ; θ > 0 , x > 0 ,
where Θ is the vector of distribution parameters, and θ is a parameter of the NLTE-X family. Various distributions were generated from the NLTE-X family, such as the exponential Fréchet distribution [5], exponential inverted Topp–Leone distribution [6], and exponential-X power family of distribution [7].
Inverted distributions have attracted many researchers’ interest. Studies have shown that inverted distributions have more flexible density and hazard function structures than non-inverted ones. In addition, the applications of inverted distributions are significant in various domains, including the biological sciences, life test problems, chemical data, and medicine. For example, In reliability analysis, the inverse Weibull distribution (IW) can accurately model the lifetime of various systems [8,9]. The inverted Kumaraswamy distribution introduced by [10] was applied to precipitation, repairable items, and vinyl chloride data. Moreover, the inverted Topp–Leone distribution was applied to the failure times of Aarset data [11].
Recently, some generalizations of the inverse distributions were studied in the literature. These include the Kumaraswamy–inverse Weibull distribution proposed by [12], the Kumaraswamy inverse exponential distribution developed by [13], the alpha power inverse Rayleigh distribution developed by [14], the exponentiated inverse Rayleigh distribution introduced by [15], and the Weibull inverted exponential distribution by [16]. The Marshall–Olkin alpha power inverse Weibull distribution by [17], the alpha-power exponentiated inverse Rayleigh distribution proposed by [18], the odd Weibull inverse Topp–Leone distribution proposed by [19], the odds generalized exponential-inverse Weibull distribution proposed by [20], the generalized inverted Kumaraswamy distribution introduced by [21], the Kumaraswamy generalized inverse Lomax distribution developed by [22], and the inverse Weibull generator of distribution proposed by [23].
Moreover, a new generalization of both the IW and generalized inverse Weibull distribution (GIW) [24], named the generalized inverse generalized Weibull (GIGW), with CDF and PDF is given by
F G I G W ( x ) = 1 1 e γ ( λ / x ) β α ; x > 0 , γ , λ , α , β > 0 ,
f G I G W ( x ) = α λ β β β γ x ( β + 1 ) e γ ( λ / x ) β ( 1 e γ ( λ / x ) β ) α 1 ; x > 0 ,
where λ is the scale parameter and α , β , and γ are the shape parameters. Some of the properties of this distribution are studied by [25]. In addition, ref. [26] has obtained some estimators of the parameter of GIGW using maximum likelihood and the Bayesian estimation methods.
The main purpose of this article is to introduce a new generalization of the IW and inverse generalized Weibull based on the NLTE-X family of distributions called the new exponential generalized inverse generalized Weibull (NEGIGW). The significance of the NEGIGW distribution and its desirable characteristics include the following:
  • The NEGIGW will improve the features and adaptability of the density and hazard rate functions, accurately capturing the behavior of several real-world phenomena. The hazard rate function of the NEGIGW exhibits a wide range of forms, including decreasing, bath-tab and upside-down bath-tab, and reversed J-shape. The density can take symmetrical and asymmetrical shapes. This will enable NEGIGW to fit a wide range of data from the engineering, medicine, and reliability fields.
  • The NEGIGW introduces new generalizations of the IW, inverse generalized Weibull (IGW), and GIGW distributions by adding new parameters, thus increasing their flexibility and improving their ability to characterize tail shapes more accurately as observed from the different shapes of the NEGIGW density and hazard functions. Therefore, this generalization will help IW, IGW, and GIGW’s inability to fit real-world data such as the lifetime of some cancer data that showed a non-monotone failure rate when it attained a maximum after a finite period and then decreased gradually.
  • The CDF and hazard rate functions, moments, and entropy of the NEGIGW are in closed forms, which are useful in analyzing complete and censored data.
  • The applications of five medical and engineering data illustrate that the NEGIGW performs better than other competing lifetime models when modeling bladder cancer, failure of the engine’s turbocharger, fatigue Fracture of Kevlar 373/epoxy and failure and service times of aircraft windshield data.
Additionally, this study extends its investigation by evaluating the performance of various estimation methods for the parameters of the NEGIGW distribution. Three prominent techniques are compared: maximum likelihood (ML), maximum product of spacing (MPS), and percentile (PC) methods. A simulation study is conducted to assess the effectiveness of these estimators across a range of sample sizes and parameter values. The statistical analysis of the simulation results will provide valuable insights into the behavior and accuracy of each method under different conditions. Moreover, the applicability of the NEGIGW distribution is explored by demonstrating its superior fit to five practical applications compared to competing distributions. This comparative analysis strengthens the case for the NEGIGW distribution as a versatile tool for modeling real-world phenomena.
This article is structured as follows: Section 2 presents the NEGIGW with some graphical representations. We derived some of the NEGIGW properties in Section 3. In Section 4, three methods of estimation are used to estimate the parameter of the NEGIGW. Section 5, demonstrates extensive simulation studies to examine the performance of the various estimators. In Section 6, five applications in various domains are investigated to examine the NEGIGW’s effectiveness. Finally, some concluding remarks are made in Section 7.

2. New Exponential Generalized Inverse Generalized Weibull Distribution

The CDF and PDF of the NEGIGW can be derived by substituting (3) and (4) in (1) and (2) as follows
G ( x ) = 1 1 e γ ( λ / x ) β α e θ 1 1 e γ ( λ / x ) β α , x > 0 , θ , γ , α , λ , β > 0 ,
g ( x ) = α λ β β γ x ( β + 1 ) e γ ( λ / x ) β ( 1 e γ ( λ / x ) β ) α 1 1 + θ 1 e γ ( λ / x ) β α e θ 1 1 e γ ( λ / x ) β α .
The survival S ( x ) function is expressed as
S ( x ) = 1 e γ ( λ / x ) β α e θ ( 1 1 e γ ( λ / x ) β α ) .
The hazard rate function (HF), which is frequently utilized in lifetime modeling as it indicates the likelihood of failure, is defined as
H F = α λ β β γ x ( β + 1 ) e γ ( λ / x ) β 1 + θ 1 e γ ( λ / x ) β α ( 1 e γ ( λ / x ) β ) .
The NEGIGW distribution’s versatility in representing real-world data is evident in its ability to generate a wide range of shapes for its PDF and HF. The density plots in Figure 1 show decreasing, right-skewed, left-skewed, and symmetrical shapes, indicating its ability to fit complicated data. The HF plots of the NEGIGW in Figure 2 exhibit diverse asymmetrical forms such as unimodal, decreasing, and upside-down bathtubs, allowing for modeling failure rates that change over time. This flexibility of NEGIGW makes it a powerful tool for researchers and analysts. Researchers can tailor the NEGIGW to capture the underlying patterns of diverse real-world phenomena. This capability can lead to more precise modeling and improved decision-making across various fields, from finance and engineering to ecology and medicine.

2.1. Linear Representation for the Density of NEGIGW

A linear representation of the PDF of NEGIGW has been developed using mathematical expansions. The derivation begins with the binomial theorem, given by
( 1 x ) n = i = 0 n ( 1 ) i n i x i .
By applying (9), the PDF of NEGIGW will become
g ( x ) = v 1 = 0 1 α λ β β γ θ 1 v 1 v 1 x ( β + 1 ) e γ ( λ / x ) β ( 1 e γ ( λ / x ) β ) α + α v 1 1 e θ 1 1 e γ ( λ / x ) β α
The exponential expansion with a positive exponent is given by
e x = i = 0 x i i ! .
Then, the PDF will be
g ( x ) = v 1 = 0 1 v 2 = 0 1 v 1 α λ β β γ θ v 1 + v 2 e θ v 2 ! x ( β + 1 ) e γ ( λ / x ) β ( 1 e γ ( λ / x ) β ) α ( v 1 + v 2 + 1 ) 1
Applying (9) again, the PDF of NEGIGW will reduced to
g ( x ) = v 3 = 0 η v 3 λ β β γ x ( β + 1 ) e ( v 3 + 1 ) γ ( λ / x ) β ,
where
η v 3 = v 1 = 0 1 v 2 = 0 1 v 1 α ( v 1 + v 2 + 1 ) 1 v 3 ( 1 ) v 3 α θ v 1 + v 2 e θ v 2 ! .

2.2. Some Special Cases of NEGIGW

Special cases of NEGIGW can be obtained as follows:
  • When γ = 1 and α = 1 , the NEGIGW reduces to the exponential Fréchet distribution with parameters β , θ , and λ presented in [5].
  • When β = 1 and γ = 1 , the NEGIGW reduces to the exponential generalized inverse exponential distribution with parameters α , θ , and λ ; (not been previously studied).
  • When β = 2 and α = 1 , the NEGIGW reduces to the exponential generalized inverse Rayleigh distribution with parameters γ , θ , and λ ; (not been previously studied).

3. Properties of the NEGIGW

In this section, some characteristic properties of the NEGIGW are derived.

3.1. Quantile Function

The uth quantile function ( 0 < u < 1 ) of X NEGIGW can be derived by inverting (5) and solving the non-linear equation by the Lambert function W [ . ] .
x u = λ 1 γ log 1 W θ e θ ( 1 u ) θ 1 α 1 β , 0 u 1 .
The median of the NEGIGW can be obtained by substituting u = 0.5 into Equation (13).
M e d i a n ( x ) = λ 1 γ log 1 W θ e θ ( 0.5 ) θ 1 α 1 β .

3.2. Moment

The r t h moment of X NEGIGW is obtained as follows
μ r = E x r = 0 x r g ( x ) d x = v 3 = 0 η v 3 λ β β γ 0 x r x ( β + 1 ) e ( v 3 + 1 ) γ ( λ / x ) β d x .
By substituting y = ( v 3 + 1 ) γ ( λ / x ) β , therefore, the r t h moment can be defined as
μ r = v 3 = 0 η v 3 λ r γ r β ( v 3 + 1 ) r β 1 Γ 1 r β , β > r ,
where η v 3 is given by (12)

3.3. Moment Generating Function

The moment-generating function (MGF) of the NEGIGW is defined as
M X ( t ) = E ( e t x ) = r = 0 t r r ! μ r .
Therefore, substituting (16) into (17), the MGF is obtained as
M X ( t ) = r = 0 t r r ! v 3 = 0 η v 3 λ r γ r β ( v 3 + 1 ) r β 1 Γ 1 r β , β > r ,
where η v 3 is given by (12).

3.4. Characteristic Function

The NEGIGW characteristic function is obtained as
ϕ x t = E ( e i t x ) = r = 0 ( i t ) r r ! v 3 = 0 η v 3 λ r γ r β ( v 3 + 1 ) r β 1 Γ 1 r β , β > r ,
where η s is given by (12).

3.5. Rényi Entropy

The Rényi entropy can be used to determine the uncertainty measurement of the random variable X. When the Rényi entropy value is high, the data’s uncertainty level increases. According to [27], the Rényi entropy, R E ( δ ) , can be expressed as
R E ( δ ) = 1 1 δ log [ g ( x ) ] δ d x
By substituting g ( x ) given in (6) into the (20) and applying the expansions (10) and (9). R E ( δ ) is presented as
[ g ( x ) ] δ = v 3 = 0 η v 3 * β δ γ δ λ δ β x δ ( β + 1 ) e γ ( δ + v 3 ) ( λ / x ) β ,
where
η v 3 * = v 1 = 0 δ v 2 = 0 δ v 1 α ( v 1 + v 2 + δ ) δ v 3 ( 1 ) v 3 α δ δ v 2 θ v 1 + v 2 e θ δ v 2 ! .
By replacing (21) in (20), and calculating the integral, the Rényi entropy of the NEGIGW can be derived as
R E ( δ ) = 1 1 δ log v 3 = 0 η v 3 * β ( δ 1 ) λ ( 1 δ ) ( δ + v 3 ) 1 β [ δ ( β + 1 ) 1 ] γ 1 β ( δ 1 ) Γ 1 β [ δ ( β + 1 ) 1 ] .

3.6. Order Statistics

Let X i : n denote the ith order statistics for NEGIGW’s random sample (RS) X 1 , X 2 , , X n . Therefore, the ith order statistics PDF, f i : n ( x ) , is presented as
f i : n ( x ) = n ! ( i 1 ) ! ( n i ) ! f ( x ) F ( x ) i 1 1 F ( x ) n i .
Using (9), the PDF of X i : n becomes
f i : n ( x ) = k = 0 n ( 1 ) k n ! ( i 1 ) ! ( n i ) ! n k f ( x ) F ( x ) k + i 1 .
Substituting (5) into (25), the PDF of X i : n
f i : n ( x ) = k = 0 n ( 1 ) k n ! ( i 1 ) ! ( n i ) ! n k f ( x ) 1 ( 1 e γ ( λ / x ) β ) α e θ 1 1 e γ ( λ / x ) β α k + i 1 ,
where f ( x ) given by (6).

4. Estimation of Parameters

This section presents three estimation methods used to estimate the parameters of the NEGIGW.

4.1. ML Estimation Method

If x 1 , , x n are NEGIGW RS of size n, the log-likelihood function () for Θ = ( θ , γ , λ , α , β ) is as follows:
Θ = n log ( α ) + n β log ( λ ) + n log ( β ) + n log ( γ ) ( β + 1 ) i = 1 n log ( x i ) γ i = 1 n λ x i β + ( α 1 ) i = 1 n log 1 e γ ( λ / x i ) β θ i = 1 n 1 1 e γ ( λ / x i ) β α + i = 1 n log 1 + θ 1 e γ ( λ / x i ) β α .
The following are the first derivatives of (27), with respect to Θ = ( θ , γ , λ , α , β )
θ = i = 1 n 1 1 e γ ( λ / x i ) β α + i = 1 n ( 1 e γ ( λ / x i ) β ) α 1 + θ 1 e γ ( λ / x i ) β α ,
γ = n γ i = 1 n ( λ / x i ) β + ( α 1 ) i = 1 n e γ ( λ / x i ) β 1 e γ ( λ / x i ) β λ x i β + θ α i = 1 n e γ ( λ / x i ) β 1 e γ ( λ / x i ) β α 1 λ x i β + θ α i = 1 n e γ ( λ / x i ) β 1 e γ ( λ / x i ) β α 1 1 + θ 1 e γ ( λ / x i ) β α λ x i β ,
λ = n β λ γ β λ β 1 i = 1 n 1 x i β + ( α 1 ) β γ λ β 1 i = 1 n e γ ( λ / x i ) β 1 e γ ( λ / x i ) β 1 x i β + θ α β γ λ β 1 i = 1 n e γ ( λ / x i ) β 1 e γ ( λ / x i ) β α 1 1 x i β + α θ γ β λ β 1 i = 1 n e γ ( λ / x i ) β 1 e γ ( λ / x i ) β α 1 1 + θ 1 e γ ( λ / x i ) β α · 1 x i β ,
α = n α + i = 1 n log 1 e γ ( λ / x i ) β + θ i = 1 n 1 e γ ( λ / x i ) β α log 1 e γ ( λ / x i ) β + θ i = 1 n 1 e γ ( λ / x i ) β α log 1 e γ ( λ / x i ) β 1 + θ 1 e γ ( λ / x i ) β α , a n d
β = n log ( λ ) + n β i = 1 n log ( x i ) γ i = 1 n λ x i β log λ x i + γ ( α 1 ) i = 1 n e γ ( λ / x i ) β 1 e γ ( λ / x i ) β λ x i β log λ x i [ 1 e x ] + θ α γ i = 1 n e γ ( λ / x i ) β 1 e γ ( λ / x i ) β α 1 λ x i β log λ x i + α θ γ i = 1 n e γ ( λ / x i ) β 1 e γ ( λ / x i ) β α 1 1 + θ 1 e γ ( λ / x i ) β α λ x i β log λ x i .
Equations (28)–(32) can be solved by numerical techniques using any optimization approach, for instance, the Newton–Raphson method.

4.2. Maximum Product of Spacing Estimation Method

The MPS method, developed by Cheng and Tong [28], estimates the parameters of the NEGIGW distribution. This method analyzes the spacings between observations in a sample. The spacings, denoted by D i , are calculated for each observation i from 1 to n + 1 in the sample of size n using the formula.
D i = G ( x i : n ) G ( x i 1 : n ) , i = 1 , 2 , , n + 1 ,
where G ( x i : n ) represents the CDF of the ith ordered observation x i : n in the sample. The MPS estimate M is then obtained by averaging the logarithm of all the spacings:
M = 1 n + 1 i = 1 n + 1 log D i .
In essence, this method maximizes the geometric mean of the observational spacings. Computing the CDF for every observation is computationally efficient for smaller datasets. However, for large samples, this can become difficult. Furthermore, the robustness of the MPS method depends on the data coming from a NEGIGW distribution, and it can be susceptible to outliers that have a large effect on spacings.

4.3. Percentile Estimation Method

Percentiles are useful in descriptive statistics and parameter estimation, see [29]. The PC method equates the sample percentile points with the corresponding population percentile points for the NEGIGW (13). The PC estimates can be obtained by minimizing (35) concerning the parameters of NEGIGW as follows:
P = i = 1 n x i : n λ 1 γ log 1 W θ e θ ( 1 u i ) θ 1 α 1 β 2 .

5. Simulation Studies

This section evaluates the performance of the various estimation methods using numerical studies. We randomly generated N = 1000 samples from NEGIGW with sizes n = 30, 100, 200, and 500 for the following three sets of parameter values:
  • S e t I : θ = 20 , γ = 0.3 , λ = 0.05 , α = 5 , β = 2
  • S e t I I : θ = 20 , γ = 0.2 , λ = 0.09 , α = 5 , β = 1.2
  • S e t I I I : θ = 20 , γ = 0.2 , λ = 0.09 , α = 5 , β = 1.7
  • S e t I V : θ = 30 , γ = 2.2 , λ = 0.2 , α = 7 and β = 1.7 .
Three estimation techniques are performed to calculate the estimation for the NEGIGW parameters using Monte Carlo simulation, following the steps below.
  • Generate a random sample from the NEGIGW distribution with size n.
  • Calculate the ML, MPS, and PC estimations for each parameter π .
  • Repeat the steps from 1 to 2, N times.
  • For each parameter, calculate the average estimate, ( π ^ ) , and mean square error (MSE), where the MSE is defined as
    M S E ( π ^ ) = var ( π ^ ) + [ Bias ( π ^ ) ] 2 = 1 N i = 1 N ( π i ^ π t r ) 2 ,
    where π ^ represents the parameter estimate, π t r represents the true value of the parameter, and
    B i a s ( π ^ ) = 1 N i = 1 N π i ^ π t r .
The R programming is used for all estimation results [30]. An analysis of the parameter estimates for the NEGIGW distribution using ML, MPS, and PC methods reveals a clear trend, as shown in Table 1, Table 2, Table 3 and Table 4. The tables display the ML, MPS, and PC estimation values and corresponding MSE. The results show that, as the sample size increases, the MSE generally decreases for all three methods, indicating improved accuracy with more data. Additionally, parameter estimates themselves tend to converge towards the true values. From Table 1, Table 2 and Table 3, the performance of ML and PC methods are similar in terms of small MSE values. However, the MPS method is considered less efficient as it has larger MSE values, especially at some parameter estimation. Out of all the estimators, the ML estimator has the lowest MSE value. ML and PC are the most reliable choices for estimating NEGIGW parameters.

6. Applications

This section demonstrates the usefulness of the NEGIGW by utilizing five real-world data in various fields. The datasets are provided below.
Data 1: Remission Periods of Bladder Cancer Patients
The data present the remission periods in months of 128 bladder cancer patients, [31]:
0.082.093.484.876.948.6613.1123.630.202.233.52
4.986.979.0213.290.402.263.575.067.099.2213.80
25.740.502.463.645.097.269.4714.2425.820.512.54
3.705.177.289.7414.7626.310.812.623.825.327.32
10.0614.7732.152.643.885.327.3910.3414.8334.260.90
2.694.185.347.5910.6615.9636.661.052.694.235.41
7.6210.7516.6243.011.192.754.265.417.6317.1246.12
1.262.834.335.497.6611.2517.1479.051.352.875.62
7.8711.6417.361.403.024.345.717.9311.7918.101.46
4.405.858.2611.9819.131.763.254.506.258.3712.02
2.023.314.516.548.5312.0320.282.023.366.7612.07
21.732.073.366.938.6512.6322.69
Data 2: Failure of Engine’s Turbocharger
The data consist of 40 observations for the time (in 103 h) of the failure of a certain kind of engine’s turbocharger, [32]:
1.63.54.85.46.06.57.07.37.78.08.4
2.03.95.05.66.16.57.17.37.88.18.4
2.64.55.15.86.36.77.37.77.98.38.5
3.04.65.36.08.78.89.0
Data 3: Failure Times of Aircraft Windshield
The data present the failure times of 84 aircraft windshields [33]:
0.0401.8662.3853.4430.3011.8762.4813.4670.3091.8992.610
3.4780.5571.9112.6253.5780.9431.9122.6323.5951.0701.914
2.6463.6991.1241.9812.6613.7791.2482.0102.6883.9241.281
2.0382.8234.0351.2812.0852.8904.1211.3032.0892.9024.167
1.4322.0972.9344.2401.4802.1352.9624.2551.5052.1542.964
4.2781.5062.1903.0004.3051.5682.1943.1034.3761.6152.223
3.1144.4491.6192.2243.1174.4851.6522.2293.1664.5701.652
2.3003.3444.6021.7572.3243.3764.663
Data 4: Service Times of Aircraft Windshield
The data present the service times of 63 aircraft windshields [33].
0.0461.4362.5920.1401.4922.6000.1501.5802.6700.2481.719
2.7170.2801.7942.8190.3131.9152.8200.3891.9202.8780.487
1.9632.9500.6221.9783.0030.9002.0533.1020.9522.0653.304
0.9962.1173.4831.0032.1373.5001.0102.1413.6221.0852.163
3.6651.0922.1833.6951.1522.2404.0151.1832.3414.6281.244
2.4354.8061.2492.4644.8811.2622.5435.140
Data 5: Fatigue Fracture of Kevlar 373/epoxy
The data represent the life of fatigue fracture of Kevlar 373/epoxy subjected to constant pressure at 90 % stress level until all had failed [34].
0.02510.08860.08910.25010.31130.34510.47630.56500.56710.6566
0.67480.67510.67530.76960.83750.83910.84250.86450.88510.9113
0.91200.98361.04831.05961.07731.17331.25701.27661.29851.3211
1.35031.35511.45951.48801.57281.57331.70831.72631.74601.7630
1.77461.82751.83751.85031.88081.88781.88811.93161.95582.0048
2.04082.09032.10932.13302.21002.24602.28782.32032.34702.3513
2.49512.52602.99113.02563.26783.40453.48463.74333.74553.9143
4.80735.40055.44355.52956.55419.0960
The total time test (TTT) plot developed by [35] is a valuable graphical tool for determining if the data are suitable for a particular distribution. Figure 3 displays the TTT plots for the fifth dataset. It can be seen that the first dataset represents a bathtub hazard rate, and the second, third, fourth, and fifth datasets represent increasing hazard rate functions. The adequacy of the five datasets for the NEGIGW is determined by comparing its fit to the following distributions with their CDFs defined as
  • The Generalized Inverse Generalized Weibull Distribution (GIGW) given in (3).
  • The Exponential Fréchet (NEXF) Distribution [5]
    F ( x ) = 1 1 e ( λ / x ) β e θ e ( λ / x ) β .
  • The Exponentiated Generalized Inverse Weibull (EGIW) Distribution [31]
    F ( x ) = [ 1 { 1 e ( λ / x ) β } α ] q .
  • The Exponentiated Weibull Exponential (EWE) Distribution [36]
    F ( x ) = [ 1 e ( j x h ) k ] z .
  • The Inverse Weibull (IW) Distribution [37]
    F ( x ) = e ( λ / x ) β .
As noted in Section 5, the ML demonstrated better results. MLs of the parameters for each distribution were computed along with their corresponding log-likelihood values. In order to assess the effectiveness of the NEGIGW, various goodness of fit (GoF) criteria were employed, namely the corrected Akaike information criterion (CAIC), Akaike information criterion (AIC), Bayesian information criterion (BIC), Hannan–Quinn information criterion (HQIC), and the Kolomogorov–Smirnov (K-S) test. The p-value corresponding to the K-S test is calculated. The optimal model is characterized by the minimum value of these statistics and the maximum p-value.
Table 5, Table 6, Table 7, Table 8 and Table 9 summarize the MLs of the parameters, the log-likelihood, and the GoF for each model. The results in Table 5, Table 6, Table 7, Table 8 and Table 9 indicate that the NEGIGW has the lowest CAIC, AIC, BIC, HQIC, and K-S measures. The NEGIGW has the highest p-values among all the fitted models. Furthermore, Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8 display the density and CDF for the NEGIGW and the competitive distributions. The histogram represents the empirical density for the data and the dot black line represents the empirical CDF for the data. The Figures demonstrate that NEGIGW best matches the actual distribution of the examined datasets. Consequently, when compared to competing distributions, the NEGIGW is the most appropriate model for the analyzed data.

7. Concluding Remarks

This article introduces a new exponential generalized inverse generalized Weibull (NEGIGW) distribution based on the NLTE-X family. The suggested NEGIGW was motivated by the idea that generalization gives more flexibility in examining practical data. The NEGIGW’s hazard rate function takes several forms, allowing it to mimic various hazard behaviors in real-world settings such as medical, biological, engineering, and other applications. Important statistical properties are derived in close form. The estimation of the NEGIGW’s parameters is obtained using three methods of estimation; ML, MPS, and PC. An extensive simulation study is performed to compare the performance of these methods. The simulation results indicate that in terms of MSEs, the ML performs better than other methods. Five applications from medical, biological, and engineering fields are used to demonstrate the usefulness of the NEGIGW. We concluded that the proposed NEGIGW fits better than other competing models. This generalization is expected to lead to further lifetimes and reliability analysis applications.

Author Contributions

Conceptualization, I.A.A. and L.A.B.; Data curation, S.F.A.; Formal analysis, I.A.A.; Investigation, S.F.A.; Methodology, L.A.B. and I.A.A.; Software, L.A.B. and S.F.A.; Validation, I.A.A. and L.A.B.; Writing—original draft, S.F.A.; Writing—review and editing, L.A.B. and I.A.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data is contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Eugene, N.; Lee, C.; Famoye, F. Beta-normal distribution and its applications. Commun. Stat.-Theory Methods 2002, 31, 497–512. [Google Scholar] [CrossRef]
  2. Jones, M. Kumaraswamy’s distribution: A beta-type distribution with some tractability advantages. Stat. Methodol. 2009, 6, 70–81. [Google Scholar] [CrossRef]
  3. Alzaatreh, A.; Lee, C.; Famoye, F. A new method for generating families of continuous distributions. Metron 2013, 71, 63–79. [Google Scholar] [CrossRef]
  4. Huo, X.; Khosa, S.K.; Ahmad, Z.; Almaspoor, Z.; Ilyas, M.; Aamir, M. A new lifetime exponential-X family of distributions with applications to reliability data. Math. Probl. Eng. 2020, 2020, 1316345. [Google Scholar] [CrossRef]
  5. Alzeley, O.; Almetwally, E.M.; Gemeay, A.M.; Alshanbari, H.M.; Hafez, E.H.; Abu-Moussa, M. Statistical inference under censored data for the new exponential-X Fréchet distribution: Simulation and application to leukemia data. Comput. Intell. Neurosci. 2021, 2021, 2167670. [Google Scholar] [CrossRef] [PubMed]
  6. Metwally, A.S.M.; Hassan, A.S.; Almetwally, E.M.; Kibria, B.G.; Almongy, H.M. Reliability analysis of the new exponential inverted Topp–Leone distribution with applications. Entropy 2021, 23, 1662. [Google Scholar] [CrossRef] [PubMed]
  7. Hassan, F.; Hafez, E.A. The exponential-X power function (NEXPF) distribution and its applications to modelling reliability. J. Algebr. Stat. 2022, 13, 741–752. [Google Scholar]
  8. Khan, M.S.; Pasha, G.; Pasha, A.H. Theoretical analysis of inverse Weibull distribution. WSEAS Trans. Math. 2008, 7, 30–38. [Google Scholar]
  9. Mudholkar, G.S.; Srivastava, D.K.; Kollia, G.D. A generalization of the Weibull distribution with application to the analysis of survival data. J. Am. Stat. Assoc. 1996, 91, 1575–1583. [Google Scholar] [CrossRef]
  10. Abd AL-Fattah, A.; El-Helbawy, A.; Al-Dayian, G. Inverted Kumaraswamy distribution: Properties and estimation. Pak. J. Stat. 2017, 33, 37–61. [Google Scholar]
  11. Hassan, A.S.; Elgarhy, M.; Ragab, R. Statistical properties and estimation of inverted Topp-Leone distribution. J. Stat. Appl. Probab. 2020, 9, 319–331. [Google Scholar]
  12. Shahbaz, M.Q.; Shahbaz, S.; Butt, N.S. The Kumaraswamy-inverse Weibull distribution. Pak. J. Stat. Oper. Res. 2012, 8, 479–489. [Google Scholar] [CrossRef]
  13. Oguntunde, P.; Adejumo, A.O.; Owoloko, E. Application of Kumaraswamy inverse exponential distribution to real lifetime data. Int. J. Appl. Math. Stat. 2017, 56, 34–47. [Google Scholar]
  14. Malik, A.; Ahmad, S. A new inverse Rayleigh distribution: Properties and application. Int. J. Sci. Res. Math. Stat. Sci. 2018, 5, 92–96. [Google Scholar] [CrossRef]
  15. Rao, G.S.; Mbwambo, S. Exponentiated inverse Rayleigh distribution and an application to coating weights of iron sheets data. J. Probab. Stat. 2019, 2019, 7519429. [Google Scholar] [CrossRef]
  16. Eghwerido, J. A new Weibull inverted exponential distribution: Properties and applications. FUPRE J. Sci. Ind. Res. (FJSIR) 2022, 6, 58–72. [Google Scholar]
  17. Basheer, A.M.; Almetwally, E.M.; Okasha, H.M. Marshall-Olkin alpha power inverse Weibull distribution: Non bayesian and bayesian estimations. J. Stat. Appl. Probab. 2021, 10, 327–345. [Google Scholar]
  18. Ali, M.; Khalil, A.; Ijaz, M.; Saeed, N. Alpha-power exponentiated inverse Rayleigh distribution and its applications to real and simulated data. PLoS ONE 2021, 16, e0245253. [Google Scholar] [CrossRef] [PubMed]
  19. Almetwally, E.M. The odd Weibull inverse Topp–leone distribution with applications to COVID-19 data. Ann. Data Sci. 2022, 9, 121–140. [Google Scholar] [CrossRef] [PubMed]
  20. Hassan, A.S.; Elsherpieny, E.A.; Mohamed, R.E. Odds generalized exponential-inverse Weibull distribution: Properties & estimation. Pak. J. Stat. Oper. Res. 2018, 14, 1–22. [Google Scholar]
  21. Iqbal, Z.; Tahir, M.M.; Riaz, N.; Ali, S.A.; Ahmad, M. Generalized inverted Kumaraswamy distribution: Properties and application. Open J. Stat. 2017, 7, 645. [Google Scholar] [CrossRef]
  22. Ogunde, A.A.; Chukwu, A.U.; Oseghale, I.O. The Kumaraswamy Generalized Inverse Lomax distribution and applications to reliability and survival data. Sci. Afr. 2023, 19, e01483. [Google Scholar] [CrossRef]
  23. Hassan, A.S.; Nassr, S.G. The Inverse Weibull generator of distributions: Properties and applications. J. Data Sci. 2018, 16, 723–742. [Google Scholar] [CrossRef]
  24. De Gusmao, F.R.; Ortega, E.M.; Cordeiro, G.M. The generalized inverse Weibull distribution. Stat. Pap. 2011, 52, 591–619. [Google Scholar] [CrossRef]
  25. Jain, K.; Singla, N.; Sharma, S.K. The generalized inverse generalized Weibull distribution and its properties. J. Probab. 2014, 2014, 736101. [Google Scholar] [CrossRef]
  26. Rao, A.K.; Pandey, H. Parameter estimation of generalized inverse generalized Weibull distribution via Bayesian approach. Eur. Sch. J. 2022, 3, 23–28. [Google Scholar]
  27. Rényi, A. On measures of information and entropy. In Proceedings of the 4th Berkeley Symposium on Mathematics, Statistics and Probability, Berkeley, CA, USA, 20–30 June 1961; Volume 1. [Google Scholar]
  28. Cheng, R.; Amin, N. Maximum product-of-spacings estimation with applications to the lognormal distribution. Math Rep. 1979, 791. [Google Scholar]
  29. Schoonjans, F.; De Bacquer, D.; Schmid, P. Estimation of population percentiles. Epidemiology 2011, 22, 750–751. [Google Scholar] [CrossRef]
  30. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2024. [Google Scholar]
  31. Elbatal, I.; Muhammed, H.Z. Exponentiated generalized inverse Weibull distribution. Appl. Math. Sci. 2014, 8, 3997–4012. [Google Scholar] [CrossRef]
  32. Xu, K.; Xie, M.; Tang, L.C.; Ho, S.L. Application of neural networks in forecasting engine systems reliability. Appl. Soft Comput. 2003, 2, 255–268. [Google Scholar] [CrossRef]
  33. Tahir, M.H.; Cordeiro, G.M.; Mansoor, M.; Zubair, M. The Weibull-Lomax distribution: Properties and applications. Hacet. J. Math. Stat. 2015, 44, 455–474. [Google Scholar] [CrossRef]
  34. Abdul-Moniem, I.; Seham, M. Transmuted gompertz distribution. Comput. Appl. Math. J. 2015, 1, 88–96. [Google Scholar]
  35. Aarset, M.V. How to identify a bathtub hazard rate. IEEE Trans. Reliab. 1987, 36, 106–108. [Google Scholar] [CrossRef]
  36. Alsolami, E.; Alsulami, D. Combining two exponentiated families to generate a new family of distributions. Symmetry 2022, 14, 1739. [Google Scholar] [CrossRef]
  37. Keller, A.Z.; Kamath, A.R.R. Alternative reliability models for mechanical systems. In Proceedings of the Third International Conference on Reliability and Maintainability, Toulouse, France, 18–22 October 1982; pp. 411–415. [Google Scholar]
Figure 1. The NEGIGW density plots.
Figure 1. The NEGIGW density plots.
Symmetry 16 01002 g001
Figure 2. The NEGIGW HF’s plots.
Figure 2. The NEGIGW HF’s plots.
Symmetry 16 01002 g002
Figure 3. TTT plots for datasets.
Figure 3. TTT plots for datasets.
Symmetry 16 01002 g003
Figure 4. The NEGIGW is compared to other distributions for the first data. (Right): CDF for all distributions. (Left): observed and expected frequencies for all distributions.
Figure 4. The NEGIGW is compared to other distributions for the first data. (Right): CDF for all distributions. (Left): observed and expected frequencies for all distributions.
Symmetry 16 01002 g004
Figure 5. The NEGIGW is compared to other distributions for the second data. (Right): CDF for all distributions. (Left): observed and expected frequencies for all distributions.
Figure 5. The NEGIGW is compared to other distributions for the second data. (Right): CDF for all distributions. (Left): observed and expected frequencies for all distributions.
Symmetry 16 01002 g005
Figure 6. The NEGIGW is compared to other distributions for the third data. (Right): CDF for all distributions. (Left): observed and expected frequencies for all distributions.
Figure 6. The NEGIGW is compared to other distributions for the third data. (Right): CDF for all distributions. (Left): observed and expected frequencies for all distributions.
Symmetry 16 01002 g006
Figure 7. The NEGIGW is compared to other distributions for the fourth data. (Right): CDF for all distributions. (Left): observed and expected frequencies for all distributions.
Figure 7. The NEGIGW is compared to other distributions for the fourth data. (Right): CDF for all distributions. (Left): observed and expected frequencies for all distributions.
Symmetry 16 01002 g007
Figure 8. The NEGIGW is compared to other distributions for the fifth data. (Right): CDF for all distributions. (Left): observed and expected frequencies for all distributions.
Figure 8. The NEGIGW is compared to other distributions for the fifth data. (Right): CDF for all distributions. (Left): observed and expected frequencies for all distributions.
Symmetry 16 01002 g008
Table 1. Estimates and MSE of NEGIGW parameters for Set I.
Table 1. Estimates and MSE of NEGIGW parameters for Set I.
Set I: θ = 20 , γ = 0.3 , λ = 0.05 , α = 5 and β = 2 .
n MLMPSPC
EstimateMSEEstimateMSEEstimateMSE
30 θ ^ 19.9997(2.59  × 10 6 )20.1601(6.7444)19.9998(6.12  × 10 5 )
γ ^ 0.2978(1.15  × 10 4 )0.2546(1.07  × 10 2 )0.2999(6.31  × 10 5 )
λ ^ 0.0511(3.57  × 10 5 )0.0591(5.45  × 10 4 )0.0491(3.72  × 10 5 )
α ^ 4.9990(3.17  × 10 5 )5.2934(1.38  × 10 1 )4.99973(3.22  × 10 5 )
β ^ 1.9860(6.21  × 10 3 )2.1186(8.49  × 10 1 )2.0001(3.14  × 10 5 )
100 θ ^ 19.9999(1.39  × 10 7 )20.0946(2.2337)19.99990(5.91  × 10 5 )
γ ^ 0.2998(4.25  × 10 5 )0.2625(0.0046)0.2996(4.25  × 10 5 )
λ ^ 0.0501(2.48  × 10 6 )0.0579(0.0001)0.0488(4.61  × 10 5 )
α ^ 4.9999(1.85  × 10 6 )5.7276(7.0071)5.0005(5.61  × 10 5 )
β ^ 1.9987(4.35  × 10 4 )1.9516(0.0950)1.9996( 8.62  × 10 5 )
200 θ ^ 20.0000(2.26  × 10 11 )20.0853(1.3915)20.0008(8.66  × 10 4 )
γ ^ 0.2999(2.46  × 10 8 )0.2681(3.10  × 10 3 )0.2996(2.57  × 10 5 )
λ ^ 0.0500(1.65  × 10 7 )0.0565(9.60  × 10 5 )0.0484(1.06  × 10 4 )
α ^ 5.0000(3.02  × 10 10 )5.6960(6.0409)5.0006(3.54  × 10 4 )
β ^ 2.0000(2.07  × 10 8 )1.9491(3.20  × 10 2 )2.0005(3.67  × 10 4 )
500 θ ^ 20.0000(1.21  × 10 14 )20.0496(9.34  × 10 1 )20.0001(3.88  × 10 5 )
γ ^ 0.2999(4.81  × 10 10 )0.2798(1.23  × 10 3 )0.2996(7.50  × 10 5 )
λ ^ 0.0499(6.23  × 10 8 )0.0555(5.94  × 10 5 )0.0480(6.73  × 10 5 )
α ^ 5.0000(1.74  × 10 13 )5.9106(4.7893)5.0000(1.76  × 10 4 )
β ^ 1.9999(1.16  × 10 10 )1.9361( 1.42  × 10 2 )1.9998(5.54  × 10 5 )
Table 2. Estimates and MSE of NEGIGW parameters for Set II.
Table 2. Estimates and MSE of NEGIGW parameters for Set II.
Set II: θ = 20 , γ = 0.2 , λ = 0.09 , α = 5 and β = 1.2 .
n MLMPSPC
EstimateMSEEstimateMSEEstimateMSE
30 θ ^ 20.0022(0.0042)23.7632(1.00  × 10 2 )20.0003(4.00  × 10 5 )
γ ^ 0.1847(0.0006)0.2033(0.0208)0.1995(1.07  × 10 4 )
λ ^ 0.0936(0.0007)0.1396(0.0130)0.0893(4.40  × 10 5 )
α ^ 5.0066(0.02501)9.6769(1.55  × 10 2 )4.9999(6.99  × 10 5 )
β ^ 1.2451(0.0274)1.3933(0.6815)1.1999(1.18  × 10 4 )
100 θ ^ 20.0004(2.40  × 10 4 )21.1927(3.62  × 10 1 )19.9996(2.40  × 10 4 )
γ ^ 0.1895(2.57  × 10 4 )0.1956(0.007)0.1996(2.80  × 10 4 )
λ ^ 0.0939(2.16  × 10 4 )0.1121(0.0035)0.0893(1.78  × 10 4 )
α ^ 5.0014(2.42  × 10 4 )7.0864(5.710  × 10 1 )5.0004(6.74  × 10 5 )
β ^ 1.2111(7.83  × 10 3 )1.2303(0.1045)1.1994(2.50  × 10 4 )
200 θ ^ 20.0001(9.07  × 10 7 )20.1055(1.58  × 10 1 )19.9999(8.37  × 10 5 )
γ ^ 0.1924(1.19  × 10 4 )0.1902(0.0063)0.1999(9.21  × 10 5 )
λ ^ 0.0934(1.17  × 10 4 )0.1081(0.0011)0.0891(7.68  × 10 5 )
α ^ 5.0005(1.10  × 10 5 )6.7276(2.66  × 10 1 )5.0001(4.89  × 10 5 )
β ^ 1.2029(3.69  × 10 3 )1.2066(0.0508)1.1998(8.30  × 10 5 )
500 θ ^ 20.0000(2.83  × 10 7 )19.7828(7.2341)20.0001(3.02  × 10 5 )
γ ^ 0.1936(7.46  × 10 5 )0.1952(2.04  × 10 3 )0.2002(7.27  × 10 5 )
λ ^ 0.0927(4.37  × 10 5 )0.1000(4.11  × 10 4 )0.0886(7.27  × 10 5 )
α ^ 5.0002(3.48  × 10 6 )5.9146(1.04  × 10 1 )4.9998(5.85  × 10 5 )
β ^ 1.2003(1.20  × 10 3 )1.1851(1.03  × 10 2 )1.1998(4.06  × 10 5 )
Table 3. Estimates and MSE of NEGIGW parameters for Set III.
Table 3. Estimates and MSE of NEGIGW parameters for Set III.
Set III: θ = 20 , γ = 0.2 , λ = 0.09 , α = 5 and β = 1.7 .
n MLMPSPC
EstimateMSEEstimateMSEEstimateMSE
30 θ ^ 20.0014(9.23  × 10 5 )22.0293(6.02  × 10 1 )20.0001(2.26  × 10 4 )
γ ^ 0.1730(3.72  × 10 4 )0.1762(0.01069)0.1994(1.72  × 10 4 )
λ ^ 0.0964(3.72  × 10 4 )0.0064(0.0057)0.0881(1.39  × 10 4 )
α ^ 5.0052(7.00  × 10 4 )7.5851(9.55  × 10 1 )4.9997(1.65  × 10 4 )
β ^ 1.7500(4.54  × 10 2 )2.0051(1.3373)1.7005(1.37  × 10 4 )
100 θ ^ 20.0003(5.14  × 10 6 )20.5078(1.83  × 10 1 )19.9996(1.29  × 10 4 )
γ ^ 0.1828(5.26  × 10 4 )0.1706(0.0058)0.1993(1.53  × 10 4 )
λ ^ 0.0952(1.43  × 10 4 )0.1093(0.0017)0.0886(1.10  × 10 4 )
α ^ 5.0012(6.53  × 10 5 )6.5891(3.43  × 10 1 )5.0007(2.23  × 10 4 )
β ^ 1.7087(1.21  × 10 2 )1.7573(0.2168)1.7002(2.35  × 10 4 )
200 θ ^ 20.0001(2.39  × 10 6 )20.2710(1.59  × 10 1 )20.0011(1.37  × 10 3 )
γ ^ 0.1872(3.32  × 10 4 )0.1788(3.78  × 10 3 )0.1991(2.52  × 10 4 )
λ ^ 0.0942(8.62  × 10 5 )0.1062(9.79  × 10 4 )0.0879(2.11  × 10 4 )
α ^ 5.0005(3.05  × 10 5 )6.4389(2.48  × 10 1 )5.0012(6.82  × 10 4 )
β ^ 1.7005(5.91  × 10 3 )1.6945(7.93  × 10 2 )1.7006(5.68  × 10 4 )
500 θ ^ 20.0000(5.19  × 10 7 )19.8140(8.4763)19.9999(1.47  × 10 4 )
γ ^ 0.1935( 1.19  × 10 4 )0.1866(2.03  × 10 3 )0.1994(2.67  × 10 4 )
λ ^ 0.0922(2.74  × 10 5 )0.1040(4.98  × 10 4 )0.0873(1.59  × 10 4 )
α ^ 5.0000(6.77  × 10 6 )6.4854(1.39  × 10 1 )5.0000(4.24  × 10 4 )
β ^ 1.6973(1.44  × 10 3 )1.6559(2.79  × 10 2 )1.6998(1.50  × 10 4 )
Table 4. Estimates and MSE of NEGIGW parameters for Set IV.
Table 4. Estimates and MSE of NEGIGW parameters for Set IV.
Set IV: θ = 30 , γ = 2.2 , λ = 0.2 , α = 7 and β = 1.7 .
n MLMPSPC
EstimateMSEEstimateMSEEstimateMSE
30 θ ^ 29.9988(5.07  × 10 3 )32.5644(6.89  × 10 1 )29.9941(8.93  × 10 2 )
γ ^ 2.1964(3.86  × 10 4 )1.1960(1.4776)2.1613(2.90  × 10 1 )
λ ^ 0.1976(2.15  × 10 4 )0.3014(3.73  × 10 2 )0.1599(1.23  × 10 1 )
α ^ 6.9989(3.30  × 10 2 )10.0888(1.54  × 10 2 )6.9908(1.01  × 10 1 )
β ^ 1.7817(6.95  × 10 2 )2.372(3.5218)1.6899(1.35  × 10 1 )
100 θ ^ 30.0003(3.12  × 10 6 )31.5468(3.80  × 10 1 )29.9949(8.78  × 10 2 )
γ ^ 2.1983(4.07  × 10 6 )1.4975(8.25  × 10 1 )2.1667(2.05  × 10 1 )
λ ^ 0.1993(6.40  × 10 5 )0.2610(1.44  × 10 2 )0.1449(1.09  × 10 1 )
α ^ 7.0013(3.89  × 10 5 )9.4165(9.63  × 10 1 )6.9984(9.66  × 10 2 )
β ^ 1.7252(1.67  × 10 2 )1.8678(3.91  × 10 1 )1.7002(1.25  × 10 1 )
200 θ ^ 30.0001(1.09  × 10 6 )30.5452(1.69  × 10 1 )29.9957(7.04  × 10 2 )
γ ^ 2.1984(3.08  × 10 6 )1.6412(5.75  × 10 1 )2.1825(9.99  × 10 2 )
λ ^ 0.1999(2.86  × 10 5 )0.23843.98  × 10 3 0.1374(2.94  × 10 2 )
α ^ 7.00050(1.50  × 10 5 )8.3229(5.07  × 10 1 )7.0067(5.17  × 10 2 )
β ^ 1.7065(6.90  × 10 3 )1.7931(1.67  × 10 1 )1.715(6.47  × 10 2 )
500 θ ^ 30.0000(3.90  × 10 7 )30.2673(5.475)29.9715(8.06  × 10 4 )
γ ^ 2.1989( 1.61  × 10 6 )1.8225(2.98  × 10 1 )2.1789(4.43  × 10 4 )
λ ^ 0.1998(1.18  × 10 5 )0.2245(1.24  × 10 3 )0.1471(2.79  × 10 3 )
α ^ 7.0003(5.70  × 10 6 )8.0011(2.47  × 10 1 )6.9577(1.78  × 10 3 )
β ^ 1.7055(2.86  × 10 3 )1.734(6.45  × 10 2 )1.6625(1.40  × 10 3 )
Table 5. Measures of ML and GoF for the first data.
Table 5. Measures of ML and GoF for the first data.
DistributionsNEGIGWGIGWNEXFEGIWEWEIW
Estimates θ ^ = 40.0485 θ ^ = 11.1467 q ^ = 10.9153 j ^ = 0.4596
γ ^ = 8.9609 γ ^ = 6.4492 k ^ = 0.5434
λ ^ = 1.9341 λ ^ = 0.7863 λ ^ = 86.1435 λ ^ = 50.5093 h ^ = 0.7251 λ ^ = 3.2582
α ^ = 31.0659 α ^ = 22.4685 α ^ = 11.6619 z ^ = 5.8631
β ^ = 0.1495 β ^ = 0.3052 β ^ = 0.3818 β ^ = 0.1951 β ^ = 0.7520
−410.9638−413.7740−417.8249−424.7263−413.7165−444.0008
CAIC832.4194835.8732841.8433857.7778835.7582892.0975
AIC831.9276835.5480841.6498857.4526835.4330892.0015
BIC846.1877846.9561850.2059868.8607846.8412897.7056
HQIC837.7216840.1832845.1262862.0878840.0682894.3191
K-S0.04950.05780.07840.09340.05440.1407
p-value0.91180.78480.41010.21380.84190.01250
Table 6. Measures of ML and GoF for the second data.
Table 6. Measures of ML and GoF for the second data.
DistributionsNEGIGWGIGWNEXFEGIWEWEIW
Estimates θ ^ = 113.6416 θ ^ = 46.7096 q ^ = 88.8453 j ^ = 3.8269
γ ^ = 10.7248 γ ^ = 49.1384 k ^ = 1.0587
λ ^ = 4.6217 λ ^ = 0.1821 λ ^ = 41.7403 λ ^ = 12.8390 h ^ = 7.6894 λ ^ = 4.6721
α ^ = 97.6858 α ^ = 71.4909 α ^ = 14.9492 z ^ = 17.8675
β ^ = 0.3634 β ^ = 0.6764 β ^ = 0.7394 β ^ = 0.3051 β ^ = 1.9445
−84.98336−87.92892−88.71552−97.82884−93.78136−101.5917
CAIC181.7314185.0007184.0977204.8005196.7056207.5079
AIC179.9667183.8578183.4310203.6577195.5627207.1836
BIC188.4111190.6134188.4977210.4132202.3182210.5613
HQIC183.0199186.3004185.2630206.1003198.0053208.4049
K-S0.11750.13210.13720.20790.12610.2438
p-value0.63790.48730.43810.06280.54810.0172
Table 7. Measures of ML and GoF for the third data.
Table 7. Measures of ML and GoF for the third data.
DistributionsNEGIGWGIGWNEXFEGIWEWEIW
Estimates θ ^ = 62.9467 θ ^ = 38.6179 q ^ = 48.1213 j ^ = 7.0647
γ ^ = 7.4028 γ ^ = 4.3172 k ^ = 0.3127
λ ^ = 4.1741 λ ^ = 2.8629 λ ^ = 46.6470 λ ^ = 27.6062 h ^ = 0.1155 λ ^ = 1.4486
α ^ = 56.5685 α ^ = 80.5273 α ^ = 17.0032 z ^ = 69.9833
β ^ = 0.2438 β ^ = 0.3792 β ^ = 0.4541 β ^ = 0.1619 β ^ = 0.8387
−138.2479−146.3092−150.8618−175.9310−163.8430−194.5367
CAIC287.2651301.1247308.0236360.3683336.1922393.2215
AIC286.4959300.6184307.7236359.8620335.6859393.0733
BIC298.6500310.3417315.0160369.5853345.4092397.9350
HQIC291.3817304.5271310.6551363.7707339.5946395.0277
K-S0.10670.14550.16160.23340.18120.3127
p-value0.29370.05690.02480.00020.00791.4533  × 10 7
Table 8. Measures of ML and GoF for the fourth data.
Table 8. Measures of ML and GoF for the fourth data.
DistributionsNEGIGWGIGWNEXFEGIWEWEIW
Estimates θ ^ = 114.3494 θ ^ = 60.2657 q ^ = 43.8506 j ^ = 7.0741
γ ^ = 8.3534 γ ^ = 3.911 k ^ = 0.2408
λ ^ = 4.4939 λ ^ = 3.1471 λ ^ = 68.1877 λ ^ = 48.6810 h ^ = 0.0186 λ ^ = 0.9307
α ^ = 97.7212 α ^ = 77.6282 α ^ = 28.6196 z ^ = 69.9850
β ^ = 0.1560 β ^ = 0.2879 β ^ = 0.3961 β ^ = 0.1765 β ^ = 0.8102
−103.7424−108.1407−115.3408−144.9830−117.5368−131.3029
CAIC218.5375224.9711237.0883298.6557243.7633266.8058
AIC217.4849224.2814236.6816297.9660243.0736266.6058
BIC228.2005232.8539243.1110306.5385251.6461270.8921
HQIC221.6994227.6530239.2103301.3376246.4452268.2916
K-S0.14060.16260.21840.39780.17520.2215
p-value0.14990.06350.00411.65  × 10 9 0.03661.49  × 10 1
Table 9. Measures of ML and GoF for the fifth data.
Table 9. Measures of ML and GoF for the fifth data.
DistributionsNEGIGWGIGWNEXFEGIWEWEIW
Estimates θ ^ = 36.4168 θ ^ = 13.5573 q ^ = 0.5567 j ^ = 3.9945
γ ^ = 1.0614 γ ^ = 3.911 k ^ = 0.5579
λ ^ = 3.4524 λ ^ = 1.0997 λ ^ = 21.1117 λ ^ = 50.4851 h ^ = 1.3137 λ ^ = 0.8207
α ^ = 31.4465 α ^ = 0.6286 α ^ = 28.4964 z ^ = 7.0722
β ^ = 0.1811 β ^ = 0.8498 β ^ = 0.4038 β ^ = 0.4141 β ^ = 0.7588
−123.6436−127.4275−131.5939−127.7306−126.4810−153.5392
CAIC258.1444263.4185269.5211264.0246261.5253311.2428
AIC257.2873262.8551269.1877263.4612260.9619311.0784
BIC268.9409272.1780276.1799272.7841270.2848315.7399
HQIC261.9446266.5810271.9821267.1871264.6878312.9414
K-S0.08450.10440.13100.11770.09340.1893
p-value0.61790.35400.13430.22390.49187.37  × 10 3
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Alsaggaf, I.A.; Aloufi, S.F.; Baharith, L.A. A New Generalization of the Inverse Generalized Weibull Distribution with Different Methods of Estimation and Applications in Medicine and Engineering. Symmetry 2024, 16, 1002. https://doi.org/10.3390/sym16081002

AMA Style

Alsaggaf IA, Aloufi SF, Baharith LA. A New Generalization of the Inverse Generalized Weibull Distribution with Different Methods of Estimation and Applications in Medicine and Engineering. Symmetry. 2024; 16(8):1002. https://doi.org/10.3390/sym16081002

Chicago/Turabian Style

Alsaggaf, Ibtesam A., Sara F. Aloufi, and Lamya A. Baharith. 2024. "A New Generalization of the Inverse Generalized Weibull Distribution with Different Methods of Estimation and Applications in Medicine and Engineering" Symmetry 16, no. 8: 1002. https://doi.org/10.3390/sym16081002

APA Style

Alsaggaf, I. A., Aloufi, S. F., & Baharith, L. A. (2024). A New Generalization of the Inverse Generalized Weibull Distribution with Different Methods of Estimation and Applications in Medicine and Engineering. Symmetry, 16(8), 1002. https://doi.org/10.3390/sym16081002

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop