Abstract
In this paper, we developed a novel superior distribution, demonstrated and derived its mathematical features, and assessed its fuzzy reliability function. The novel distribution has numerous advantages, including the fact that its CDf and PDf have a closed shape, making it particularly relevant in many domains of data science. We used both conventional and Bayesian approaches to make various sorts of estimations. A simulation research was carried out to investigate the performance of the classical and Bayesian estimators. Finally, we fitted a COVID-19 mortality real data set to the suggested distribution in order to compare its efficiency to that of its rivals.
1. Introduction
Nowadays we face a huge waves of viruses such as COVID-19, which attracted a great attention in the last two years. However, as statisticians we must play an important role in understanding and modelling the COVID-19 infections, so we had to make a statistical model that is capable of fitting and modelling them, whether it is continuous or discrete random variable. In this paper, we made a lot of effort to choose an appropriate model which is superior in fitting the COVID-19 infections in many countries. At last, we settled to introduce a new superior distribution as mixture of other ones to overcome the defects of the baseline distribution.
Over the years and decades, modelling for lifetime distributions has garnered considerable attention, and their popularity has increased as the relevance of modelling phenomena and pandemics has grown. Researchers in distribution theory choose to model the data whether by adding a new parameter that increases the flexibility of the distribution of interest or by creating a new distribution family. Modelling is extremely useful in a wide variety of sectors, including commerce, economics, dependability, and health research. See [1] for more reading. In the next paragraph we will discuss in details the new proposed distribution.
We developed a statistical model for COVID-19 mortality data in Saudi Arabia in this study; for further information about the other models, refer to Kumar [2], Khakharia et al. [3], Shafiq et al. [4], Sindhu et al. [5], Wang [6], Lalmuanawma et al. [7], and Bullock et al. [8].
The new proposed model is a combination between the Kumaraswamy distribution’s (Ku) and the modified Kies distribution. The CDF of Ku distribution and PDF with parameters were put forth by [9] in the following manner:
Many authors worked on extensions of the Kies distribution as an example; see Kumar and Dharmaja [10] who worked on the exponentiated Kies distribution. Also, Dey et al. [11] worked on the modified Kies distribution when the sample they worked on censored data specially under type II progressive censored sample. Al-Babtain et al. [12] developed a different distribution family based on the modified Kies (MK) family. Let us say that is the baseline CDF depending on a parameter vector ; (3) defines the CDF of the MK family. Kumar and Dharmaja examined the characteristics of the exponentiated Kies distribution as we can see in reference [10]. One of the authors who have created a MK distribution for a type II progressive censored sample; in addition, the estimation of the distribution parameters has been obtained (see Dey et al. [11] for more information). Referring to Al-Babtainet al. [12], if is the baseline CDF , then is the MK family’s CDF that may be expressed as follows:where is the parameter vector and is parameter vector of .
The PDF of (3) is as follows:
Using the distributions above, the two-parameter modified Kies Kumaraswamy (MKKu) is produced in this paper. The MKKu has a number of interesting properties. The suggested MKKu distribution has an extremely adaptable PDF since it may be positive skewed, negative skewed, or symmetric, which allows for more flexibility for tail. It is capable of simulating declining, increasing, bathtub, and reverse-J hazard rates. Also, the suggested distribution has a precise closed-form CDF and is extremely easy to manage. These advantages make the distribution an attractive choice for applications in a variety of fields, including biological life testing, dependability, and actuarial data. This article showed one real-world data application and demonstrated via modelling that the novel distribution is a great rival to many well-known and conventional distributions with the parameters of the scale and shape such as type II power Topp-Leone inverse exponential (Bantan et al. [13]), Topp-Leone generalized exponential (Kunjiratanachot et al. [14]), Kumaraswamy, beta, Gompertz Lomax (Oguntunde et al. [15]), alpha power inverted TL (Ibrahim et al. [16]), and MK exponential (Al-Babtain et al. [12]) distributions.
The paper is arranged and organized as follows. In Section 2, the MKKu distribution is obtained. In Section 3, some mathematical characteristics of the MKKu distribution are discussed. In Section 4, we establish an MKKu distribution estimate technique. In Section 5, a fuzzy reliability function was obtained. In Section 6, we present the simulation study for the estimation techniques and concluded our major findings in this section. In Section 7, we acquired applications of real-world data analysis. Section 8 summarises and concludes the article. Future research will include the development of an extension for the bivariate MK inverted Topp-Leone based on copula as in [17]. We propose to investigate a unique application of the MKITL distribution using a censored sample; see [18] for further information.
2. MKKu Distribution
By employing the Kumaraswamy distribution with these parameters and the CDF specified (for 0 < x < 1), (1) and (2) are used. We construct the CDF for MKKu distribution by putting the Ku distribution’s CDF into (3).
The corresponding PDF iswhere is parameter vector .
The survival function of the MKKu distribution is represented by
The function of hazard rate (HR) of the MKKu distribution is given by
Figures 1 and 2 provide specific plots of the MKKu distribution for the values indicated for and in the equation of the MKKu distribution. It can be seen by the diagrams provided in Figure 2 that the MKKu distribution HR function may be increased, lowered, and shaped like a bathtub. In contrast to a Ku distribution, which is a bad model for the data and phenomena that indicate growing, decreasing forms, and failure rates of the bathtub, the MKKu distribution is more adaptable in terms of evaluating lifespan data than its rival; this is one of its advantages.


3. The Suggested Distribution Statistical Properties
3.1. A General Expansion of the MKKu Density
In this article, we provide a linear representation for the MK family and utilise it to construct a suitable linear representation for the MKKu distribution.
The last equation of the MKKu distribution may be reformulated using the PDF and CDF of the Ku distribution.
Equation (10) denotes the Ku density with parameter .
3.2. Quantile for the MKKu Distribution
In this part of the paper, we drive the formula of the MKKu distribution’s quantile function; for example, let us say that is derived by inverting (5) as follows:
The first quartile, denoted by the symbol , the second quartile, denoted by the symbol , and the third quartile, denoted by the symbol , are produced by putting U = 0.25, 05, 0.75, accordingly in (11).
3.3. The Proposed Distribution Related Moments
The rth of the distribution can be find by using (10) as
4. Estimation Techniques for Distribution Parameters
Using Bayesian and non-Bayesian estimate approaches, we are able to solve the estimation issue for the MKKu distribution parameters. The maximum likelihood estimators (MLE), maximum product of spacings, Anderson–Darling, Cramér–von-Misse, and Bayesian approaches based on the squared error loss function (SELF) are used. These non-Bayesian methods have been discussed to estimate parameters for different models as per [19–21], while for Bayesian estimation, see [22–24].
4.1. The MLE Estimator
Suppose is a random sample of any size, then likelihood function of the proposed distribution can be written as follows:
For such MKKu distribution, the log-likelihood formula is represented by
We will differentiate with respect to every parameter in the distribution
We can evaluate the estimates of the three parameters , , and by determining the greatest value for (14) and then by computing its first derivative in terms of , , and . Using the Newton-Rapshon approach, the R program may be used to maximize the log-likelihood and yield the MLE.
4.2. Maximum Product of Spacings (MPS) Method
MPS approach was established by Cheng and Amin [25] and has been employed in censored applications as [26–28]. Assume . Given a random sample of size , the MKKu distribution’s uniform spacing may be written as indicate the uniform spacings, , and . The MPS estimators (MPSE) of the MKKu parameters can be found by maximizingwith respect to , and .
The MPSE of the MKKu parameters can found by solving
4.3. Cramer–von-Mises (CVME) Method
This method may be calculated using the difference between both the CDF and empirical distribution function estimations by Luceno [29]. The CVME values for the MKKu parameters are produced by minimizing the next function in terms of , and :
Further, the CVME follow by solving the next nonlinear equations
Solve these equations to get the estimates (20)–(22).
4.4. Anderson–Darling Method
Another sort of minimal distance estimators is the (ADE) Anderson–Darling estimators, minimizing the ADE of the MKKu parameters with respect to , and . Additionally, the ADE may be produced by solving nonlinear equations, which yields the ADE of the MKKu parameters estimates.
4.5. Bayesian Estimation Method
This section contains Bayesian estimations of , and . We assume that , and have gamma priors; also, the gamma prior density was first used. The following are the gamma priors for distribution parameters:
As we know the two priors are independent of each other, so the joint prior of and is obtained as
To get acceptable and superior values for the independent joint prior’s hyperparameters, we may employ the MLE method’s estimate and variance-covariance matrix. The calculated hyperparameters may be represented by equating the mean and variance of gamma priors:where is the number of iteration.
By multiplying (13) by (25) and making some simplifications, the posterior distribution is formed as follows:
As , the Bayesian estimations of the distribution parameter using the squared error loss function (SELF) is
It is obvious that SELF estimations of in (28) is considered as the division of more than one integration over each other. As the multiple integrals are hard to deal with, so we used Mathematica 11 to apply an approximation technique, which is very useful in solving these kinds of integration. So in this context, we used the Markov Chain Monte Carlo (MCMC) technique to get an approximate value of integrals. The important methods of the MCMC technique is the Metropolis-Hastings (MH) algorithm. The following approach is used to produce random samples from conditional posterior densities of the MKKu distribution:
5. Fuzzy Reliability
The concepts of dependability and HR function are probabilities that describe the life time (a random variable) from the start of a failure to the point where we need to adjust the manufacturing process. The use of dependability and HR functions has been broadened now that fuzzy factors have been included to these functions. The lifespan probability of the system’s components is treated as accurate numbers in classic reliability models. In the real world, however, this accuracy of system lifetimes is not true since the values of system parameters acquired via experiment, faulty measurement, human judgement, or guess are all subject to some uncertainty [30]. As a result, incorporating the idea of fuzziness will be more acceptable. As a result, while studying system behaviour or addressing system reliability, it will be more suitable to include the idea of fuzziness. The membership functions are utilised to describe the level of fuzziness associated with the lifespan data or system parameters in order to generate the fuzzy parameters [31]. The fuzzy set theory may be used to explain the real world in a realistic and practical way. Zadeh [32] was the one who introduced it.
5.1. Basic Concepts in the Theory of Fuzzy Sets
The data type utilised to estimate the parameters, as well as the trustworthiness of the probability distribution, is critical to the correctness of the findings we acquire. As a result, the data type has to be given. Fuzzy data are one of these sorts of data, and it is one of the new and essential developments in statistics since many occurrences in our actual world have no precise limits. They, too, lack precision in their measurements. As a result, we will go over some fundamental notions in the theory of fuzzy sets [31].
Let be a continuous random variable that represents a system’s failure time (component). Chen and Pham [33] may then calculate the fuzzy dependability using the fuzzy probability formula:where is a membership function that indicates the degree to which each element of a given universe belongs to the fuzzy set.
Now, assume that is
The lifespan for can be calculated using the computational approach of the function of fuzzy numbers and corresponds to a given value of as follows (Chen and Pham [33]):
then
Thus, the fuzzy reliability values can be obtained for all values of , and we consider fuzzy reliability of MKKu distribution depends on the definition of fuzzy reliability as per the following equations:
If ,
If has value,where , and if ,
6. Simulation Results
In this part, we estimated the MKKu parameters using two different approaches and compared their performance using a simulation study. We explore several sample sizes for various parameter choices of , and . We produce random samples of from the MKKu distribution. We derive the relative bias (RB) values of the estimates, the mean square error (MSE), and the accompanying fuzzy reliability for each estimate as follows:
The effectiveness of different estimators is quantified in terms of RB, MSE, and confidence interval, with the most effective methods having MSE values near to zero. The R software is used to acquire the simulation results. Tables 1 and 2 show the estimates RB, MSE, and confidence interval for MLE and Bayesian.
6.1. Concluding Remarks on the Simulation’s Findings
(1)All estimation methods work admirably and have a very tiny MSE, and their mean values trend to the original values of the parameters.(2)The discrepancies in the results of all estimators are extremely minor, both in terms of the MSE and the average of the parameters.(3)The fuzzy reliability with high value of cut has the largest value with the largest MSE.(4)In median sample size, we note the MPS is better than MLE.(5)We noted that the CVM estimation method is better than MLE, MPS, and AD estimation method in Table 2.(6)We noted that AD estimation method is better than MLE, MPS, and AD estimation method in Table 1.(7)We noted that the Bayesian estimation method is better than another estimation methods.7. Analyses of Data Used as an Application to the Proposed Distribution
This section will demonstrate the MKKu distribution’s potential for real-world data set. The MKKu distribution is compared to many different models, including Kumaraswamy, Beta, and Gompertz Lomax (GL) (Oguntunde et al. [15]), Topp-Leone generalized exponential (TLGE) [14], Type II Power Topp-Leone invers exponential (TIIPTLIE) [13], Kumaraswamy-Kumaraswamy (KuKu) [34], alpha power inverted TL (APITL) [16], Toop-Leone Weibull Lomax (TLWL) [35], and Marshall-Olkin Kumaraswamy (MOKu) [36] distributions.
Table 3 contains the values of the AIC, BIC, CAIC, HQIC, Cramer-von Mises statistics (WS), Anderson–Darling statistics (AS), and Kolmogorov-Smirnov statistics (KSS), as well as their associated P-values, for all models fitted using a real data set. Additionally, these tables provide the mean and standard error (SE) of the parameters for the investigated models. Figure 3 illustrates the fit empirical, histogram, and PP-plot for the MKKu distribution for Saudi Arabia’s COVID-19 data.

The data contains 40 days of COVID-19 data from Saudi Arabia, from 22 July to 30 August 2021 (see the link [37]). These results illustrate the drought mortality rate. The data are as follows: 0.2375 0.2962 0.2167 0.2752 0.2353 0.2347 0.1951 0.2140 0.2329 0.2711 0.2126 0.2314 0.1924 0.2113 0.2683 0.2487 0.2674 0.1716 0.2666 0.2091 0.2278 0.1706 0.2271 0.1890 0.2077 0.2452 0.1319 0.2259 0.1504 0.1879 0.1689 0.2063 0.2249 0.1686 0.1310 0.1497 0.1309 0.1495 0.1121 0.1120. The results of application of Saudi Arabia are shown in Table 3 and Figure 3.
Table 3 shows MLE and different measures for parameter estimates of MKKu distribution for COVID-19 data of Saudi Arabia. Table 4 shows different estimators by using different estimation methods. Figure 4 shows convergence plots of MCMC for parameter estimates of MKKu distribution. Figure 5 shows density of posterior marginal with histogram of MCMC results for parameter estimates of MKKu distribution.


7.1. Concluding Remarks on the Real Data Application
When compared to all other models employed to fit the COVID-19 data, the MKKu distribution has the greatest P-value and the smallest distance of Kolmogorov-Smirnov (KS), W, and A values in Table 3. So we can say that the proposed distribution is superior among all its competitors as it outperforms all its rivals in fitting COVID-19 data.
8. Additional Points
We propose a novel MKKu distribution generalization, in this paper. We have investigated its statistical properties and determined a linear representation for its PDF that effectively determined the linear representations of the PDF and CDF, moments, moment generating function, and fuzzy reliability function. Numerous classical and Bayesian estimation methods were considered to obtain point estimates for the unknown MKKu parameters , , . The simulation study was done by using the R package, to see and compare the performance of various estimation methods. The MCMC method was used to accomplish this, and we concluded that the Bayesian approach outperforms all the classical approaches considered. We used a COVID-19 data set from the KSA. It was determined that MKKu fitted the data more effectively than most other competing distributions.
9. Future Work
In future, we will work on different kinds of fuzzy reliability function under different kinds of censored data, not a complete sample to study the performance of the estimators under censoring scheme; also we will update the work to find the optimal censoring scheme. We will also go forward to apply accelerated life test on the censored data under different censoring schemes.
Data Availability
All data are available in the paper.
Conflicts of Interest
The authors declare that they have no conflicts of interest.