Previous Article in Journal
Fixed Points of Some Asymptotically Regular Multivalued Mappings Satisfying a Kannan-Type Condition

Article

# A New Inverted Topp-Leone Distribution: Applications to the COVID-19 Mortality Rate in Two Different Countries

by 2 and
1
Faculty of Business Administration, Delta University of Science and Technology, El-Mansoura 11152, Egypt
2
Department of Statistics, University of Tabuk, Tabuk 71491, Saudi Arabia
3
Department of Mathematics, Faculty of Science, Helwan University, Helwan 11795, Egypt
*
Authors to whom correspondence should be addressed.
Received: 30 December 2020 / Revised: 9 February 2021 / Accepted: 17 February 2021 / Published: 26 February 2021

## Abstract

This paper aims to find a statistical model for the COVID-19 spread in the United Kingdom and Canada. We used an efficient and superior model for fitting the COVID 19 mortality rates in these countries by specifying an optimal statistical model. A new lifetime distribution with two-parameter is introduced by a combination of inverted Topp-Leone distribution and modified Kies family to produce the modified Kies inverted Topp-Leone (MKITL) distribution, which covers a lot of application that both the traditional inverted Topp-Leone and the modified Kies provide poor fitting for them. This new distribution has many valuable properties as simple linear representation, hazard rate function, and moment function. We made several methods of estimations as maximum likelihood estimation, least squares estimators, weighted least-squares estimators, maximum product spacing, Cram$e´$r-von Mises estimators, and Anderson-Darling estimators methods are applied to estimate the unknown parameters of MKITL distribution. A numerical result of the Monte Carlo simulation is obtained to assess the use of estimation methods. also, we applied different data sets to the new distribution to assess its performance in modeling data.

## 1. Introduction

Modeling for lifetime distributions have attracted great attention over the years and decades, and their interest has grown over time because the importance of modeling phenomena and pandemics is very important. Researchers in distribution theory prefer to make modeling for data either by introducing a new parameter to make the distribution of interest more flexible or perhaps by producing a new distribution family. Modeling is very interesting in numerous fields, such as industry, engineering, reliability, and medical research, see (Anake et al. [1]) for more reading. The inverted distributions have a great importance due to their applicability in many areas such as biological sciences, life test problems, medical, etc. The density and hazard ratio forms of the inverted distributions show a different structure from the other non-inverted distributions. Many authors used the inverted distributions in many applications, the reader can refer to Abd AL-Fattah et al. [2], Barco et al. [3], Hassan and Abd-Allah [4], Muhammed [5], Chesneau et al. [6], Usman and ul Haq [7], Eferhonore et al. [8] and others. Also, in this paper, we made a statistical modeling for COVID-19 mortality data in two different countries, for more reading about this, see Kumar [9], Khakharia et al. [10], Wang [11], Lalmuanawma et al. [12] and Bullock et al. [13].
The cumulative distribution function (CDF) and probability density function (PDF) of the inverted Topp-Leone distribution (ITL) distribution with shape parameter $δ > 0$, have been proposed by Hassan et al. [14] as follows:
$G ( z ; δ ) = 1 − ( 1 + 2 z ) δ ( 1 + z ) 2 δ ; z ≥ 0 , δ > 0$
and,
$g ( z ; δ ) = 2 δ z ( 1 + z ) − 2 δ − 1 ( 1 + 2 z ) δ − 1 ; z , δ > 0$
Kumar and Dharmaja [15] studied the properties of the exponentiated Kies distribution. The product moments of a modified Kies distribution under type II progressive censored sample, as well as an estimate of the distribution parameters, have been derived, see Dey et al. [16]. A new family of distributions based on the modified Kies (MK) distribution family was presented by Al-Babtain et al. [17]. If $G ( z ; δ )$ is the baseline CDF depending on a parameter vector $δ$, then the CDF of the MK family is defined by
$F ( z ; Θ ) = 1 − e − G ( z ; δ ) 1 − G ( z ; δ ) α , z > 0 , α > 0 ,$
where $Θ$ is parameter vector $( α , δ )$. The PDF of (3) is as below
$f ( z ; Θ ) = α g ( z ; δ ) G ( z ; δ ) α − 1 1 − G ( z ; δ ) α + 1 e − G ( z ; δ ) 1 − G ( z ; δ ) α , z > 0 , α > 0 ,$
Referring to the distributions above the two-parameter modified Kies inverted Topp-Leone (MKITL), which has many attractive characteristics, is obtained in this article. The proposed MKITL distribution has a very versatile PDF because it can form either positive skewed, negative skewed, and also symmetric, enabling further flexibility of tail. It can model diminishing, rising, bathtub, and reverse-J hazard rates. One more advantage of the proposed distribution is that it also has an exact closed-form CDF and very convenient to handle. These merits make the distribution an elegant candidate for use in varied areas such as biomedical life testing, reliability, and actuarial data. This article introduced two real data applications and deduced from modeling the results that the new distribution is perfect competitor with some traditional and common distributions with scale and shape parameters such as ITL, Marshall-Olkin exponential, modified Kies exponential, reverse Weibull, reverse exponential, and reverse Rayleigh distributions.
In future studies, we intend to make a new extension bivariate modified Kies inverted Topp-Leone based on copula such as and Muhammed and Almetwally [18], Almetwally and Muhammed [19], El-Morshedy et al. [20], and Kim et al. [21]. We intend to discuss a new application based on a censored sample for MKITL distribution see for more reading, Almetwally et al. [22] and Aslam et al. [23].
The rest of this paper is ordered as follows: In Section 2, we obtain the MKITL distribution. In Section 3, we discuss some of the mathematical properties of the MKITL distribution. In Section 4, we obtain an estimation methods MKITL distribution. In Section 5, we obtain simulation results of MKITL distribution. In Section 6, two applications of real data analysis were obtained. The paper is summarized and concluded in Section 7.

## 2. MKITL Distribution

By using the ITL distribution with the parameter of $δ$ and the CDF given (for $z > 0$) by using Equations (1) and (2). We define the CDF for the distribution MKITL, by substituting in Equation (3) the CDF of the ITL distribution.
$F ( z ; Θ ) = 1 − e − ( 1 + z ) 2 δ ( 1 + 2 z ) − δ − 1 α , z > 0 , α , δ > 0 ,$
and the corresponding PDF is
$f ( z ; Θ ) = 2 α δ z ( 1 + z ) 2 δ − ( 1 + 2 z ) δ α − 1 ( 1 + 2 z ) α δ + 1 ( 1 + z ) 1 − 2 δ e − ( 1 + z ) 2 δ ( 1 + 2 z ) − δ − 1 α , z > 0 , α , δ > 0 .$
where $Θ$ is parameter vector $( α , δ )$. The hazard rate (HR) function of the MKITL distribution is shown as
$h ( z ; Θ ) = 2 α δ z ( 1 + z ) 2 δ − ( 1 + 2 z ) δ α − 1 ( 1 + 2 z ) α δ + 1 ( 1 + z ) 1 − 2 δ .$
Figure 1 and Figure 2 demonstrate certain plots of the MKITL distribution for the values specified for $α$ and $δ$. The diagrams shown in Figure 2 demonstrate that the MKITL distribution HR function can be increased, decreased, and shaped a bathtub. One of the merits of MKITL distribution over an ITL distribution is that the latter is very poor model for data and phenomena that shows increasing, declining shapes, failure rates of the bathtub so it becomes more versatile to evaluate lifetime data than its competitor.

## 3. All Statistical Properties of the Proposed Distribution of MKITL

#### 3.1. Linear Representation for the MKITL Distribution

We introduce a linear representation in this subsection for the MK family and use it to provide a useful linear representation for the MKITL distribution. A mixture representation of the MK family can be provided as follows,
$f ( z ; Θ ) = α ∑ j = 0 ∞ ( − 1 ) j j ! g ( z ; δ ) G ( z ; δ ) α ( j + 1 ) − 1 1 − G ( z ; δ ) α ( j + 1 ) + 1 .$
The last equation of the MKITL distribution can be rewritten with the ITL distribution PDF and CDF.
$f ( z ; Θ ) = ζ j , k , h 2 δ ( h + 1 ) z ( 1 + z ) − 2 δ ( h + 1 ) − 1 ( 1 + 2 z ) δ ( h + 1 ) − 1 ,$
where $ζ j , k , h = ∑ j , k , h = 0 ∞ α ( − 1 ) j + k + h ( h + 1 ) j ! α ( j + 1 ) + k k α ( j + 1 ) + k − 1 h$. Equation (9) denotes the ITL density with parameter $δ ( h + 1 )$.

#### 3.2. Quantile for the MKITL Distribution

The quantile function of the MKITL distribution, say $z = F ( z , Θ ) − 1 ( Q )$ is derived by inverting (5) as follows:
$z = 1 − 1 + − l n ( 1 − Q ) 1 α − 1 δ − 0.5 − 1 − 1 ; 0 < Q < 1$
In particular, the first quartile, say Q1, the second quartile, say Q2, and the third quartile, say Q3 are obtained by setting $Q = 0.25 , 0.5 , 0.75$, respectively, in (10).

#### 3.3. Moments for the MKITL Distribution

The $r t h$ moment of z comes from Equation (9) as
$μ r ´ = E ( z r ) = ζ j , k , h 2 δ ( h + 1 ) β ( j + r + 2 , δ ( h + 1 ) − r )$
The $r t h$ incomplete moment of z can be obtained from (9) as
$ω r ( t ) = ∫ 0 t z r f ( z , Θ ) d z = ζ j , k , h 2 δ ( h + 1 ) ∑ b = 0 ∞ δ ( h + 1 ) − 1 b B ( b + r + 2 , δ ( h + 1 ) − r , t 1 + t )$
where $B ( b + r + 2 , δ ( h + 1 ) − r , t 1 + t )$ is the incomplete beta function.

## 4. Methods of Estimation for the Distribution Parameters

We managed to make a solution for the estimation problem of the MKITL distribution parameters using six different estimation methods. These methods are the maximum likelihood estimators (MLE), least squares estimators (LSEs), maximum product of spacing estimators (MPSE), weighted least-squares estimators (WLSE), Cram$e ´$r-von Mises estimators (CVME) and Anderson-Darling estimators (ADE).

#### 4.1. Maximum Likelihood Estimators

Let $z 1 , . . . , z n$ be a random sample from the MKITL distribution with parameters $α$ and $δ$. the log-likelihood function for MKITL distribution is given by
$l ( Θ ) = n log ( 2 ) + log ( α ) + l o g ( δ ) + ∑ i = 1 n l o g ( z i ) − ( 1 − 2 δ ) ∑ i = 1 n log ( 1 + z i ) + ( α − 1 ) ∑ i = 1 n log ( 1 + z i ) 2 δ − ( 1 + 2 z i ) δ − ( α δ + 1 ) ∑ i = 1 n log ( 1 + 2 z i ) − ∑ i = 1 n ( 1 + z i ) 2 δ ( 1 + 2 z i ) − δ − 1 α$
The partial derivatives of $l ( Θ )$ with respect to the model parameters $α$ and $δ$ are
$∂ l ( Θ ) ∂ α = n α + ∑ i = 1 n log ( 1 + z i ) 2 δ − ( 1 + 2 z i ) δ − α ∑ i = 1 n log ( 1 + 2 z i ) − ∑ i = 1 n ( 1 + z i ) 2 δ ( 1 + 2 z i ) − δ − 1 α log ( 1 + z i ) 2 δ ( 1 + 2 z i ) − δ − 1$
and
$∂ l ( Θ ) ∂ δ = n δ + 2 ∑ i = 1 n log ( 1 + z i ) − α ∑ i = 1 n log ( 1 + 2 z i ) + ( α − 1 ) ∑ i = 1 n 2 ( 1 + z i ) 2 δ log ( 1 + z i ) − ( 1 + 2 z i ) δ log ( 1 + 2 z i ) ( 1 + z i ) 2 δ − ( 1 + 2 z i ) δ − α ∑ i = 1 n ( 1 + z i ) 2 δ ( 1 + 2 z i ) δ − 1 α − 1 ( 1 + z i ) 2 1 + 2 z i δ log ( 1 + z i ) 2 1 + 2 z i$
The MLE of the distribution parameters $α$ and $δ$ can be calculated by finding the maximum value for Equation (13) by finding the first derivative for it with respect to $α$ and $δ$. The R software can be used to optimise the log-likelihood and obtain the MLE by using the Newton-Rapshon method.

#### 4.2. Least Squares and Weighted Least Squares Methods

Here we introduce two methods of estimation that we can use in this paper; these methods are the least squares (LS) and weighted least square (WLS) methods. They are very important common methods for estimating the parameters of different distributions. Suppose that we have a random sample say $z ( 1 ) < . . . < z ( n )$ from the MKITL distribution with the two parameters $α$ and $δ$. Now we can calculate both the LS estimators (LSE) and WLS estimators (WLSE) of the MKITL distribution parameters $α$ and $δ$ by making the following equation as minimum as possible
$V ( Θ ) = ∑ i = 1 n υ i 1 − e − ( 1 + z i ) 2 δ ( 1 + 2 z i ) − δ − 1 α − i n + 1 2$
with respect to $α$ and $δ$, where $υ i = 1$ in case of LSE and $υ i = ( n + 1 ) 2 ( n + 2 ) [ i ( n − i + 1 ) ]$ in case of WLSE. Furthermore, the LSE and WLSE follow by solving the nonlinear equations
$∂ V ( Θ ) ∂ α = 2 ∑ i = 1 n υ i 1 − e − ( 1 + z i ) 2 δ ( 1 + 2 z i ) − δ − 1 α − i n + 1 e − ( 1 + z i ) 2 δ ( 1 + 2 z i ) − δ − 1 α ( 1 + z i ) 2 δ ( 1 + 2 z i ) − δ − 1 α log ( 1 + z i ) 2 δ ( 1 + 2 z i ) − δ − 1 ,$
and
$∂ V ( Θ ) ∂ δ = 2 α ∑ i = 1 n υ i 1 − e − ( 1 + z i ) 2 δ ( 1 + 2 z i ) − δ − 1 α − i n + 1 e − ( 1 + z i ) 2 δ ( 1 + 2 z i ) − δ − 1 α ( 1 + z i ) 2 δ ( 1 + 2 z i ) δ − 1 α − 1 ( 1 + z i ) 2 1 + 2 z i δ log ( 1 + z i ) 2 1 + 2 z i .$

#### 4.3. Maxzimum Product of Spacings Method

In this sub section, we will use an alternative method for the MLE method which is called the maximum product of spacings (MPS), this method is used to estimate the parameters of continuous univariate models, this method was introduced by Cheng and Amin [24] and used for censored application by Singh et al. [25], Almetwally and Almongy [26], Basu et al. [27], Almetwally et al. [28], El-Sherpieny et al. [29] and Alshenawy et al. [30] and [31]. Let $z ( 1 ) < . . . < z ( n )$, be a random sample of size n, the uniform spacing of the MKITL distribution can be defined by
$D i ( Θ ) = F ( z ( i ) , Θ ) − F ( z ( i − 1 ) , Θ ) ; i = 1 , . . . , n + 1$
where $D i ( Θ )$ denotes to the uniform spacings, $F ( z ( 0 ) , Θ ) = 0$, $F ( z ( n + 1 ) , Θ ) = 1$ and $∑ i = 1 n + 1 D i ( Θ ) = 1$.
$G ( Θ ) = 1 n + 1 ∑ i = 1 n + 1 log e − 1 + z ( i − 1 ) 2 1 + 2 z ( i − 1 ) δ − 1 α − e − 1 + z ( i ) 2 1 + 2 z ( i ) δ − 1 α$
The MPS estimators (MPSE) of the MKITL parameters can be obtained by maximizing Equation (20) with respect to $α$ and $δ$.
$∂ G ( Θ ) ∂ α = 1 n + 1 ∑ i = 1 n + 1 e − 1 + z ( i ) 2 1 + 2 z ( i ) δ − 1 α 1 + z ( i ) 2 1 + 2 z ( i ) δ − 1 α log 1 + z ( i ) 2 1 + 2 z ( i ) δ − 1 e − 1 + z ( i − 1 ) 2 1 + 2 z ( i − 1 ) δ − 1 α − e − 1 + z ( i ) 2 1 + 2 z ( i ) δ − 1 α − e − 1 + z ( i − 1 ) 2 1 + 2 z ( i − 1 ) δ − 1 α 1 + z ( i − 1 ) 2 1 + 2 z ( i − 1 ) δ − 1 α log 1 + z ( i − 1 ) 2 1 + 2 z ( i − 1 ) δ − 1 e − 1 + z ( i − 1 ) 2 1 + 2 z ( i − 1 ) δ − 1 α − e − 1 + z ( i ) 2 1 + 2 z ( i ) δ − 1 α ,$
$∂ G ( Θ ) ∂ δ = α n + 1 ∑ i = 1 n + 1 e − 1 + z ( i ) 2 1 + 2 z ( i ) δ − 1 α ( 1 + z ( i ) ) 2 δ ( 1 + 2 z ( i ) ) δ − 1 α − 1 ( 1 + z ( i ) ) 2 1 + 2 z ( i ) δ log ( 1 + z ( i ) ) 2 1 + 2 z ( i ) e − 1 + z ( i − 1 ) 2 1 + 2 z ( i − 1 ) δ − 1 α − e − 1 + z ( i ) 2 1 + 2 z ( i ) δ − 1 α − e − 1 + z ( i − 1 ) 2 1 + 2 z ( i − 1 ) δ − 1 α ( 1 + z ( i − 1 ) ) 2 δ ( 1 + 2 z ( i − 1 ) ) δ − 1 α − 1 ( 1 + z ( i − 1 ) ) 2 1 + 2 z ( i − 1 ) δ log ( 1 + z ( i − 1 ) ) 2 1 + 2 z ( i − 1 ) e − 1 + z ( i − 1 ) 2 1 + 2 z ( i − 1 ) δ − 1 α − e − 1 + z ( i ) 2 1 + 2 z ( i ) δ − 1 α .$
Furthermore, the MPSE of the MKITL parameters can also be obtained by solving Equations (21) and (22), using the R program.

#### 4.4. Cram$e ´$r-von-Mises Method

The Cram$e ´$r-von-Mises estimators (CVME), can be calculated depending on the difference between both the CDF and the empirical distribution function estimates. For more information see (Luceno [32]).
$C ( Θ ) = 1 12 + ∑ i = 1 n 1 − e − 1 + z ( i ) 2 1 + 2 z ( i ) δ − 1 α − 2 i − 1 2 n 2 .$
We can find the estimates by The CVME method for the two parameters of MKITL, which are $α$ and $δ$ by finding the minimum value for Equation (23) with respect to $α$ and $δ$.
$∂ C ( Θ ) ∂ α = 2 ∑ i = 1 n 1 − e − 1 + z ( i ) 2 1 + 2 z ( i ) δ − 1 α − 2 i − 1 2 n e − 1 + z ( i ) 2 1 + 2 z ( i ) δ − 1 α log 1 + z ( i ) 2 1 + 2 z ( i ) δ − 1 ,$
and
$∂ C ( Θ ) ∂ δ = α 2 ∑ i = 1 n 1 − e − 1 + z ( i ) 2 1 + 2 z ( i ) δ − 1 α − 2 i − 1 2 n e − 1 + z ( i ) 2 1 + 2 z ( i ) δ − 1 α 1 + z ( i ) 2 1 + 2 z ( i ) δ − 1 α − 1 1 + z ( i ) 2 1 + 2 z ( i ) δ log 1 + z ( i ) 2 1 + 2 z ( i ) .$
Further, the CVME found by solving the nonlinear Equations (24) and (25).

#### 4.5. Anderson-Darling Method

Here we introduce another minimum distance estimation method. This method is called the Anderson-Darling estimation (ADE).
$A ( Θ ) = − n − 1 n ∑ i = 1 n ( 2 i − 1 ) log 1 − e − 1 + z ( i ) 2 1 + 2 z ( i ) δ − 1 α − 1 + z ( n + 1 − i ) 2 1 + 2 z ( n + 1 − i ) δ − 1 α 2$
We can manipulate the ADE to estimate the parameters of the MKITL distribution by making the Equation (26) at minimum value with respect to $α$ and $δ$.

## 5. Simulation Results

In this section, we made an estimation by using six methods of the MKITL parameters, and compared between them to assess the performance of each method by a simulation study. We consider different sample sizes $n = 25 , 70 , 150$ for different parameters values $α = ( 0.6 , 1.7 , 3 )$ and $δ = ( 0.6 , 1.7 , 3 )$. We generate $N = 10 , 000$ random samples from MKITL distribution. For each estimate, we obtain the average values of the estimates (mean) and their corresponding mean squares error (MSE).
The performance of different estimators is evaluated in terms of MSE, i.e., the most efficient method of estimating will be those whose MSE values are closer to zero. Simulation results are obtained via the R program. Table 1, Table 2 and Table 3 show the estimates mean and MSE for MLE, LS, WLS, MPS, CVM, and AD. Moreover, the mean using all estimation methods tends towards the true parameter values, as the sample size increases in most cases.

#### Concluded Observation on the Simulation Results

• When $α = 0.6 , 1.7 , 3$, and as the value of $δ$ increases the MSE of the two parameters increases in most cases, while the mean values of both parameters tends to the initial values.
• For a fixed value of $α$, and $δ$, by increasing the sample size, the mean values of the parameters tend to the initial values, and MSE decreases.
• In most cases, the MPS was the best method for estimating the parameter referring to the mean values of the parameters and MSE as in Table 1, Table 2 and Table 3, while the AD was the best efficient method for estimating the parameters referring to the mean values of the parameters and MSE in Table 3 when $δ = 0.6$.
• All estimators perform very well and provides very small MSE and the mean value of estimates for all estimators tend to the initial value of the parameters.
• The differences between all estimators values are very small, referring to the MSE values and the mean value of the parameters.

## 6. Application of Real Data Analysis

This section is devoted to illustrating the potentiality of the MKITL distribution for two real data sets. MKITL distribution is compared with other competitive models, namely: ITL (Hassan et al. [14]), the Marshall-Olkin exponential (MOEx) [Marshall and Olkin [33], modified Kies exponential (MKEx) (Al-Babtain et al. [17]), inverse Weibull (IW), inverse exponential (IE) and inverse Rayleigh (IR) (Ahmad et al. [34]).
Table 4 and Table 5 provide values of Cram$e ´$r-von Mises (W*), Anderson-Darling (A*) and Kolmogorov- Smirnov (KS) statistic along with its P-value for the all models fitted based on Two real data sets. In addition, these tables contain the MLE and standard errors (SE) of the parameters for the considered models. In Table 4 and Table 5 the MKITL distribution has the highest P-value and the lowest distance of Kolmogorov-Smirnov(KS), W*, and A* value when it compares with all other models used here to fit the COVID-19. Figure 3 and Figure 4 show the fit empirical, histogram, QQ-plot, and PP-plot for the MKITL distribution for COVID-19 data of the United Kingdom and Canada.

#### 6.1. Application 1

Firstly: The data represents a COVID-19 data belong to Canada of 36 days, from 10 April to 15 May 2020 see the link [https://covid19.who.int/].
These data formed of drought mortality rate. The data are as follows: 3.1091 3.3825 3.1444 3.2135 2.4946 3.5146 4.9274 3.3769 6.8686 3.0914 4.9378 3.1091 3.2823 3.8594 4.0480 4.1685 3.6426 3.2110 2.8636 3.2218 2.9078 3.6346 2.7957 4.2781 4.2202 1.5157 2.6029 3.3592 2.8349 3.1348 2.5261 1.5806 2.7704 2.1901 2.4141 1.9048.

#### 6.2. Application 2

Secondly: The data represents a COVID-19 data belong to The United Kingdom of 24 days, from 15 October to 7 November 2020 [https://covid19.who.int/]. These data formed of drought mortality rate. The data are as follows: 0.2240 0.2189 0.2105 0.2266 0.0987 0.1147 0.3353 0.2563 0.2466 0.2847 0.2150 0.1821 0.1200 0.4206 0.3456 0.3045 0.2903 0.3377 0.1639 0.1350 0.3866 0.4678 0.3515 0.3232.

#### 6.3. Concluding Remarks on the Two Application

• Referring to data set one we can see that MKITL provides the highest P-value, and the lowest W*, A* and lowest KS distance.
• From Figure 3 , we can deduce that MKITL was the best model for fitting the first data set.
• Referring to data set two, we can see that MKITL provides the highest P-value, and the lowest A*, W* and KS distance.
• From Figure 4, we can deduce that MKITL was the best model for fitting the second data set.
• Referring to Table 4 we can see that IE, IR, ITL and MOEx distribution provides poor fitting for the first data set.
• Referring to Table 5 we can see that IE, MOEx distribution provides poor fitting for the second data set.
• We can conclude that from both applications that the MKITL distribution provides the best fitting among all its competitive distributions, which gives it superiority in fitting this kind of mortality rate for COVID-19 data.

## 7. Summary

In this paper, we propose a new two-parameter model, called the modified Kies inverted Topp-Leone distribution, which can be denoted as MKITL distribution. The MKITL distribution is more flexible to analyze lifetime data than other common distributions. Survival function, hazard function, linear representation, quantile, and moments of the MKITL distribution are provided. We compare MLE, LSE, MPSE, WLSE, CVME, and ADE methods, and we conclude the alternative methods of MLE are better than the MLE method. We provide two applications on COVID 19 mortality rate and we proved that the MKITL distribution is the best model among all its competitors to fit this kind of data. The parameter estimation of MKITL distribution is derived by MLE, LSE, MPSE, WLSE, CVME, and ADE. The methods of estimation are used to estimate the model parameters, and simulation results are provided to assess the model performance. Two real-life data proposed model provides a consistently better fit than the MKITL, IE, IW, MOEx, MKEx, and IR distributions.

## Author Contributions

All authors contribute equally to this paper. All authors have read and agreed to the published version of the manuscript.

## Funding

This research received no external funding.

## Conflicts of Interest

The authors declare no conflict of interest.

## References

1. Anake, T.A.; Oguntunde, P.E.; Odetunmibi, O.A. On a fractional beta- distribution. Int. J. Math. Comput. 2015, 26, 26–34. [Google Scholar]
2. Abd AL-Fattah, A.M.; El-Helbawy, A.A.; Al-Dayian, G.R. Inverted Kumaraswamy Distribution: Properties and Estimation. Pak. J. Stat. 2017, 33, 37–61. [Google Scholar]
3. Barco, K.V.P.; Mazucheli, J.; Janeiro, V. The inverse power Lindley distribution. Commun. Stat.-Simul. Comput. 2017, 46, 6308–6323. [Google Scholar] [CrossRef]
4. Hassan, A.S.; Abd-Allah, M. On the inverse power Lomax distribution. Ann. Data Sci. 2019, 6, 259–278. [Google Scholar] [CrossRef]
5. Muhammed, H.Z. On the inverted Topp Leone distribution. Int. J. Reliab. Appl. 2019, 20, 17–28. [Google Scholar]
6. Chesneau, C.; Tomy, L.; Gillariose, J.; Jamal, F. The inverted modified Lindley distribution. J. Stat. Theory Pract. 2020, 14, 1–17. [Google Scholar] [CrossRef]
7. Usman, R.M.; ul Haq, M.A. The Marshall-Olkin extended inverted Kumaraswamy distribution: Theory and applications. J. King Saud Univ.-Sci. 2020, 32, 356–365. [Google Scholar] [CrossRef]
8. Eferhonore, E.E.; Thomas, J.; Zelibe, S.C. Theoretical Analysis of the Weibull Alpha Power Inverted Exponential Distribution: Properties and Applications. Gazi Univ. J. Sci. 2020, 33, 265–277. [Google Scholar]
9. Kumar, S. Monitoring Novel Corona Virus (COVID-19) Infections in India by Cluster Analysis. Ann. Data Sci. 2020, 7, 417–425. [Google Scholar] [CrossRef]
10. Khakharia, A.; Shah, V.; Jain, S.; Shah, J.; Tiwari, A.; Daphal, P.; Mehendale, N. Outbreak Prediction of COVID-19 for Dense and Populated Countries Using Machine Learning. Ann. Data Sci. 2020, 8, 1–19. [Google Scholar] [CrossRef]
11. Wang, Y.z.J. A call for caution in extrapolating chest CT sensitivity for COVID-19 derived from hospital data to patients among general population. Quant. Imaging Med. Surg. 2020, 10, 798. [Google Scholar] [CrossRef]
12. Lalmuanawma, S.; Hussain, J.; Chhakchhuak, L. Applications of machine learning and artificial intelligence for Covid-19 (SARS-CoV-2) pandemic: A review. Chaos Solitons Fractals 2020, 139, 110059. [Google Scholar] [CrossRef] [PubMed]
13. Bullock, J.; Pham, K.H.; Lam, C.S.N.; Luengo-Oroz, M. Mapping the landscape of artificial intelligence applications against COVID-19. arXiv 2020, arXiv:2003.11336. [Google Scholar]
14. Hassan, A.S.; Elgarhy, M.; Ragab, R. Statistical Properties and Estimation of Inverted Topp-Leone Distribution. J. Stat. Appl. Probab. 2020. forthcoming. [Google Scholar]
15. Kumar, C.S.; Dharmaja, S.H.S. The exponentiated reduced Kies distribution: Properties and applications. Commun. Stat.-Theory Methods 2017, 46, 8778–8790. [Google Scholar] [CrossRef]
16. Dey, S.; Nassar, M.; Kumar, D. Moments and estimation of reduced Kies distribution based on progressive type-II right censored order statistics. Hacet. J. Math. Stat. 2019, 48, 332–350. [Google Scholar] [CrossRef]
17. Al-Babtain, A.A.; Shakhatreh, M.K.; Nassar, M.; Afify, A.Z. A New Modified Kies Family: Properties, Estimation Under Complete and type-II Censored Samples, and Engineering Applications. Mathematics 2020, 8, 1345. [Google Scholar] [CrossRef]
18. Muhammed, H.Z.; Almetwally, E.M. Bayesian and Non-Bayesian Estimation for the Bivariate Inverse Weibull Distribution Under Progressive Type-II Censoring. Ann. Data Sci. 2020. [Google Scholar] [CrossRef]
19. Almetwally, E.M.; Muhammed, H.Z. On a bivariate Fréchet distribution. J. Stat Appl Probab. 2020, 9, 1–21. [Google Scholar]
20. El-Morshedy, M.; Alhussain, Z.A.; Atta, D.; Almetwally, E.M.; Eliwa, M.S. Bivariate Burr X generator of distributions: Properties and estimation methods with applications to complete and type-II censored samples. Mathematics. 2020, 8, 264. [Google Scholar] [CrossRef]
21. Kim, J.M.; Ju, H.; Jung, Y. Copula Approach for Developing a Biomarker Panel for Prediction of Dengue Hemorrhagic Fever. Ann. Data Sci. 2020, 7, 697–712. [Google Scholar] [CrossRef]
22. Almetwally, E.M.; Almongy, H.M.; Rastogi, M.K.; Ibrahim, M. Maximum Product Spacing Estimation of Weibull Distribution Under Adaptive Type-II Progressive Censoring Schemes. Ann. Data Sci. 2020, 7, 257–279. [Google Scholar] [CrossRef]
23. Aslam, M.; Yousaf, R.; Ali, S. Bayesian Estimation of Transmuted Pareto Distribution for Complete and Censored Data. Ann. Data Sci. 2020, 7, 663–695. [Google Scholar] [CrossRef]
24. Cheng, R.C.H.; Amin, N.A.K. Estimating parameters in continuous univariate distributions with a shifted origin. J. R. Stat. Soc. Ser. B (Methodol.) 1983, 45, 394–403. [Google Scholar] [CrossRef]
25. Singh, R.K.; Singh, S.K.; Singh, U. Maximum product spacings method for the estimation of parameters of generalized inverted exponential distribution under Progressive Type II Censoring. J. Stat. Manag. Syst. 2016, 19, 219–245. [Google Scholar]
26. Almetwally, E.M.; Almongy, H.M. Maximum product spacing and Bayesian method for parameter estimation for generalized power Weibull distribution under censoring scheme. J. Data Sci. 2019, 17, 407–444. [Google Scholar]
27. Basu, S.; Singh, S.K.; Singh, U. Estimation of Inverse Lindley Distribution Using Product of Spacings Function for Hybrid Censored Data. Methodol. Comput. Appl. Probab. 2019, 21, 1377–1394. [Google Scholar] [CrossRef]
28. Almetwally, E.M.; Almongy, H.M.; ElSherpieny, E.A. Adaptive type-II progressive censoring schemes based on maximum product spacing with application of generalized Rayleigh distribution. J. Data Sci. 2019, 17, 802–831. [Google Scholar]
29. El-Sherpieny, E.S.A.; Almetwally, E.M.; Muhammed, H.Z. Progressive Type-II hybrid censored schemes based on maximum product spacing with application to Power Lomax distribution. Phys. A Stat. Mech. Its Appl. 2020, 553, 124251. [Google Scholar] [CrossRef]
30. Alshenawy, R.; Al-Alwan, A.; Almetwally, E.M.; Afify, A.Z.; Almongy, H.M. Progressive Type-II Censoring Schemes of Extended Odd Weibull Exponential Distribution with Applications in Medicine and Engineering. Mathematics 2020, 8, 1679. [Google Scholar] [CrossRef]
31. Alshenawy, R.; Sabry, M.A.; Almetwally, E.M.; Almongy, H.M. Product Spacing of Stress–Strength under Progressive Hybrid Censored for Exponentiated-Gumbel Distribution. Comput. Mater. Contin. 2021, 66, 2973–2995. [Google Scholar] [CrossRef]
32. Luceno, A. Fitting the generalized Pareto distribution to data using maximum goodness of fit estimators. Comput. Stat. Data Anal. 2006, 51, 904–917. [Google Scholar] [CrossRef]
33. Marshall, A.W.; Olkin, I. A new method for adding a parameter to a family of distributions with application to the exponential and Weibull families. Biometrika 1997, 84, 641–652. [Google Scholar] [CrossRef]
34. Ahmad, A.; Ahmad, S.P.; Ahmed, A. Transmuted Inverse Rayleigh distribution: A generalization of the Inverse Rayleigh distribution. Math. Theory Model. 2014, 4, 90–98. [Google Scholar]
Figure 1. PDF of the MKITL distribution for certain parameter values.
Figure 1. PDF of the MKITL distribution for certain parameter values.
Figure 2. HF of the MKITL distribution for certain parameter values.
Figure 2. HF of the MKITL distribution for certain parameter values.
Figure 3. Cumulative function and empirical CDF, histogram, Q-Q plot and P-P plot for the MKITL distribution for COVID-19 data of Canada.
Figure 3. Cumulative function and empirical CDF, histogram, Q-Q plot and P-P plot for the MKITL distribution for COVID-19 data of Canada.
Figure 4. Cumulative function and empirical CDF, histogram, Q-Q plot and P-P plot for the MKITL distribution for COVID-19 data of The United Kingdom.
Figure 4. Cumulative function and empirical CDF, histogram, Q-Q plot and P-P plot for the MKITL distribution for COVID-19 data of The United Kingdom.
Table 1. Mean and MSE of the MLE, LS, MPS, WLS, CVM and AD estimate for MKITL distribution when $α = 0.6$.
Table 1. Mean and MSE of the MLE, LS, MPS, WLS, CVM and AD estimate for MKITL distribution when $α = 0.6$.
$δ$$n$ MeanMSEMeanMSEMeanMSEMeanMSEMeanMSEMeanMSE
0.625$α$0.56930.01670.60210.01800.56920.01310.60750.01640.63760.02230.61200.0152
$δ$0.61140.02080.61910.02800.61130.01890.61940.02520.63310.03100.61960.0226
70$α$0.58390.00470.60240.00570.58390.00440.60550.00490.61450.00620.60580.0048
$δ$0.59900.00650.60460.00830.59900.00620.60460.00730.60910.00860.60500.0072
150$α$0.58920.00220.60080.00290.58910.00210.60280.00240.60640.00300.60230.0024
$δ$0.59880.00280.60220.00380.59870.00280.60270.00340.60430.00390.60260.0033
1.725$α$0.57770.01750.61780.02030.57790.01250.62270.01880.65420.02600.62480.0170
$δ$1.72360.15681.75460.21151.72290.13911.75220.18831.79380.23501.75250.1714
70$α$0.58090.00510.60270.00620.58090.00480.60640.00560.61490.00680.60500.0052
$δ$1.71250.05501.72870.07811.71300.05141.72940.06781.74180.08131.72760.0646
150$α$0.58970.00220.60120.00270.58970.00220.60310.00240.60670.00280.60280.0023
$δ$1.71020.02261.72380.03041.71020.02151.72390.02651.72980.03111.72400.0264
325$α$0.57310.01800.60620.01960.57310.01350.61080.01800.64240.02450.61640.0168
$δ$3.04670.43463.10180.60143.04660.34973.11710.61543.17280.68213.11980.5300
70$α$0.58110.00480.59820.00610.58110.00450.60180.00530.61020.00660.60180.0051
$δ$3.00470.15543.03840.21083.00470.13463.03840.18503.06160.21943.03960.1788
150$α$0.59240.00220.60160.00290.59230.00210.60450.00250.60720.00300.60400.0024
$δ$2.99880.06873.02430.08922.99890.06503.02490.07833.03480.09093.02380.0771
Table 2. Mean and MSE of the MLE, LS, MPS, WLS, CVM and AD estimate for MKITL distribution when $α = 1.7$.
Table 2. Mean and MSE of the MLE, LS, MPS, WLS, CVM and AD estimate for MKITL distribution when $α = 1.7$.
$δ$$n$ MeanMSEMeanMSEMeanMSEMeanMSEMeanMSEMeanMSE
0.625$α$1.58560.08621.65910.08901.58590.07911.68310.08331.75910.11191.70540.0786
$δ$0.60370.00290.60470.00320.60360.00280.60590.00300.60940.00340.60720.0030
70$α$1.64830.03281.68980.03381.64820.03061.70880.03301.72200.03511.71000.0320
$δ$0.60040.00100.60120.00120.60030.00100.60170.00110.60280.00120.60190.0011
150$α$1.67060.01461.70230.01821.67060.01441.70980.01541.71850.01891.71010.0158
$δ$0.59980.00050.60050.00060.59980.00050.60070.00050.60120.00060.60080.0005
1.725$α$1.61160.11451.70110.12761.61170.08961.71880.11521.80550.15821.73700.1032
$δ$1.70860.02371.71050.02521.70860.02271.71390.02401.72380.02661.71730.0239
70$α$1.63710.03071.68670.03851.63730.03081.70010.03371.72210.04071.70170.0324
$δ$1.70280.00751.70420.00861.70270.00741.70580.00801.70890.00881.70670.0080
150$α$1.66710.01401.70320.01951.66700.01401.70920.01611.71980.02041.70830.0156
$δ$1.70030.00351.70190.00401.70020.00341.70260.00381.70400.00411.70280.0037
325$α$1.61100.11581.70900.14311.61130.08971.72560.12671.81460.17981.73910.1109
$δ$2.99870.06493.00100.07312.99850.06233.00650.06903.02470.07663.01190.0677
70$α$1.64070.03061.70400.04451.64070.03011.71280.03641.74000.04851.71160.0338
$δ$3.00270.02643.00740.03003.00250.02593.00990.02823.01560.03063.01040.0280
150$α$1.66270.01371.69820.01871.66280.01421.70400.01591.71470.01941.70330.0152
$δ$3.00050.01183.00280.01303.00030.01173.00400.01223.00660.01313.00420.0122
Table 3. Mean and MSE of the MLE, LS, MPS, WLS, CVM and AD estimate for MKITL distribution when $α = 3$.
Table 3. Mean and MSE of the MLE, LS, MPS, WLS, CVM and AD estimate for MKITL distribution when $α = 3$.
$δ$$n$ MeanMSEMeanMSEMeanMSEMeanMSEMeanMSEMeanMSE
0.625$α$2.83370.14082.89690.10582.83390.16062.93090.09852.96070.07142.96180.0585
$δ$0.59880.00090.59910.00110.59880.00090.59990.00100.60100.00110.60050.0010
70$α$2.91550.04962.97820.02952.91560.05062.99620.02593.01490.03442.98650.0424
$δ$0.59820.00030.59910.00040.59820.00030.59940.00030.59990.00040.59960.0003
150$α$2.94430.02752.98530.00822.94400.02383.00120.01292.99540.00832.99490.0090
$δ$0.59940.00010.60010.00020.59940.00010.60020.00020.60040.00020.60030.0002
1.725$α$2.81130.35292.91960.26942.81130.29872.97330.27913.07590.28443.01070.2637
$δ$1.70340.00801.70380.00851.70340.00771.70570.00811.71130.00881.70850.0080
70$α$2.89620.09672.99260.11392.89660.09703.00820.09873.05290.12033.01670.1010
$δ$1.70050.00271.70220.00301.70050.00271.70280.00291.70480.00311.70320.0028
150$α$2.94700.04452.99550.04572.94680.04353.02170.04523.02130.04543.01750.0429
$δ$1.69740.00131.69700.00141.69740.00131.69810.00131.69810.00141.69820.0013
325$α$2.83510.32452.99180.32902.83470.25963.02730.31483.17300.40193.06260.2887
$δ$2.99310.02292.99670.02512.99330.02242.99960.02393.00990.02583.00310.0236
70$α$2.90730.09463.00500.12012.90720.08963.02790.10753.06840.13033.02750.1001
$δ$3.00170.00793.00270.00863.00180.00783.00450.00813.00740.00883.00510.0080
150$α$2.93510.03902.99250.05302.93490.04053.00760.04473.02180.05463.00550.0435
$δ$2.99670.00402.99810.00442.99680.00402.99910.00413.00030.00442.99900.0041
Table 4. MLE, SE, KS test and P-values for COVID-19 data of Canada.
Table 4. MLE, SE, KS test and P-values for COVID-19 data of Canada.
$α$$δ$KSp-ValueW *A *
MKITLMLE3.35550.72860.14540.43230.10870.6473
SE0.40460.0275
MKExMLE2.28550.18640.16950.25230.22991.3306
SE0.27050.0103
IWMLE3.169123.40530.17370.22740.25481.5284
SE0.36568.1037
IRMLE8.2339 0.27900.00740.17701.0691
SE1.3723
IEMLE3.0078 0.42840.00000.12890.7751
SE0.5013
ITLMLE1.1523 0.43900.00000.15050.8600
SE0.1920
MOExMLE83.08361.44110.98900.00000.38942.2258
SE39.15210.1530
Table 5. MLE, SE, KS test and P-values for COVID-19 data of The United Kingdom.
Table 5. MLE, SE, KS test and P-values for COVID-19 data of The United Kingdom.
$α$$δ$KSp-ValueW*A*
MKITLMLE1.359912.79310.09360.97150.02440.1762
SE0.23191.4051
MKEXMLE2.17552.32410.09560.96580.02540.1831
SE0.35980.1636
ITLMLE21.7398 0.23710.11360.04030.2776
SE4.4376
IRMLE0.0399 0.16660.46850.15410.9476
SE0.0081
IWMLE2.30350.02290.18690.32920.16981.0330
SE0.33500.0149
IEMLE0.2200 0.37520.00150.10520.6731
SE0.0449
MOEXMLE50.893815.56590.98170.00000.11960.7700
SE36.25382.5549
 Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.