# Comparison of COVID-19 Pandemic Dynamics in Asian Countries with Statistical Modeling

Authors: Min Zuo; Saima K. Khosa; Zubair Ahmad; Zahra Almaspoor
Journal: Computational and Mathematical Methods in Medicine
 (2020)

---

## Abstract

In the current scenario, the outbreak of a pandemic disease COVID-19 is of great interest. A broad statistical analysis of this event is still to come, but it is immediately needed to evaluate the disease dynamics in order to arrange the appropriate quarantine activities, to estimate the required number of places in hospitals, the level of individual protection, the rate of isolation of infected persons, and among others. In this article, we provide a convenient method of data comparison that can be helpful for both the governmental and private organizations. Up to date, facts and figures of the total the confirmed cases, daily confirmed cases, total deaths, and daily deaths that have been reported in the Asian countries are provided. Furthermore, a statistical model is suggested to provide a best description of the COVID-19 total death data in the Asian countries.


---

## Body

## 1. Introduction

Coronavirus disease (COVID-19) is an infectious disease caused by a newly discovered coronavirus. The name “coronavirus” is derived from the Latin corona, meaning “crown” or “wreath.” The name refers to the characteristic appearance of virions (the infective form of the virus) by electron microscopy.Coronaviruses were first discovered in the 1930s when an acute respiratory infection of domesticated chickens was caused by infectious bronchitis virus (IBV). Later, in the 1940s, two more animal coronaviruses, mouse hepatitis virus (MHV) and transmissible gastroenteritis virus (TGEV), were isolated [1]. For the first time, human coronaviruses were discovered in the 1960s [2]. The earliest ones studied were from human embryonic tracheal organ cultures obtained from the respiratory tract of an adult with a common cold, which were later named human coronavirus 229E and human coronavirus OC43 [3]. Other human coronaviruses have since been identified, including SARS-CoV in 2003, HCoV NL63 in 2004, HKU1 in 2005, and MERS-CoV in 2012 ([4]). Most of these have involved serious respiratory tract infections.Recently, a new type of coronaviruses observed in a place called Wuhan city of China, which has a well-known seafood wholesale market, where a large number of people come to sell or buy live seafood. On 31 December 2019, the Wuhan Municipal Health Commission (WMHC) reported a bunch of 27 pneumonia cases of unknown aetiology. Later, on 11 January 2020, the World Health Organization (WHO) named this novel coronavirus as SARS-CoV-2, the virus causing COVID-19, see [5].The SARS-CoV, MERS-CoV, and COVID-2019 viruses are highly pathogenic Betacoronaviruses and responsible for causing a respiratory and gastrointestinal syndrome. The average incubation period for coronavirus infection is 5 days, with an interval that can reach up to 16 days. The transmissibility of patients infected with SARSCoV is on average 7 days after the onset of symptoms. However, preliminary data from COVID-19 suggests that transmission may occur, even without the appearance of signs and symptoms; seehttps://en.wikipedia.org/wiki/Coronavirus_disease_2019.


## 2. Detail and Comparison of COVID-19 Cases in the Asian Countries

Asia is one of the most affected region due to COVID-19. In this section, we provide the detailed information and comparison of the total cases, total deaths, total recovered, and active cases in Asian countries. The detail description of the total cases, total deaths, total recovered, and active cases in the Asian countries up to 8th April 2020, are provided in Tables1 and 2. For details, we refer to https://www.worldometers.info/coronavirus/#countries. Note that the graphical visualization of total cases, total deaths, total recovered, and active cases of the COVID-19 of the Asian countries are displayed in Appendix A. We provide a very simple method for comparison which is not only limited to the Asian countries but it can also be applied for every country to analyze the impact of the disease.Table 1
Detailed report and comparison of COVID-19 cases in the Asian countries.
CountriesTotal casesTotal deathsTotal recoveredActive casesAfghanistan (AF)4441429401Azerbaijan (AZ)822863751Bahrain (BH)8215467349Bangladesh (BD)2182033165Bhutan (BT)5023Brunei (BN)13519143Cambodia (KH)11706354China (CN)818023333772791190Cyprus (CY)526952465Hong Kong (HG)9614264693India (IN)57491785065065Indonesia (ID)29562402222494Iran (IR)6458639932981230781Iraq (IQ)120269452681Israel (IL)9409728018531Japan (JP)4257936223542Jordan (JO)3586150202Kazakhstan (KZ)718754657Kuwait (KW)8551111743Kyrgyzstan (KG)270433233Laos (LA)150015Lebanon (LB)5751962494Macao (MO)4501035Malaysia (MY)41196514872567Maldives (MV)190136Mongolia (MN)160412Myanmar (MM)223019Nepal (NP)9018Oman (OM)419272345Pakistan (PK)4196604673669Philippines (PH)3870182963592Qatar (QA)221061782026Saudi Arabia (SA)3122416312450Table 2
Detailed report and comparison of COVID-19 cases in the Asian countries.
CountriesTotal casesTotal deathsTotal recoveredActive casesSingapore (SG)162364061211South Korea (SK)1038420067763408Sri Lanka (LK)189744138Palestine (PS)263144218Syria (SY)192314Taiwan (TW)379567307Thailand (TH)2369308881451Timor-Leste (TL)1001Turkey (TR)38226812184635568United Arab Emirates (AE)2659122392408Uzbekistan (UZ)545330512Vietnam (VN)2510126125


## 3. Proposed Family of Statistical Models

In the practice of big data sciences, particularly in statistical theory, there has be an increased interest in defining new statistical models or new families of statistical models to provide a better description of the problems under consideration; see [6, 7]. For more details, we refer to [8].Often, adding extra parameter(s) gives more flexibility to a class of distribution functions, improves the characteristics, and provides better fits to the real-life data than the other modified models. But, unfortunately, on the other hand, the reparametrization problem arises. To avoid such problems and provide a better description of real phenomena of nature, we further carry this branch of statistical theory and propose a new class of statistical models. The proposed class of distributions may be called a new flexible extended-X (NFE-X) class of distributions.Letpt be the density of a random variable T∈a1,a2 for −∞≤a1<a2<∞ and let KFx;ξ be a function of Fx;ξ of a random variable X. The cumulative distribution function (cdf) of the T-X family of distributions [9] is given by
(1)Gx=∫a1KFx;ξptdt,x∈ℝ,where KFx;ξ fulfills some certain conditions, see [9]. The density function corresponding to (1) is
(2)gx=∂∂xKFx;ξpKFx;ξ,x∈ℝ.Ifpt=1−e−t,t≥0, and setting KFx;ξ=−log1−Fx;ξ2/eFx;ξ2 in (1), we get the cdf of the proposed class of distributions. The random variable X is said to have a NFE-X class of distributions, if the cumulative distribution function (cdf) of X, denoted by Gx;ξ is given by
(3)Gx;ξ=1−1−Fx;ξ2eFx;ξ2,x,ξ∈ℝ.The density function corresponding to (3) is
(4)gx;ξ=2fx;ξFx;ξeFx;ξ22−Fx;ξ2,x∈ℝ.One of the most prominent motivations of the proposed approach is to introduce a new class of distributions without adding additional parameter results in avoiding rescaling problems. The next section offers, a special submodel of the proposed class called a new flexible extended-Weibull (NFE-Weibull) distribution and investigates the graphical behaviour of its density function.


## 4. Submodel Description

This section offers a special submodel of the NFE-X class of distributions. Let Fx;ξ be the distribution function of the Weibull model given by Fx;ξ=1−e−ηxθ,x≥0,η,θ>0, where ξ=η,θ. Then, the cdf of the NFE-Weibull has the expression given by
(5)Gx;ξ=1−1−1−e−ηxθ2e1−e−ηxθ2,x≥0,η,θ>0,with density function
(6)gx;ξ=2ηθxθ−1e−ηxθ1−e−ηxθe1−e−ηxθ22−1−e−ηxθ2,x>0.For different values of the model parameters, plots of the density function of the NFE-Weibull model are sketched in Figure1.Figure 1
Different plots for the density function of the NFE-Weibull distribution.


## 5. Mathematical Properties

In this section, some mathematical and statistical properties of the NFE-Weibull distribution derived are discussed.


### 5.1. Quantile Function

The quantile function of the NFE-X family is the function Qu;ξ that satisfies the nonlinear equation
(7)QGu;ξ;ξ=u,u∈0,1.By using (3) in (7), after some algebraic manipulation, we get
(8)QGu;ξ;ξ=G−1u=F−1u,where t is the solution of log1−u+Fx;ξ2−log1−Fx;ξ2.


### 5.2. Moments

Here, we derive some of the moments for the NFE-X family. For the sake of simplicity we omit the dependency of gx;ξ and Gx;ξ on the parameter vector ξ. The density (4) can be represented as follows:
(9)gx=2fx;ξFx;ξeFx;ξ22−Fx;ξ2,gx=∑i=01−1i24−i/21ifx;ξFx;ξ2i+1eFx;ξθ,gx=∑i=01∑j=0∞−1i+j24−i/2j!1ifx;ξFx;ξ2i+j+1,Using the pdf and cdf of the Weibull distribution in (9), we get
(10)gx=∑i=01∑j=0∞∑k=02i+j+1−1i+j+k24−i/2j!1i2i+j+1kτk,η,θ,where τk,η,θ=ηθxθ−1e−ηk+1xθ. For any positive integer r, the rth moment of the NFE-Weibull distribution is given by
(11)μr′=EXr=∫0∞xrgxdx.On using (10) in (11), we get the rth moment of the NFE-Weibull distribution.Forr=1,2,3,4 we get the first four moments of the NFE-X distributions. The effects of the shape parameters on the skewness and kurtosis can be detected on the moments. Based on moments, we obtain skewness and kurtosis measures of the NFE-Weibull distribution. The skewness of the NFE-Weibull distribution is obtained as using the following expression:
(12)Skewness=μ3μ23/2,where μ2 and μ3 are the second and third moments of the random variable X with pdf (6). Furthermore, the kurtosis of X is derived as
(13)Kurtosis=μ4μ22,where μ4 is the fourth moment of X. These measures are less sensitive to outliers. Plots for the mean, variance, skewness, and kurtosis of the NFE-Weibull distribution are displayed in Figures 2 and 3.Figure 2
Plots of mean, variance, skewness, and kurtosis of the NFE-Weibull distribution.Figure 3
Plots of mean, variance, skewness, and kurtosis of the NFE-Weibull distribution.


### 5.3. On Other Means and Moments

Witht>0, the following result proposes an expansion of the primitive
(14)∫0txrgxdx=∑i=01∑j=0∞∑k=02i+j+1−1i+j+k24−i/2j!1i2i+j+1kτr,k,η,θ,where τr,k,η,θ=∫0tηθxr+θ−1e−ηk+1xθ.Several crucial conditional moments can be obtained using the integral∫0txrgxdx for various values of r. The most useful of them are presented below. For any t>0,
(i)
Therth conditional moments of X is given by,
(15)EXr∣X>t=11−Gt∫t+∞xrgxdx=11−GtEXr+∫0txrgxdx(ii)
Therth reversed moments of X is given by
(16)EXr∣X≤t=1Gt∫0txrgxdx(iii)
The mean deviations ofX about the mean, say μ is given by
(17)δ=EX−μ=2μGμ−2∫0μxgxdxwhere μ=EX.(iv)
The mean deviations ofX about the median, say M is given by
(18)τ=EX−M=μ−2∫0MxgxdxThe residual life parameters can be also determined usingEXr and ∫0txrgx;θ,ηdx for several values of r. In particular,(v)
The mean residual life is defined as(19)Kt=EX−t∣X>t=1StEX−∫0txgxdx−tand the variance residual life is given by
(20)Vt=VarX−t∣X>t=1StEX2−∫0tx2gxdx−t2−2tKt−Kt2(vi)
The mean reversed residual life is defined as(21)Lt=Et−X∣X≤t=t−1Gt∫0txgxdxand the variance reversed residual life is defined as
(22)Wt=Vart−X∣X≤t=1Gt∫0tx2gxdx+2tLt−Lt2−t2


## 5.1. Quantile Function

The quantile function of the NFE-X family is the function Qu;ξ that satisfies the nonlinear equation
(7)QGu;ξ;ξ=u,u∈0,1.By using (3) in (7), after some algebraic manipulation, we get
(8)QGu;ξ;ξ=G−1u=F−1u,where t is the solution of log1−u+Fx;ξ2−log1−Fx;ξ2.


## 5.2. Moments

Here, we derive some of the moments for the NFE-X family. For the sake of simplicity we omit the dependency of gx;ξ and Gx;ξ on the parameter vector ξ. The density (4) can be represented as follows:
(9)gx=2fx;ξFx;ξeFx;ξ22−Fx;ξ2,gx=∑i=01−1i24−i/21ifx;ξFx;ξ2i+1eFx;ξθ,gx=∑i=01∑j=0∞−1i+j24−i/2j!1ifx;ξFx;ξ2i+j+1,Using the pdf and cdf of the Weibull distribution in (9), we get
(10)gx=∑i=01∑j=0∞∑k=02i+j+1−1i+j+k24−i/2j!1i2i+j+1kτk,η,θ,where τk,η,θ=ηθxθ−1e−ηk+1xθ. For any positive integer r, the rth moment of the NFE-Weibull distribution is given by
(11)μr′=EXr=∫0∞xrgxdx.On using (10) in (11), we get the rth moment of the NFE-Weibull distribution.Forr=1,2,3,4 we get the first four moments of the NFE-X distributions. The effects of the shape parameters on the skewness and kurtosis can be detected on the moments. Based on moments, we obtain skewness and kurtosis measures of the NFE-Weibull distribution. The skewness of the NFE-Weibull distribution is obtained as using the following expression:
(12)Skewness=μ3μ23/2,where μ2 and μ3 are the second and third moments of the random variable X with pdf (6). Furthermore, the kurtosis of X is derived as
(13)Kurtosis=μ4μ22,where μ4 is the fourth moment of X. These measures are less sensitive to outliers. Plots for the mean, variance, skewness, and kurtosis of the NFE-Weibull distribution are displayed in Figures 2 and 3.Figure 2
Plots of mean, variance, skewness, and kurtosis of the NFE-Weibull distribution.Figure 3
Plots of mean, variance, skewness, and kurtosis of the NFE-Weibull distribution.


## 5.3. On Other Means and Moments

Witht>0, the following result proposes an expansion of the primitive
(14)∫0txrgxdx=∑i=01∑j=0∞∑k=02i+j+1−1i+j+k24−i/2j!1i2i+j+1kτr,k,η,θ,where τr,k,η,θ=∫0tηθxr+θ−1e−ηk+1xθ.Several crucial conditional moments can be obtained using the integral∫0txrgxdx for various values of r. The most useful of them are presented below. For any t>0,
(i)
Therth conditional moments of X is given by,
(15)EXr∣X>t=11−Gt∫t+∞xrgxdx=11−GtEXr+∫0txrgxdx(ii)
Therth reversed moments of X is given by
(16)EXr∣X≤t=1Gt∫0txrgxdx(iii)
The mean deviations ofX about the mean, say μ is given by
(17)δ=EX−μ=2μGμ−2∫0μxgxdxwhere μ=EX.(iv)
The mean deviations ofX about the median, say M is given by
(18)τ=EX−M=μ−2∫0MxgxdxThe residual life parameters can be also determined usingEXr and ∫0txrgx;θ,ηdx for several values of r. In particular,(v)
The mean residual life is defined as(19)Kt=EX−t∣X>t=1StEX−∫0txgxdx−tand the variance residual life is given by
(20)Vt=VarX−t∣X>t=1StEX2−∫0tx2gxdx−t2−2tKt−Kt2(vi)
The mean reversed residual life is defined as(21)Lt=Et−X∣X≤t=t−1Gt∫0txgxdxand the variance reversed residual life is defined as
(22)Wt=Vart−X∣X≤t=1Gt∫0tx2gxdx+2tLt−Lt2−t2


## 6. Maximum Likelihood Estimation and Monte Carlo Simulation

The section deals with the estimation of the model parameters and Monte Carlo simulation to assess the performance of the estimators.


### 6.1. Maximum Likelihood Estimation

The maximum likelihood estimation procedure is the commonly employed method of estimating the model parameters. The estimators that are obtained based on this procedure enjoy desirable asymptotic properties, and therefore, they are often utilized to obtain confidence intervals (CI) and test of statistical hypotheses. Suppose thatx1,x2,⋯,xn be the observed values of a random sample of size n obtained from (4). The corresponding log-likelihood function can be expressed as
(23)ℓξ=nlog2+∑i=0nlogfxi;ξ+∑i=0nlogFxi;ξ−∑i=0nFxi;ξ2+∑i=0nlog2−Fxi;ξ2.The log-likelihood function can be maximized either directly or by solving the nonlinear likelihood function obtained by differentiating. The first-order partial derivative of the log-likelihood function with respect toξ is given by
(24)∂∂ξℓξ=∑i=0n∂fxi;ξ/∂ξfxi;ξ+∑i=0n∂Fxi;ξ/∂ξFxi;ξ−2∑i=0nFxi;ξ∂Fxi;ξ∂ξ−2∑i=0nFxi;ξ∂Fxi;ξ/∂ξ2−Fxi;ξ2.Setting∂/∂ξℓξ equal to zero and solving numerically yields the maximum likelihood estimators (MLEs) of ξ=θ,η. An optimization software such as the R function optim or nlminb can be used to find ξ^ that minimizes the negative log-likelihood function (i.e., maximizes the log-likelihood function). Although the specification of the derivatives is optional in these R functions, fast and rapid convergence may be achieved if the expressions for the negative log-likelihood function are provided. In our implementation (R codes are given in Appendix B), we use optim() R-function with the argument method = “SANN” to obtain the MLEs.


### 6.2. Monte Carlo Simulation

A numerical investigation is established to examine the behaviour of MLEs for the NFE-Weibull model. For different sample sizes, measures like biases, absolute biases, and mean square errors (MSEs) are calculated to evaluate the performance of the estimators.(i)
We generate 500 from NFE-Weibull distribution of sizes;n=25,50,⋯,500(ii)
An optimization algorithm requires a set of initial values for the parameters. Certain values of the model parametersθ,η are chosen as Set1:θ=0.5,η=1; Set2:θ=1.2,η=1; and Set3:θ=0.8,η=1.2(iii)
MLEs of the parametersθ and η are calculated for each n and for all sets(iv)
Calculate the biases, absolute biases, and MSE for eachnThe simulation results are displayed in Figures4–6.Figure 4
Plots of MLEs, MSEs, biases, and absolute biases forθ=0.5 and η=1.Figure 5
Plots of MLEs, MSEs, biases, and absolute biases forθ=1.2 and η=1.Figure 6
Plots of MLEs, MSEs, biases, and absolute biases forθ=0.8 and η=1.2.


## 6.1. Maximum Likelihood Estimation

The maximum likelihood estimation procedure is the commonly employed method of estimating the model parameters. The estimators that are obtained based on this procedure enjoy desirable asymptotic properties, and therefore, they are often utilized to obtain confidence intervals (CI) and test of statistical hypotheses. Suppose thatx1,x2,⋯,xn be the observed values of a random sample of size n obtained from (4). The corresponding log-likelihood function can be expressed as
(23)ℓξ=nlog2+∑i=0nlogfxi;ξ+∑i=0nlogFxi;ξ−∑i=0nFxi;ξ2+∑i=0nlog2−Fxi;ξ2.The log-likelihood function can be maximized either directly or by solving the nonlinear likelihood function obtained by differentiating. The first-order partial derivative of the log-likelihood function with respect toξ is given by
(24)∂∂ξℓξ=∑i=0n∂fxi;ξ/∂ξfxi;ξ+∑i=0n∂Fxi;ξ/∂ξFxi;ξ−2∑i=0nFxi;ξ∂Fxi;ξ∂ξ−2∑i=0nFxi;ξ∂Fxi;ξ/∂ξ2−Fxi;ξ2.Setting∂/∂ξℓξ equal to zero and solving numerically yields the maximum likelihood estimators (MLEs) of ξ=θ,η. An optimization software such as the R function optim or nlminb can be used to find ξ^ that minimizes the negative log-likelihood function (i.e., maximizes the log-likelihood function). Although the specification of the derivatives is optional in these R functions, fast and rapid convergence may be achieved if the expressions for the negative log-likelihood function are provided. In our implementation (R codes are given in Appendix B), we use optim() R-function with the argument method = “SANN” to obtain the MLEs.


## 6.2. Monte Carlo Simulation

A numerical investigation is established to examine the behaviour of MLEs for the NFE-Weibull model. For different sample sizes, measures like biases, absolute biases, and mean square errors (MSEs) are calculated to evaluate the performance of the estimators.(i)
We generate 500 from NFE-Weibull distribution of sizes;n=25,50,⋯,500(ii)
An optimization algorithm requires a set of initial values for the parameters. Certain values of the model parametersθ,η are chosen as Set1:θ=0.5,η=1; Set2:θ=1.2,η=1; and Set3:θ=0.8,η=1.2(iii)
MLEs of the parametersθ and η are calculated for each n and for all sets(iv)
Calculate the biases, absolute biases, and MSE for eachnThe simulation results are displayed in Figures4–6.Figure 4
Plots of MLEs, MSEs, biases, and absolute biases forθ=0.5 and η=1.Figure 5
Plots of MLEs, MSEs, biases, and absolute biases forθ=1.2 and η=1.Figure 6
Plots of MLEs, MSEs, biases, and absolute biases forθ=0.8 and η=1.2.


## 7. Modeling COVID-19 Total Deaths of the Asian Countries

We mentioned earlier that a broad statistical analysis of the events that occurred due to COVID-19 is still to come. But, now it is immediately needed to propose a suitable model to provide a better description of the COVID-19 total death data to estimate the required number of places in hospitals, the level of individual protection, the rate of isolation of infected persons, etc. In this section, we model the COVID-19 total deaths that have occurred in the Asian countries up to April 8, 2020. The NFE-Weibull distribution applied to this dataset in comparison with the other well-known distributions such as the two-parameter Weibull, three-parameter Marshall-Olkin Weibull (MOW), and exponentiated Weibull (EW) distributions. It is important to emphasize that the EW distribution is a popular model for analyzing data in the applied areas, particularly in medical sciences, see [10]. The MOW distribution is another nonnested model and offers the characteristics of the Weibull and gamma distributions, see [11]. The cdfs of the competing distributions are as follows:
(1)
Weibull distribution(25)Gx;θ,η=1−e−ηxθ,x≥0,θ,η>0(2)
EW distribution(26)Gx;a,θ,η=1−e−ηxθa,x≥0,a,θ,η>0(3)
MOW distribution(27)Gx;θ,η,σ=1−e−ηxθ1−1−σ1−e−ηxθ,x≥0,θ,η,σ>0Selection of an appropriate approximation model is desirable to assign some preference to the alternatives. Therefore, we consider certain analytical measures in order to verify which distribution fits better the considered data. These analytical measures include (i) four discrimination measures such as the Akaike information criterion (AIC), Bayesian information criterion (BIC), Hannan-Quinn information criterion (HQIC), and consistent Akaike information criterion (CAIC) and (ii) three other goodness-of-fit measures including the Anderson Darling (AD) test statistic, Cramer-Von-Messes (CM) test statistic, and Kolmogorov-Smirnov (KS) test statistics with correspondingp values. A model with lowest values for these statistics is considered a best candidate model. The formulae for these measures can be explored in [12]For the COVID-19 total death data of the Asian countries, the estimates with the standard error (in parentheses) of the model parameters are provided in Table3. The analytical measures of the NFE-Weibull and other considered models are provided in Tables 4 and 5.Table 3
Estimated parameters along with standard errors (in parenthesis) of the fitted models.
DistributionθηaσNFE-Weibull0.314 (0.0392)0.281 (0.0567)Weibull0.424 (0.0484)0.178 (0.5544)MOW0.660 (0.0442)0.014 (0.4690)0.113 (1.8756)EW0.173 (0.0323)1.794 (0.5113)13.375 (7.858)Table 4
Discrimination measures of the NFE-Weibull model and other competing models.
DistributionAICBICCAICHQICNFE-Weibull386.178389.400386.531387.314Weibull394.615397.837394.968395.751MOW389.535392.367388.262389.238EW392.503395.336392.231393.207Table 5
Goodness-of-fit measures of the NFE-Weibull model and other competing models.
DistributionCMADKSp valueNFE-Weibull0.1881.1650.1520.353Weibull0.2841.7860.1930.276MOW0.2091.2060.1710.329EW0.2321.4990.1860.303As we see, the results (Tables4 and 5) show that the NFE-Weibull distribution has smaller values of the analytical measures and the maximum p value reveals that the proposed model provides better fit than the other considered competitors. Hence, the proposed model can be used as a best candidate model for modeling the COVID-19 total death data of the Asian countries. In support of the results provided in Tables 4 and 5, the estimated cdfs of the fitted distributions are plotted in Figure 7, whereas the Kaplan-Meier survival plots of the proposed and other fitted distributions are presented in Figure 8. From Figures 7 and 8, it is clear that the proposed model fit the estimated cdf and survival function very closely than the other competitors.Figure 7
Plots of the estimated cdfs of the NFE-Weibull and other competitive distributions for the COVID-19 total death data of the Asian countries.Figure 8
The Kaplan-Meier survival plots of the NFE-Weibull and other competitive distributions for the COVID-19 total death data of the Asian countries.


## 8. Concluding Remarks

The COVID-19 is one among the most deadly viruses that has greatly affected daily life affairs. The government and a number of other organizations should be interested to provide bases for comparison and to provide a better description of the data under consideration to get reliable estimates of the parameters of interest. In this article, a brief comparison of the COVID-19 events such as total cases, total deaths, total recovered, and active cases of the Asian countries are provided. Such clear cut comparison should be helpful to facilitate the COVID-19 affected peoples. Furthermore, a new class of statistical models is introduced. Some mathematical properties of the proposed class are derived. The maximum likelihood estimators of the model parameters are obtained. Finally, a special submodel of the proposed class called a new flexible extended Weibull distribution is studied in detail. The flexibility provided by the proposed model could be very useful in adequately describing the total death data in the Asian countries due to the COVID-19. We observed that the proposed model may provide a close fit to the COVID-19 total death data.



---

*Source: 4296806-2020-06-28.xml*