# A Flexible Extension of Pareto Distribution: Properties and Applications

Authors: Huda M. Alshanbari; Abd Al-Aziz Hosni El-Bagoury; Ahmed M. Gemeay; E. H. Hafez; Ahmed Sedky Eldeeb
Journal: Computational Intelligence and Neuroscience
 (2021)

---

## Abstract

This paper introduced a relatively new mixture distribution that results from a mixture of Fréchet–Weibull and Pareto distributions. Some properties of the new statistical model were derived, such as moments with their related measures, moment generating function, mean residual life function, and mean deviation. Furthermore , different estimation methods were introduced for determining the unknown parameters of the proposed model. Finally, we introduced three real data sets which were applied to our distribution and compared them with other well-known statistical competitive models to show the superiority of our model for fitting the three real data sets, and we can clearly see that our distribution outperforms its competitors. Also, to verify our results, we carried out the existence and uniqueness test to the log-likelihood to determine whether the roots are global maximum or not.


---

## Body

## 1. Introduction

Modeling new phenomena is very important in the field of big data and data science. There are many ways for modeling and representing data. One of these ways is the statistical modeling for real data sets. Statistical modeling is very important in real-life sciences, as many applications and phenomena appear every era of time, so the continuous need for new distribution grows larger. As we know, many of the phenomena that arise nowadays need modeling, but, unfortunately, the traditional distribution could not model them. So, sometimes researchers turn to add new parameters, may be two parameters, to overcome these deficiencies in modeling new phonemes. But there is a new way to overcome all the deficiencies in the traditional distribution.This method is formulated by making a mixture from two or three distributions to formulate a new superior that can model all the data that the traditional ones failed to model. Many authors worked on the Pareto distribution; see [1], where the authors worked on the Pareto-IV distribution and estimated its parameters under accelerated life test, when the items were under type-II censored sample. Also, as an example for authors that worked on the Fréchet distribution, see [2], where the authors estimated the parameters of the Fréchet distribution under type -II censoring scheme using classical and Bayesian estimation methods.Mixture distribution may appropriately be utilized for specific data set where various subsets of the entire data set have various properties that can best be demonstrated independently. They can be more mathematically manageable because the individual mixture components deal with that more nicely compared to the overall mixture density. Applications of the mixture of distributions play an important role in reliability theory, insurance risk theory, and the oil industry. Willmot [3] presented the asymptotic tail behavior of Poisson mixtures with applications. Giudici et al. [4] made a novel methodology, dependent on mixtures of the product of Dirichlet process priors, which gave a formal inferential device to think about the logical influence of each covariate.Without characterizing the system, Bucar et al. [5] demonstrated that the reliability of this system could be approximated by utilizing a finite Weibull mixture distribution. Nakhi and Kalla [6] discussed the mixture of hyper-Poisson distribution with mixing a generalized gamma distribution and hyper-Poisson distribution generalized gamma mixtures.Panjer and Willmot [7] discover the estimator of the scale parameter in mixture models and the inadmissibility of the unusual estimator set up by displaying better estimators. They used these outcomes in mixtures of normal distributions and mixtures of exponential distributions. Karim et al. [8] introduced Rayleigh mixture distribution with various weight functions, and two correlated Rayleigh random variables have been determined.By presuming that the random variableX has a mixture of distributions if at least one parameter of the distribution of X is also a random variable. Let gx;θ be probability density function (PDF) of X, where θ is a parameter of the distribution of X. If θ is a random variable, then X has a mixture of distributions. The PDF of X is defined as(1)fx=∫θgx;θhθdθ,θis a continuous random variable,∑θgx;θhθ,θis a discrete random variable.Extreme point distributions have developed as one of the most important statistical fields for the applied sciences. Techniques of extreme point are also becoming heavily utilized in many other fields. Extreme point analyses often involve estimate of the likelihood of occurrences which are more extreme than any previously recorded event. Fréchet and Weibull distributions are the most important models for extreme values, and many statisticians have studied these models in many published papers according to their importance in many fields such as earthquakes, floods, engineering, physics, quality control, and medicine. For more information about Fréchet and Weibull distributions, see [1, 2, 9, 10].So, the main concern of this research is that we derive a mixture distribution called Fréchet–Weibull mixture Pareto distribution (FWMPD) from mixing Fréchet–Weibull distribution with Pareto distribution. This new mixture has a lot of significant advantages, which are very flexible and versatile. This distribution can model skewed and symmetric as well as asymmetric data. Now we will introduce the concept that we based our proposed distribution on.In order to make the paper easier for the reader, we sectioned and written the paper as follows: In Section2, we introduce the proposed distribution and the steps to formulate it. In Section 3, we deduce some of the statistical properties of the proposed distribution mathematically. In Section 4, we introduce eight different classical methods for estimating unknown parameters of the proposed model. For more about different kinds of classical methods of estimation, see [11–14] and [15]. In Section 5, we introduce three real data sets as an application to assess the performance of the distribution and to show its efficiency for fitting different real data sets. In Section 6, we introduce the conclusions illustrated from the paper along with the major findings.


## 2. The Mixture of Fréchet–Weibull and Pareto Distributions

The formulation of the new mixture model is presented in this part of the paper . The PDF and the cumulative distribution function (CDF) of the Fréchet–Weibull distribution [16] (X>0) is represented as follows:(2)fx=αkβαλαkx−1−αkexp−βαλxαk,Fx=exp−βαλxαk,where α and k are shape parameters and λ and β are scale parameters.The PDF and the CDF for the Pareto random variableX≥b are, respectively, given by(3)fx=abaxa+1,Fx=1−bxa,where b is a scale parameter and a is a shape parameter.If a random variableX follows Fréchet–Weibull distribution and by taking one of its four parameters (λ) as a random variable following Pareto distribution, then it is said to have FWMPD when its PDF and CDF are, respectively, defined as follows:(4)fx=abaβa/kx−a−1Γ1−akα,βαbxkα,x,a,b,α,β,k>0,(5)Fx=e−βαb/xαk−bxaβa/kΓ1−akα,βαbxkα,where a,α, and k are shape parameters, b and β are scale parameters, and Γ1−a/kα,βαb/xkα is upper incomplete gamma function.


### 2.1. Survival and Hazard Functions

The characteristics dependent on the reliability function and its correlated functions are very useful to study the example of any lifetime phenomenon. The survival function [Sx], hazard function [hx], and reverse hazard function [rx] of FWMPD are defined as follows:(6)Sx=1−e−βαb/xαk+bxaβa/kΓ1−akα,βαbxkα,hx=ab/xaβa/kΓ1−a/kα,βαb/xkαxb/xaβa/kΓ1−a/kα,βαb/xkα−e−βαb/xαk+1,rx=ab/xaβa/kΓ1−a/kα,βαb/xkαxe−βαb/xak−b/xaβa/kΓ1−a/kα,βαb/xkα,where Γ1−a/kα,βαb/xkα is upper incomplete gamma function.


### 2.2. Asymptotic Behavior

This section contains studies on the behaviors of PDF, CDF, andSx of FWMPD at x=0 and x=∞, respectively, as follows:(7)limx⟶0fx=abaβa/klimx⟶0x−a−1limx⟶0Γ1−akα,βαbxkα=abaβa/k×0×0=0,limx⟶∞fx=abaβa/klimx⟶∞x−a−1limx⟶∞Γ1−akα,βαbxkα=abaβa/k×0×Γ1−akα=0,limx⟶0Fx=limx⟶0e−βαb/xαk−βa/klimx⟶0bxaΓ1−akα,βαbxkα=0−βa/k×0=0,limx⟶∞Fx=limx⟶∞e−βαb/xαk−βa/klimx⟶∞bxaΓ1−akα,βαbxkα=1−βa/k×0×Γ1−akα=1,and since Fx+Sx=1, we have(8)limx⟶0Sx=1,limx⟶∞Sx=0.


### 2.3. Impact of Changing Parameters Values

In this section, we display the impact of changing parameters values on drawing PDF, CDF,Sx, and hx of FWMPD, which are graphed and plotted in Figures 1–4.Figure1(a) explains how the behavior of PDF of FWMPD is affected by increasing the value of parameter k, where α=0.5, β=3, a=1.5, and b=0.75, and Figure 1(b) explains how its behavior is affected by increasing the value of parameter β, where α=1, k=3, a=1.5, and b=2.Figure 1
The effects of parametersk and β on the PDF of FWMPD.
(a)(b)Figure2(a) shows how the behavior of CDF is changed when the significance increasing happened of the parameter a, as we can see this effect very clearly from the graph, where α=2, β=3, k=0.75, and b=0.5, and Figure 2(b) shows how the behavior of CDF is affected by increasing the value of parameter b, where α=0.5, β=0.75, a=1.5, and k=3.Figure 2
The effects of parametersa and b on the CDF of FWMPD.
(a)(b)Figure3(a) shows how the behavior of Sxis changed when the significance increasing happened of the parameter α, as we can see this effect very clearly from the graph, where a=2, β=0.5, k=0.75, and b=1 are still fixed, and Figure 3(b) shows how the behavior of Sx is affected with changing the value of parameter β, where a=5, α=3, k=1.5, and b=0.5.Figure 3
The effects of parametersα and β on Sx of FWMPD.
(a)(b)Figure4(a) shows how the behavior of hx is changed when the significanceincreasing happened of the parameter k, as we can see this effect very clearly from the graph, where a=4, β=0.5, α=1, and b=0.75, and Figure 4(b) shows how the behavior of hx is affected by the change of parameter β, where a=0.75, k=2, α=3, and b=0.5.Figure 4
The effects of parametersk and β on hx of FWMPD.
(a)(b)


## 2.1. Survival and Hazard Functions

The characteristics dependent on the reliability function and its correlated functions are very useful to study the example of any lifetime phenomenon. The survival function [Sx], hazard function [hx], and reverse hazard function [rx] of FWMPD are defined as follows:(6)Sx=1−e−βαb/xαk+bxaβa/kΓ1−akα,βαbxkα,hx=ab/xaβa/kΓ1−a/kα,βαb/xkαxb/xaβa/kΓ1−a/kα,βαb/xkα−e−βαb/xαk+1,rx=ab/xaβa/kΓ1−a/kα,βαb/xkαxe−βαb/xak−b/xaβa/kΓ1−a/kα,βαb/xkα,where Γ1−a/kα,βαb/xkα is upper incomplete gamma function.


## 2.2. Asymptotic Behavior

This section contains studies on the behaviors of PDF, CDF, andSx of FWMPD at x=0 and x=∞, respectively, as follows:(7)limx⟶0fx=abaβa/klimx⟶0x−a−1limx⟶0Γ1−akα,βαbxkα=abaβa/k×0×0=0,limx⟶∞fx=abaβa/klimx⟶∞x−a−1limx⟶∞Γ1−akα,βαbxkα=abaβa/k×0×Γ1−akα=0,limx⟶0Fx=limx⟶0e−βαb/xαk−βa/klimx⟶0bxaΓ1−akα,βαbxkα=0−βa/k×0=0,limx⟶∞Fx=limx⟶∞e−βαb/xαk−βa/klimx⟶∞bxaΓ1−akα,βαbxkα=1−βa/k×0×Γ1−akα=1,and since Fx+Sx=1, we have(8)limx⟶0Sx=1,limx⟶∞Sx=0.


## 2.3. Impact of Changing Parameters Values

In this section, we display the impact of changing parameters values on drawing PDF, CDF,Sx, and hx of FWMPD, which are graphed and plotted in Figures 1–4.Figure1(a) explains how the behavior of PDF of FWMPD is affected by increasing the value of parameter k, where α=0.5, β=3, a=1.5, and b=0.75, and Figure 1(b) explains how its behavior is affected by increasing the value of parameter β, where α=1, k=3, a=1.5, and b=2.Figure 1
The effects of parametersk and β on the PDF of FWMPD.
(a)(b)Figure2(a) shows how the behavior of CDF is changed when the significance increasing happened of the parameter a, as we can see this effect very clearly from the graph, where α=2, β=3, k=0.75, and b=0.5, and Figure 2(b) shows how the behavior of CDF is affected by increasing the value of parameter b, where α=0.5, β=0.75, a=1.5, and k=3.Figure 2
The effects of parametersa and b on the CDF of FWMPD.
(a)(b)Figure3(a) shows how the behavior of Sxis changed when the significance increasing happened of the parameter α, as we can see this effect very clearly from the graph, where a=2, β=0.5, k=0.75, and b=1 are still fixed, and Figure 3(b) shows how the behavior of Sx is affected with changing the value of parameter β, where a=5, α=3, k=1.5, and b=0.5.Figure 3
The effects of parametersα and β on Sx of FWMPD.
(a)(b)Figure4(a) shows how the behavior of hx is changed when the significanceincreasing happened of the parameter k, as we can see this effect very clearly from the graph, where a=4, β=0.5, α=1, and b=0.75, and Figure 4(b) shows how the behavior of hx is affected by the change of parameter β, where a=0.75, k=2, α=3, and b=0.5.Figure 4
The effects of parametersk and β on hx of FWMPD.
(a)(b)


## 3. Statistical Properties

In this part of the paper, we introduce the mathematical properties for the proposed distribution. These properties are the moments, moment generating function, mean residual life function, and the mean deviation of the proposed distribution.


### 3.1. Moments

In this subsection, we present therth moments of the proposed distribution. Now let μr′ be the rth about the origin of FWMPD and it is defined as follows:(9)μr′=∫x=0∞xrfxdx=abrβr/kΓ1−r/kαa−r,a>r.By settingr = 1, 2, 3, and 4, we can get so easily the first four moments by assigning FWMPD, respectively. Therefore, the mean and variance of FWMPD are given by(10)μ1′=μ=abβ1/kΓ1−1/kαa−1,σ2=ab2β2/ka−12Γ1−2/kα−a−2aΓ1−1/kα2a−2a−12,respectively, and, by using the moments about the origin, we can determine the first four central moments about the mean of FWMPD, which are given by the following relations:(11)μ1=μ1′−μ=0,μ2=μ2′−9μ1′2=ab2β2/ka−12Γ1−2/kα−a−2aΓ1−1/kα2a−2a−12,μ3=μ3′−3μ2′μ1′+2μ1′3=ab3β3/kaΓ1−1/kα2aΓ1−1/kα2−3a−12Γ1−2/kα/a−2a−13+Γ1−3/kαa−3,μ4=μ4′−4μ3′μ1′+6μ2′μ1′2−3μ1′4=ab4β4/k×aΓ1−1/kα−3a2Γ1−1/kα3−4a−13Γ1−3/kα/a−3+6aa−12Γ1−1/kαΓ1−2/kα/a−2a−14+Γ1−4/kαa−4,respectively, which will be used to determine coefficients of skewness, kurtosis, and variation, respectively, as follows:(12)β1=μ32μ23=a−2a×a−2a−13/a−3Γ1−3/kα−aΓ1−1/kα3a−12Γ1−2/kα−2a−2aΓ1−1/kα22aa−12Γ1−2/kα−a−2aΓ1−1/kα23,β2=μ4μ22=a−22a−14aΓ1−1/kα−3a2Γ1−1/kα3−4a−13Γ1−3/kα/a−3+6aa−12Γ1−1/kαΓ1−2/kα/a−2/a−14+Γ1−4/kα/a−4aa−12Γ1−2/kα−a−2aΓ1−1/kα22,CV=σμ×100=aa−12Γ1−2/kα−a−2aΓ1−1/kα2/a−2aΓ1−1/kα.


### 3.2. Moment Generating Function

The moment generating function of FWMPD is given by(13)Mt=abaλ−a−1∑m=0∞tmλmβm/kΓ1−m/kαm!,a>m,and its characteristic function is given by(14)ϕt=abaλ−a−1∑m=0∞itmλmβm/kΓ1−m/kαm!,a>m.


### 3.3. Mean Residual Life Function

The mean residual life function of a continuous random variableX and survival function SX following FWMPD is given by(15)μx=EX−x|X>x=1Sx∫x∞Sudu=1Sx∫x∞ufudu−x=abaβa/ka−1Sxx1−aΓ1−akα,βαbxkα+β1−a/kb1−aγ1−1kα,βαbxkα−x,a>1,where Γ1−a/kα,βαb/xkα and γ1−1/kα,βαb/xkα denote upper and lower in complete gamma function, respectively. We can notice that(16)μ0=abaβa/ka−1S00+β1−a/kb1−aΓ1−1kα=abβ1/kΓ1−1/kαa−1=Ex,which is an important property for μx.


### 3.4. Mean Deviation

The mean deviation about the mean for FWMPD is given by(17)MD=∫x=0∞x−μfxdx=2μFμ−1+2∫μ∞xfxdx=2μFμ−1+2abaβa/ka−1μ1−aΓ1−akα,βαbμkα+β1−a/kb1−aγ1−1kα,βαbμkα,a>1,where Γ1−a/kα,βαb/μkα and γ1−1k/α,βαb/μkα indicate upper and lower incomplete gamma function, respectively, and, by changing μ with any measure of central tendency, we can find its mean deviation.


## 3.1. Moments

In this subsection, we present therth moments of the proposed distribution. Now let μr′ be the rth about the origin of FWMPD and it is defined as follows:(9)μr′=∫x=0∞xrfxdx=abrβr/kΓ1−r/kαa−r,a>r.By settingr = 1, 2, 3, and 4, we can get so easily the first four moments by assigning FWMPD, respectively. Therefore, the mean and variance of FWMPD are given by(10)μ1′=μ=abβ1/kΓ1−1/kαa−1,σ2=ab2β2/ka−12Γ1−2/kα−a−2aΓ1−1/kα2a−2a−12,respectively, and, by using the moments about the origin, we can determine the first four central moments about the mean of FWMPD, which are given by the following relations:(11)μ1=μ1′−μ=0,μ2=μ2′−9μ1′2=ab2β2/ka−12Γ1−2/kα−a−2aΓ1−1/kα2a−2a−12,μ3=μ3′−3μ2′μ1′+2μ1′3=ab3β3/kaΓ1−1/kα2aΓ1−1/kα2−3a−12Γ1−2/kα/a−2a−13+Γ1−3/kαa−3,μ4=μ4′−4μ3′μ1′+6μ2′μ1′2−3μ1′4=ab4β4/k×aΓ1−1/kα−3a2Γ1−1/kα3−4a−13Γ1−3/kα/a−3+6aa−12Γ1−1/kαΓ1−2/kα/a−2a−14+Γ1−4/kαa−4,respectively, which will be used to determine coefficients of skewness, kurtosis, and variation, respectively, as follows:(12)β1=μ32μ23=a−2a×a−2a−13/a−3Γ1−3/kα−aΓ1−1/kα3a−12Γ1−2/kα−2a−2aΓ1−1/kα22aa−12Γ1−2/kα−a−2aΓ1−1/kα23,β2=μ4μ22=a−22a−14aΓ1−1/kα−3a2Γ1−1/kα3−4a−13Γ1−3/kα/a−3+6aa−12Γ1−1/kαΓ1−2/kα/a−2/a−14+Γ1−4/kα/a−4aa−12Γ1−2/kα−a−2aΓ1−1/kα22,CV=σμ×100=aa−12Γ1−2/kα−a−2aΓ1−1/kα2/a−2aΓ1−1/kα.


## 3.2. Moment Generating Function

The moment generating function of FWMPD is given by(13)Mt=abaλ−a−1∑m=0∞tmλmβm/kΓ1−m/kαm!,a>m,and its characteristic function is given by(14)ϕt=abaλ−a−1∑m=0∞itmλmβm/kΓ1−m/kαm!,a>m.


## 3.3. Mean Residual Life Function

The mean residual life function of a continuous random variableX and survival function SX following FWMPD is given by(15)μx=EX−x|X>x=1Sx∫x∞Sudu=1Sx∫x∞ufudu−x=abaβa/ka−1Sxx1−aΓ1−akα,βαbxkα+β1−a/kb1−aγ1−1kα,βαbxkα−x,a>1,where Γ1−a/kα,βαb/xkα and γ1−1/kα,βαb/xkα denote upper and lower in complete gamma function, respectively. We can notice that(16)μ0=abaβa/ka−1S00+β1−a/kb1−aΓ1−1kα=abβ1/kΓ1−1/kαa−1=Ex,which is an important property for μx.


## 3.4. Mean Deviation

The mean deviation about the mean for FWMPD is given by(17)MD=∫x=0∞x−μfxdx=2μFμ−1+2∫μ∞xfxdx=2μFμ−1+2abaβa/ka−1μ1−aΓ1−akα,βαbμkα+β1−a/kb1−aγ1−1kα,βαbμkα,a>1,where Γ1−a/kα,βαb/μkα and γ1−1k/α,βαb/μkα indicate upper and lower incomplete gamma function, respectively, and, by changing μ with any measure of central tendency, we can find its mean deviation.


## 4. Classical Methods of Estimation

This section discusses the conventional techniques for estimating the suggested model parametersθ=a,b,α,β,k⊤ by eight different classical estimation methods. Many papers discussed these methods (for more information, see [17–21]). Determining the estimated parameters in explicit form is mathematically complicated, so these estimates will be obtained numerically by using Wolfram Mathematica software version 12.0.


### 4.1. Classical Methods for the Complete Sample

In this subsection, we introduce eight methods of estimation which were used for estimating the parameters of the proposed distribution.


#### 4.1.1. Maximum Likelihood Estimates (MLEs)

Letx1,…,xn is a randomized sample having a size n from the PDF (3). So, the log-likelihood function for θ is as follows:(18)ℓ=nalogb+nalogβk+nloga−a+1∑i=1nlogxi+∑i=1nlogΓ1−akα,βαbxikα.The MLEs ofθ can be obtained by maximizing ℓ.


#### 4.1.2. Ordinary Least-Squares Estimates (OLSEs)

Letx1:n,…,xn:n be the corresponding order statistics. The OLSEs for the distribution parameters can be easily obtained by making the following equation at minimum value, by using any mathematical software, may be MATHEMATICA 12 or any advanced program.(19)O=∑i=1nFxi:n−in+12=∑i=1ne−βαb/xi:nαk−bxi:naβa/kΓ1−akα,βαbxi:nkα−in+12.


#### 4.1.3. Weighted Least-Squares Estimates (WLSEs)

By minimizing the following equation, the WLSEs of proposed model parameters can be computed:(20)W=∑i=1nn+12n+2in−i+1Fxi:n−in+12=∑i=1nn+12n+2in−i+1e−βαb/xi:nαk−bxi:naβa/kΓ1−akα,βαbxi:nkα−in+12.


#### 4.1.4. Anderson–Darling Estimates (ADEs), Right-Tail Anderson–Darling Estimates (RTADEs), and Left-Tail Anderson–Darling Estimates (LTADEs)

The ADEs for the distribution parameters can be easily obtained by making the following equation at minimum value, by using any mathematical software, may be MATHEMATICA 12 or any advanced program.(21)A=−n−1n∑i=1n2i−1logFxi:n+logSxi:n.The RTADEs for the distribution parameters can be easily obtained by making the following equation at minimum value, by using any mathematical software, may be MATHEMATICA 12 or any advanced program.(22)R=n2−2∑i=1nFxi:n−1n∑i=1n2i−1logSxi:n.The LTADEs for the distribution parameters can be easily obtaine by making the following equation at minimum value, by using any mathematical software, may be MATHEMATICA 12 or any advanced program.(23)L=−32n+2∑i=1nFxi:n−1n∑i=1n2i−1logFxi:n.


#### 4.1.5. Cramér–von Mises Estimates (CVMEs) and Maximum Product of Spacing Estimates (MPSEs)

The CVMEs are determined by minimizing(24)CV=112n+∑i=1nFxi:n−2i−12n2=112n+∑i=1ne−βαb/xi:nαk−bxi:naβa/kΓ1−akα,βαbxi:nkα−2i−12n2.The MPSEs are determining by maximizing the following equation:(25)G=1n+1∑i=1n+1logDi,where Di=Fxi−Fxi−1, Fx0=0, Fxn+1=1, and ∑i=1n+1Di=1.


## 4.1. Classical Methods for the Complete Sample

In this subsection, we introduce eight methods of estimation which were used for estimating the parameters of the proposed distribution.


### 4.1.1. Maximum Likelihood Estimates (MLEs)

Letx1,…,xn is a randomized sample having a size n from the PDF (3). So, the log-likelihood function for θ is as follows:(18)ℓ=nalogb+nalogβk+nloga−a+1∑i=1nlogxi+∑i=1nlogΓ1−akα,βαbxikα.The MLEs ofθ can be obtained by maximizing ℓ.


### 4.1.2. Ordinary Least-Squares Estimates (OLSEs)

Letx1:n,…,xn:n be the corresponding order statistics. The OLSEs for the distribution parameters can be easily obtained by making the following equation at minimum value, by using any mathematical software, may be MATHEMATICA 12 or any advanced program.(19)O=∑i=1nFxi:n−in+12=∑i=1ne−βαb/xi:nαk−bxi:naβa/kΓ1−akα,βαbxi:nkα−in+12.


### 4.1.3. Weighted Least-Squares Estimates (WLSEs)

By minimizing the following equation, the WLSEs of proposed model parameters can be computed:(20)W=∑i=1nn+12n+2in−i+1Fxi:n−in+12=∑i=1nn+12n+2in−i+1e−βαb/xi:nαk−bxi:naβa/kΓ1−akα,βαbxi:nkα−in+12.


### 4.1.4. Anderson–Darling Estimates (ADEs), Right-Tail Anderson–Darling Estimates (RTADEs), and Left-Tail Anderson–Darling Estimates (LTADEs)

The ADEs for the distribution parameters can be easily obtained by making the following equation at minimum value, by using any mathematical software, may be MATHEMATICA 12 or any advanced program.(21)A=−n−1n∑i=1n2i−1logFxi:n+logSxi:n.The RTADEs for the distribution parameters can be easily obtained by making the following equation at minimum value, by using any mathematical software, may be MATHEMATICA 12 or any advanced program.(22)R=n2−2∑i=1nFxi:n−1n∑i=1n2i−1logSxi:n.The LTADEs for the distribution parameters can be easily obtaine by making the following equation at minimum value, by using any mathematical software, may be MATHEMATICA 12 or any advanced program.(23)L=−32n+2∑i=1nFxi:n−1n∑i=1n2i−1logFxi:n.


### 4.1.5. Cramér–von Mises Estimates (CVMEs) and Maximum Product of Spacing Estimates (MPSEs)

The CVMEs are determined by minimizing(24)CV=112n+∑i=1nFxi:n−2i−12n2=112n+∑i=1ne−βαb/xi:nαk−bxi:naβa/kΓ1−akα,βαbxi:nkα−2i−12n2.The MPSEs are determining by maximizing the following equation:(25)G=1n+1∑i=1n+1logDi,where Di=Fxi−Fxi−1, Fx0=0, Fxn+1=1, and ∑i=1n+1Di=1.


## 4.1.1. Maximum Likelihood Estimates (MLEs)

Letx1,…,xn is a randomized sample having a size n from the PDF (3). So, the log-likelihood function for θ is as follows:(18)ℓ=nalogb+nalogβk+nloga−a+1∑i=1nlogxi+∑i=1nlogΓ1−akα,βαbxikα.The MLEs ofθ can be obtained by maximizing ℓ.


## 4.1.2. Ordinary Least-Squares Estimates (OLSEs)

Letx1:n,…,xn:n be the corresponding order statistics. The OLSEs for the distribution parameters can be easily obtained by making the following equation at minimum value, by using any mathematical software, may be MATHEMATICA 12 or any advanced program.(19)O=∑i=1nFxi:n−in+12=∑i=1ne−βαb/xi:nαk−bxi:naβa/kΓ1−akα,βαbxi:nkα−in+12.


## 4.1.3. Weighted Least-Squares Estimates (WLSEs)

By minimizing the following equation, the WLSEs of proposed model parameters can be computed:(20)W=∑i=1nn+12n+2in−i+1Fxi:n−in+12=∑i=1nn+12n+2in−i+1e−βαb/xi:nαk−bxi:naβa/kΓ1−akα,βαbxi:nkα−in+12.


## 4.1.4. Anderson–Darling Estimates (ADEs), Right-Tail Anderson–Darling Estimates (RTADEs), and Left-Tail Anderson–Darling Estimates (LTADEs)

The ADEs for the distribution parameters can be easily obtained by making the following equation at minimum value, by using any mathematical software, may be MATHEMATICA 12 or any advanced program.(21)A=−n−1n∑i=1n2i−1logFxi:n+logSxi:n.The RTADEs for the distribution parameters can be easily obtained by making the following equation at minimum value, by using any mathematical software, may be MATHEMATICA 12 or any advanced program.(22)R=n2−2∑i=1nFxi:n−1n∑i=1n2i−1logSxi:n.The LTADEs for the distribution parameters can be easily obtaine by making the following equation at minimum value, by using any mathematical software, may be MATHEMATICA 12 or any advanced program.(23)L=−32n+2∑i=1nFxi:n−1n∑i=1n2i−1logFxi:n.


## 4.1.5. Cramér–von Mises Estimates (CVMEs) and Maximum Product of Spacing Estimates (MPSEs)

The CVMEs are determined by minimizing(24)CV=112n+∑i=1nFxi:n−2i−12n2=112n+∑i=1ne−βαb/xi:nαk−bxi:naβa/kΓ1−akα,βαbxi:nkα−2i−12n2.The MPSEs are determining by maximizing the following equation:(25)G=1n+1∑i=1n+1logDi,where Di=Fxi−Fxi−1, Fx0=0, Fxn+1=1, and ∑i=1n+1Di=1.


## 5. Modeling Real Data Sets

This section discusses the flexibility of the proposed model for fitting three real-world data sets and compares it with other well-known competing models. The three analyzed data sets are used to show the flexibility of FWMED as we used very common distribution for comparison that they are known by their flexibility such as Fréchet–Weibull mixture exponential distribution (FWMED) [22], Weibull distribution (WD), exponential distribution (ED), gamma distribution (GD), and inverse Pareto distribution (IPD) [23].The competing distributions are compared using goodness-of-fit measures, including Anderson–Darling (AD), Cramér–von Mises (CM), and Kolmogorov–Smirnov (KS) with itsp value (KS-p value).To evaluate the validity of competing models, the MLEs method is used for esitimatingthe parameters of the competing models, and the analytical measurements are generated using the Wolfram Mathematica version 12 program.


### 5.1. Data Set I

This real data set is for the relief times of twenty patients taking a acertain kind of medecine called analgesic. These data were introduced by Clark and Gross [24], page 105, and they are given in Table 1.Table 1
Numerical values of the first data set.
1.11.41.31.71.91.81.62.21.72.74.11.81.51.21.431.72.31.62Table2 provides the analytical measures along with MLEs. The fitted PDF, CDF, SF, and P-P plots of the FWMPD model for the first data set are depicted in Figure 5. The results in Table 2 show that the FWMPD distribution is the best fit one compared to other models that are comparable for the first data set. Figure 6 provides profile-likelihood plots of the FWMPD parameters for the first real data set. These plots illustrate the unimodality of profile-likelihood functions for all estimated parameters. Table 3 presents the values of estimates, negative log-likelihood function, CM, AD, KS, and KSP of the proposed model for the eight different estimation methods. Figure 7 displays P-P plots for the proposed model by using different estimation methods along with fitted PDFs by results of these methods.Table 2
Goodness-of-fit measures and estimates of FWMPD and other fitted models for the first data set.
ModelADCMKSKSPEstimatesFWMPD0.1599420.0275550.1019770.985411α^=8.26658β^=1.9737k^=0.62992a^=5.07903b^=0.45048FWMED0.2731640.0343230.1107370.966914α^=1.35077a^=3.04669λ^=2.11748k^=4.3601ED4.60350.9629670.4395130.00088α^=0.526317ExED0.310460.0476590.1343140.863396α^=36.6832a^=2.23524WD1.083540.183430.184970.50057a^=2.78703b^=2.12998GD0.599020.102510.1734060.584477α^=9.66948λ^=0.19649IPD4.802070.987220.387230.00496α^=139428θ^=0.000012Figure 5
The fitted FWMPD PDF, CDF, SF, and P-P plots for the first data set. (a–d) FWMPD.
(a)(b)(c)(d)Figure 6
Plots of the profile-likelihood functions for the proposed model’s MLEs of the first data set.
(a)(b)(c)(d)(e)Table 3
The estimates and log-likelihood function of the proposed distribution parameters along with goodness-of-fit measures for the first data set by different estimation methods.
α^β^k^a^b^−LADCMKSKSPMLEs8.266581.97370.6299295.079030.45048315.39860.1599420.02755570.1019770.985411ADEs3.664310.4278791.436924.914132.3864715.40120.1602710.02768450.1000690.988169CVMEs3.99030.3904471.406095.436432.6352815.5430.1959290.02524850.09329220.99498MPSEs2.705651.17131.564854.82451.1777515.5430.1959290.02524850.09329220.99498OLSEs2.623020.4416881.960494.786271.9929615.41570.1630910.02962340.1018260.985646RTADEs2.601320.8473941.523643.422671.3902116.8960.6056710.1091020.1498260.760347WLSEs2.472070.5272181.174672.868482.0373519.81851.39160.2462370.2036420.378108LTADEs3.465081.464021.455745.033631.0209115.40740.1615810.02899180.0971620.991618Figure 7
The probability-probability (P-P) plot and the fitted PDFs of the proposed model for the first data set.
(a)(b)


### 5.2. Data Set II

This data set was taken from McCool, and it represents the fatigue lifetime in hours for 10 bearings of certain types. It was studied by Wu and Wong [25], and it is given in Table 4.Table 4
The second data set's statistical analysis.
152.7172.0172.5173.3193.0204.7216.5234.9262.6422.6Table5 provides the analytical measures along with ML estimates. The fitted PDF, CDF, SF, and P-P plots of the FWMPD model for the second data set are depicted in Figure 8. The results in Table 5 show that the FWMPD distribution is the best fit compared to other models that are comparable for the second data set. Figure 9 provides the profile-likelihood plots of the FWMPD parameters for the second real data set. These plots illustrate the unimodality of profile-likelihood functions for all estimated parameters. Table 6 presents the values of estimates, negative log-likelihood function, CM, AD, KS, and KSP of the proposed model for the eight different estimation methods. Figure 10 displays the P-P plots for the proposed model by using different estimation methods along with fitted PDFs by results of these methods.Table 5
Goodness-of-fit measures and estimates of FWMPD and other fitted models for the second data set.
ModelADCMKSKSPEstimatesFWMPD0.2031510.02816030.1420740.987635α^=22.8864β^=12.5295k^=0.773245a^=3.5865b^=5.87626FWMED2.70670.5500770.4491940.0353535α^=1.13686a^=0.000914λ^=0.107839k^=0.823594ED2.589770.5507840.4996950.0135584α^=0.004535ExED0.4515390.0528460.1850380.883325α^=88.952a^=0.023466WD0.9482930.1547080.2193760.721648a^=2.93592b^=246.409GD0.6815320.0987720.1853980.881913α^=11.5632λ^=19.0674IPD2.725480.5654040.4399310.0416838α^=114313θ^=0.0017835Figure 8
The fitted FWMPD PDF, CDF, SF, and P-P plots for the second data set. (a–d) FWMPD.
(a)(b)(c)(d)Figure 9
Plots of the profile-likelihood functions for the proposed model’s MLEs of the second data set.
(a)(b)(c)(d)(e)Table 6
The estimates and log-likelihood function of the proposed distribution parameters along with goodness-of-fit measures for the second data set by different estimation methods.
α^β^k^a^b^−LADCMKSKSPMLEs22.890112.52830.7731623.586455.8748752.70570.2031530.02816060.1420780.987631ADEs10.348240.95921.291553.623648.6746552.78030.1898190.02757190.1446040.984983CVMEs8.4469939.08441.239634.195938.0954952.86880.1993110.02693410.1467730.982409MPSEs10.924226.40121.189442.832519.4167952.86880.1993110.02693410.1467730.982409OLSEs22.890112.52830.7750233.586515.874852.72260.2062660.03043520.15965270.960749RTADEs12.742946.32211.395243.471639.8543152.71040.1989780.02817690.1433020.986394WLSEs8.1239146.45111.444753.1889310.489253.10350.229020.03449340.1666360.944053LTADEs20.961443.25790.46692974.490990.049164652.91210.208540.0270480.1491970.979187Figure 10
The probability-probability (P-P) plot and the fitted PDFs of the proposed model for the second data set.
(a)(b)


### 5.3. Data Set III

This data set indicates the survival times for head neck cancer of 45 patients; we consider this data set as a complete one. For more details about the data, see [26]. This data set is given in Table 7. Table 8 provides the analytical measures along with ML estimates. The fitted PDF, CDF, SF, and P-P plots of the FWMPD model for the third data set are depicted in Figure 11. The results in Table 8 show that the FWMPD distribution is the best fit compared to other models that are comparable for the third data set. Figure 12 provides profile-likelihood plots of the FWMPD parameters for the third real data set. These plots illustrate the unimodality of profile-likelihood functions for all estimated parameters. Table 9 presents the values of estimates, negative log-likelihood function, CM, AD, KS, and KSP of the proposed model for the eight different estimation methods. Figure 13 displays P-P plots for the proposed model by using different estimation methods along with fitted PDFs by results of these methods.Table 7
Numerical values of the third data set.
3784929411011211912713013314014615515916917317919419520924928131933943246951952854761363372575981710921245133115571642177117761897202321462297Table 8
Goodness-of-fit measures and estimates of FWMPD and other fitted models for the third data set.
ModelADCMKSKSPEstimatesFWMPD0.6393630.0846630.0966990.794056α^=10.7787β^=6.99619k^=0.123512a^=1.32604b^=0.000016FWMED1.014880.1521550.1312860.419929α^=1.55186a^=0.00109λ^=3.91053k^=0.976671ED1.482530.2411530.1655970.169403α^=0.001564ExED1.583440.2624420.1730830.134882α^=1.03827a^=0.001603WD1.389130.2202130.1578020.212421a^=0.979504b^=632.973GD1.557820.2571970.1712390.142806α^=1.02661λ^=622.634IPD0.651750.08696020.0984140.776073α^=45.269θ^=4.91826Figure 11
The fitted FWMPD PDF, CDF, SF, and P-P plots for the third data set. (a–d) FWMPD.
(a)(b)(c)(d)Figure 12
Plots of the profile-likelihood functions for the proposed model MLEs of the third data set.
(a)(b)(c)(d)(e)Table 9
The estimates and log-likelihood function of the proposed distribution parameters along with goodness-of-fit measures for the third data set by different estimation methods.
α^β^k^a^b^−LADCMKSKSPMLEs10.77876.996190.1235121.326040.00001628333.5080.6706160.092083570.09119920.848307ADEs3.771385.374460.374461.070621.14398333.8240.6351780.0791010.09533210.808067CVMEs6.32314.022130.2743730.823390.546817334.7570.6772780.07517010.1112280.633673MPSEs3.039654.291190.4129761.218813.029334.7570.6772780.07517010.1112280.633673OLSEs4.870734.121840.3581810.7843521.60477334.9770.714750.07659750.117670.561555RTADEs4.346445.328080.296791.261540.40316333.5890.6376970.08328880.09701360.790793WLSEs3.7713875.374460.374461.070621.14398333.8240.6351780.079100970.09533160.808072LTADEs2.9938112.56280.573330.8656171.09981334.5780.6576120.07639290.1051590.702176Figure 13
The probability-probability (P-P) plot and the fitted PDFs of the proposed model for the third data set.
(a)(b)


### 5.4. Concluding Remarks on the Results of the Real Data Sets

(1)
Regarding the data sets in Tables1, 4, and, 7, we applied three data sets to the proposed distribution, and we deduced that the distribution outperforms all its competitors.(2)
Referring to the values of the KS, AD, and CM, we can deduce that the proposed distribution has the least measures, and this assures its superiority.(3)
Referring to theP values of the distribution, we can deduce that the proposed distribution has the highest value, and this assures its superiority.(4)
Figures6, 9, and 12 provide profile-likelihood plots of the FWMPD parameters for the three real data sets, respectively. These plots illustrate that the estimated parameters give a maximum value of the log-likelihood function, and these estimates are global maximum estimates.


## 5.1. Data Set I

This real data set is for the relief times of twenty patients taking a acertain kind of medecine called analgesic. These data were introduced by Clark and Gross [24], page 105, and they are given in Table 1.Table 1
Numerical values of the first data set.
1.11.41.31.71.91.81.62.21.72.74.11.81.51.21.431.72.31.62Table2 provides the analytical measures along with MLEs. The fitted PDF, CDF, SF, and P-P plots of the FWMPD model for the first data set are depicted in Figure 5. The results in Table 2 show that the FWMPD distribution is the best fit one compared to other models that are comparable for the first data set. Figure 6 provides profile-likelihood plots of the FWMPD parameters for the first real data set. These plots illustrate the unimodality of profile-likelihood functions for all estimated parameters. Table 3 presents the values of estimates, negative log-likelihood function, CM, AD, KS, and KSP of the proposed model for the eight different estimation methods. Figure 7 displays P-P plots for the proposed model by using different estimation methods along with fitted PDFs by results of these methods.Table 2
Goodness-of-fit measures and estimates of FWMPD and other fitted models for the first data set.
ModelADCMKSKSPEstimatesFWMPD0.1599420.0275550.1019770.985411α^=8.26658β^=1.9737k^=0.62992a^=5.07903b^=0.45048FWMED0.2731640.0343230.1107370.966914α^=1.35077a^=3.04669λ^=2.11748k^=4.3601ED4.60350.9629670.4395130.00088α^=0.526317ExED0.310460.0476590.1343140.863396α^=36.6832a^=2.23524WD1.083540.183430.184970.50057a^=2.78703b^=2.12998GD0.599020.102510.1734060.584477α^=9.66948λ^=0.19649IPD4.802070.987220.387230.00496α^=139428θ^=0.000012Figure 5
The fitted FWMPD PDF, CDF, SF, and P-P plots for the first data set. (a–d) FWMPD.
(a)(b)(c)(d)Figure 6
Plots of the profile-likelihood functions for the proposed model’s MLEs of the first data set.
(a)(b)(c)(d)(e)Table 3
The estimates and log-likelihood function of the proposed distribution parameters along with goodness-of-fit measures for the first data set by different estimation methods.
α^β^k^a^b^−LADCMKSKSPMLEs8.266581.97370.6299295.079030.45048315.39860.1599420.02755570.1019770.985411ADEs3.664310.4278791.436924.914132.3864715.40120.1602710.02768450.1000690.988169CVMEs3.99030.3904471.406095.436432.6352815.5430.1959290.02524850.09329220.99498MPSEs2.705651.17131.564854.82451.1777515.5430.1959290.02524850.09329220.99498OLSEs2.623020.4416881.960494.786271.9929615.41570.1630910.02962340.1018260.985646RTADEs2.601320.8473941.523643.422671.3902116.8960.6056710.1091020.1498260.760347WLSEs2.472070.5272181.174672.868482.0373519.81851.39160.2462370.2036420.378108LTADEs3.465081.464021.455745.033631.0209115.40740.1615810.02899180.0971620.991618Figure 7
The probability-probability (P-P) plot and the fitted PDFs of the proposed model for the first data set.
(a)(b)


## 5.2. Data Set II

This data set was taken from McCool, and it represents the fatigue lifetime in hours for 10 bearings of certain types. It was studied by Wu and Wong [25], and it is given in Table 4.Table 4
The second data set's statistical analysis.
152.7172.0172.5173.3193.0204.7216.5234.9262.6422.6Table5 provides the analytical measures along with ML estimates. The fitted PDF, CDF, SF, and P-P plots of the FWMPD model for the second data set are depicted in Figure 8. The results in Table 5 show that the FWMPD distribution is the best fit compared to other models that are comparable for the second data set. Figure 9 provides the profile-likelihood plots of the FWMPD parameters for the second real data set. These plots illustrate the unimodality of profile-likelihood functions for all estimated parameters. Table 6 presents the values of estimates, negative log-likelihood function, CM, AD, KS, and KSP of the proposed model for the eight different estimation methods. Figure 10 displays the P-P plots for the proposed model by using different estimation methods along with fitted PDFs by results of these methods.Table 5
Goodness-of-fit measures and estimates of FWMPD and other fitted models for the second data set.
ModelADCMKSKSPEstimatesFWMPD0.2031510.02816030.1420740.987635α^=22.8864β^=12.5295k^=0.773245a^=3.5865b^=5.87626FWMED2.70670.5500770.4491940.0353535α^=1.13686a^=0.000914λ^=0.107839k^=0.823594ED2.589770.5507840.4996950.0135584α^=0.004535ExED0.4515390.0528460.1850380.883325α^=88.952a^=0.023466WD0.9482930.1547080.2193760.721648a^=2.93592b^=246.409GD0.6815320.0987720.1853980.881913α^=11.5632λ^=19.0674IPD2.725480.5654040.4399310.0416838α^=114313θ^=0.0017835Figure 8
The fitted FWMPD PDF, CDF, SF, and P-P plots for the second data set. (a–d) FWMPD.
(a)(b)(c)(d)Figure 9
Plots of the profile-likelihood functions for the proposed model’s MLEs of the second data set.
(a)(b)(c)(d)(e)Table 6
The estimates and log-likelihood function of the proposed distribution parameters along with goodness-of-fit measures for the second data set by different estimation methods.
α^β^k^a^b^−LADCMKSKSPMLEs22.890112.52830.7731623.586455.8748752.70570.2031530.02816060.1420780.987631ADEs10.348240.95921.291553.623648.6746552.78030.1898190.02757190.1446040.984983CVMEs8.4469939.08441.239634.195938.0954952.86880.1993110.02693410.1467730.982409MPSEs10.924226.40121.189442.832519.4167952.86880.1993110.02693410.1467730.982409OLSEs22.890112.52830.7750233.586515.874852.72260.2062660.03043520.15965270.960749RTADEs12.742946.32211.395243.471639.8543152.71040.1989780.02817690.1433020.986394WLSEs8.1239146.45111.444753.1889310.489253.10350.229020.03449340.1666360.944053LTADEs20.961443.25790.46692974.490990.049164652.91210.208540.0270480.1491970.979187Figure 10
The probability-probability (P-P) plot and the fitted PDFs of the proposed model for the second data set.
(a)(b)


## 5.3. Data Set III

This data set indicates the survival times for head neck cancer of 45 patients; we consider this data set as a complete one. For more details about the data, see [26]. This data set is given in Table 7. Table 8 provides the analytical measures along with ML estimates. The fitted PDF, CDF, SF, and P-P plots of the FWMPD model for the third data set are depicted in Figure 11. The results in Table 8 show that the FWMPD distribution is the best fit compared to other models that are comparable for the third data set. Figure 12 provides profile-likelihood plots of the FWMPD parameters for the third real data set. These plots illustrate the unimodality of profile-likelihood functions for all estimated parameters. Table 9 presents the values of estimates, negative log-likelihood function, CM, AD, KS, and KSP of the proposed model for the eight different estimation methods. Figure 13 displays P-P plots for the proposed model by using different estimation methods along with fitted PDFs by results of these methods.Table 7
Numerical values of the third data set.
3784929411011211912713013314014615515916917317919419520924928131933943246951952854761363372575981710921245133115571642177117761897202321462297Table 8
Goodness-of-fit measures and estimates of FWMPD and other fitted models for the third data set.
ModelADCMKSKSPEstimatesFWMPD0.6393630.0846630.0966990.794056α^=10.7787β^=6.99619k^=0.123512a^=1.32604b^=0.000016FWMED1.014880.1521550.1312860.419929α^=1.55186a^=0.00109λ^=3.91053k^=0.976671ED1.482530.2411530.1655970.169403α^=0.001564ExED1.583440.2624420.1730830.134882α^=1.03827a^=0.001603WD1.389130.2202130.1578020.212421a^=0.979504b^=632.973GD1.557820.2571970.1712390.142806α^=1.02661λ^=622.634IPD0.651750.08696020.0984140.776073α^=45.269θ^=4.91826Figure 11
The fitted FWMPD PDF, CDF, SF, and P-P plots for the third data set. (a–d) FWMPD.
(a)(b)(c)(d)Figure 12
Plots of the profile-likelihood functions for the proposed model MLEs of the third data set.
(a)(b)(c)(d)(e)Table 9
The estimates and log-likelihood function of the proposed distribution parameters along with goodness-of-fit measures for the third data set by different estimation methods.
α^β^k^a^b^−LADCMKSKSPMLEs10.77876.996190.1235121.326040.00001628333.5080.6706160.092083570.09119920.848307ADEs3.771385.374460.374461.070621.14398333.8240.6351780.0791010.09533210.808067CVMEs6.32314.022130.2743730.823390.546817334.7570.6772780.07517010.1112280.633673MPSEs3.039654.291190.4129761.218813.029334.7570.6772780.07517010.1112280.633673OLSEs4.870734.121840.3581810.7843521.60477334.9770.714750.07659750.117670.561555RTADEs4.346445.328080.296791.261540.40316333.5890.6376970.08328880.09701360.790793WLSEs3.7713875.374460.374461.070621.14398333.8240.6351780.079100970.09533160.808072LTADEs2.9938112.56280.573330.8656171.09981334.5780.6576120.07639290.1051590.702176Figure 13
The probability-probability (P-P) plot and the fitted PDFs of the proposed model for the third data set.
(a)(b)


## 5.4. Concluding Remarks on the Results of the Real Data Sets

(1)
Regarding the data sets in Tables1, 4, and, 7, we applied three data sets to the proposed distribution, and we deduced that the distribution outperforms all its competitors.(2)
Referring to the values of the KS, AD, and CM, we can deduce that the proposed distribution has the least measures, and this assures its superiority.(3)
Referring to theP values of the distribution, we can deduce that the proposed distribution has the highest value, and this assures its superiority.(4)
Figures6, 9, and 12 provide profile-likelihood plots of the FWMPD parameters for the three real data sets, respectively. These plots illustrate that the estimated parameters give a maximum value of the log-likelihood function, and these estimates are global maximum estimates.


## 6. Conclusion and Major Findings

In this article, we introduced a new mixture of distribution FWMPD, and we estimated its parameters by the classical methods of estimation: the maximum likelihood estimation and 7 other methods. We introduced its mathematical properties and graphed its PDF and CDF to study its behavior under different values of estimates. Last but not least, we made an application on the proposed distribution to assure its superiority compared to its competitors. We evaluated its KS, AD, CM, andP value, and we deduced that it has the lowest values for KS, AD, and CM and the greatest values for P values which make it a better candidate among all its competitors. Also, to make sure that the roots for the MLEs for the proposed distribution give a maximum value, we graphed Figures 6, 9, and 12 for the profile-likelihood function of the proposed model with its parameters for the three real data sets, respectively. These plots illustrate the unimodality of profile-likelihood functions for all estimated parameters. We expect that the presented model will find a broader range of applications in fields like engineering, survival and lifespan data, meteorology, hydrology, and economics.



---

*Source: 9819200-2021-08-17.xml*