The Modi Exponentiated Exponential Distribution

In this study, a new four-parameter distribution called the Modi Exponen-tiated Exponential distribution was proposed and studied. The new distribution has three shape and one scale parameters. Its mathematical and statistical properties were investigated. The parameters of the new model were estimated using the method of Maximum Likelihood Estimation. Monte Carlo simulation was used to evaluate the performance of the MLEs through average bias and RMSE. The flexibility and goodness-of-fit of the proposed distribution were demonstrated by applying it to two real data sets and comparing it with some existing distributions.


Introduction
Probability distributions are essential tools in a variety of disciplines as they provide statistical interpretations that help make sense of data.However, the existing distributions such as the exponential distribution have limitations making them insufficient in modeling a variety of data.
The exponential distribution is a continuous probability distribution commonly used in statistics and probability theory.It is frequently used to describe the time between events in a process where events occur independently at a fixed average rate [1].
To overcome these limitations, many authors have proposed modifications of existing distributions.Gupta and Kundu [2] modified the exponential distribution by adding a shape parameter.However, their model is not flexible enough to control both skewness and kurtosis and also accommodate non-monotonic hazard rate shapes.
In the theory of probability distributions, choosing a particular probability distribution to represent real-life phenomena may be motivated by two factors such as tractability and flexibility [3].Although a probability distribution's tractability may be advantageous in theory since it is simple to use, particularly when simulating random samples, practitioners and other stakeholders may be more interested in a distribution's flexibility.In reality, it is preferable to employ probability distributions that best fit the given data set rather than transforming the data set because doing so may compromise the data set's original features.Due to this, numerous attempts have been made in recent years to ensure that the current standard theoretical distributions are improved and extended [4] [5] [6].Marshall and Olkin [7] proposed and studied a new method for adding a parameter to a family of distributions.Margaretha, et al. [8], in a conversation about employing Bayesian methodology to estimate parameters for the Exponentiated Exponential distribution in the context of left-censored data , in their work, Mahdavi and Kundu [9] introduced a new method for generating distributions and applied it to the exponential distribution.This generator was reported to control the skewness in the distribution which is not normal distributed.Many Researchers have made significant advancements in developing and extending various distributions.
Singh, et al. [10] and Niyoyunguruza, et al. [11] used the Marshall-Olkin generator method to extend distributions.They subsequently conducted parameter estimation using Maximum Likelihood Estimation (MLE) technique.By applying these extended distributions to a range of datasets, they demonstrated that these new distributions provided a better fit to the data when compared to the baseline distributions.Yahaya and Ieren [12] proposed the odd generalized exponential Gumbel distribution (OGEGD) for modeling lifetime data, Salem and selim [13] proposed Generalized Weibull-Exponential Distribution (GWED), the power Exponentiated family was proposed by Modi [14], Uwadi, et al. [15] introduced and studied the Exponentiated Gumbel Exponential (EGuE) distribution, and studied its mathematical and statistical properties.Modi, et al. [16], in their work proposed and studied a new family of distributions called Modi Exponential distribution, and applied it to two Real Datasets.In this study, we introduce the Modi Exponentiated Exponential distribution, which comprises four parameters.This distribution can be effectively used for fitting and analyzing data in various fields.
This paper is structured as follows: Section 2 provides a definition of the Modi generator, while Section 3 introduces the Exponentiated Exponential distribution.In Section 4, we present the Modi Exponentiated Exponential distribution, and develop its cumulative distribution function (CDF), probability density function (PDF), hazard rate function, and survival function.We also derive some statistical properties and mathematical properties of the proposed distribution.The estimation of the model parameters is given in Section 5. Simulation of study, Application to real data set, and Conclusion were given respectively in Sections 6, 7, and 8.

Modi Family
Modi, et al. [16] proposed and studied Modi family of distributions which is flexible and can be used to model a wide range of phenomena in various fields, including engineering, economics, and finance.The CDF, F(x) and PDF, f(x) of the Modi family are respectively, given by: for all 0 x > , where 0 α > , 0 β > and S(x) is the CDF of the existing distribu- tion, and s(x) is the PDF of the existing distribution.

Exponentiated Exponential Distribution
The practice of adding a new parameter to an existing family of distribution functions is a common one in statistical distribution theory.
Gupta and Kundu [2] proposed Exponentiated Exponential (EE) distribution.The CDF and PDF of EE are given respectively by: ( ) ( ) where , , 0 where λ is scale parameter, and δ is a shape parameter.

Modi Exponentiated Exponential (MEE) Distribution
In the field of statistical distribution theory, it is common to enhance the flexibility of a class of distribution functions by introducing an additional parameter.This practice can be very useful for data analysis, as it allows for greater versatility in modeling various types of data.

Probability Density Function and Hazard Rate Function
The PDF plays a crucial role in the modeling and analysis of continuous random variables.
The shapes of hazard rate function for the MEE are depicted in Figure 2. It demonstrates that the hazard rate function can take various shapes such as increasing, decreasing, constant, inverted bathtub, etc. Journal of Data Analysis and Information Processing 1 e log 1 1 ; ; ; ; where ( ) is the inverse of the distribution function of MEE and 0 1 u ≤ ≤ .
To obtain the values of the lower quartile, median, and upper quartile, one can use the quantile function by replacing u with 1/4, 1/2, and 3/4, respectively, giving: The lower quartile as: The median as: The upper quartile as:

Skewness and Kurtosis
Mathematically, the Moors Kurtosis and Galton skewness of the MEE distribution are stated as: where the quartiles and octiles value is indicated by Q(.) Table 1 provides information about how quantiles of the MEE distribution vary for different combinations of parameter values.It can be used to look up specific quantile values based on the desired parameters, which can be useful for statistical analysis and modeling.

The r th Moments of MEE Distribution
Calculation of the moments of a distribution is crucial for statistical analysis, especially in practical applications.Moments are used to determine various statistical measures such as measures of central tendency, dispersion, and shape.
The mathematical expression for the r th moment is given by ( ) where ( ) By substituting (4.3) in (4.10)we get By using binomial expansion, we By using gamma function the integral becomes: Comparing with this, we can identify t as x, p as (r + 1), and q as j λ .Finally, Table 2 illustrates the adaptability of the MEE distribution concerning its mean and variance.The Coefficient of Skewness (CS) values indicate its potential to exhibit right, or near symmetrical skewness.Similarly, the Coefficient of Kurtosis (CK) values suggest that the MEE distribution can display mesokurtic, leptokurtic, or platykurtic traits.These attributes highlight the MEE distribution's versatility, making it an attractive choice for modeling needs.The Probability Density Function (PDF) of the j th order statistic, with 1 j n ≤ ≤ , can be expressed as follows: by substituting Equations (4.3) and (4.2) in (5.8), the PDF of the j th order statistic for the MEEE can be expressed as follows: ( The PDF of smallest order statistic for the MEE occurs when the value of j is 1, and is given by 1 1 e 1 1 1 e e 1 1 e 1 1 e The PDF of largest order statistic for the MEE can be found when the index j is equal to n, and its PDF is given by 1 e e 1 1 e

Maximum Likelihood Estimation (MLE)
Maximum Likelihood Estimation (MLE) is a method of finding the values of the parameters that maximize the likelihood function.If we have n values , , , n x x x  that are randomly selected from the MEE distribution, the log-likelihood function denoted by ( ) l ϕ is given by: 1 e e log log 1 e 2 0 1 e Based on our observations, it is evident that Equations 5.4 -5.7 can not be solved analytically, meaning there is no direct mathematical solution for them.Consequently, we need to resort to a numerical optimization technique to find their solutions.In this study, the Broyden-Fletcher-Goldfarb-Shannon (BFGS) algorithm was employed to estimate the parameters of the MEE distribution.

Simulation Study
This section describes a Monte Carlo simulation study conducted to examine the behavior of Maximum Likelihood Estimators (MLEs) for the parameters of the MEE model.A simulation study was conducted to examine the accuracy of the Maximum Likelihood Estimators (MLEs) for four model parameters ( α , β , δ , and λ ) in terms of their average biases (ABs) and root mean squared errors (RMSEs).In order to obtain random samples from the MEE, the inverse of the CDF presented in Equation 4.5 was utilized.To accomplish this, we generated 1000 replications for each of the sample sizes The average biases were computed by: and the root mean square errors by: where η is the parameter in question, while ˆi η is its estimated value at the i th replication at each sample size, and R is the total number of replication.
The following table presents the MLEs, ABs, and RMSEs values for different sample sizes, corresponding to the parameters, λ , δ , α , β .Table 3 presents the simulation results for Maximum Likelihood Estimators (MLEs), Average Biases (ABs), and Root Mean Square Errors (RMSEs) for the parameter values in Set I and set II. Analyzing Table 3 in set I, it can be observed that as the sample size increases, the MLEs approach the true parameter values for the MEE distribution.Generally, the ABs, and RMSEs of the parameter estimators decrease as the sample size increases.
Table 3, on the other hand, displays the simulation results for the parameter values in Set II.
Indeed, the simulation results for the MEE distribution in Table 3 for set II demonstrate a similar pattern.As the sample size increases, the MLEs become closer to the true values of the parameters.Additionally, the Average Biases (ABs), and Root Mean Square Errors (RMSEs) for the parameter estimators tend to decrease, indicating improved accuracy, and precision as the sample size increases.
Sensitivity Analysis methodology was used to determine parameter values, and parameter estimations were carried out through Monte Carlo simulation using R software.

Application to Real Data Set
In this section, we fitted the Modi Exponentiated Exponential distribution to two real data sets and observed its flexibility compared to other well-known distributions.The analysis was conducted using R software.We calculated values of various information criteria such as Akaike Information Criterion (AIC), Hannan Quin Information Criterion (HQIC), Bayesian Information Criterion (BIC), and Consistent Akaike Information Criterion (CAIC).Additionally, we performed Kolmogorov-Smirnov (K-S), Cramér-von Mises test (W*), and Anderson-Darling test (A*) to assess the goodness of fit for the considered distributions.The distribution with the highest log-likelihood and the highest p-value for the K-S test and the lowest AIC, BIC, HQIC, CAIC, W*, A*, K-S was considered the best.
The PDFs of the existing distributions compared with the Modi Exponentiated Exponential distribution are presented in Table 4.
Data set I: This data set presented in [17]    Table 5 presents the summary characteristics of data set I. The data is skewed to the right, as indicated by a skewness coefficient of 3.325, and it exhibits significant tailing in its distribution, with a kurtosis coefficient of 16.15.
The parameter known as the Kurtosis coefficient was obtained from the data set using R software through descriptive statistics.
In Figure 3, the histogram of the data displays a right-skewed and the presence of outliers is effectively revealed by the box plot.
In Figure 4, the TTT plot of the data set, it can be observed that the hazard rate function is an inverted bathtub shape while the violin plot highlights that the majority of values are concentrated around the median.
Table 6 and Table 7 provide the AIC, HQIC, BIC, and CAIC values, along with the K-S, W*, and A* tests.Based on these results, the MEE model emerges as the most favorable choice because it has the lowest values for AIC, HQIC, BIC, CAIC, W*, K-S and A*, indicating better goodness-of-fit.Additionally, it exhibits the highest p-value for the K-S statistic and log-likelihood function value, further supporting its superiority.
Figure 5 illustrates a plot of fitted densities, comparing the MEE distribution to its sub-models using the bladder cancer data set.The plot reveals that the MEE distribution demonstrates a favorable and encouraging fit when compared to the existing distributions.
Table 8 presents the summary characteristics of data set II.The data is skewed to the right, as indicated by a skewness coefficient of 1.47277, and it can be   platykurtic based on kurtosis coeficient of 2.54029.The parameter known as the Kurtosis coefficient was obtained from the data set using descriptive statistics in the R software.
In the Figure 6, the TTT plot of the data set shows that the hazard rate function is on the rise, indicating an increasing shape and the histogram illustrates right skewed of the data.In the Figure 7, the violin plot emphasizes that most values are centered around the median.The box plot effectively identifies the presence of outliers in the data.
The MEE model is considered the best model based on the information provided in Table 9 and Table 10.This is because it has the lowest values for AIC, HQIC, BIC, W*, A* and CAIC, indicating better model fit.Additionally, it has the highest p-value for the K-S statistic and log-likelihood function value, further supporting its superiority compared to other models.
Figure 8 illustrates a plot of fitted densities, comparing the MEE distribution to its sub-models using the waiting time data set.The plot reveals that the MEE distribution demonstrates a favorable and encouraging fit when compared to the existing distributions.

Conclusion
In this paper, we introduced a new four-parameter model called the Modi Exponentiated Exponential (MEE) distribution, and applied to two real data sets.We have examined the mathematical and statistical properties of this proposed distribution.We derived expressions for its r th moment, survival function, hazard rate function, cumulative distribution, and quantile function.Furthermore, through various plots, we have observed that the MEE distribution exhibits different shapes, indicating its versatility in fitting data sets with diverse distributions.We have also obtained the Probability Density Function (PDF) of its minimum and maximum order statistics.Journal of Data Analysis and Information Processing   To estimate the parameters of the MEE distribution, we employed the method of maximum likelihood estimation.Monte Carlo simulation was used to assess the performance of MLEs.The study observed that MLEs demonstrate good accuracy and consistent estimating of model parameters.As the sample size increases, MLEs tend to approach the true values of the parameters, as indicated by the decreasing ABs.Additionally, RMSEs also decrease with increasing sample size.Our analysis demonstrates that MEE distribution outperforms the existing distributions, in modeling the two data sets considered in this study.

Figure 1
Figure 1 depicts several potential forms of the PDF for the MEE distribution across different parameter values.The MEE PDF has the potential to exhibit symmetry or unimodal as well as a reversed-J shape.Parameter values were derived through the technique of Sensitivity Analysis, a valuable technique that enhances the understanding of input-variable-to-modeloutput connections.It offers insights into uncertainties and empowers decision-makers to make well-informed and resilient choices.The hazard rate function is obtained by using Equations (4.2) and (4.3).

Figure 1 .
Figure 1.Plot of PDF of the MEE for various values of , , δ λ β and α .

Figure 2 .
Figure 2. Plot of hazard rate function of the MEE for various values of , , δ λ β and α .

3 )
Differentiating 5.3 partially with respect to each parameter and equating to zero gives 50,75,100,125, 200, 275,600,625,650 n = , and 700 n = using Equation 4.5 for various combinations of parameter values ( α , β , λ , δ ).The parameter values were provided in two different sets.

Figure 3 .
Figure 3. Histogram and box plots of bladder cancer data set.

Figure 4 .
Figure 4. TTT and violin plots of bladder cancer data set.

Figure 6 .
Figure 6.TTT and histogram plots of waiting time data set.

Figure 7 .
Figure 7. Violin and box plots of waiting data set.

Figure 8 .
Figure 8. Histogram and fitted densities of the waiting time data set for different distributions.

1. Quantile Function The
quantile function is crucial in generating random samples from a specified distribution.Journal of Data Analysis and Information Processing

Table 1 .
Quantiles of the MEE distribution for some parameter values.

Table 2 .
The first five moments, skewness, and kurtosis of the MEE distribution for different parameter values.

Table 3 .
The outcomes of a Monte Carlo simulation, investigation for the parameters in set I and set II.

Table 4 .
The existing distributions compared with Modi Exponentiated Exponential distribution.

Table 5 .
Comprehensive overview of the bladder cancer data set: Descriptive analysis.
(a) Histogram of cancer data set (b) Box plot of cancer data set

Table 6 .
Maximum Likelihood Estimates and goodness-of-fit tests for data set I.

Table 7 .
A summary of the results from the information criteria analysis conducted on data set I.
Figure 5.The densities of the bladder cancer data set estimated using different distribution models.Journal of Data Analysis and Information Processing

Table 8 .
Investigation of customer waiting duration for bank services: A descriptive analysis.

Table 9 .
Maximum likelihood estimation and goodness-of-fit analysis for data set II.

Table 10 .
A summary of the results from the information criteria analysis conducted on data set II. Journal of Data Analysis and Information Processing