# Modelling the COVID-19 Mortality Rate with a New Versatile Modification of the Log-Logistic Distribution

Authors: Abdisalam Hassan Muse; Ahlam H. Tolba; Eman Fayad; Ola A. Abu Ali; M. Nagy; M. Yusuf
Journal: Computational Intelligence and Neuroscience
 (2021)

---

## Abstract

The goal of this paper is to develop an optimal statistical model to analyze COVID-19 data in order to model and analyze the COVID-19 mortality rates in Somalia. Combining the log-logistic distribution and the tangent function yields the flexible extension log-logistic tangent (LLT) distribution, a new two-parameter distribution. This new distribution has a number of excellent statistical and mathematical properties, including a simple failure rate function, reliability function, and cumulative distribution function. Maximum likelihood estimation (MLE) is used to estimate the unknown parameters of the proposed distribution. A numerical and visual result of the Monte Carlo simulation is obtained to evaluate the use of the MLE method. In addition, the LLT model is compared to the well-known two-parameter, three-parameter, and four-parameter competitors. Gompertz, log-logistic, kappa, exponentiated log-logistic, Marshall–Olkin log-logistic, Kumaraswamy log-logistic, and beta log-logistic are among the competing models. Different goodness-of-fit measures are used to determine whether the LLT distribution is more useful than the competing models in COVID-19 data of mortality rate analysis.


---

## Body

## 1. Introduction

Models are at the heart of almost all statistical work. A statistical model is a classification of probability distributions. The distribution family can be parametric, semiparametric, or nonparametric. Parametric models produce more efficient estimates with lower standard errors than nonparametric and semiparametric models [1], more specifically, if the distributional assumption is correct. In general, probability distributions have been widely used to model lifetime data in a variety of fields, particularly biomedical sciences and engineering. Because of the variability of the data, the statistical models chosen have a significant impact on the quality of the modelling in terms of providing the best description of the phenomenon under consideration.The generated family of distributions has a large influence on the quality of statistical analysis procedures, and much effort has gone into developing new statistical models. There are, however, a number of significant issues with actual data that do not fit into any of the commonly used mathematical models. As a result, the technique of extending a family of distributions by introducing new parameters is acknowledged in the statistical literature.The log-logistic distribution is an excellent choice for analyzing data with unimodal or decreasing failure rates. However, in a number of situations where data behave monotonically, such as increasing failure rates, or nonmonotonically, such as bathtub- or modified bathtub-shaped failure rates, the log-logistic model is not a good candidate model to use [2–5]. As a result, new extensions and generalizations to existing models are required for accurate and precise data modelling. As a result, numerous statistical techniques are intended to discover new modifications to classical models in order to achieve a better fit to the data of interest. Other techniques, on the other hand, such as the Sin-G family, Cos-G family, Tan-G family, and Sec-G family, provide a versatile generalization of the existing probability distribution without adding any extra parameters [6–10]. The techniques that are effectively used to extend the classical distributions without adding extra parameters are termed as “new trigonometric classes of probability distribution” [11].A random variableX is said to have a log-logistic distribution with shape parameter β>0 and scale parameter α>0,denoted by X ∼ LLog (α,β), if its cumulative distribution function (cdf) is defined by the following equation:(1)Fx;φ=x/αβ1+x/αβ,x≥0,φ>0.The probability density function (pdf) is given by(2)fx;φ=β/αx/αβ−11+x/αβ2,x≥0,φ>0.The reliability (survival) function is given by(3)Sx;φ=1−Fx;φ=x/α−β1+x/α−β,x≥0,φ>0.The failure (hazard) rate function is given by(4)hx;φ=fx;φSx;φ=fx;φ1−Fx;φ=β/αx/αβ−11+x/αβ,x≥0,φ>0.The reversed hazard rate function (also known as the retro function) is given by(5)rx;φ=fx;φFx;φβ/αx/αβ−1/1+x/αβ2x/αβ/1+x/αβ=θ/αx/α−11+x/αβ;x≥0,φ>0.The cumulative hazard rate function is given by(6)Hx=−logSx;ϕ=−logx/α−β1+x/α−β,x≥0,ϕ>0,where ϕ=α,β′ is a vector of parameters.It is well understood that the classical log-logistic distribution fails to capture the accurate phenomenon under investigation in many cases. As a result, several modifications have been proposed and researched. Inducting one or more parameters to the classical log-logistic distribution yields a modified form of the log-logistic distribution. When compared to the classic log-logistic distribution, several of these modified distributions have been found to be more flexible and adaptable of modelling real-life data. An up-to-date survey of recent modifications of the log-logistic distribution can be found in [12].Most modifications of the log-logistic model in the statistical literature have been derived by adding extra parameters to control the shape of the skewness (or asymmetry) and the kurtosis of the distribution; see the exponentiated LL distribution [13], beta LL distribution [14], gamma LL distribution [15], Marshall–Olkin LL distribution [16], transmuted LL distribution [17], cubic transmuted LL distribution [18], McDonald LL distribution [19], and alpha power transformed LL distribution [20, 21]. Other generalizations and modifications of the log-logistic distribution developed recently can be seen in [22, 23].The majority of techniques for extending the classical LL model produce a heavy-tailed distribution. Unfortunately, the abovementioned generalizations of the log-logistic distribution have some limitations. For instance, (i) adding extra parameter(s) to the distribution enhances flexibility, but even so, such practices usually result in reparameterization issues; (ii) the number of model parameters is increased, causing difficulty in estimating the model's parameters; (iii) some extending techniques reduce the tractability of the cdf, making manual computation of statistical properties more difficult; and (iv) other generalization techniques complicate the pdf, resulting in computational issues; incorporating new extra parameters into existing models increases flexibility, which is a desirable feature. On the other hand, it makes inferences more difficult [24, 25]. As a result, in this study, we modified the classical log-logistic distribution with a two-parameter continuous distribution referred to the log-logistic tangent distribution (LLT).The primary goal of this current study is it to introduce and investigate a versatile modification of the log-logistic distribution using the Tan-G family generator method. As previously stated, the proposed two-parameter distribution may be more flexible than other popular log-logistic model extensions. It would eventually become clear that the proposed distribution family accommodates both monotone and nonmonotone hazard rates.The rest of the paper is organized as follows. The proposed family is discussed in Section2. Section 3 presents the proposed distribution. Some mathematical and statistical properties of the proposed distribution are discussed in Section 4. Section 5 presents the estimation of the unknown parameters of the model. Section 6 discusses a Monte Carlo simulation for the proposed model. In Section 7, an application to a real-life dataset is presented and discussed. The comparison of some of the parametric probability distributions and their submodels is presented. Finally, concluding remarks and a work summary are presented in Section 8.


## 2. The Proposed Family

Several continuous probability distributions involving trigonometric functions have been developed in the statistics literature by many researchers, with the tangent distribution being notable due to its wide range of applications in many real-life datasets from various disciplines. A useful survey in this context can be found in [11, 24–27]. In this study, we concentrate on the trigonometric classes of continuous probability distributions, which are described by a cdf concerning trigonometric functions (tangent, cosine, sine, secant, and different mixtures of them). The Tan-G family developed in [11, 26–29] is the fundamental work.As indicated by its name, Souza et al. [28] introduced a new method for extending the classical probability distributions, leading to greater flexibility in analyzing and modelling different data types. They looked at a parent distribution, which is an arbitrary continuous probability distribution with a cdf and a corresponding pdf. The cumulative distribution function (cdf) is defined by the associated tangent family of distributions.(7)FTanx=∫0π/4Gxsec2tdt.Equation (7) can be written as(8)FTanx=tanπ4Gx,where Gx is the cdf of the parent (baseline) distribution, and if Gx has pdf gx, the pdf of the class is given by(9)fTanx=π4gxsec2π4Gx.The survival function is obtained by(10)STanx=1−tanπ4Gx.The failure rate function is given by(11)hTanx=fTanxSTanx=π/4gxsec2π/4Gx1−tanπ/4Gx.The reversed hazard rate function is given by(12)rTanx=fTanxFTanx=π/4gxsec2π/4Gxtanπ/4Gx.Also, the cumulative hazard rate function is given by(13)HTanx=−logSTanx=−log1−tanπ4Gx.Recently, Souza et al. [28] proposed and studied the Burr-XII tangent distribution which serves as a potential lifetime model. Ampadu [29] introduced and studied the Weibull tangent distribution with applications to health science data. Hence, the Tan-G family is a special family for extending the classical well-known lifetime distributions without adding an extra parameter.


## 3. The Log-Logistic Tangent Distribution

In this section, the LLT distribution has been examined, considering thatGxis the cdf of the log-logistic distribution.The distribution function (cdf) of the log-logistic tangent distribution, forx>0, can be expressed as(14)FTan−LLx;ϕ=tanπ4x/αβ1+x/αβ.The corresponding probability density function to the abovementioned cdf is given by(15)fTan−LLx;ϕ=π4β/αx/αβ−11+x/αβ2sec2π4x/αβ1+x/αβ.The survival function is expressed as follows:(16)STan−LLx;ϕ=1−tanπ4x/αβ1+x/αβ.The hazard rate function is obtained by(17)hTan−LLx;ϕ=π/4β/αx/αβ−1/1+x/αβ2sec2π/4x/αβ/1+x/αβ1−tanπ/4x/αβ/1+x/αβ.The reversed hazard rate function is expressed as follows:(18)rTan−LLx;ϕ=π/4β/αx/αβ−1/1+x/αβ2sec2π/4x/αβ/1+x/αβtanπ/4x/αβ/1+x/αβ.The cumulative hazard function can be given as follows:(19)HTan−LLx=−logSTan−LLx=−log1−tanπ4x/αβ1+x/αβ.Withx≥0,α,β>0,andφ=α,β' is the vector parameter in all of the abovementioned equations, respectively.Some possible shapes of the pdf, cdf, hazard function, survival function, and the reversed hazard functions of the log-logistic tangent distribution are displayed in Figures1–7.Figure 1
Visual display plots of the cdf of the LLT distribution.Figure 2
Visual display plots of the pdf of the LLT distribution.Figure 3
Plots of the survival function of the LLT distribution.Figure 4
Unimodal hazard rate of the LLT distribution.Figure 5
Plots of the LLT increasing hazard rate.Figure 6
Plots of the LLT decreasing hazard rate.Figure 7
Plots of the retro hazard function of the LLT distribution.The hazard rate function of the proposed distribution can accommodate for both monotone (increasing and decreasing) and nonmonotone (i.e., unimodal) hazard rates, as it can be seen in Figures4–6, and Figure 7 represents the reversed hazard function of the LLT distribution.


## 4. Some Statistical Properties of the Proposed Distribution

In this section, we use numerical examples to derive some mathematical properties of the log-logistic tangent distribution, such as the quantile function, skewness and kurtosis, moments, and residual and reverse residual life.


### 4.1. Quantile Function

For this model, the quantile function of the LLT distribution is used in theoretical aspects of distribution theory such as statistical simulations and applications. To generate random samples, the simulation algorithm used a quantile function by following the steps in Algorithm1.LetX be the Tan-G-distributed random variable. The quantiles can be utilized to obtain data of the distribution according to(20)x=Qu=F−1u=G−14πarctanu.Algorithm 1: Random generator for the Tan-G class.
(1)
Generateu∼U0,1(2)
SpecifyG−1.(3)
Obtain an outcome ofX by X=QuThe LLT distribution's quantile function is given by(21)QLLTu=α4/πarctanu1−4/πarctanu1β.The lower quartile (Q1), median (Q2), and upper quartile (Q3) of LLT distribution can be derived from the equation of the quantile function by setting u=1/4,1/2,and3/4, respectively.The quantiles of the LLT distribution for some parameter values are shown in Table1.Table 1
The quantiles of the LLT distribution for some values of the parameters.
Quantilesα,β(0.5, 1.5)(0.2, 0.9)(0.3, 0.8)(0.1, 1.2)(0.05, 0.9)0.10.13820.02350.02690.02010.00590.20.24150.05950.07670.04030.01490.30.35180.11130.15510.06440.02780.40.47970.18670.27760.09500.04670.50.63790.30010.47370.13560.07510.60.84730.48180.80650.19330.12040.71.15180.80361.43430.28380.20090.81.66881.49082.87450.45110.37270.92.89563.73518.07810.89830.9338


### 4.2. Skewness and Kurtosis

Some of the properties of the continuous distribution can be studied through its asymmetry and kurtosis. The mathematical form of the Moors Kurtosis and Galton asymmetry (or skewness) of the LLT model with two parameters is defined by the following relationship:(22)SK=Q3/4+Q1/4−2Q2/4Q3/4−Q1/4,KM=Q7/8+Q3/8−Q5/8−Q1/8Q6/8−Q2/8,where Q describes different quartile values.The equations above can be solved as functions of the LLT quantile function. These measures have the advantage of being less sensitive in the presence of outliers and existing even when the distribution is devoid of moments.


### 4.3. Moments

Moments are essential in statistical modelling, particularly in role in applications. The LLT distribution’srth moment is defined as(23)μr′=∫−∞∞xrfx;α,βdx.In fact, we have(24)μr′=∫−∞∞xrπ4β/αx/αβ−11+x/αβ2sec2π4x/αβ1+x/αβdx.The first five moments followed by the standard deviation, coefficient of variance, skewness, and kurtosis for some parameter values are shown in Table2.Table 2
First five moments, SD, CV, asymmetry, and kurtosis of the LLT distribution for different parameter values.
Momentsα,β(1.1, 0.)(1.2, 0.9)(0.3, 0.8)(0.1, 1.2)(0.2, 0.4)μ′10.13830.14450.17270.17720.0920μ′20.08040.08620.09140.07440.0471μ′30.05620.06080.06110.04400.0314μ′40.04310.04690.04570.03060.0236μ′50.03490.03810.03640.02340.0188SD0.24760.25550.24820.20740.1965CV1.79061.76851.43771.17092.1388CS1.85341.77021.57281.74082.6374CK5.35404.97624.51955.69509.5119


### 4.4. Residual Life and Reverse Residual Life

The residual lifetime function (rlf) has broader uses in risk management and survival analysis. The rlf of the LLT random variable (r.v.) can be expressed as follows:(25)Rtx=Sx+tSt,Rtx=1−tanπ/4x+t/αβ/1+x+t/αβ1−tanπ/4x/αβ/1+x/αβ.Furthermore, the reverse residual lifetime function (rrlf) of the LLT r.v. can be obtained as follows:(26)R^tx=Sx+tSt,R^tx=1−tanπ/4x−t/αβ/1+x−t/αβ1−tanπ/4x/αβ/1+x/αβ.


## 4.1. Quantile Function

For this model, the quantile function of the LLT distribution is used in theoretical aspects of distribution theory such as statistical simulations and applications. To generate random samples, the simulation algorithm used a quantile function by following the steps in Algorithm1.LetX be the Tan-G-distributed random variable. The quantiles can be utilized to obtain data of the distribution according to(20)x=Qu=F−1u=G−14πarctanu.Algorithm 1: Random generator for the Tan-G class.
(1)
Generateu∼U0,1(2)
SpecifyG−1.(3)
Obtain an outcome ofX by X=QuThe LLT distribution's quantile function is given by(21)QLLTu=α4/πarctanu1−4/πarctanu1β.The lower quartile (Q1), median (Q2), and upper quartile (Q3) of LLT distribution can be derived from the equation of the quantile function by setting u=1/4,1/2,and3/4, respectively.The quantiles of the LLT distribution for some parameter values are shown in Table1.Table 1
The quantiles of the LLT distribution for some values of the parameters.
Quantilesα,β(0.5, 1.5)(0.2, 0.9)(0.3, 0.8)(0.1, 1.2)(0.05, 0.9)0.10.13820.02350.02690.02010.00590.20.24150.05950.07670.04030.01490.30.35180.11130.15510.06440.02780.40.47970.18670.27760.09500.04670.50.63790.30010.47370.13560.07510.60.84730.48180.80650.19330.12040.71.15180.80361.43430.28380.20090.81.66881.49082.87450.45110.37270.92.89563.73518.07810.89830.9338


## 4.2. Skewness and Kurtosis

Some of the properties of the continuous distribution can be studied through its asymmetry and kurtosis. The mathematical form of the Moors Kurtosis and Galton asymmetry (or skewness) of the LLT model with two parameters is defined by the following relationship:(22)SK=Q3/4+Q1/4−2Q2/4Q3/4−Q1/4,KM=Q7/8+Q3/8−Q5/8−Q1/8Q6/8−Q2/8,where Q describes different quartile values.The equations above can be solved as functions of the LLT quantile function. These measures have the advantage of being less sensitive in the presence of outliers and existing even when the distribution is devoid of moments.


## 4.3. Moments

Moments are essential in statistical modelling, particularly in role in applications. The LLT distribution’srth moment is defined as(23)μr′=∫−∞∞xrfx;α,βdx.In fact, we have(24)μr′=∫−∞∞xrπ4β/αx/αβ−11+x/αβ2sec2π4x/αβ1+x/αβdx.The first five moments followed by the standard deviation, coefficient of variance, skewness, and kurtosis for some parameter values are shown in Table2.Table 2
First five moments, SD, CV, asymmetry, and kurtosis of the LLT distribution for different parameter values.
Momentsα,β(1.1, 0.)(1.2, 0.9)(0.3, 0.8)(0.1, 1.2)(0.2, 0.4)μ′10.13830.14450.17270.17720.0920μ′20.08040.08620.09140.07440.0471μ′30.05620.06080.06110.04400.0314μ′40.04310.04690.04570.03060.0236μ′50.03490.03810.03640.02340.0188SD0.24760.25550.24820.20740.1965CV1.79061.76851.43771.17092.1388CS1.85341.77021.57281.74082.6374CK5.35404.97624.51955.69509.5119


## 4.4. Residual Life and Reverse Residual Life

The residual lifetime function (rlf) has broader uses in risk management and survival analysis. The rlf of the LLT random variable (r.v.) can be expressed as follows:(25)Rtx=Sx+tSt,Rtx=1−tanπ/4x+t/αβ/1+x+t/αβ1−tanπ/4x/αβ/1+x/αβ.Furthermore, the reverse residual lifetime function (rrlf) of the LLT r.v. can be obtained as follows:(26)R^tx=Sx+tSt,R^tx=1−tanπ/4x−t/αβ/1+x−t/αβ1−tanπ/4x/αβ/1+x/αβ.


## 5. Estimation of the Parameters

The maximum likelihood approach is used in this section to estimate the unknown parameters of the log-logistic tangent distribution based on a complete sample. Letx=x1,x2,…,xnT represent n independent random variables drawn from the LLT distribution. The sample’s likelihood function is defined as(27)L=∏i=1nfxi,α,β,Lx;α,β=∏i=1nπ4β/αx/αβ−11+x/αβ2sec2π4x/αβ1+x/αβ.The log-likelihood function can be written as follows:(28)ℓ=nlogπ4+nlogβ−nβlogα+β−1∑i=1nlogxi−2∑i=1nlog1+xiαβ+2∑i=1nlogsecπ4xiαβ1+xiαβ−1.By taking the 1st derivatives of the log-likelihood function in equation (28) with respect to α and β parameters and setting the result to zero, we get(29)∂ℓ∂α=−nβα+2βα∑i=1nxiαβ1+xiαβ−1−π2∑i=1nxiαββα1+xiαβ−2tanπ4xiαβ1+xiαβ−1,∂ℓ∂β=nβ−nlogα+∑i=1nlogxi−2∑i=1nxiαβlogxiα1+xiαβ−1+π2∑i=1nxiαβlogxiαxiα1+xiαβ−2tanπ4xiαβ1+xiαβ−1.It is worth noting that the MLE’sα^andβ^ of αandβ,respectively, can be achieved by equating the outcomes to zero and solving the system of nonlinear equations numerically. The well-known theory on MLE can be used under some standard regularity conditions, ensuring nice asymptotic properties [30].


## 6. Simulation Study

In this section, we use Monte Carlo (MC) simulation to evaluate the effectiveness of the Maximum Likelihood Estimation (MLE) method for estimating the log-logistic tangent distribution parameters. The simulation study is carried out to investigate the average bias (AB), mean square error (MSE), and root mean square error (RMSE) for the parameters of the proposed model. The simulation experiment was conducted by running a number of simulations with varying sample sizes and parameter values. The quantile function given in Eq. (21) was used to generate random samples for the LLT. The MC simulation study was iterated 750 times with sample sizes n=25,50,…,750 and the parameter scenarios α=2.0andβ=1.5 in set I and α=3.0andβ=2.0 in set II.The MLEs are ascertained for each item of simulated data, say (α^,β^) for i=1,2,…,750, and the AB, MSEs, and RMSEs of the parameters were calculated by(30)AB=1N∑i=1Nθ^−θ,MSE=1N∑i=1Nθ^−θ2,RMSE=1N∑i=1Nθ^−θ2,where θ=αandβ.Table3 shows the AB, MSE, and RMSE values of the parameters for various sample sizes. Based on these results, we draw the conclusion that MLEs do a good job of estimating parameters and that the estimates seem to be reasonably constant and closer to the true values for these sample sizes. In addition, Table 4 shows that the RMSE decreases as the sample size increases, as expected. Furthermore, as sample size increases, so does the AB. As a result, even if the sample size is small, the MLEs and their asymptotic results can be used to calculate confidence intervals for model parameters.Table 3
MC simulation outcomes for the proposed distribution using the values for MLEs, ABs, MSEs, and RMSEs.
ParametersIIInMLEABMSERMSEMLEABMSERMSEα252.0590.0590.2730.5223.0600.0600.3680.607502.0300.0300.1170.3413.0220.0220.1420.3781002.0100.0100.0660.2563.0170.0170.0860.2942002.0080.0080.0300.1723.0130.0130.0410.2024002.0040.0040.0150.1243.0080.0080.0190.1376002.0030.0030.0100.1003.0030.0030.0120.1087002.0020.0020.0070.0863.0020.0020.0090.098β251.5840.0840.0790.2822.1240.1240.1560.395501.5410.0410.0350.1872.0600.0600.0650.2551001.5180.0180.0170.1302.0230.0230.0290.1732001.5130.0130.0080.0892.0110.0110.0140.1204001.5050.0050.0040.0622.0070.0070.0070.0856001.5020.0020.0030.0502.0040.0040.0040.0687001.5010.0010.0020.0442.0030.0030.0030.058The simulation results for the aforementioned measures are depicted in Figures8 and 9. These plots demonstrate that increasing the sample size n reduces the estimated biases. Furthermore, as the sample size n is increased, the estimated MSEs and RMSEs decay toward zero. These findings demonstrate the MLEs' efficiency as well as their consistency.Figure 8
MLEs, MSEs, biases, and RMSE plots of the proposed distribution for set I values.Figure 9
MLEs, MSEs, biases, and RMSE plots of the proposed distribution for set II values.


## 7. Applications to the COVID-19 Dataset

The log-logistic tangent distribution derivation's main interest is its application in data analysis goals, which makes it valuable in different disciplines, especially those associated with mortality data analysis. A number of proposed distributions for COVID-19 datasets have recently been proposed; for more information on these, please see [31–39].The log-logistic tangent distribution was applied to a real-life COVID-19 mortality rate dataset from Somalia for demonstration purposes, and its performance was compared to the performance of other fitted models such as the log-logistic, Weibull, Gompertz, and kappa distributions. For the log-likelihood, the most common information criteria including Akaike Information Criterion (AIC), Hannan–Quinn Information Criterion (HQIC), Corrected Akaike Information Criterion (CAIC), and Bayesian Information Criterion (BIC) values were used to select the most appropriate. The Anderson–Darling (A∗) statistic, the Cramer–von Mises (W∗) distance value, and the Kolmogorov–Smirnov (K–S) statistic, as well as the corresponding p value, are all recorded.(1)
Log-logistic:(31)fx=β/αx/αβ−11+x/αβ2,forx≥0.(2)
Gompertz distribution:(32)fx=αβexpβ+αx−βeαx,forx≥0.(3)
Kappa distribution:(33)fx=αβα+xβα−α+1/α,forx≥0.See also the other four recently modified log-logistic distributions with three parameters and four parameters [14, 16, 40, 41].The AIC is(34)AIC=2k−2l.The BIC is(35)BIC=klnn−2l.The CAIC is(36)CAIC=2nkn−k−1−2l.The HQIC is(37)HQIC=2klnlnn−2l,where the log-likelihood function is represented by l, the sample size by n, and the number of model parameters by k. The following goodness-of-fit measures are being considered.The Anderson–Darling (A∗) test statistic is given by(38)A∗=−n−1n∑i=1n2l−1×lnGXi+ln1−GXn−i+1.The Cramer–von Mises (W∗) test statistics is given by(39)W∗=112n+∑i=1n2i−12n+GXi2,where xi is the ith observation in the sample and n is the sample size. When the data are sorted in ascending order, xi is calculated.The best model has the lowest AIC, CAIC, BIC, and HQIC, a well as the lowestA∗, W∗, and K–S values. In addition, the best model is chosen as the one with the highest log-likelihood function value, and p values for the K–S statistics are applied to compare the competitive distributions.


### 7.1. Data I : Somalia COVID-19 Mortality Rate Dataset

The dataset contains COVID-19 mortality rate from Somalia during the time between 1st March 2021 to 20th April 2021 (see https://covid19.who.int/). These data are made up of a rough mortality rate. The data frame contains 51 observations, which are formed by using the daily cumulative cases (DCCs) and daily new deaths (DNDs). The data are as follows:3.008, 2.963, 4.762, 5.263, 11.382, 5.330, 10.471, 2.857, 4.274, 8.633, 5.882, 10.280, 5.330, 8.730, 7.377, 8.696, 8.257, 4.294, 5.330, 6.769, 11.627, 8.547, 9.302, 7.742, 6.173, 6.015, 5.330, 9.770, 5.330, 4.294, 5.330, 8.120, 7.547, 7.563, 5.330, 6.875, 8.800, 11.429, 5.330, 9.898, 5.330, 5.330, 9.630, 3.750, 5.330, 11.808, 5.330, 5.330, 5.330, 5.785, and 3.265.


### 7.2. Exploratory Data Analysis

The primary goal of data analysis is to extract information from the data. We used five different techniques to perform exploratory data analysis in this work: (1) data descriptive statistics; (2) Box plot; (3) TTT plot; (4) histogram; and (5) time-series plot. Figure10 displays the Box, TTT, histogram, and times-series plots for our dataset.Figure 10
Box plot, TTT plot, histogram, and time-series plot for dataset I.The total-time-on-test (TTT) plot is a graphical method for determining the failure rate function's shape. There is qualitative information about the shape of the failure rate function in many real-world applications that can aid in the selection of a particular distribution. The TTT plot for our dataset in this work is shown, and it shows an increasing failure rate shape.Table4 shows us the descriptive statistics of the COVID-19 data for mortality rate between 1st March to 20 April 2021 by computing specific aspects of the data (central tendency and spread).Table 4
Descriptive measures of data I.
NMin.MedianModeVarianceAsymmetryKurtosisMeanMax.CVGeo. mean512.8575.8825.55.9620.450−0.7326.79311.810.3596.363The proposed LLT distribution has the lowest AIC, CAIC, BIC, and HQIC values and the highest log-likelihood values in Table5, so it is selected as the most suitable model among the competing distributions considered in this work.Table 5
Likelihood and the information criteria for data I.
DistributionℓAICBICCAICHQICLLT−116.839236.677240.541237.927238.154LL−118.198238.997241.860239.246239.473Gompertz−120.495244.991248.854245.8241246.467Kappa−127.196258.391262.255258.641259.867Expo LL−116.995239.991245.786240.501242.205Beta LL−115.559239.119246.846239.988242.071Kumw LL−115.970239.940247.667240.810242.892Mo-LL−116.998239.996245.792240.507242.212Table6 shows the parameter estimate and p value for the Cramer–von Mises (W∗), Anderson–Darling (A∗), and Kolmogorov–Smirnov (K–S) tests for all competing distributions using abovementioned dataset I. According to Table 6, the proposed LLT distribution has the lowest A∗, W∗, and K–S tests, as well as the highest p value. As a result, among the competing distributions considered in this study, the proposed LLT distribution is chosen as the most appropriate model.Table 6
MLE model parameter estimators and goodness-of-fit tests for the COVID-19 data.
DistributionsEstimatesW∗A∗K−SPvalueLLT (δ,θ)δ = 5.904θ  = 4.8140.1991.0670.169 (0.206)llog (α, β)α = 6.388β  = 4.6330.2031.0890.169 (0.109)Gompertz (α, β)α = 0.365β  = 0.0210.2211.2660.181 (0.016)Kappa (α, β)α = 58.227β  = 11.0710.3221.8560.267 (0.002)Exp LL (α, β, γ)α = 1.050β  = 4.561γ = 6.2850.2101.0850.197 (0.113)Beta LL (α, β, γ, δ)α = 5.431β  = 20.368γ = 1.360δ  = 17.7050.2141.0750.178 (0.104)Kum LL (α, β, γ, δ)α = 2.319β  = 5.749γ = 2.182δ  = 7.9030.2011.0710.161 (0.101)Mo-LL (α, β, γ)α = 6.064β  = 4.633γ = 4.3300.2031.0890.169 (0.109)The estimated pdf for the competing models is shown in Figure11, the estimated cdf for the competing models is shown in Figure 12, and the estimated cdf, Kaplan–Meier, pdf, and PP plots for the proposed model are shown in Figure 13.Figure 11
Some fitted pdf’s of the fitted distributions in dataset I.Figure 12
Some fitted cdf’s of the fitted distributions in dataset I.Figure 13
The estimated cdf, Kaplan–Meier, estimated pdf, and PP plots of the log-logistic tangent distribution for dataset I.


## 7.1. Data I : Somalia COVID-19 Mortality Rate Dataset

The dataset contains COVID-19 mortality rate from Somalia during the time between 1st March 2021 to 20th April 2021 (see https://covid19.who.int/). These data are made up of a rough mortality rate. The data frame contains 51 observations, which are formed by using the daily cumulative cases (DCCs) and daily new deaths (DNDs). The data are as follows:3.008, 2.963, 4.762, 5.263, 11.382, 5.330, 10.471, 2.857, 4.274, 8.633, 5.882, 10.280, 5.330, 8.730, 7.377, 8.696, 8.257, 4.294, 5.330, 6.769, 11.627, 8.547, 9.302, 7.742, 6.173, 6.015, 5.330, 9.770, 5.330, 4.294, 5.330, 8.120, 7.547, 7.563, 5.330, 6.875, 8.800, 11.429, 5.330, 9.898, 5.330, 5.330, 9.630, 3.750, 5.330, 11.808, 5.330, 5.330, 5.330, 5.785, and 3.265.


## 7.2. Exploratory Data Analysis

The primary goal of data analysis is to extract information from the data. We used five different techniques to perform exploratory data analysis in this work: (1) data descriptive statistics; (2) Box plot; (3) TTT plot; (4) histogram; and (5) time-series plot. Figure10 displays the Box, TTT, histogram, and times-series plots for our dataset.Figure 10
Box plot, TTT plot, histogram, and time-series plot for dataset I.The total-time-on-test (TTT) plot is a graphical method for determining the failure rate function's shape. There is qualitative information about the shape of the failure rate function in many real-world applications that can aid in the selection of a particular distribution. The TTT plot for our dataset in this work is shown, and it shows an increasing failure rate shape.Table4 shows us the descriptive statistics of the COVID-19 data for mortality rate between 1st March to 20 April 2021 by computing specific aspects of the data (central tendency and spread).Table 4
Descriptive measures of data I.
NMin.MedianModeVarianceAsymmetryKurtosisMeanMax.CVGeo. mean512.8575.8825.55.9620.450−0.7326.79311.810.3596.363The proposed LLT distribution has the lowest AIC, CAIC, BIC, and HQIC values and the highest log-likelihood values in Table5, so it is selected as the most suitable model among the competing distributions considered in this work.Table 5
Likelihood and the information criteria for data I.
DistributionℓAICBICCAICHQICLLT−116.839236.677240.541237.927238.154LL−118.198238.997241.860239.246239.473Gompertz−120.495244.991248.854245.8241246.467Kappa−127.196258.391262.255258.641259.867Expo LL−116.995239.991245.786240.501242.205Beta LL−115.559239.119246.846239.988242.071Kumw LL−115.970239.940247.667240.810242.892Mo-LL−116.998239.996245.792240.507242.212Table6 shows the parameter estimate and p value for the Cramer–von Mises (W∗), Anderson–Darling (A∗), and Kolmogorov–Smirnov (K–S) tests for all competing distributions using abovementioned dataset I. According to Table 6, the proposed LLT distribution has the lowest A∗, W∗, and K–S tests, as well as the highest p value. As a result, among the competing distributions considered in this study, the proposed LLT distribution is chosen as the most appropriate model.Table 6
MLE model parameter estimators and goodness-of-fit tests for the COVID-19 data.
DistributionsEstimatesW∗A∗K−SPvalueLLT (δ,θ)δ = 5.904θ  = 4.8140.1991.0670.169 (0.206)llog (α, β)α = 6.388β  = 4.6330.2031.0890.169 (0.109)Gompertz (α, β)α = 0.365β  = 0.0210.2211.2660.181 (0.016)Kappa (α, β)α = 58.227β  = 11.0710.3221.8560.267 (0.002)Exp LL (α, β, γ)α = 1.050β  = 4.561γ = 6.2850.2101.0850.197 (0.113)Beta LL (α, β, γ, δ)α = 5.431β  = 20.368γ = 1.360δ  = 17.7050.2141.0750.178 (0.104)Kum LL (α, β, γ, δ)α = 2.319β  = 5.749γ = 2.182δ  = 7.9030.2011.0710.161 (0.101)Mo-LL (α, β, γ)α = 6.064β  = 4.633γ = 4.3300.2031.0890.169 (0.109)The estimated pdf for the competing models is shown in Figure11, the estimated cdf for the competing models is shown in Figure 12, and the estimated cdf, Kaplan–Meier, pdf, and PP plots for the proposed model are shown in Figure 13.Figure 11
Some fitted pdf’s of the fitted distributions in dataset I.Figure 12
Some fitted cdf’s of the fitted distributions in dataset I.Figure 13
The estimated cdf, Kaplan–Meier, estimated pdf, and PP plots of the log-logistic tangent distribution for dataset I.


## 8. Conclusions

The two-parameter log-logistic model has found widespread use in statistical sciences, especially survival analysis, biomedical sciences, engineering, and actuarial sciences. In this paper, we proposed a new probability distribution that derives from a combination of log-logistic distribution and a Tan-G trigonometric class of distribution. This model was named log-logistic tangent distribution “LLT distribution” and has been successfully extended in this work. Reliability function, pdf, cdf, failure rate function, reversed failure rate function, and the cumulative failure rate function were carefully derived, and expressions for the basic statistical properties for the log-logistic tangent distribution have been developed. The log-logistic tangent distribution was applied to the COVID-19 mortality rate dataset and provided a better fit than the known two-parameter, three-parameter, and four-parameter competitors. Gompertz, log-logistic, kappa, exponentiated log-logistic, Marshall–Olkin log-logistic, Kumaraswamy log-logistic, and beta log-logistic based on goodness-of-fit tests, selection criteria such as AIC, BIC, CAIC, and HQIC values, and the log-likelihood value were applied. We, therefore, conclude that the log-logistic tangent distribution is the most adequate model among the ones considered, as well as a very competitive model for describing different datasets in different areas of application.



---

*Source: 8640794-2021-11-13.xml*