Next Article in Journal
Third-Order Tensor Decorrelation Based on 3D FO-HKLT with Adaptive Directional Vectorization
Previous Article in Journal
Differential Transform Method for Axisymmetric Vibration Analysis of Circular Sandwich Plates with Viscoelastic Core
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Khalil New Generalized Weibull Distribution Based on Ranked Samples: Estimation, Mathematical Properties, and Application to COVID-19 Data

Department of Statistics and Operations Research, Faculty of Science, King Saud University, P.O. Box 2455, Riyadh 11451, Saudi Arabia
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(5), 853; https://doi.org/10.3390/sym14050853
Submission received: 21 March 2022 / Revised: 13 April 2022 / Accepted: 17 April 2022 / Published: 20 April 2022

Abstract

:
In this paper, a five-parameter distribution, Khalil’s new generalized Weibull distribution, is defined and studied in detail. Some mathematical and statistical functions are studied. The effects of shape parameters on skewness and kurtosis are studied. Extensions for density and distribution functions are provided. Estimation of the intended model parameters based on ranked samples is investigated. The behavior of the maximum likelihood estimators is examined using a Monte Carlo simulation. In order to predict unique symmetric and asymmetric patterns and illustrate the applicability and potential of the intended distribution, a COVID-19 dataset is analyzed. The goodness-of-fit results of the new generalized Weibull model of Khalil are compared with some other models. Finally, we make some concluding remarks.

1. Introduction

In order to achieve the goals that the researcher has set for his or her research, the researcher must collect data and information about that research. Since some data may be difficult to obtain or unacceptable due to cost, labor, and time constraints, he must choose the sampling method that will guarantee that he will achieve the research objectives with the least amount of time and expense—or, in many practical cases, it is not possible to obtain true measurements of the variables of interest, which is expensive and therefore a waste of time. To solve these problems, it is necessary to use a sampling method that ensures that time, effort, and cost are reduced in obtaining data.
The Ranked Set Sampling (RSS) method was proposed by McIntyre [1] as an inexpensive and effective method for estimating pasture yield. Although the RSS method is not parametric, many authors have used it to estimate parameters for many distributions and have demonstrated that its estimates are more effective than estimates based on simple random sampling and the same sample size. More detailed information can be found in [2,3,4,5]. Therefore, the aim of this work was to estimate the parameters of the Khalil new generalized Family-Weibull Distribution (KHGWD) considering ordered groups.
We now briefly introduce the RSS strategy used in the supplements. Consider an absolutely continuous random variable χ with the cumulative distribution function (CDF) and the probability density function (PDF). Then, a simple random sample of size n derived from the random variable χ is denoted by χ = { χ i : n , i = 1 , , n } . Suppose further that a random sample of size n 2 is selected and randomly divided into n groups of equal size. Then, RSS is observed according to the following pattern:
(1) χ ( 1 : n ) 1 ̲ χ ( 2 : n ) 1 χ ( n : n ) 1 χ 1 , 1 = χ ( 1 : n ) 1
(2) χ ( 1 : n ) 2 χ ( 2 : n ) 2 ̲ χ ( n : n ) 2 χ 2 , 2 = χ ( 2 : n ) 2
. . . . .
. . . . .
. . . . .
(n) χ ( 1 : n ) n χ ( 2 : n ) n χ ( n : n ) r ̲ χ n , n = χ ( n : n ) n
The RSS vector of observations is given by χ j , i , i = 1 , , n , j = 1 , , n , where χ j , i is the statistic of order ith in the group jth based on a given simple random sample of size n. Then, as is well known, the PDF of χ j , i is
f j , i ( x ) = n ! ( i 1 ) ! ( n i ) ! f ( x ) F ( x ) i 1 1 F ( x ) n i , j , i = 1 , , n .
For more details, see Arnold et al. [6].
Next, Najma et al. [7] proposes a new method to extend the family of lifetime distributions. The method is called Khalil new generalized family of distributions. For any baseline CDF G ( x ) and a PDF g ( x ) , the CDF and the PDF of the new generalized family of Khalil distributions are respectively given by
F ( x ; α , β ) = e α G ( x ) β 1 e α 1 , x R ; α , β > 0 ,
and
f ( x ; α , β ) = α β e G ( x ) β α g ( x ) G ( x ) β 1 1 e α , x R ; α , β > 0 ,
where α and β are the scale and shape parameters, respectively.
In 1939, Swedish scientist Waloddi Weibull established the Weibull distribution in a study of the breaking strength of instruments. The Weibull distribution is one of the most commonly used failure models. The Weibull distribution is used to simulate many probabilistic applications; this is due to its unique symmetric and asymmetric patterns. The distribution has several desirable properties, acceptable physical interpretations, and the ability to fit the failure rates of various systems, whether those rates are high, low, or constant.
Let X be a random variable (R.V.) that follows the Weibull distribution with three parameters ( λ , μ , ν ) , then its CDF, denoted by G ( x ; λ , μ , ν ) , is given by
G ( x ; λ , μ , ν ) = 1 e x ν μ λ , x ν ; λ , μ , ν > 0 .
Here, ν > 0 is the location parameter, x ν , λ > 0 , and μ > 0 are the shape and the scale parameter, respectively. The corresponding PDF, denoted by g ( x ; λ , μ , ν ) , is given by
g ( x ; λ , μ , ν ) = λ μ x ν μ λ 1 e x ν μ λ , x ν ; λ , μ , ν > 0 .

2. The Khalil New Generalized Family-Weibull Distribution (KHGWD)

Take G ( x ) and g ( x ) in Equations (2) and (3) as G ( x ; λ , μ , ν ) of Equation (4) and g ( x ; λ , μ , ν ) of Equation (5), respectively. The CDF and the PDF of KHGWD are respectively given by
F ( x ; α , β , λ , μ , ν ) = 1 e α 1 e x ν μ λ β 1 e α ,
and
f ( x ; α , β , λ , μ , ν ) = α β λ x ν μ λ 1 1 e x ν μ λ β 1 e x ν μ λ α 1 e x ν μ λ β μ 1 e α ,
where x ν , α > 0 , and μ > 0 are two scale parameters, β > 0 and λ > 0 are two shape parameters, and ν > 0 is a location parameter.
The survival function and hazard rate function of time t via KHGWD are, respectively, given be
S ( t ; α , β , λ , μ , ν ) = 1 1 e α 1 e t ν μ λ β 1 e α ,
and
H ( t ; α , β , λ , μ , ν ) = α β λ t ν μ λ 1 1 e t ν μ λ β 1 e t ν μ λ α 1 e t ν μ λ β μ 1 e α 1 1 e α 1 e t ν μ λ β 1 e α .
In what follows, an R.V. X with the KHGWD (7) is denoted by X KHGWD ( α , β , λ , μ , ν ) . Figure 1, Figure 2 and Figure 3 show various survival functions (SF), PDF, and hazard rate functions (HRF) for the KHGWD ( α , β , λ , μ , ν ) . The right panel for KHGWD(1.10, 2.80, 0.50, 1.40, 0.20), KHGWD(0.30, 0.80, 0.60, 2.90, 0.12), KHGWD(9.02, 6.90, 0.50, 1.90, 0.20), KHGWD(1.2, 3.90, 0.65, 0.90, 0.20), and KHGWD(0.20, 3.90, 0.30, 1.10, 0.01), while the left panel for KHGWD(5.10, 0.50, 6.50, 1.40, 1.20) and KHGWD(0.10, 0.50, 6.50, 1.40, 1.20), KHGWD(5.10, 0.80, 6.50, 1.40, 1.20), KHGWD(5.10, 0.50, 2.80, 1.40, 1.20), and KHGWD(5.10, 0.50, 6.50, 2.40, 1.20).
The model is sometimes very flexible. It seems to approach the bell curve with some torsion, as seen in the right panel of Figure 2. At other times, it seems to have strong tails, as seen in the left panel of Figure 2, which depends on the particular values of the parameters. The left panel of Figure 2 also shows that the proposed model has heavy tails when the parameters are increased. Based on the behavior of the proposed model shown in Figure 2, it is a good candidate for modeling semi-normal data (right part) and data with heavy tails (left part) in various financial, industrial, medical, and global epidemiological applications which have the same behavior. However, from the plots in Figure 3, it is clear that the KHGWD has unimodal and increasing failure rate functions. The unimodal and increasing failure rate functions are another superiority of the proposed model along with the heavy-tailed behavior. Therefore, the proposed model is suitable for modeling COVID-19.

3. Mathematical Properties

This section is devoted to deriving some mathematical properties of the KHGWD ( α , β , λ , μ , ν ) .

3.1. Quantile Function

To generate random variables by Monte Carlo simulation, the quantile function of the distribution is required. Assuming p U n i f o r m ( 0 , 1 ) , we solved the following equation for the quantile function Λ ( p ) :
p = 1 e α ( F ( Λ ( p ) ) ) β 1 e α .
Letting y = F ( Λ ( p ) ) , we have
p = 1 e α ( y ) β 1 e α .
By solving for y, we have
F ( Λ ( p ) ) = log 1 p 1 e α α 1 β .
Thus,
Λ ( p ) = F 1 { log 1 p 1 e α α 1 β } ,
where F 1 is the quantile of the baseline distribution KHGWD ( α , β , λ , μ , ν ) . Inverting F ( x ) = p in (6), we can write
Λ KHGWD ( x , α , β , λ , μ , ν ) = ν + μ log 1 log 1 1 e α x α 1 β 1 λ .
By setting x as a uniform R.V. in the unit interval (0, 1), we can also use (14) for simulating KHGWD ( α , β , λ , μ , ν ) R.V.s. Figure 4 plots different quantile functions for the KHGWD ( α , β , λ , μ , ν ) , and particularly for (1) β = 2.3 ; λ = 0.2 , μ = 0.2 , ν = 3.2 . (2) α = 4.3 , λ = 0.2 , μ = 0.2 , ν = 3.2 . (3) α = 4.3 , β = 0.3 , μ = 0.2 , ν = 3.2 . (4) α = 4.3 , β = 0.3 , λ = 0.2 , ν = 3.2 . (5) α = 0.3 , β = 0.3 , λ = 1.2 , μ = 0.2 , respectively. Some numerical values of the quantile measure are provided in Table 1.
In addition, the effects of shape parameters on skewness and kurtosis can be determined using quantile measures. We obtain skewness and kurtosis measures of KHGWD ( α , β , λ , μ , ν ) . The skewness (SK) (see Bowley [8]) of X is given by
SK ( α , β , λ , μ , ν ) = 2 Λ 1 / 2 Λ 3 / 4 Λ 1 / 4 Λ 1 / 4 Λ 3 / 4 ,
and the Kurtosis (K) (see Moor [9]) is given by
K ( α , β , λ , μ , ν ) = Λ 1 / 8 Λ 3 / 8 + Λ 5 / 8 Λ 7 / 8 Λ 2 / 8 Λ 6 / 8 .
Figure 5 shows SK ( α , β , 0.2 , 0.2 , 3.2 ) and K ( α , β , 0.2 , 0.2 , 3.2 ) for different values of α and β . Figure 6, on the other hand, shows SK ( 1.3 , 5.3 , λ , μ , 3.2 ) and K ( 1.3 , 5.3 , λ , μ , 3.2 ) for different values of λ and μ .

3.2. The Expansion for KHGWD Density Function

Using the general binomial and the power series expansion, we have
1 e x ν μ λ β 1 = κ = 0 β 1 κ ( 1 ) κ e κ x ν μ λ = κ = 0 β 1 κ ( 1 ) κ h = 0 ( 1 ) h h ! x ν μ λ h .
Next, we can also write
e x ν μ λ α 1 e x ν μ λ β = ρ = 0 ( 1 ) ρ ρ ! x ν μ λ + α 1 e x ν μ λ β ρ = ρ = 0 ( 1 ) ρ ρ ! q = 0 ρ q x ν μ λ q α ρ q 1 e x ν μ λ β ( ρ q ) ,
Using (15) and (16), the PDF of the expanded KHGWD ( α , β , λ , μ , ν ) is given by
f ( x ; α , β , λ , μ , ν ) = α β λ μ 1 e α κ = 0 β 1 κ h = 0 ρ = 0 q = 0 α ρ q × τ = 0 δ = 0 Ψ κ , h , ρ , q , τ , δ x ν μ λ ( h + q + δ + 1 ) 1 ,
where
Ψ κ , h , ρ , q , τ , δ = ρ q β ( ρ q ) τ ( 1 ) κ + h + ρ + τ + δ τ δ h ! ρ ! δ ! .

3.3. The Expansion for the KHGWD Distribution Function

The CDF (15) of the expanded KHGWD ( α , β , λ , μ , ν ) can be written as
F ( x ; α , β , λ , μ , ν ) = 1 h = 0 ( 1 ) h h ! α h 1 e x ν μ λ h β 1 e α 1 = 1 h = 0 α h τ = 0 h β τ ρ = 0 Ψ h , τ , ρ x ν μ λ ρ 1 e α 1 ,
where
Ψ h , τ , ρ = ( 1 ) h + τ + ρ τ ρ h ! ρ ! .

4. Estimation of the Parameters Based on the Ranked Set Samples

Let { X j , i , i = 1 , , m , j = 1 , , r } be the observed sample, where X j , i is the ith order statistic in the jth group. The likelihood function based on the KHGWD ( α , β , λ , μ , ν ) , is given by
l = m ! α β 1 e α r m j = 1 r i = 1 m e α G j , i β e α G j , i β 1 e α 1 i 1 1 e α G j , i β 1 e α 1 m i g j , i G j , i β 1 ( i 1 ) ! ( m i ) ! ,
where G j , i = G j , i ( x j , i , λ , μ , ν ) and g j , i = g j , i ( x j , i , λ , μ , ν ) are the baseline CDF (6) and PDF (7), respectively. The corresponding log-likelihood function is given by
L = K 1 + r m log [ α β ] n log 1 e α j = 1 r i = 1 m α G j , i β + j = 1 r i = 1 m ( i 1 ) log e α G j , i β 1 e α 1 + ( m i ) log 1 e α G j , i β 1 e α 1 + j = 1 r i = 1 m log g j , i + ( β 1 ) log G j , i ,
where K 1 is constant. The first partial derivatives of log-likelihood (20) with respect to α , β , λ , μ , ν , respectively, are given by
L α = n 1 α e α 1 e α j = 1 r i = 1 m G j , i β + j = 1 r i = 1 m ( i 1 ) e α G j , i β 1 ( m i ) e α e α G j , i β e α e α G j , i β 1 e α 1 e α G j , i β G j , i β ,
L β = n β + α j = 1 r i = 1 m G j , i β 1 log G j , i 1 + m i e α e α G j , i β i 1 e α G j , i β 1 e α G j , i β ,
L λ = j = 1 r i = 1 m g j , i ( λ ) g j , i + G j , i ( λ ) β 1 G j , i α β G j , i β 1 α β j = 1 r i = 1 m G j , i β 1 G j , i ( λ ) e α G j , i β i 1 e α G j , i β 1 + m i e α e α G j , i β ,
L μ = j = 1 r i = 1 m g j , i ( μ ) g j , i + G j , i ( μ ) β 1 G j , i α β G j , i β 1 α β j = 1 r i = 1 m G j , i β 1 G j , i ( μ ) e α G j , i β i 1 e α G j , i β 1 + m i e α e α G j , i β ,
and
L ν = j = 1 r i = 1 m g j , i ( ν ) g j , i + G j , i ( ν ) β 1 G j , i α β G j , i β 1 α β j = 1 r i = 1 m G j , i β 1 G j , i ( ν ) e α G j , i β i 1 e α G j , i β 1 + m i e α e α G j , i β ,
where
G j , i ( λ ) = G j , i λ = x j , i ν μ λ e x j , i ν μ λ log x j , i ν μ ,
G j , i ( μ ) = G j , i μ = λ μ x j , i ν μ λ e x j , i ν μ λ ,
G j , i ( ν ) = G j , i ν = λ μ x j , i ν μ λ 1 e x j , i ν μ λ ,
g j , i ( λ ) = g j , i λ = μ x j , i ν μ λ 1 1 + λ log x j , i ν μ x j , i ν μ λ 1 e x j , i ν μ λ ,
g j , i ( μ ) = g j , i μ = λ μ 2 x j , i ν μ λ 1 x j , i ν μ λ 1 e x j , i ν μ λ ,
and
g j , i ( ν ) = g j , i ν = λ μ 2 x j , i ν μ λ 2 1 + λ x j , i ν μ λ 1 e x j , i ν μ λ .
The maximum likelihood estimators α ^ M L , β ^ M L , λ ^ M L , μ ^ M L , and ν ^ M L of the KHGWD parameters are the solutions of the nonlinear Equations (21)–(25) for r = m = n . Characteristically, they can be solved by fixed-point iteration methods or the Newtonian approach.

5. Monte Carlo Simulation Study

This section is concerned with evaluating the performance of the maximum likelihood estimators of KHGWD ( α , β , λ , μ , ν ) through a Monte Carlo simulation study. The R program through the optimum function can used to compute the simulation and application results. The Weibull location parameter ν was set to zero to simplify the calculations.
The simulation of the KHGWD model is performed for two parameter sets:
  • Set 1: α = 3.2 , β = 1.5 , λ = 3.8 , μ = 2.1 , ν = 0 .
  • Set 2: α = 4.3 , β = 1.6 , λ = 3.5 , μ = 1.8 , ν = 0 .
The simulation study is performed as follows.
  • Random samples of sizes n = 25 , 100 , 225 , 400 , 625 , 900 are generated from KHGWD and are randomly divided into r groups of equal size m, where r = m = 5 , 10 , 15 , 20 , 25 , and 30, respectively.
  • The model parameters have been estimated via the maximum likelihood method.
  • Five-thousand repetitions are made to calculate these estimators’ biases, absolute biases, and mean square errors (MSEs).
  • The formulas for obtaining the estimate, biases, and MSEs are given respectively, by
    α ^ = 1 5000 i = 1 5000 α ^ i ,
    Bias α ^ = 1 5000 i = 1 5000 α ^ i α ,
    and
    MSE α ^ = 1 5000 i = 1 5000 α ^ i α 2 ,
    respectively.
  • Step (4) is also repeated for the parameters β , λ , and μ .
The simulation results of the KHGWD for Set 1 and Set 2 are presented in Table 2 and Table 3, respectively. Figure 7 and Figure 8 displays graphically the results provided in Table 2 and Table 3, respectively.
Figure 7 and Figure 8 illustrate the simulation results for the above measures. These plots show that increasing the sample size n leads to a reduction in the estimated biases. In addition, increasing the sample size n leads to a decrease in the estimated MSEs, which approaches zero as n increases. These results demonstrate both the efficiency and consistency properties of the MLEs.

6. COVID-19 Data

The COVID-19 dataset of Hong and Li [10] is considered an application of the KHGWD model and is randomly divided into seven groups of similar seven observations. Table 4 shows the descriptive statistics of the proposed COVID-19 data. The boxplot and Q-Q plot are shown in Figure 9. By performing the maximum likelihood estimators of KHGWD ( α , β , λ , μ , ν ) by Monte Carlo simulation study for ν = 0 , the estimated parameters are α ^ = 2.4222994 , β ^ = 0.1311977 , λ ^ = 5.3593414 , and μ ^ = 20.1015298 . Figure 10: Plots of the fitted PDF of KHGWD and CDF. Figure 11: PP plot and Kaplan–Meier survival function of KHGWD.
Table 5 compares the KHGWD based on one-sample Kolmogorov–Smirnov test. The goodness-of-fit results of the KHGWD model are compared with some other models, including the Khalil generalized exponential distribution (KHGEXP), the Khalil generalized gamma distribution (KHGGamma), and Weibull distribution. The distribution functions of these competitive distributions are given by
(1) KHGEXP distribution:
F ( x ; α , β , μ ) = 1 e α 1 e x μ β 1 e α , x 0 , α , β , μ > 0 .
(2) KHGGamma distribution:
F ( x ; α , β , λ , μ ) = e α 1 Γ λ γ ( λ , μ x ) β 1 e α 1 , x 0 , α , β , λ , μ > 0 ,
where γ ( λ , μ x ) is the lower incomplete gamma function.
(3) Weibull distribution:
F ( x ; λ , μ ) = 1 e x μ λ , x 0 ; λ , μ > 0 .
Based on the results presented in Table 5, we see that KHGWD is a good competitor among the competing models for modeling the COVID-19 data.

7. Conclusions

We introduced the new generalized Weibull distribution of Khalil with five parameters. The model has a high degree of flexibility to fit the appropriate data. The model is unimodal, has a strong tail-heavy behavior, and exhibits increasing failure rate functions. We have studied some mathematical and statistical functions and the effects of shape parameters on skewness and kurtosis. We have provided density and distribution functions in an extended form. Based on ranked samples, the maximum likelihood estimators for the intended model parameters and a Monte Carlo simulation study are provided. A COVID-19 dataset is analyzed to illustrate the applicability and potential of the intended distribution. The biases and mean square errors decrease as the sample size increases. It is clear that the proposed model fits the estimated PDF and CDF plots well. The COVID-19 dataset has a strong deflection that is skewed to the right in the boxplot. The proposed model fits the Kaplan–Meier survival plot very well. The new generalized Weibull distribution of Khalil, based on the Kolmogorov–Smirnov one-sample test, provides a better fit than other competing models. The results show that the model is considered ideal for modeling COVID-19.

Author Contributions

Formal analysis, W.E.; Funding acquisition, Y.A.T.; Methodology, W.E.; Software, W.E.; Supervision, W.E.; Writing—original draft, W.E.; Writing—review and editing, Y.A.T. All authors have read and agreed to the published version of the manuscript.

Funding

The study was funded by Researchers Supporting Project number (RSP2022R488), King Saud University, Riyadh, Saudi Arabia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are available on request ([email protected]).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. McIntyre, G.A. A method for unbiased selective sampling, using ranked sets. Aust. J. Agric. Res. 1952, 3, 385–390. [Google Scholar] [CrossRef]
  2. Frey, J. A note on Fisher information and imperfect ranked-set sampling. Commun. Stat.-Theory Methods 2014, 43, 2726–2733. [Google Scholar] [CrossRef]
  3. Park, S.; Lim, J. On the effect of imperfect ranking on the amount of Fisher information in ranked set samples. Commun. Stat.-Theory Methods 2012, 413, 3608–3620. [Google Scholar] [CrossRef]
  4. Chen, Z.; Bai, Z.; Sinha, B. Ranked Set Sampling: Theory and Applications; Springer: New York, NY, USA, 2013. [Google Scholar]
  5. Tahmasebi, S.; Longobardi, M.; Kazemi, M.R.; Alizadeh, M. Cumulative Tsallis entropy for maximum ranked set sampling with unequal samples. Phys. A Stat. Mech. Appl. 2020, 556, 124763. [Google Scholar] [CrossRef]
  6. Arnold, B.C.; Balakrishnan, N.; Nagaraja, H.N. A First Course in Order Statistics; John Wiley & Sons: New York, NY, USA, 1992. [Google Scholar]
  7. Najma, S.; Alamgir, K.; Wali, K.; Habib, S.; Pijitra, J.; Thammarat, P. A Novel Generalized Family of Distributions for Engineering and Life Sciences Data Applications. Math. Probl. Eng. 2021, 2021, 9949999. [Google Scholar]
  8. Bowley, A.L. Elements of Statistics, 4th ed.; Charles Scribner’s Sons: New York, NY, USA, 1920. [Google Scholar]
  9. Moors, J.J.A. The meaning of kurtosis: Darlington re-examined. Am. Stat. 1986, 40, 283–284. [Google Scholar]
  10. Hong, H.G.; Li, Y. Estimation of time-varying reproduction numbers underlying epidemiological processes: A new statistical tool for the COVID-19 pandemic. PLoS ONE 2020, 15, e0236464. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Different SF for the KHGWD ( α , β , λ , μ , ν ) .
Figure 1. Different SF for the KHGWD ( α , β , λ , μ , ν ) .
Symmetry 14 00853 g001
Figure 2. Different PDF for the KHGWD ( α , β , λ , μ , ν ) .
Figure 2. Different PDF for the KHGWD ( α , β , λ , μ , ν ) .
Symmetry 14 00853 g002
Figure 3. Different HRF for the KHGWD ( α , β , λ , μ , ν ) .
Figure 3. Different HRF for the KHGWD ( α , β , λ , μ , ν ) .
Symmetry 14 00853 g003
Figure 4. Different quantile functions for the KHGWD ( α , β , λ , μ , ν ) .
Figure 4. Different quantile functions for the KHGWD ( α , β , λ , μ , ν ) .
Symmetry 14 00853 g004
Figure 5. Plots for the SK ( α , β , 0.2 , 0.2 , 3.2 ) and K ( α , β , 0.2 , 0.2 , 3.2 ) .
Figure 5. Plots for the SK ( α , β , 0.2 , 0.2 , 3.2 ) and K ( α , β , 0.2 , 0.2 , 3.2 ) .
Symmetry 14 00853 g005
Figure 6. Plots for the SK ( 1.3 , 5.3 , λ , μ , 3.2 ) and K ( 1.3 , 5.3 , λ , μ , 3.2 ) .
Figure 6. Plots for the SK ( 1.3 , 5.3 , λ , μ , 3.2 ) and K ( 1.3 , 5.3 , λ , μ , 3.2 ) .
Symmetry 14 00853 g006
Figure 7. Plot of the simulation results of the KHGWD for Set 1.
Figure 7. Plot of the simulation results of the KHGWD for Set 1.
Symmetry 14 00853 g007
Figure 8. Plot of the simulation results of the KHGWD for Set 2.
Figure 8. Plot of the simulation results of the KHGWD for Set 2.
Symmetry 14 00853 g008
Figure 9. The boxplot and Q-Q plot of the COVID-19 data.
Figure 9. The boxplot and Q-Q plot of the COVID-19 data.
Symmetry 14 00853 g009
Figure 10. The fitted PDF and CDF of the KHGWD.
Figure 10. The fitted PDF and CDF of the KHGWD.
Symmetry 14 00853 g010
Figure 11. The PP plot and the Kaplan–Meier survival function of the KHGWD.
Figure 11. The PP plot and the Kaplan–Meier survival function of the KHGWD.
Symmetry 14 00853 g011
Table 1. Some quantile values for α = 1.3 , β = 5.3 , λ = 0.2 , μ = 0.2 , and ν = 3.2 .
Table 1. Some quantile values for α = 1.3 , β = 5.3 , λ = 0.2 , μ = 0.2 , and ν = 3.2 .
x Λ ( x )
0.13.30459
0.23.54091
0.33.97569
0.44.74677
0.56.13344
0.68.75958
0.714.2517
0.828.0657
0.980.4795
Table 2. Simulation results of the KHGWD for Set 1.
Table 2. Simulation results of the KHGWD for Set 1.
n α ^ β ^ λ ^ μ ^
2.072.5514.0441.805
254.0874.5982.9460.255
−1.131.0510.244−0.295
2.5862.1043.7871.917
1003.32.0761.6920.165
−0.6140.604−0.013−0.183
2.8381.93.6751.973
2252.9491.011.0190.12
−0.3610.4−0.124−0.127
2.9881.7923.6452.007
4002.5230.540.6810.089
−0.2120.292−0.155−0.093
3.0531.7223.6452.028
6252.1990.3170.4690.069
−0.1470.222−0.155−0.072
3.1151.7023.642.04
9002.0530.2710.3970.061
−0.0850.202−0.16−0.06
The first, second, and third rows represent the estimate, MSE, and Bias, respectively.
Table 3. Simulation results of the KHGWD for Set 2.
Table 3. Simulation results of the KHGWD for Set 2.
n α ^ β ^ λ ^ μ ^
3.4762.5144.3961.609
256.1694.5515.5540.26
−0.8240.9140.895−0.1907
3.1822.1823.7991.606
1004.9082.2722.6020.155
−1.1180.5820.299−0.194
3.5111.9413.6141.663
2253.8121.0721.3240.095
−0.7890.3410.114−0.137
3.7121.8443.5121.693
4003.2570.5780.7680.066
−0.5880.2440.012−0.107
3.8831.7993.4531.716
6252.6440.3610.4710.047
−0.4170.199−0.047−0.084
4.0541.7653.4431.737
9002.2740.270.3860.035
−0.2460.165−0.057−0.063
The first, second, and third rows represent the estimate, MSE, and Bias, respectively.
Table 4. Descriptive statistics of the COVID-19 data.
Table 4. Descriptive statistics of the COVID-19 data.
Min.1st Qu.MedianMean3rd Qu.Max
0.0540.7043.0794.7876.74320.083
Table 5. Relative quality of the KHGWD vs. competing models.
Table 5. Relative quality of the KHGWD vs. competing models.
Modelp-Value
KHGWD0.4347
KHGEXP0.3776
KHGGamma0.4056
Weibull0.3776
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Emam, W.; Tashkandy, Y.A. Khalil New Generalized Weibull Distribution Based on Ranked Samples: Estimation, Mathematical Properties, and Application to COVID-19 Data. Symmetry 2022, 14, 853. https://doi.org/10.3390/sym14050853

AMA Style

Emam W, Tashkandy YA. Khalil New Generalized Weibull Distribution Based on Ranked Samples: Estimation, Mathematical Properties, and Application to COVID-19 Data. Symmetry. 2022; 14(5):853. https://doi.org/10.3390/sym14050853

Chicago/Turabian Style

Emam, Walid, and Yusra A. Tashkandy. 2022. "Khalil New Generalized Weibull Distribution Based on Ranked Samples: Estimation, Mathematical Properties, and Application to COVID-19 Data" Symmetry 14, no. 5: 853. https://doi.org/10.3390/sym14050853

APA Style

Emam, W., & Tashkandy, Y. A. (2022). Khalil New Generalized Weibull Distribution Based on Ranked Samples: Estimation, Mathematical Properties, and Application to COVID-19 Data. Symmetry, 14(5), 853. https://doi.org/10.3390/sym14050853

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