Next Article in Journal
Optimization Strategy for Modal Test Measurement Points of Large-Span Steel Beams Based on Improved Particle Swarm Optimization Algorithm with Random Weights
Next Article in Special Issue
Characterisation of Physiological Responses to Odours in Autism Spectrum Disorders: A Preliminary Study
Previous Article in Journal
In-Depth Evaluation of Association between Crash and Hand Arthritis via Naturalistic Driving Study
Previous Article in Special Issue
A Genetic Programming Approach for Economic Forecasting with Survey Expectations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Power Zeghdoudi Distribution: Properties, Estimation, and Applications to Real Right-Censored Data

1
Laboratory of Probability and Statistics, University of Badji Mokhtar, Annaba 23000, Algeria
2
Department of Mathematics, Faculty of Science, Al al-Bayt University, Mafraq 25113, Jordan
3
Department of Mathematical Sciences, Umm Al-Qura University, Makkah 24382, Saudi Arabia
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(23), 12081; https://doi.org/10.3390/app122312081
Submission received: 9 October 2022 / Revised: 19 November 2022 / Accepted: 21 November 2022 / Published: 25 November 2022
(This article belongs to the Special Issue Applied Biostatistics & Statistical Computing)

Abstract

:
A new two-parameter power Zeghdoudi distribution (PZD) is suggested as a modification of the Zeghdoudi distribution using the power transformation method. As a result, the PZD may have increasing, decreasing, and unimodal probability density function and decreasing mean residual life function. In addition, other properties are presented, such as moments, order statistics, reliability measures, Bonferroni and Lorenz curves, Gini index, stochastic ordering, mean and median deviations, and quantile function. Following this, a section is devoted to the related model parameters which are estimated using the maximum likelihood estimation method, the weighted least squares and least squares methods, the maximum product of spacing method, the Cramer–von Mises method, and the right-tail and left-tail Anderson–Darling methods, and the Nikulin–Rao–Robson test statistic is considered. A simulation study is conducted to assess these methods and to investigate the distribution properties with right-censored data. The applicability of the proposed model is studied based on three real data sets of failure times, bladder cancer patients, and glass fiber data with a comparison with such competitors as the gamma, xgamma, Lomax, Darna, power Darna, power Lindley, and exponentiated power Lindley models. According to several established criteria, the comparative findings are overwhelmingly favorable to the suggested model.

1. Introduction

Ref. [1] investigated the characteristics of the Lindley distribution (LD) in great depth. According to [2], there are some circumstances in which the LD may not be appropriate from a theoretical or applied point of view. Due to this, ref. [2] introduced the power Lindley (PL) distribution, a more adaptable distribution, using a power transformation method. Ref. [2] showed that based on the values of the parameters, the hazard function of the PL distribution can rise, decrease, or decrease–increase–decrease. Some of the other distributions that have been suggested utilizing the power transformation method are the power Shanker distribution by [3], the power Ishita and Aradhana distributions by [4], and the power Rama and Garima distributions by [5,6]. The power Ishita distribution is proposed by [7], the power Sujatha distribution is considered by [8], and the power Prakaamy distribution is proposed by [9]. In addition, ref. [10] suggested the power half-logistic distribution, ref. [11] introduced the power binomial exponential distribution, ref. [12] proposed the power length-biased Suja distribution, ref. [13] discussed a novel three-parameter distribution known as the power Darna distribution (PDD), ref. [14] offered the power XLindley distribution, and ref. [15] proposed the power-modified Kies-exponential distribution. All of the power-transformed distributions are shown to be more adaptable than their respective baseline distributions, and more effective for analyzing complicated data structures in various spheres of life.
Ref. [16] provided a novel lifetime distribution with one parameter for modeling lifetime data from engineering and medical science, known as Zeghdοudi distribution (ZD), with a cumulative distribution function (cdf) given by
F ( y ; θ ) = 1 ( 1 + θ 2 y 2 + θ ( θ + 2 ) y 2 + θ ) e θ y , y > 0 , θ > 0 ,
and a survival function defined as
S ( y ; θ ) = 1 F ( y ; θ ) = ( 1 + θ 2 y 2 + θ ( θ + 2 ) y 2 + θ ) e θ y .
The probability density function (pdf) of the ZD is
f ( y ; θ ) = θ 3 2 + θ y ( 1 + y ) e θ y , y > 0 , θ > 0 ,
and its hazard function has the form
h ( y ; θ ) = f ( y , θ ) 1 F ( y , θ ) = θ 3 y ( 1 + y ) 2 + θ + θ 2 y 2 + θ ( θ + 2 ) y .
The ZD has the mean E ( Y ) = 2 ( θ + 3 ) θ ( θ + 2 ) , variance V a r ( Y ) = 2 ( θ 2 + 6 θ + 6 ) θ 2 ( 2 + θ ) 2 , coefficient of skewness S K = 24 ( θ + 5 ) ( θ + 2 ) 2 ( 2 θ 2 + 12 θ + 12 ) 3 2 , and coefficient of kurtosis K u = 30 ( θ + 6 ) ( θ + 2 ) 3 ( θ 2 + 6 θ + 6 ) 6 .
According to [16], the maximum likelihood (MLE) and method of moment (MOM) estimators of the distribution parameter have the same forms as
θ ^ MOM = θ MLE = 1 y ( 1 y ¯ + y ¯ 2 + 4 y ¯ + 1 )  
where y ¯ = 1 n i = 1 n y i .
The accuracy of the parametric distributions is tested using several statistical methods. Numerous fresh initiatives are investigated to address these problems. Based on the well-known Kaplan and Meyer estimators, refs. [17,18] recommended a new form of the Chi-squared test for randomly censored actual data. Then, for the accelerated loss distributions, ref. [19] assumed certain nonparametric adjustments for the Kolmogorov-Smirnov statistic (KSS), Anderson–Darling statistic (ADS), and Cramer–von Mises statistic (KVMS). A different Chi-squared GOF test statistic is recently presented and demonstrated by [20] for the correct censored data. For distributional validation using right-censored approaches, Bagdonavicius and Nikulin’s (B-N) is the most current Chi-squared GOF test statistic in use.
The B-N GOF test and a new modified Chi-square GOF test are both utilized in this research to provide a censored maximum likelihood estimation for the right-censored validation. The results are checked for accuracy using a modified GOF statistical test on the right-censored actual data set. A modified GOF test that is based on the censored maximum likelihood estimators on the initial data is able to recover the information loss, whereas the grouped data have a Chi-square distribution. The modified criteria tests’ components are attained. Finally, the filtered technique is demonstrated with real-world data set applications for validation.
To the best of our knowledge, this paper is the first to suggest the power Zeghdoudi distribution, which offers more precision and flexibility for fitting actual data than some of the competitors considered in this study. The following sections make up the article: The density and cumulative distribution functions of the PZD are provided in Section 2 with graphical representations. The statistical properties of the PZD and related measures are described in Section 3. Various estimation techniques are presented in Section 4 and some simulations are given in Section 5. Applications of three actual real data sets in engineering and medical sciences are provided in Section 6. Closing observations and recommendations for further works are formulated in Section 7.

2. The Power Zeghdοudi Distribution

In this section, based on the same power transformation considered by [2], X = Y 1 β , a random variable X is said to have a power Zeghdοudi distribution with parameters   θ   and   β   if its survival function is given by
S P Z D ( x ; θ , β ) = P ( X > x ) = 1 F P Z D ( x ; θ , β ) = e x β θ ( 2 + θ + x β θ ( 2 + θ + x β θ ) ) 2 + θ , x > 0 , θ > 0 , β > 0
Figure 1 involves some plots of the survival function of the PZD for some distribution parameters. We can note that the survival function of the PZD has many shapes, such as decreasing–constant ( θ = 10 ,   β = 2 ) and decreasing ( θ = 2.5 ,   β = 0.9 ), which are heavy-tailed to the right.
Hence, from (3) the corresponding cdf of PZD is given by
F P Z D ( x ; θ , β ) = 1 ( 1 + θ x β ( θ + 2 + θ x β ) θ + 2 ) e   θ x β , x > 0 , θ > 0 , β > 0 .
It is clear that for β = 1 , the PZD reduces to the ZD. Plots of the PZD cdf for some selections of   β   and   θ   are presented in Figure 2. It is noted that the cdf shapes of the PZD vary with the distribution parameter values.
The pdf of the PZD can be obtained by taking the derivative of F P Z D ( x ; θ , β ) as
f P Z D ( x ; β , θ ) = β θ 3 2 + θ x 2 β 1 ( 1 + x β ) e θ x β , x > 0 , θ > 0 , β > 0 .
Different PZD pdf plots are given in Figure 3 with various parameters. The suggested distribution pdf is asymmetric, unimodal, semi-symmetric ( θ = 8 ,   β = 5 ), increasing–decreasing ( θ = 16 ,   β = 1.5 ), skewed to the left ( θ = 0.2 ,   β = 5 ), and skewed to the right ( θ = 5 ,   β = 0.8 ), according to the parameter’s values.

The Mode PZD

The logarithm of the PZD pdf given in Equation (4) is
l o g [ f P Z D ( x ; β , θ ) ] = l o g ( β θ 3 2 + θ ) + ( 2 β 1 ) log ( x ) + log ( 1 + x β ) θ x β .
The first derivative of this equation with respect to x is
l o g [ f P Z D ( x ; β , θ ) ] x = 2 β + x β [ ( β ( θ + θ x β 3 ) + 1 ) ] 1 x ( x β + 1 )
Let l o g [ f P Z D ( x ; β , θ ) ] x = 0 to obtain the mode of the distribution as
x 0 = ( 1 + 3 β β θ + 4 ( 1 2 β ) β θ + ( 1 3 β + β θ ) 2 2 β θ ) 1 β .

3. Statistical Properties

In this section, some of the main statistical properties of the PZD are derived and some simulations are provided.

3.1. Moments

The rth moment about the origin of the PZD is defined as
E ( X P Z D r ) = β ( θ + 2 ) + r β ( θ + 2 ) θ r β Γ ( r β + 2 ) , r = 1 , 2 , 3 , ...
and hence the mean of the distribution is
E ( X P Z D ) = β ( θ + 2 ) + 1 β ( θ + 2 ) θ 1 β Γ ( 1 β + 2 )
The inverse moment of the PZD is
E ( 1 X P Z D r ) = β ( θ + 2 ) r β ( θ + 2 ) θ r / β Γ ( 2 r β ) , r = 1 , 2 , 3 , ...
While kurtosis of a random variable measures the degree of tail heaviness, skewness of a random variable measures the degree of long-tail. The coefficient of skewness ( ƛ 1 ) , coefficient of variation ( ƛ 2 ) , and coefficient of kurtosis ( ƛ 3 ) of the PZD are given, respectively, by
ƛ 1 = β 2 ( θ + 2 ) 2 [ β ( θ + 2 ) + 3 ] Γ ( 2 + 3 β ) + 2 [ β ( θ + 2 ) + 1 ] 3 Γ ( 2 + 1 β ) 3 3 β ( θ + 2 ) [ β ( θ + 2 ) + 1 ] [ β ( θ + 2 ) + 2 ] Γ ( 2 + 2 β ) Γ ( 2 + 1 β ) ( β ( θ + 2 ) [ β ( θ + 2 ) + 2 ] Γ ( 2 + 2 β ) [ β ( θ + 2 ) + 1 ] 2 Γ ( 2 + 1 β ) 2 ) 3 / 2
ƛ 2 = β ( θ + 2 ) [ β ( θ + 2 ) + 2 ] Γ ( 2 + 2 β ) [ β ( θ + 2 ) + 1 ] 2 Γ ( 2 + 1 β ) 2 [ β ( θ + 2 ) + 1 ] Γ ( 2 + 1 β )
ƛ 3 = A 1 A 2 A 3 ( [ β ( θ + 2 ) + 1 ] 2 Γ ( 2 + 1 β ) 2 β ( θ + 2 ) [ β ( θ + 2 ) + 2 ] Γ ( 2 + 2 β ) ) 2
where
A 1 = β 3 ( θ + 2 ) 3 [ β ( θ + 2 ) + 4 ] Γ ( 2 + 4 β ) ,
A 2 = 4 β 2 ( θ + 2 ) 2 [ β ( θ + 2 ) + 1 ] [ β ( θ + 2 ) + 3 ] Γ ( 2 + 3 β ) Γ ( 2 + 1 β ) ,
A 3 = 3 [ β ( θ + 2 ) + 1 ] 2 Γ ( 2 + 1 β ) 2 { [ β ( θ + 2 ) + 1 ] 2 Γ ( 2 + 1 β ) 2 + 2 β ( θ + 2 ) [ β ( θ + 2 ) + 2 ] Γ ( 2 + 2 β ) } .
Table 1 gives the mean ( μ X ) , standard deviation ( σ X ) , ƛ 1 , ƛ 2 , and ƛ 3 of the PZD for various selections of the parameters.
The results in Table 1 indicate that the distribution is skewed to the right for these parameter value selections, with relatively small values of skewness. In addition, it is noted that the mean and variance values are decreasing as one parameter’s values are increasing for a fixed value of the other parameter.

3.2. Order Statistics

Let X 1 , X 2 , ... , X m be a random sample selected from PZD with cdf (4) and pdf (5) and X ( 1 ) ,   X ( 2 ) , ... , X ( m ) be its order statistics; then, the pdf of the ith ( i = 1 , 2 , ... , m ) order statistic, X ( i ) , is given by
f ( i ) ( x ; θ , β ) = m ! ( i 1 ) ! ( m i ) ! [ F ( x ; θ , β ) ] i 1 [ S ( x ; θ , β ) ] m i f ( x ; θ , β ) = β θ 3 m ! x 2 β 1 ( x β + 1 ) e θ i x β [ θ x β ( θ + θ x β + 2 ) + θ + 2 θ + 2 ] i [ θ + 2 e θ x β ( θ x β ( θ + θ x β + 2 ) + θ + 2 ) ] m i Γ ( i ) Γ ( 1 i + m ) ( θ + 2 ) m [ θ + θ x β ( θ + θ x β + 2 ) + 2 ]
For i = 1 and m in (12), the pdfs of the smallest and largest order statistics, respectively, are given by
f ( 1 ) ( x ; θ , β ) = β θ 3 m x 2 β 1 ( x β + 1 ) e θ x β { θ + 2 [ θ x β ( θ + θ x β + 2 ) + θ + 2 ] e θ x β } m 1 ( θ + 2 ) m ,
and
f ( m ) ( x ; θ , β ) = β θ 3 m x 2 β 1 ( x β + 1 ) e m θ x β [ θ x β ( θ + θ x β + 2 ) + θ + 2 ] m 1 ( θ + 2 ) m .
Figure 4 represents the three-dimensional plots of the minimum and maximum pdfs with m = 5 , β = 2 and 1 θ 4 .

3.3. Reliability Measures

The reliability and hazard functions of any lifetime equipment are the fundamental tools for analyzing aging and associated features. The conditional probability of failure, assuming it has lasted up to time x 0 , is represented by the hazard function. The hazard, reversed hazard rate, and odds functions of the PZD are given, respectively, by
H P Z D ( x ; θ , β ) = f P Z D ( x ; θ , β ) S P Z D ( x ; θ , β ) = β θ 3 x 2 β 1 ( x β + 1 ) θ + θ x β ( θ + θ x β + 2 ) + 2 ,
R H P Z D ( x ; θ , β ) = f P Z D ( x ; θ , β ) F P Z D ( x ; θ , β ) = β θ 3 x 2 β 1 ( x β + 1 ) ( θ + 2 ) e θ x β ( θ + 2 + θ x β ( θ + θ x β + 2 ) ) ,
O P Z D ( x ; θ , β ) = F P Z D ( x ; θ , β ) S P Z D ( x ; θ , β ) = e x β θ ( 2 + θ ) 2 + θ + x β θ ( 2 + θ + x β θ ) 1 .
Plots of the hazard, reversed hazard rate, and odds functions of the PZD are illustrated in Figure 5, Figure 6, and Figure 7, respectively.
Mean residual life (MRL) is crucial in demography, reliability, and survival analysis. The mean residual life function of PZD is defined as
M R L P Z D ( x ; θ , β ) = E ( X x | X > x ) = 1 1 F P Z D ( x ; θ , β ) x ( 1 F P Z D ( t ; θ , β ) ) d t = e θ x β [ ( ( θ + 2 ) Γ ( 1 + 1 β , x β θ ) + Γ ( 2 + 1 β , x β θ ) + ( θ + 2 ) Γ ( 1 β , x β θ ) ) ] β θ 1 / β [ θ + θ x β ( θ + θ x β + 2 ) + 2 ] .
MRL plots of the PZD for some particular parameters are offered in Figure 8. This figure reveals that the MRL for PZD is decreasing for these values with different behavior decreasing–constant ( θ = 8 ,   β = 0.6 ) and decreasing ( θ = 2 ,   β = 3 ) .
Remark 1. 
It is to be noted that M R L P Z D ( 0 ; θ , β ) = μ = β ( θ + 2 ) + 1 β ( θ + 2 ) θ 1 / β Γ ( 2 + 1 β ) .

3.4. Mean and Median Deviations

The mean deviations from the mean and median, which are both essential measurements of the scatter in the population, are defined as
𝒦 μ = 0 μ ( μ x ) f ( x ) d x + μ ( x μ ) f ( x ) d x = 2 μ F ( μ ) 2 0 μ x f ( x ) d x = 2 μ β ( θ + 2 ) 2 ( ( θ μ ) 1 / β [ θ Γ ( 2 + 1 β , θ μ β ) + Γ ( 3 + 1 β , θ μ β ) θ Γ ( 2 + 1 β ) Γ ( 3 + 1 β ) ] [ θ ( θ μ β + θ + 2 ) μ β + θ + 2 ] e θ μ β + θ + 2 )
and
𝒦 m = 0 m ( m x ) f ( x ) d x + m ( x m ) f ( x ) d x = μ 2 0 m x f ( x ) d x = 2 β [ θ Γ ( 2 + 1 β , m β θ ) + Γ ( 3 + 1 β , m β θ ) ] [ β ( θ + 2 ) + 1 ] Γ ( 2 + 1 β ) β ( θ + 2 ) θ 1 / β ,
where μ = β ( θ + 2 ) + 1 β ( θ + 2 ) θ 1 / β Γ ( 2 + 1 β ) and m is the population median.

3.5. Quantile Function and Stochastic Ordering

In this section, the quantile function and stochastic ordering are presented and discussed.

3.5.1. Quantile Function

Let X f P Z D ( x ; β , θ ) given in (4). The quantile function, say Q ( p ) , 0 < p < 1 defined by F ( Q ( p ) ) = p , is the root of the equation
ln ( θ + 2 + θ Q ( p ) β ( θ + 2 + θ Q ( p ) β ) ) θ Q ( p ) β = ln ( 1 p ) + ln ( θ + 2 )
Table 2 presents some quantile values for some choices of p where the real and positive roots are assigned. It can be noted that the Q ( p ) values are increasing in p.

3.5.2. Stochastic Ordering

To define the stochastic ordering, a random variable X is smaller than a random variable Y in the following:
(1)
Mean residual life order ( X M R L O Y ) , if m X ( t ) m Y ( t ) for all t.
(2)
Likelihood ratio order ( X L R O Y ) , if f X ( t ) f Y ( t ) decreases for all t.
(3)
Hazard rate order ( X H R O Y ) , if h X ( t ) h Y ( t ) for all t.
(4)
Stochastic order ( X S O Y ) , if F X ( t ) F Y ( t ) for all t.
Remark 2. 
X L R O Y X H R O X S O Y Y X M R L O Y .
Theorem 1. 
Let X P Z D f X P Z D ( t ; β , θ ) and Y P Z D f Y P Z D ( t ; α , υ ) be two random variables. If ( β > α , θ , υ ), or ( β = α with θ > υ ), then X L R O Y , and thus X H R O Y , X M R L O Y and X S O Y .
Proof. 
Let X P Z D f X P Z D ( x ; β , θ ) and Y P Z D f Y P Z D ( x ; α , υ ) . Therefore,
f X P Z D ( t ; β , θ ) f Y P Z D ( t ; α , υ ) = β θ 3 2 + θ t 2 β 1 ( 1 + t β ) e θ t β α υ 3 2 + υ t 2 α 1 ( 1 + t α ) e υ t α . = β θ 3 ( 2 + υ ) α υ 3 ( 2 + θ ) t 2 β 1 ( 1 + t β ) t 2 α 1 ( 1 + t α ) e ( θ t β υ t α ) .
For simplification, we have
ln f X P Z D ( t ; β , θ ) f Y P Z D ( t ; α , υ ) = ln ( β θ 3 ( 2 + υ ) α υ 3 ( 2 + θ ) ) + ln ( t 2 β 2 α ) + ln ( 1 + t β 1 + t α ) ( θ t β υ t α )
Thus,
t ln f X P Z D ( t ; β , θ ) f Y P Z D ( t ; α , υ ) = β ( 3 1 1 + t β t β θ ) + α ( 3 + 1 1 + t α + t α υ ) t
To this end, if ( β > α , θ , υ ), or ( β = α with θ > υ ), then t ln ( f X P Z D ( t ; β , θ ) f Y P Z D ( t ; α , υ ) ) < 0 . This shows that X L R O Y , and according to Remark 2, the theorem is proved. □

3.6. Some Curves and Gini Index

Assume that X is a continuous random variable with a non-negative, twice-differentiable cumulative distribution function. The Bonferroni curve of the random variable X is defined as
B ( p ) = 1 p μ   ( μ q x f ( x ) d x ) = [ β ( θ + 2 ) + 1 ] Γ ( 2 + 1 β ) β [ θ Γ ( 2 + 1 β , q β θ ) + Γ ( 3 + 1 β , q β θ ) ] p [ β ( θ + 2 ) + 1 ] Γ ( 2 + 1 β ) ,  
  q = F 1 ( p ) and p ( 0 , 1 ] , and the Lorenz curve is given as
L ( p ) = 1 μ   ( μ q x f ( x ) d x ) = 1 β [ θ Γ ( 2 + 1 β , q β θ ) + Γ ( 3 + 1 β , q β θ ) ] [ β ( θ + 2 ) + 1 ] Γ ( 2 + 1 β ) ,   q = F 1 ( p ) .
The Gini index of the PZD can be expressed as
G I = 1 1 μ   0 ( 1 F ( x ) ) 2 d x = 1 2 1 β 4 ( β + 1 ) { β [ 2 β ( 8 β ( θ + 2 ) 2 + 2 θ ( θ + 8 ) + 27 ) + 4 θ + 13 ] + 1 } Γ ( 1 β ) β 4 ( θ + 2 ) [ β ( θ + 2 ) + 1 ] Γ ( 2 + 1 β )

4. Methods of Estimation and Test

As is well known, the properties of maximum likelihood estimators are not always verified for small samples. Therefore, in recent years, classical and new estimation methods have been developed. One purpose of this work is to investigate different methods to estimate the unknown distribution parameters of the new model, such as maximum likelihood estimation (MLE), ordinary least square (LS), weighted least square (WLS), maximum product of spacing (MPS), Cramer–von Mises estimation (CME), and Anderson–Darling estimation (ADE) methods based on the empirical function distribution. In addition, the Nikulin–Rao–Robson test is considered in this section.

4.1. The MLE

Because of their properties, the maximum likelihood estimators are preferred in providing the values of the unknown parameters. Considering x 1 , x 2 , ... x n to be a random sample following the PZD with parameters   Ψ = ( θ , β ) , the log-likelihood function becomes
log L = n ln ( β ) + 3 n ln ( θ ) n ln ( 2 + θ ) + ( 2 β 1 ) i = 1 n ln ( x i ) + i = 1 n ln ( 1 + x i β ) θ i = 1 n x i β
The maximum likelihood estimators θ ^ and β ^ of the unknown parameters θ and β , respectively, can be obtained from the following nonlinear score equations as
log L θ = 0 3 n θ n 2 + θ i = 1 n x i 2 = 0 .
Hence,
θ ^ ( β ^ ) = 4 n ( i = 1 n x i β ^ ) + ( i = 1 n x i β ^ ) 2 + n 2 i = 1 n x i β ^ + n i = 1 n x i β ^ .
In addition,
log L β = 0   n β + 2 i = 1 n ln ( x i ) + i = 1 n x i β ln ( x i ) 1 + x i β θ i = 1 n x i β ln ( x i ) = 0 .
The simultaneous solution of the above likelihood equations constitutes the MLEs of θ and β . However, analytical solutions are not feasible. Therefore, using R Core Team [21] software, we approximated the solution to the aforementioned normal equations. In addition, we have
2 log L β 2 = θ i = 1 n x i β log 2 ( x i ) + i = 1 n ( x i β log 2 ( x i ) x i β + 1 x i 2 β log 2 ( x i ) ( x i β + 1 ) 2 ) n β 2 < 0 ,
and
2 log L β 2 = θ i = 1 n x i β log 2 ( x i ) + i = 1 n ( x i β log 2 ( x i ) x i β + 1 x i 2 β log 2 ( x i ) ( x i β + 1 ) 2 ) n β 2 < 0 ,
2 log L θ 2 = n ( θ + 2 ) 2 3 n θ 2 < 0 .

4.2. Methods of LS and WLS

As the explicit forms of the maximum likelihood estimators cannot be obtained every time, some other methods are developed to overcome this problem. The LS and the WLS are well-known methods used for estimating unknown parameters [22]. Here, we consider these two approaches for estimating the unknown parameters of the P Z D . Let x 1 , x 2 ,   , x n be the ordered observations obtained from a sample of size n from the P Z D . By calculating the minimum of the function
S ( Ψ ) = i = 1 n η i { 1 ( 1 + θ x i β ( θ + 2 + θ x i β ) 2 + θ ) e θ x i β i n + 1 } 2
with respect to θ and β , respectively, the LS estimates θ ^ L S E and β ^ L S E can be obtained by setting η i = 1 , while we can obtain the WLS estimates θ ^ W L S and β ^ W L S by setting η i = ( n + 1 ) 2 ( n + 2 ) i ( n i + 1 ) . The following equations can be solved to yield these estimates as well as
S ( Ψ ) θ = i = 1 n η i { 1 ( 1 + θ x i β ( θ + 2 + θ x i β ) 2 + θ ) e θ x i β i n + 1 } φ 1 ( x i , Ψ ) = 0 ,
S ( Ψ ) β = i = 1 n η i { 1 ( 1 + θ x i β ( θ + 2 + θ x i β ) 2 + θ ) e θ x i β i n + 1 } φ 2 ( x i , Ψ ) = 0 ,
where
φ 1 ( x i , Ψ ) = ( θ x i β + 2 x i β + θ + 3 ) θ 2 x i 2 β e θ x i β 2 + θ ,   and   φ 2 ( x i , Ψ ) = ( 1 + x i β ) θ 3 x i 2 β ln ( x i ) e θ x i β 2 + θ .

4.3. Method of Maximum Product of Spacing

According to [23], and based on the differences among the values of the cdf at successive data points, the maximum product of spacing is an estimation method as interesting as that of MLE. Moreover, ref. [24] concluded that the MPS estimator method outperforms all the other estimator methods.
The uniform spacings can be defined using a random sample of size n from a distribution with cdf F ( x i , Ψ )   as follows:
D i ( Ψ ) = F ( x i , Ψ ) F ( x i 1 , Ψ ) , i = 1 , 2 , ... , n
where F ( x 0 , Ψ ) = 0 and F ( x n + 1 , Ψ ) = 1 . The MPS estimates denoted by   θ ^ M P S   and   β ^ M P S   maximize the function M ( Ψ ) with respect to the unknown parameters θ and β as
M ( Ψ ) = 1 n + 1 i = 1 n + 1 log D i ( Ψ ) = 1 n + 1 i = 1 n + 1 log { e θ x i 1 β ( 1 + θ x i 1 β ( θ + 2 + θ x i 1 β ) 2 + θ ) e θ x i β ( 1 + θ x i β ( θ + 2 + θ x i β ) 2 + θ ) } ,
or by solving the subsequent equations
M ( Ψ ) θ = 1 n + 1 i = 1 n + 1 φ 1 ( x i 1 , Ψ ) φ 1 ( x i , Ψ ) D i = 0 ,  
M ( Ψ ) β = 1 n + 1 i = 1 n + 1 φ 2 ( x i 1 , Ψ ) φ 2 ( x i , Ψ ) D i = 0 ,
where φ 1 ( x i , Ψ ) and φ 2 ( x i , Ψ ) are given earlier. For more details on this method, see [25].

4.4. Cramer–Von Mises and Anderson–Darling Methods

In recent years, some authors have used classical GOF statistics such as Cramer–von Mises and Anderson–Darling statistics to derive the estimators of the unknown parameters. Ref. [26] used these methods for the binomial-exponential 2 distributions, while [27] evaluated the process capability index for normal distribution. Aside from these methods and using this approach, we propose the use of the Pearson Chi-square statistic X 2 where the estimators noted C S M M are obtained by calculating the minimum of X 2 with respect to the unknown parameters.
Ref. [28] and later [29] showed that the Cramer–von Mises estimates based on the distance of Cramer–von Mises GOF statistics are the least biased compared to the other estimators. For calculation purposes, ref. [29] gives the formula C ( Ψ ) for this statistic, where the estimators ensure its minimum concerning the unknown parameters, where
C ( Ψ ) = 1 12 n + i = 1 n { 1 ( 1 + θ x i β ( θ + 2 + θ x i β ) 2 + θ ) e θ x i β 2 i 1 2 n } 2 .
They also can be gained as the solution of the subsequent equations:
C ( Ψ ) θ = i = 1 n { 1 ( 1 + θ x i β ( θ + 2 + θ x i β ) 2 + θ ) e θ x i β 2 i 1 2 n } φ 1 ( x i , Ψ ) = 0
C ( Ψ ) β = i = 1 n { 1 ( 1 + θ x i β ( θ + 2 + θ x i β ) 2 + θ ) e θ x i β 2 i 1 2 n } φ 2 ( x i , Ψ ) = 0 .
The Anderson–Darling estimation method is based on the classical Anderson and Darling   ( A D )   goodness-of-fit statistic test proposed by [30]. This statistic is essentially used to fit data to a theoretical hypothesized model F 0 ( x ) . When the parameters are unknown, they can be estimated by θ ^ A D E , β ^ A D E which minimize the A D statistic given in the form
A D E ( Ψ ) = n n 1 i = 1 n ( 2 i 1 ) { log F ( x i ) + log ( 1 F ( x i + n + 1 ) ) } ,
which is equivalent to canceling the first derivatives of this function with respect to β and θ . For more details on this method, see [29]. In their paper, ref. [31] showed that the A D E estimation method gives the most efficient estimators. As is known, the Anderson–Darling statistic gives more weight to the tails of the distribution, so for right- or left-tailed distributions, the A D E right tail ( A D E R T ) and A D E left tail ( A D E L T ) can be considered for estimating the unknown parameters, as in [27,31,32].
The A D E R T estimators of θ and β denoted by θ ^ A D E ( R T ) , β ^ A D E ( R T ) , respectively, are found by minimizing the equation given by
A D E ( R T ) ( Ψ ) = n 2 2 i = 1 n log F ( x i ) 1 n i = 1 n ( 2 i 1 ) log { 1 F ( x i + n + 1 ) }
or
θ [ A D E ( R T ) ( Ψ ) ] = 0 , β [ A D E ( R T ) ( Ψ ) ] = 0
Similarly, the A D E L T estimators denoted by θ ^ A D E ( L T ) , β ^ A D E ( L T ) are acquired by minimizing the equation
A D E ( L T ) ( Ψ ) = 3 2 n + 2 i = 1 n F ( x i ) 1 n i = 1 n ( 2 i 1 ) log { F ( x i ) }
or
θ [ A D E ( L T ) ( Ψ ) ] = 0 , β [ A D E ( L T ) ( Ψ ) ] = 0
with respect to the unknown distribution parameters. For more about the Anderson–Darling estimators, see [27,31,32]. Conversely, ref. [33] demonstrated that the A D E minimum-distance estimation method is a highly competitive method compared to the ML estimation. In addition, ref. [24] confirmed the superiority of this method.

4.5. The Nikulin–Rao–Robson Test

The probability distribution used to describe any phenomenon is very important in this analysis. Since the twentieth century, researchers have continued to develop techniques to validate the different models. In addition to the classical procedures, we considered in this work a criteria test statistic based on the Nikulin–Rao–Robson statistic (NRR) to fit data to the PZD . The great interest of the NRR statistic is that it recovers all the information lost while grouping data. Based on the MLE on initial data, this statistic Y 2 introduced by [34,35] is a modification of the well-known Chi-square Pearson statistic T 2 . For testing the null hypothesis H 0 , given by
H 0 = P ( X i x ) = F ( x , Ψ ) , x 0 .  
According to the above, the sample X 1 , X 2 , ,   X n belongs to the parametric family F ( x , Ψ ) with unknown parameters, and Y 2 consists in grouping data into r equiprobable intervals I 1 , I 2 , ... , I r , where
I j = ] a j 1 , a j ] ; I i I j = ϕ i j ; j = 1 r I j = R 1 ,
such as
p j = a j 1 a j f ( x , Ψ ) d x = 1 r , j = 1 , 2 , ... , r .
If υ = ( υ 1 , υ 2 , ... , υ r ) T represents the number of observed X i grouping into these intervals I j , and the vector T n ( Ψ ) is given by
T n ( Ψ ^ ) = ( v 1 n p 1 ( Ψ ^ ) n p 1 ( Ψ ^ ) , v 2 n p 2 ( Ψ ^ ) n p 2 ( Ψ ^ ) , , v r n p r ( Ψ ^ ) n p r ( Ψ ^ ) ) t .
The NRR statistic Y 2 is defined as
Y n 2 ( Ψ ^ ) = T n 2 ( Ψ ^ ) + 1 n l T ( Ψ ^ ) ( I ( Ψ ^ ) J ( Ψ ^ ) ) 1 l ( Ψ ^ ) ,
where I ( Ψ ^ ) and J ( Ψ ^ ) are the estimated information matrices on non-grouped and grouped data, respectively, and Ψ ^ is the vector of the ML estimators on initial data. The components of the vector l ( Ψ ^ ) = ( l s ( Ψ ^ ) ) 1 × s t and the matrix J ( Ψ ^ ) are
l s ( Ψ ^ ) = i = 1 r ν j p j p j ( Ψ ^ ) Ψ ^ s   ,  
J ( Ψ ^ ) = ( 1 p j p j ( Ψ ^ ) Ψ ^ k ) r × s t ( 1 p j p j ( Ψ ^ ) Ψ ^ k ) r × s ,   j = 1 , 2 , , r   and   k = 1 , 2 , , s ,
where s represents the number of parameters. The distribution of Y 2 ( Ψ ^ ) is a Chi-square with r 1 degrees of freedom. To construct the test statistic Y 2 corresponding to the PZD ( Ψ ^ ) we calculated the ML estimators Ψ ^ = ( β ^ , θ ^ ) t , the limit intervals a j , and the derivatives p j ( Ψ ^ ) Ψ ^ k for our model.
p j ( Ψ ^ ) θ ^ = φ 1 ( x j , Ψ ^ ) φ 1 ( x j 1 , Ψ ^ )   , p j ( Ψ ^ ) β ^ = φ 2 ( x j , Ψ ^ ) φ 2 ( x j 1 , Ψ ^ ) .
Then, after computing the estimated information matrices I ( Ψ ^ ) and J ( Ψ ^ ) , we obtained the statistic Y 2 ( Ψ ^ ) for the P Z D distribution. This method is used to construct GOF test statistics for some generalized models; for examples see [36,37,38,39].

5. Simulations

In this section, we investigate the performance of the estimation procedures used in this paper. To that end, we generated 10 , 000 samples with different sizes ( n = 10 ,   n = 25 ,   n = 35 ,   n = 75 ,   n = 200 , and n = 400 ) from the P D Z model, and then computed the different estimators and their mean square errors (MSE) (values in brackets) for the unknown parameters using R statistical software when θ = 5 and β = 0.5 . The results are given by sample size in Table 3.
Based on Table 2, it can be noted that the bias and MSE values of different estimators decrease as the sample size increases. For example, when θ = 5 with n = 10 and n = 400 , the MLE values are θ ^ = 4.8841 and θ ^ = 4.9989 , respectively. In addition, for the same case of θ = 5 , n = 10 , and n = 400 , the values of the LSE are θ ^ = 5.0759 and θ ^ = 5.0167 , respectively. In general, we can say that the LSE method performs better than other methods up to n 35 , while the MLE is the best one for n 75 in terms of the MSE.

6. Applications

The analysis of three real data sets is considered to show the usefulness of the proposed distribution in modeling different phenomena and the performances of the methods used to determine the unknown parameters. In addition to the traditional model selection techniques, we calculated the Y 2 statistic to prove that the proposed model fits some data sets better than some commonly used alternative distributions, such as:
-
The xgamma distribution suggested by [40] with pdf given by
f ( x ; θ ) = θ 2 1 + θ ( 1 + θ 2 x 2 ) e θ x , x > 0 , θ > 0 .
-
Gamma distribution with pdf given by
f ( x ; α , β ) = 1 β α Γ ( α ) x α 1 e x β , x > 0 , α , β > 0 .
-
Lomax distribution with pdf given by
f ( x ; α , λ ) = α λ ( 1 + x λ ) ( α + 1 ) , x > 0 , θ , β > 0 .
-
Darna distribution proposed by [41] with pdf given by
f D D ( y ; θ , α ) = θ 2 α 2 + θ 2 ( 2 α + θ 4 y 2 2 α 3 ) e θ y α ; y > 0 , θ > 0 , α > 0 ,
-
Power Darna distribution introduced by [13] with pdf given by
f P D D ( x ; α , β , θ ) = β θ 2 α 3 ( 2 α 2 + θ 2 ) x β 1 ( 4 α 4 + θ 4 x 2 β ) e θ x β α , x > 0 , θ , β , α > 0 .
-
Power Lindley distribution proposed by [2] with pdf defined as
f P L D ( x ; β , θ ) = β θ 2 1 + θ x β 1 ( 1 + x β ) e θ x β , x > 0 , θ , β > 0 .
-
Exponentiated power Lindley distribution (EPLD) offered by [42] with pdf given by
f E P L D ( x ; β , θ ) = α β θ 2 x β 1 1 + θ ( 1 + x β ) e θ x β [ 1 ( 1 + θ x β 1 + θ ) e θ x β ] α 1 , x > 0 , θ , β > 0 .

6.1. First Data: Failure Times of Devices

The first real data set refers to the failure times of 50 items put under a life test. It is generally known that this data set displays bathtub behavior of the hazard rate function [43]. This data set has been analyzed by many authors. The data are 0.1, 1.0, 1.0, 1.0, 2.0, 3.0, 82.0, 83.0, 84.0, 6.0, 7.0, 11.0, 36.0, 40.0, 45.0, 12.0, 18.0, 18.0, 21.0, 18.0, 0.2, 1.0, 1.0, 18.0, 18.0, 32.0, 45.0, 47.0, 63.0, 63.0, 67.0, 67.0, 50.0, 55.0, 60.0, 67.0, 67.0, 72.0, 85.0, 85.0, 85.0, 85.0, 75.0, 79.0, 82.0, 84.0, 84.0, 85.0, 86.0, 86.0. Figure 9 displays the density, total time on test transform (TTT), probability–probability (P–P) plot, and quantile–quantile (Q–Q) plot of the PZD for the failure times data.
Table 4 represents the values of the parameter estimators for the hypothesized distribution P Z D obtained by the different methods, and Kolmogorov–Smirnov (KS) and the corresponding p-values. Because n = 50 > 30 , it is established by the simulation study that the M L E methods give the best results.
As the maximum likelihood parameter estimators for the competing distributions are required to calculate their corresponding NRR test statistic Y 2 , the MLE values are calculated and given in Table 5.
To distinguish between the proposed model and its alternatives, we used classical criteria for model selection, and the NRR statistics with the corresponding values are summarized in Table 6.
From Table 6, we can see that the smallest values of the different criteria for model selection are obtained for the PZD, which confirms that the proposed PZD model describes these data better than all the alternatives.

6.2. Second Data: Remission Times

The second case addressed the lengths of remission times (in months) of 128 bladder cancer patients who are randomly selected, as reported by [44]. The data observations are 0.08, 4.87, 6.94, 8.66, 2.09, 3.48, 13.11, 23.63, 0.20, 13.80, 25.74, 0.50, 2.23, 3.52, 4.98, 6.97, 9.02, 3.88, 5.32, 7.39, 10.34, 13.29, 0.40, 2.26, 9.22, 2.46, 3.64, 5.09, 7.26, 0.51, 2.54, 14.76, 26.31, 0.81, 3.70, 5.17, 7.28, 9.47, 14.24, 25.82, 9.74, 2.62, 3.82, 5.32, 2.69, 4.23, 5.41, 7.62, 7.32, 10.06, 14.77, 32.15, 2.64, 14.83, 34.26, 4.33, 5.49, 7.66, 0.90, 5.34, 7.59, 10.66, 15.96, 2.69, 4.18, 36.66, 12.05, 10.75, 16.62, 43.01, 5.41, 7.63, 11.25, 17.14, 79.05, 17.12, 46.12, 4.40, 5.85, 8.26, 1.26, 2.83, 1.35, 2.87, 5.62, 1.19, 2.75, 4.26, 7.87, 11.64, 4.34, 5.71, 7.93, 11.79, 17.36, 1.40, 3.02, 18.10, 1.46, 11.98, 19.13, 12.02, 2.02, 3.31, 1.76, 3.25, 12.03, 20.28, 2.02, 4.50, 6.25, 8.37, 4.51, 6.54, 8.53, 3.36, 6.76, 8.65, 12.63, 22.69, 12.07, 21.73, 3.57, 5.06, 7.09, 2.07, 3.36, 6.93. Figure 10 represents the density, TTT, P–P, and Q–Q plots of the PZD for the 128 bladder cancer patients’ data. Various parameter estimators based on cancer patients’ data are summarized in Table 7.
As the ML parameter estimators for the competing distributions are needed to compute their corresponding NRR test statistic Y 2 , the MLE values are calculated and given in Table 8.
To show that this model can also describe this type of data, we used all classical criteria for model selection and the modified Chi-square NRR. All of the corresponding values are computed and given in Table 9.
From Table 9, we can see that the smallest values of the different criteria for model selection are obtained for the PZD, which confirms that the proposed model describes these data better than all the alternatives considered here.

6.3. Third Data: Strengths of 1.5 cm Glass Fibers

We then considered the data set representing the strengths of 1.5 cm glass fibers. The data set is also analyzed by [45,46]. The data are 0.55, 0.93, 1.25, 1.36, 1.49, 1.52, 1.58, 1.61, 1.64, 1.68, 1.73, 1.81, 2.00, 0.74, 1.04, 1.27, 1.39, 1.49, 1.53, 1.59, 1.61, 1.66, 1.68, 1.76, 1.82, 2.01, 0.77, 1.11, 1.28, 1.42, 1.50, 1.54, 1.60, 1.62, 1.66, 1.69, 1.76, 1.84, 2.24, 0.81, 1.13, 1.29, 1.48, 1.50, 1.55, 1.61, 1.62, 1.66, 1.70, 1.77, 1.84, 0.84, 1.24, 1.30, 1.48, 1.51, 1.55, 1.61, 1.63, 1.67, 1.70, 1.78, 1.89. Figure 11 represents the density, TTT, P–P, and Q–Q plots of the PZD for the glass fiber data.
Table 10 represents the values of the parameter estimators for the hypothesized distribution PZD obtained by the different methods and the corresponding p-values.
As the maximum likelihood parameter estimators for the competing distributions are needed to compute their corresponding NRR test statistic Y2, the MLE values are calculated and given in Table 11.
To show that this model can also describe this type of data, we used all classical criteria for model selection and the modified chi-square NRR. All the corresponding values are computed and given in Table 12.

7. Conclusions and Suggestions

The power transformation approach is used in this research to generalize the Zeghdoudi distribution of one parameter to a power Zeghdoudi distribution of two parameters. Some statistical properties, including moments, reliability functions, order statistics, Bonferroni and Lorenz curves, Gini index, stochastic ordering, quantile function, and mean and median deviations are derived. Numerous estimation techniques are considered, such as the maximum likelihood estimation, least squares, weighted least squares, maximum product of spacing methods, the Cramer–von Mises method, and the right-tail and left-tail Anderson–Darling estimation methods. A simulation study is carried out to assess these estimation approaches and determine the distribution properties. Three actual data sets are examined and fitted and the Anderson–Darling and Nikulin–Rao–Robson test statistics are considered for the novel PZD to investigate the distribution’s applicability in different fields compared with the gamma, xgamma, Lomax, Darna, power Darna, power Lindley, and exponentiated power Lindley models. It is found that PZD often performs better than other distributions in fitting the three data sets considered in this study. In addition, the MSE values of different proposed estimators decrease for increased values of the sample size. Future work could include bivariate extensions, transmutation extensions, discrete extensions, and additional regression models. In addition, the PZD parameters can be estimated using the ranked set sampling method and some of its modifications. See [47] for stratified median ranked set sampling and [48] for double percentile ranked set samples.

Author Contributions

Conceptualization, A.I.A.-O.; Methodology, K.A. and A.I.A.-O.; Software, K.A. and A.I.A.-O.; Formal analysis, K.A.; Investigation, K.A., A.I.A.-O. and R.A.; Data curation, K.A.; Writing—original draft, K.A., A.I.A.-O. and R.A.; Writing—review & editing, A.I.A.-O. and R.A.; Funding acquisition, R.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data used in this study are available within the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ghitany, M.E.; Atieh, B.; Nadarajah, S. Lindley distribution and its application. Math. Comput. Simul. 2008, 78, 493–506. [Google Scholar] [CrossRef]
  2. Ghitany, M.; Al-Mutairi, D.K.; Balakrishnan, N.; Al-Enezi, L. Power Lindley distribution and associated inference. Comput. Stat. Data Anal. 2013, 64, 20–33. [Google Scholar] [CrossRef]
  3. Shanker, R.; Shukla, K.K. Power Shanker distribution and its application. Turk. Klin. Biyoistatistik 2017, 9, 175–187. [Google Scholar] [CrossRef]
  4. Shanker, R.; Shukla, K.K. A two-parameter power Aradhana distribution with properties and application. Indian Soc. Ind. Appl. Math. 2018, 9, 210–220. [Google Scholar] [CrossRef]
  5. Abebe, B.; Tesfay, M.; Eyob, T.; Shanker, R. A two-parameter power Garima distribution with properties and applications. Annal. Biostat. Biomed. Appl. 2019, 1, 515. [Google Scholar]
  6. Abebe, B.; Tesfay, M.; Eyob, T.; Shanker, R. A two-parameter power Rama distribution with properties and applications. Biom. Biostat. Int. J. 2019, 8, 6–11. [Google Scholar]
  7. Shukla, K.K.; Shanker, R. Power Ishita distribution and its application to model lifetime data. Stat. Transit. New Ser. 2018, 19, 135–148. [Google Scholar] [CrossRef] [Green Version]
  8. Shanker, R.; Shukla, K.K. A two-parameter power Sujatha distribution with properties and application. Int. J. Math. Stat. 2019, 20, 11–22. [Google Scholar]
  9. Shukla, K.K.; Shanker, R. Power Prakaamy distribution and its applications. Int. J. Comput. Theor. Stat. 2020, 7, 25–36. [Google Scholar]
  10. Krishnarani, S. On a power transformation of half-logistic distribution. J. Probab. Stat. 2016, 2016, 2084236. [Google Scholar] [CrossRef] [Green Version]
  11. Habibi, M.; Asgharzadeh, A. Power binomial exponential distribution: Modeling, simulation and application. Commun. Stat.-Simul. Comput. 2018, 47, 3042–3061. [Google Scholar] [CrossRef]
  12. Al-Omari, A.; Alhyasat, K.; Ibrahim, K.; Abu Bakar, M. Power length-biased Suja distribution: Properties and application. Electron. J. Appl. Stat. Anal. 2019, 12, 429–452. [Google Scholar]
  13. Al-Omari, A.I.; Aidi, K.; Alsultan, R. Power Darna distribution with right censoring: Estimation, testing, and applications. Appl. Sci. 2022, 12, 8272. [Google Scholar] [CrossRef]
  14. Meriem, B.; Gemeay, A.M.; Almetwally, E.M.; Halim, Z.; Alshawarbeh, E.; Abdulrahman, A.T.; Hussam, E. The Power XLindley Distribution: Statistical Inference, Fuzzy Reliability, and COVID-19 Application. J. Funct. Spaces 2022, 2022, 9094078. [Google Scholar] [CrossRef]
  15. Afify, A.Z.; Gemeay, A.M.; Alfaer, N.M.; Cordeiro, G.M.; Hafez, E.H. Power-Modified Kies-Exponential Distribution: Properties, Classical and Bayesian Inference with an Application to Engineering Data. Entropy 2022, 24, 883. [Google Scholar] [CrossRef]
  16. Messaadia, H.; Zeghdoudi, H. Zeghdoudi distribution and its applications. Int. J. Comput. Sci. Math. 2018, 9, 58–65. [Google Scholar] [CrossRef]
  17. Habib, M.; Thomas, D. Chi-square goodness-if-fit tests for randomly censored data. Ann. Stat. 1986, 14, 759–765. [Google Scholar] [CrossRef]
  18. Hollander, M.; Pena, E.A. A chi-squared goodness-of-fit test for randomly censored data. J. Am. Stat. Assoc. 1992, 87, 458–463. [Google Scholar] [CrossRef]
  19. Galanova, N.; Lemeshko, B.Y.; Chimitova, E. Using nonparametric goodness-of-fit tests to validate accelerated failure time models. Optoelectron. Instrum. Data Process. 2012, 48, 580–592. [Google Scholar] [CrossRef]
  20. Bagdonavicius, V.; Nikulin, M. Chi-squared goodness-of-fit test for right censored data. Int. J. Appl. Math. Stat. 2011, 24, 30–50. [Google Scholar]
  21. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020. [Google Scholar]
  22. Swain, J.J.; Venkatraman, S.; Wilson, J.R. Least-squares estimation of distribution functions in Johnson’s translation system. J. Stat. Comput. Simul. 1988, 29, 271–297. [Google Scholar] [CrossRef]
  23. Cheng, R.C.H.; Amin, N.A.K. Estimating parameters in continuous univariate distribution with a shifted origin. J. R. Stat. Soc. 1983, 45, 394–403. [Google Scholar] [CrossRef]
  24. Al-Mofleh, H.; Afify, A.Z. A generalization of Ramos-Louzada distribution: Properties and estimation, properties and estimation. arXiv 2019, arXiv:1912.08799v1. [Google Scholar]
  25. Cheng, R.C.H.; Stephens, M.A. A goodness-of-fit test using Moran’s statistic with estimated parameters. Biometrika 1989, 76, 385–392. [Google Scholar] [CrossRef]
  26. Bakouch, H.S.; Dey, S.; Ramos, P.L.; Louzada, F. Binomial-exponential 2 distribution: Different estimation methods with weather applications. TEMA 2017, 18, 233–251. [Google Scholar] [CrossRef] [Green Version]
  27. Saha, M.; Dey, S.; Yadav, A.S.; Ali, S. Confidence intervals of the index Cpk for normally distributed quality characteristics using classical and Bayesian methods of estimation. Braz. J. Probab. Stat. 2021, 35, 138–157. [Google Scholar] [CrossRef]
  28. MacDonald, P.D.M. Comment on an estimation procedure for mixtures of distributions” by choi and bulgren. J. R. Stat. Society. Ser. B (Methodol.) 1971, 33, 326–329. [Google Scholar]
  29. Boos, D.D. Minimum distance estimators for location and goodness of fit. J. Am. Stat. Assoc. 1981, 76, 663–667. [Google Scholar] [CrossRef]
  30. Anderson, T.W.; Darling, D.A. Asymptotic Theory of Certain “Goodness of Fit” Criteria Based on Stochastic Processes. Ann. Math. Statist. 1952, 23, 193–212. [Google Scholar] [CrossRef]
  31. Rodrigues, G.C.; Louzada, F.; Ramos, P.L. Poisson–exponential distribution. different methods of estimation. J. Appl. Stat. 2018, 45, 128–144. [Google Scholar] [CrossRef]
  32. ZeinEldin, R.A.; Ahsan ul Haq, M.; Hashmi, S.; Elsehety, M. Alpha power transformed inverse Lomax distribution with different methods of estimation and applications. Complexity 2020, 2020, 1860813. [Google Scholar] [CrossRef] [Green Version]
  33. Akgül, F.G. Comparison of the estimation methods for the parameters of exponentiated reduced Kies distribution. J. Nat. Appl. Sci. 2018, 22, 1209–1216. [Google Scholar]
  34. Nikulin, M.S. Chi-square test for continuous distributions with location and scale parameters. Teor. Veroyatnostei Primen. 1973, 18, 583–591. [Google Scholar]
  35. Rao, K.C.; Robson, D.S. A Chi-square statistic for goodness-of-fit tests within the exponential family. Commun. Stat. 1974, 3, 1139–1153. [Google Scholar] [CrossRef]
  36. Goual, H.; Seddik-Ameur, N. Chi-squared type test for the AFT-generalized inverse Weibull distribution. Commun. Stat.-Theory Methods 2014, 43, 2605–2617. [Google Scholar] [CrossRef]
  37. Seddik-Ameur, N.; Aidi, K. Chi-square tests for generalized exponential distributions with censored data. Electron. J. Appl. Stat. Anal. 2016, 9, 371–384. [Google Scholar]
  38. Tilbi, D.; Seddik-Ameur, N. Chi-squared goodness-of-fit tests for the generalized Rayleigh distribution. J. Stat. Theory Pract. 2017, 11, 594–603. [Google Scholar] [CrossRef]
  39. Treidi, W.; Seddik-Ameur, N. NRR statistic for the extension Weibull distribution. Glob. J. Pure Appl. Math. 2016, 12, 2809–2818. [Google Scholar]
  40. Sen, S.; Maiti, S.S.; Chandra, N. The Xgamma distribution: Statistical properties and application. J. Mod. Appl. Stat. Methods 2016, 15, 774–788. [Google Scholar] [CrossRef] [Green Version]
  41. Al-Omari, A.I.; Shraa, D. Darna distribution: Properties and application. Electron. J. Appl. Stat. Anal. 2019, 12, 520–541. [Google Scholar]
  42. Ashour, S.K.; Eltehiwy, M.A. Exponentiated power Lindley distribution. J. Adv. Res. 2015, 6, 895–905. [Google Scholar] [CrossRef] [Green Version]
  43. Aarset, M.V. How to identify a bathtub hazard rate. IEEE Trans. Reliab. 1987, 36, 106–108. [Google Scholar] [CrossRef]
  44. Lee, E.T.; Wang, J.W. Statistical Methods for Survival Data Analysis, 3rd ed.; Wiley: New York, NY, USA, 2003. [Google Scholar]
  45. Bourguignon, M.; Silva, R.B.; Cordeiro, G.M. The Weibull G family of probability distributions. J. Data Sci. 2014, 12, 53–68. [Google Scholar] [CrossRef] [PubMed]
  46. Smith, R.L.; Naylor, J.C. A Comparison of maximum likelihood and Bayesian estimators for the three-parameter Weibull distribution. Appl. Stat. 1987, 36, 358–360. [Google Scholar] [CrossRef]
  47. Ibrahim, K.; Syam, M.; Al-Omari, A.I. Estimating the population mean using stratified median ranked set sampling. Appl. Math. Sci. 2010, 4, 2341–2354. [Google Scholar]
  48. Jemain, A.A.; Al-Omari, A.I. Double percentile ranked set samples for estimating the population mean. Adv. Appl. Stat. 2006, 6, 261–276. [Google Scholar]
Figure 1. Plots of the PZD survival function for different parameter values.
Figure 1. Plots of the PZD survival function for different parameter values.
Applsci 12 12081 g001
Figure 2. Plots of the cdf of PZD for selected values of β and θ .
Figure 2. Plots of the cdf of PZD for selected values of β and θ .
Applsci 12 12081 g002
Figure 3. Pdf plots of the PZD with various values of parameters.
Figure 3. Pdf plots of the PZD with various values of parameters.
Applsci 12 12081 g003
Figure 4. Three-dimensional plots of maximum and minimum pdf of the PZD for m = 5 ,   β = 2 .
Figure 4. Three-dimensional plots of maximum and minimum pdf of the PZD for m = 5 ,   β = 2 .
Applsci 12 12081 g004
Figure 5. The hazard function plots of PZD for various values of parameters.
Figure 5. The hazard function plots of PZD for various values of parameters.
Applsci 12 12081 g005
Figure 6. The RHF plots of PZD for various values of parameters.
Figure 6. The RHF plots of PZD for various values of parameters.
Applsci 12 12081 g006
Figure 7. The odds function plots of PZD for selected values of parameters.
Figure 7. The odds function plots of PZD for selected values of parameters.
Applsci 12 12081 g007
Figure 8. The MRL plots of PZD for selected values of θ ,   β .
Figure 8. The MRL plots of PZD for selected values of θ ,   β .
Applsci 12 12081 g008
Figure 9. The density, TTT, P–P, and Q–Q plots of the PZD for the failure times data.
Figure 9. The density, TTT, P–P, and Q–Q plots of the PZD for the failure times data.
Applsci 12 12081 g009aApplsci 12 12081 g009b
Figure 10. The density, TTT, P–P, and Q–Q plots of the PZD for bladder cancer patients’ data.
Figure 10. The density, TTT, P–P, and Q–Q plots of the PZD for bladder cancer patients’ data.
Applsci 12 12081 g010
Figure 11. The density, TTT, P–P, and Q–Q plots of the PZD with glass fiber data.
Figure 11. The density, TTT, P–P, and Q–Q plots of the PZD with glass fiber data.
Applsci 12 12081 g011
Table 1. Values of μ X , σ X , ƛ 1 , ƛ 2 , ƛ 3 for the PZD for selected parameters.
Table 1. Values of μ X , σ X , ƛ 1 , ƛ 2 , ƛ 3 for the PZD for selected parameters.
θ μ X σ X ƛ 1 ƛ 2 ƛ 3 β μ X σ X ƛ 1 ƛ 2 ƛ 3
β = 3 θ = 2
30.8810.2110.2390.0032.8564.01.0120.1810.179−0.1952.977
3.10.8700.2080.2390.0042.8554.11.0120.1770.175−0.2092.990
3.20.8600.2060.2400.0062.8544.21.0110.1730.171−0.2223.002
3.30.8500.2040.2400.0072.8534.31.0100.1690.167−0.2343.015
3.40.8410.2020.2400.0092.8534.41.0100.1650.163−0.2463.028
3.50.8320.2000.2410.0102.8524.51.0090.1610.160−0.2583.041
3.60.8230.1980.2410.0112.8524.61.0090.1580.157−0.2693.054
3.70.8150.1960.2410.0132.8514.71.0080.1550.153−0.2803.067
3.80.8070.1950.2410.0142.8514.81.0080.1520.150−0.2903.080
3.90.7990.1930.2410.0152.8514.91.0080.1490.147−0.3003.093
Table 2. Values for some selected values of p.
Table 2. Values for some selected values of p.
p θ = 2 ,   β = 4 θ = 3 ,   β = 3 θ = 7 ,   β = 5 θ = 0.5 ,   β = 0.8
0.10.774930.607310.611842.14411
0.20.860980.700590.667583.43875
0.250.893110.736380.688274.05443
0.300.921660.768610.706564.67193
0.400.972440.826870.739045.95421
0.501.018940.881250.768627.36366
0.601.064520.935440.797468.99532
0.701.112240.993090.8274911.0074
0.751.138201.024840.8437612.2477
0.801.166751.060030.8615813.7405
0.901.240051.151750.9070218.2684
0.951.298681.226430.9430122.7118
0.991.404361.363851.0070732.9481
Table 3. Estimators’ values of θ = 5 and β = 0.5 with different estimation methods.
Table 3. Estimators’ values of θ = 5 and β = 0.5 with different estimation methods.
θ = 5 β = 0.5 θ = 5 β = 0.5
n = 10 n = 25
MLE4.8841 (0.0122)0.5899 (0.0124)4.8993 (0.0105)0.5776 (0.0106)
WLS4.8973 (0.0120)0.5832 (0.0117)4.9176 (0.0103)0.5719 (0.0099)
MPS4.8812 (0.0136)0.6118 (0.0147)4.8935 (0.0119)0.5995 (0.0129)
CMEs5.0735 (0.0126)0.3099 (0.0131)5.0632 (0.0109)0.4023 (0.0113)
LSE5.0759 (0.0098)0.4244 (0.0093)5.0607 (0.0086)0.4355 (0.0085)
ADE5.1405 (0.0158)0.3768 (0.0159)5.1282 (0.0141)0.3928 (0.0141)
ADEL-R5.1477 (0.0167)0.3709 (0.0171)5.1354 (0.0150)0.3851 (0.0153)
ADEL-T5.1656 (0.0175)0.3653 (0.0177)5.1533 (0.0158)0.3782 (0.0159)
n = 35 n = 75
MLE4.9146 (0.0090)0.5611 (0.0089)4.9644 (0.0040)0.5348 (0.0037)
WLS4.9391 (0.0088)0.5567 (0.0082)4.9271 (0.0079)0.5445 (0.0070)
MPS4.9069 (0.0104)0.5832 (0.0112).4.9224 (0.0090)0.5567 (0.0077)
CMEs5.0517 (0.0094)0.4152 (0.0096)5.0364 (0.0072)0.4274 (0.0061)
LSE5.0454 (0.0071)0.4507 (0.0068)5.0329 (0.0055)0.4623 (0.0052)
ADE5.1068 (0.0127)0.4011 (0.0124)5.0852 (0.0111)0.4143 (0.0093)
ADEL-R5.1143 (0.0136)0.3980 (0.0136)5.0924 (0.0120)0.4112 (0.0099)
ADEL-T5.1319 (0.0144)0.3914 (0.0142)5.1103 (0.0125)0.4056 (0.0116)
n = 200 n = 400
MLE4.9869 (0.0027)0.5123 (0.0024)4.9989 (0.0012)0.5003 (0.0009)
WLS4.9586 (0.0061)0.5183 (0.0053)4.9813 (0.0043)0.5054 (0.0035)
MPS4.9437 (0.0072)0.5304 (0.0060)4.9763 (0.0056)0.5078 (0.0042)
CMEs5.0195 (0.0053)0.5412 (0.0044)5.0096 (0.0038)0.4834 (0.0026)
LSE5.0167 (0.0039)0.4785 (0.0035)5.0035 (0.0023)0.4908 (0.0015)
ADE5.0527 (0.0093)0.4408 (0.0076)5.0234 (0.0077)0.4733 (0.0058)
ADEL-R5.0633 (0.0102)0.4377 (0.0081)5.0302 (0.0086)0.4702 (0.0063)
ADEL-T5.0741 (0.0110)0.4321 (0.0094)5.0416 (0.0092)0.4646 (0.0076)
Table 4. Different parameter estimators for failure times data.
Table 4. Different parameter estimators for failure times data.
Method β θ KSp-Value
MLE0.67634.67770.17260.9496
WLS0.36462.31620.19760.9284
MPS0.23463.46770.20110.9202
CMEs0.36795.66260.19380.9334
LSE0.59764.76470.17840.9412
ADE0.66966.61330.21340.9177
ADEL-R0.80107.67270.22120.9086
ADEL-T0.92628.76890.22760.9010
Table 5. ML parameter estimates for the alternative distributions based on failure times data.
Table 5. ML parameter estimates for the alternative distributions based on failure times data.
DistributionsEstimates
PZD β ^ = 0.6763 , θ ^ = 4.6777
ZD θ ^ = 1.7826
Lomax α ^ = 2.3467 ,   λ ^ = 1.0934
Xgamma θ ^ = 1.3476
Gamma α ^ = 1.7266 ,   β ^ = 0.0839
Darna α ^ = 1.6346 ,   θ ^ = 2.7682
Power Lindley α ^ = 1.6235 ,   β ^ = 2.0145 ,
Power Darna α ^ = 2.856 ,   β ^ = 3.1624 ,   θ ^ = 0.0935
EP LD α ^ = 8.6321 ,   β ^ = 0.8562 ,   θ ^ = 3.4152
Table 6. Values of criteria statistics for model selection based on failure times data.
Table 6. Values of criteria statistics for model selection based on failure times data.
DistributionsY2WAK-S (p-Value)-NLLAICCAICBIC
PZD8.46840.29860.68440.1726 (0.9496)213.8102431.6204431.8757435.4444
ZD10.33440.41440.74680.2633 (0.9264)236.7156475.4312475.5145477.3432
Lomax10.98160.52820.76770.4476 (0.9103)253.4021510.8042511.0595514.6282
Xgamma12.63260.56330.94360.3466 (0.9098)241.0956484.1312484.2145486.1032
Gamma11.33920.44960.84340.5233 (0.9126)262.3364528.6728528.9281532.4968
Darna9.22040.36420.69390.1833 (0.9312)221.9302447.8604448.1157451.6844
Power Lindley12.85960.56980.94450.5325 (0.8965)268.9236 541.8472542.1025545.6710
Power Darna8.76590.32450.68620.1756 (0.9423)221.3546 448.7092449.2309454.4450
EP LD13.02510.57020.95120.5366 (0.8902)272.3256 550.6512551.1729556.3870
Table 7. Different parameter estimators for cancer patients’ data.
Table 7. Different parameter estimators for cancer patients’ data.
Method β θ KSp-Value
MLE0.43425.12670.03140.8434
WLS0.67646.36160.03860.8383
MPS0.69897.77680.04970.8198
CMEs0.42134.96740.04780.8266
LSE0.34174.16250.04120.8317
ADE0.78958.46160.05190.8078
ADEL-R0.76188.77860.05810.7984
ADEL-T0.82169.76330.05530.8010
Table 8. ML parameter estimates for the alternative distributions based on cancer patients’ data.
Table 8. ML parameter estimates for the alternative distributions based on cancer patients’ data.
DistributionsEstimates
PZD β ^ = 0.4342 , θ ^ = 5.1267
ZD θ ^ = 2.3461
Lomax α ^ = 1.6798 ,   λ ^ = 0.9436
Xgamma θ ^ = 1.9826
Gamma α ^ = 2.3639 ,   β ^ = 0.1364
Darna α ^ = 0.9671 ,   θ ^ = 1.3478
Power Lindley α ^ = 1.8963 ,   β ^ = 3.0596 ,
Power Darna α ^ = 3.0614 ,   β ^ = 2.9356 ,   θ ^ = 0.1796
EP LD α ^ = 5.6947 ,   β ^ = 0.4972 ,   θ ^ = 1.7693
Table 9. Values of criteria statistics for model selection based on cancer patients’ data.
Table 9. Values of criteria statistics for model selection based on cancer patients’ data.
DistributionsY2AWK-S (p-Value)-NLLAICBIC
PZD9.89060.67820.09340.0314 (0.8434)346.3315696.6630700.4870
ZD12.76480.73410.15980.0346 (0.7946) 352.9865707.9730709.8850
Darna10.96360.69140.12410.0324 (0.8126)349.6614703.3228707.1468
Xgamma13.92460.76330.17390.0377 (0.7767)354.6978711.3956713.3076
Gamma11.82460.71460.14330.0334 (0.8098)351.2365706.4730710.2970
Lomax13.01460.74240.16460.0356 (0.7833)353.6748711.3496715.1736
Power Lindley12.5362 0.81250.17620.0374 (0.7614)356.2350716.4700722.1740
Power Darna10.0236 0.68250.11430.0320 (0.8234)348.6320703.2640711.8200
EP LD14.5236 0.91250.17260.0402 (0.7522)362.2530730.5060739.0620
Table 10. Different parameter estimators with glass fiber data.
Table 10. Different parameter estimators with glass fiber data.
Method β θ KSp-Value
MLE0.3152 3.51620.04220.8496
WLS0.2365 5.02150.04850.8377
MPS0.5326 4.99360.05020.8222
CMEs0.4152 3.51620.05640.8143
LSE0.3465 4.03210.04460.8435
ADE0.2526 4.51360.05340.8102
ADEL-R0.0236 2.63510.05990.7933
ADEL-T0.1635 3.51620.05780.8089
Table 11. ML parameter estimates for the alternative distributions with glass fiber data.
Table 11. ML parameter estimates for the alternative distributions with glass fiber data.
DistributionsEstimates
PZD β ^ = 0.3152 , θ ^ = 3.5162
ZD θ ^ = 3.5614
Lomax α ^ = 2.1635 ,   λ ^ = 0.7263
Xgamma θ ^ = 1.5312
Gamma α ^ = 3.1259 ,   β ^ = 0.0985
Darna α ^ = 0.8132 ,   θ ^ = 1.7452
Power Lindley   α ^ = 1.3416 ,   β ^ = 1.9236
Power Darna α ^ = 1.285 ,   β ^ = 2.6351 ,   θ ^ = 0.0632
EP LD α ^ = 6.2451 ,   β ^ = 0.7325 ,   θ ^ = 4.201
Table 12. Values of criteria statistics for model selection with glass fiber data.
Table 12. Values of criteria statistics for model selection with glass fiber data.
DistributionsY2AWK-S (p-Value)-NLLAICBIC
PZD7.2365 0.56230.07850.0422 (0.8496)89.623183.246187.532
ZD8.6235 0.67480.09460.0512 (0.8348)93.516189.032191.175
Darna9.1532 0.71460.10160.0533 (0.8273)94.623193.246197.532
Xgamma12.130 0.81890.12030.0862 (0.7869)99.889201.778203.921
Gamma11.846 0.81330.11760.0785 (0.7902)98.495200.99205.276
Lomax11.526 0.80410.11350.0674 (0.7945)97.648199.296203.582
Power Lindley9.8462 0.76880.10780.0584 (0.8135)95.784195.568199.854
Power Darna7.5263 0.63250.08230.0476 (0.8412)91.253188.506194.935
EP LD10.236 0.79140.11060.0624 (0.8053)96.667199.334205.763
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Aidi, K.; Al-Omari, A.I.; Alsultan, R. The Power Zeghdoudi Distribution: Properties, Estimation, and Applications to Real Right-Censored Data. Appl. Sci. 2022, 12, 12081. https://doi.org/10.3390/app122312081

AMA Style

Aidi K, Al-Omari AI, Alsultan R. The Power Zeghdoudi Distribution: Properties, Estimation, and Applications to Real Right-Censored Data. Applied Sciences. 2022; 12(23):12081. https://doi.org/10.3390/app122312081

Chicago/Turabian Style

Aidi, Khaoula, Amer Ibrahim Al-Omari, and Rehab Alsultan. 2022. "The Power Zeghdoudi Distribution: Properties, Estimation, and Applications to Real Right-Censored Data" Applied Sciences 12, no. 23: 12081. https://doi.org/10.3390/app122312081

APA Style

Aidi, K., Al-Omari, A. I., & Alsultan, R. (2022). The Power Zeghdoudi Distribution: Properties, Estimation, and Applications to Real Right-Censored Data. Applied Sciences, 12(23), 12081. https://doi.org/10.3390/app122312081

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