Parametric Modeling Approach to Covid-19 Pandemic Data

Nofiu Idowu Badmus^1*, Olanrewaju Faweya², Sikiru Ajibade Ige³

¹Department of Statistics, University of Lagos, Akoka, Nigeria.
²Department of Statistics, Ekiti State University, Ado-Ekiti, Ekiti State, Nigeria.
³Department of Mathematics, Yaba College of Technology, Lagos, Nigeria.

DOI: 10.4236/ojs.2023.131005                

Abstract

The problem of skewness is common among clinical trials and survival data, which has been the research focus derivation and proposition of different flexible distributions. Thus, a new distribution called Extended Rayleigh Lomax distribution is constructed from Rayleigh Lomax distribution to capture the excessiveness of some survival data. We derive the new distribution by using beta logit function proposed by Jones (2004). Some statistical properties of the distribution such as density, cumulative density, reliability rate, hazard rate, reverse hazard rate, moment generating and likelihood functions; skewness, kurtosis and coefficient of variation are obtained. We also performed the expected estimation of model parameters by maximum likelihood; goodness of fit and model selection criteria, including Anderson Darling, CramerVon Misses, Kolmogorov Smirnov (KS), Akaike Information, Bayesian Information, and Consistent Akaike Information Criterion is employed to select the better distribution from those models considered in the work. The results from the statistics criteria show that the intended distribution performs well and has a good representation of the States in Nigeria’s Covid-19 death cases data than other competing models.

Keywords

Anderson Darling, Cramer-von Mises, Covid-19, Kolmogorov Smirnov, Link Function, Survival Analysis

1. Introduction

In survival analysis, problems are encountered in the analysis of clinical data because distributions proposed are not flexible enough to follow the movement of the data to give accurate results. In light of this, there is a need to develop a more flexible parametric model using Covid-19 data for example. In recent times, there was the outbreak of the third wave of Covid-19 pandemic called Delta Variant after the second wave generating a global outcry. Many researches/works have been done by several researchers since the breakup of the pandemic in December 2019 from various fields, such as Medicine, Statistics, Economics, etc., with different ideas, models, methods and approaches in their respective works. These include Badmus et al. (2020), Dey et al. (2020), WHO, (2020), Yoo, (2020) amongst others [1] [2] [3] [4] . Moharraza et al. [5] in their work developed a simple model to assess the spread of global Covid-19 pandemic and tested the validity of their model using the cases of the coronavirus disease compiled by the WHO. Meanwhile, the results from the analysis supported the validity of their developed model.

Zheng and Bobasera [6] compared and contrast the Spanish flu and Covid-19 pandemic using an approach called a “powerful approach” to quantify the degree of mean distance between them and reported according to their results that vaccines from the progress in science and technology have reduced the death probability in some countries. E.g., about 0.0001 United Kingdom about 0.001 in Italy, the United States, Canada, and San Marino considered in their work

The main goal of the extensions of the proposed model is to enhance the robustness and flexibility of the classical model. Also, motivating factors, including the model, are a family of beta generalized models in such that when equating one or more parameter(s) to one (1), the model becomes the baseline distribution. In addition, the model has capability to increase and decrease shapes due to heavily skewed and heavy tail probability distribution that plays an important role in modeling skewed data like clinical data.

Most clinical data are always skewed, thus a new distribution is constructed and generated from a parent distribution called Rayleigh Lomax (RL) distribution by Kawsar et al. [7] is generated using beta link function introduced by Jones [8] . This is expected to have different shapes for the survival and hazard rate functions. More parameters are added to the parent distribution, and the flexibility and capability of the distribution to model real life data are established. There is still room for further study from the topic such as “Log-Beta Raleigh Lomax model, Modified Raleigh Lomax distribution and mixture of Rayleigh, Exponentiated Rayleigh Lomax, Lehmann Type II Rayleigh Lomax, Lomax and exponential distribution in modeling real life phenomenon”.

The paper is arranged in the following order. Section 2 contains material and methods which involve the properties of the proposed Extended Rayleigh Lomax distribution. In Section 3, we obtain the moments and generating function. Section 4 includes an estimation of parameters, analysis of secondary data obtained from Covid-19 weekly report in Nigeria, results and discussion. Section 5 includes based on the findings.

2. Material and Methods

There are several methods in literature which have been used by many researchers. In this study, we consider beta logit function introduced by Jones [8] , which can jointly convolute two or more distributions.

2.1. Properties of Extended Rayleigh Lomax (ERL) Distribution

Density Function

The density function of the above distribution is obtained using the beta link function given as:

$g\left(z,p,q\right)=\frac{c\left(z\right){\left[C\left(z\right)\right]}^{p-1}{\left[1-C\left(z\right)\right]}^{q-1}}{B\left(p,q\right)}$ (1)

$z,p,q,\theta ,\lambda ,\beta \simeq ERL\left(z,p,q,\theta ,\lambda ,\beta >0\right)$

$C\left(z\right)=1-{\text{e}}^{-\frac{\beta }{2}{\left(\frac{\theta }{\theta +z}\right)}^{-2\lambda }}$

and

$c\left(z\right)=\frac{\beta \lambda }{\theta }{\left(\frac{\theta }{\theta +z}\right)}^{-2\lambda +1}{\text{e}}^{-\frac{\beta }{2}{\left(\frac{\theta }{\theta +z}\right)}^{-2\lambda }}$ $x\ge -\theta$ and $\theta ,\lambda ,\beta >0$

where, C(z) and c(z) are cdf and pdf of the parent distribution respectively, a and b are additional shape parameters to the parent distribution.

$g\left(z,p,q,\theta ,\lambda ,\beta \right)=\frac{1}{B\left(p,q\right)}{\left[C\left(z\right)\right]}^{p-1}{\left(1-C\left(z\right)\right)}^{q-1}c\left(z\right)<z,p,q,\theta ,\lambda ,\beta >0$ (2)

where, θ, λ and β are initial scale and shape parameters, while a and b are new parameters called shape, introduced to the distribution. Then (2) becomes extended Rayleigh Lomax (ERL) Distribution and it plots are shown in Figure 1.

Some new and existing distributions from the extended Rayleigh Lomax distribution are mentioned below:

1) When $p=1$ in (2), we have the Lehmann Type II Rayleigh Lomax Distribution (New)

2) If $q=1$ in (2), we get the pdf of exponentiated Rayleigh Lomax distribution (New)

3) When $\beta =\text{1}$ in (2), this consists of beta Lomax distribution (New)

4) When $\frac{\beta }{\theta }=\beta$ , $\lambda =1$ and $\lambda =1$ in (2), it then becomes beta Rayleigh distribution (New)

5) If $p=q=1$ in (2), it yields Rayleigh Lomax distribution which is the parent distribution [7] .

6) When $q=\beta =1$ in (2), this consists exponential Lomax distribution by El-Bassiouny et al. [9] .

7) If $p=q=\lambda =\theta =1$ in (2), it yields Rayleigh Lomax distribution which is the parent distribution [10] .

2.2. Distribution Function of BRL Distribution

The associative distribution function cdf in (2) is given as

Figure 1. Depicts the: (1) pdf, plot of density function; (b) cdf, plot of distribution function; (c) survival, plot of survival rate function; and (d) hazard function in plots of ERL distribution, plot of hazard rate function.

$\begin{array}{c}P\left(Z\le z\right)={\displaystyle {\int }_0^{z}f\left(c\right)\text{d}c}=G\left(z,p,q,\theta ,\lambda ,\beta \right)\\ ={\displaystyle {\int }_0^z\frac{1}{B\left(p,q\right)}{\left[C\left(z\right)\right]}^{p-1}{\left(1-C\left(z\right)\right)}^{q-1}c\left(z\right)\text{d}c}\end{array}$ (3)

We set

$k\left(x\right)=1-{\text{e}}^{-\frac{\beta }{2}{\left(\frac{\theta }{\theta +x}\right)}^{-2\lambda }}$ (4)

$\frac{\text{d}c}{\text{d}z}=\frac{\beta \lambda }{\theta }{\left(\frac{\theta }{\theta +z}\right)}^{-2\lambda +1}{\text{e}}^{-\frac{\beta }{2}{\left(\frac{\theta }{\theta +z}\right)}^{-2\lambda }}\left[1-{\text{e}}^{-\frac{\beta }{2}{\left(\frac{\theta }{\theta +z}\right)}^{-2\lambda }}\right]$

$\frac{\text{d}c}{\text{d}z}=\frac{\beta \lambda }{\theta }{\left(\frac{\theta }{\theta +z}\right)}^{-2\lambda +1}{\text{e}}^{-\frac{\beta }{2}{\left(\frac{\theta }{\theta +z}\right)}^{-2\lambda }}-{\left[{\text{e}}^{-\frac{\beta }{2}{\left(\frac{\theta }{\theta +z}\right)}^{-2\lambda }}\right]}^2$

Also

$\text{d}z=\frac{\text{d}c}{\frac{\beta \lambda }{\theta }{\left(\frac{\theta }{\theta +z}\right)}^{-2\lambda +1}{\text{e}}^{-\frac{\beta }{2}{\left(\frac{\theta }{\theta +z}\right)}^{-2\lambda }}\left[1-{\text{e}}^{-\frac{\beta }{2}{\left(\frac{\theta }{\theta +z}\right)}^{-2\lambda }}\right]}$

Putting dx in Equation (2), we realize:

$g\left(z,p,q,\theta ,\lambda ,\beta \right)=\frac{1}{B\left(p,q\right)}{\displaystyle {\int }_0^{\infty }{\left[1-{\text{e}}^{-\frac{\beta }{2}{\left(\frac{\theta }{\theta +z}\right)}^{-2\lambda }}\right]}^{p-1}{\left[{\text{e}}^{-\frac{\beta }{2}{\left(\frac{\theta }{\theta +z}\right)}^{-2\lambda }}\right]}^{q-1}\text{d}c}$ (5)

and c in Equation (4) becomes

$c\left(z\right)=\frac{\text{d}c\left(z\right)}{\text{d}z}=\frac{\beta \lambda }{\theta }{\left(\frac{\theta }{\theta +z}\right)}^{-2\lambda +1}{\text{e}}^{-\frac{\beta }{2}{\left(\frac{\theta }{\theta +z}\right)}^{-2\lambda }}$

Equation (5) can be expressed as

$g\left(z,p,q,\theta ,\lambda ,\beta \right)=\frac{1}{B\left(p,q\right)}{C}^{a-1}{\left(1-C\right)}^{b-1}\frac{\text{d}c}{\text{d}z}$ (6)

Now, putting (2) in (6), we get

$\begin{array}{c}G\left(z,p,q,\theta ,\lambda ,\beta \right)=P\left(Z\le z\right)={\displaystyle {\int }_0^z\frac{1}{B\left(p,q\right)}{C}^{p-1}{\left(1-C\right)}^{q-1}\frac{\text{d}c}{\text{d}z}}\\ ={\displaystyle {\int }_0^z\frac{1}{B\left(p,q\right)}{C}^{p-1}{\left(1-C\right)}^{q-1}\frac{\text{d}c}{\text{d}z}}\end{array}$

where, $B\left(z,p,q\right)={\displaystyle {\int }_0^z{C}^{p-1}{\left(1-C\right)}^{q-1}\text{d}c}$ and it is called the incomplete beta function.

$G\left(z,p,q,\theta ,\lambda ,\beta \right)={\displaystyle {\int }_0^z\frac{C^{p-1}{\left(1-C\right)}^{q-1}\text{d}c}{B\left(p,q\right)}}=\frac{B\left(c,p,q\right)}{B\left(p,q\right)}$ (7)

(7) becomes the cumulative distribution function of ERL distribution.

2.3. The Survival Rate/Reliability Function

The reliability function of BRL distribution is given by

$Sur\left(z,p,q,\theta ,\lambda ,\beta \right)=1-G\left(z,p,q,\theta ,\lambda ,\beta \right)=1-{\displaystyle {\int }_0^{z}f\left(c\right)\text{d}c}$

$1-{\displaystyle {\int }_0^z\frac{C^{p-1}{\left(1-C\right)}^{q-1}\text{d}c}{B\left(p,q\right)}}=1-\frac{B(c,p,q)}{B\left(p,q\right)}$

$Sur\left(z,p,q,\theta ,\lambda ,\beta \right)=\frac{B\left[\left(p,q\right)-\left(z,p,q\right)\right]}{B\left(p,q\right)}$ (8)

2.4. The Hazard Rate Function

$\begin{array}{c}haz\left(z,p,q,\theta ,\lambda ,\beta \right)=\frac{g\left(z,p,q,\theta ,\lambda ,\beta \right)}{1-G\left(z,p,q,\theta ,\lambda ,\beta \right)}=\frac{\frac{C^{p-1}{\left(1-C\right)}^{q-1}c}{B\left(p,q\right)}}{\frac{B\left[\left(p,q\right)-\left(z,p,q\right)\right]}{B\left(p,q\right)}}\\ =\frac{C^{p-1}{\left(1-C\right)}^{q-1}c}{B\left[\left(p,q\right)-\left(z,p,q\right)\right]}\end{array}$ (9)

2.5. The Reversed Hazard Rate Function

$\begin{array}{c}Rhaz\left(z,p,q,\theta ,\lambda ,\beta \right)=\frac{g\left(z,p,q,\theta ,\lambda ,\beta \right)}{G\left(z,p,q,\theta ,\lambda ,\beta \right)}=\frac{\frac{C^{p-1}{\left(1-C\right)}^{q-1}c}{B\left(c,p,q\right)}}{B\left(p,q\right)}\\ =\frac{B\left(p,q\right)C^{p-1}{\left(1-C\right)}^{q-1}c}{B\left(c,p,q\right)}\end{array}$ (10)

2.6. Testing the Trueness of the PDF of ERL Distribution

The ERL distribution is a probability density function with the use of:

${\displaystyle {\int }_0^{\infty }{g}_{ERL}\left(z\right)\text{d}z}=1$ (11)

Jones (2004) in his generalized beta distribution of first kind is given by:

$g_z\left(z,p,q,u\right)={\left[B\left(p,q\right)\right]}^{-u}{\left[C\left(y\right)\right]}^{pu-1}{\left[1-C{\left(y\right)}^u\right]}^{q-1}K\left(y\right)0<y<1$

where $p,q$ and $u>0$ , therefore differentiating (p) above, we obtain

$g_{ERLD}\left(z\right)={\left[B\left(p,q\right)\right]}^{-u}{G}^{pu-1}{\left(1-G^u\right)}^{q-1}\frac{\text{d}G}{\text{d}y}$

${\displaystyle {\int }_{-\infty }^{\infty }{g}_{ERLD}\left(z\right)}={\displaystyle {\int }_{-\infty }^{\infty }\frac{u}{B\left(p,q\right)}{G}^{pu-1}{\left(1-G^u\right)}^{q-1}\text{d}G}$

Putting $M=G^u$ , then differentiating M with respect to G

$\frac{\text{d}M}{\text{d}G}=U{G}^{u-1}$

$\text{d}G=\frac{\text{d}M}{U{G}^{u-1}}$

$G=M^{\frac{1}{u}}$

$\begin{array}{c}{\displaystyle {\int }_{-\infty }^{\infty }{g}_{ERLD}\left(z\right)\text{d}z}={\displaystyle {\int }_0^1\frac{u}{B\left(p,q\right)}{\left(M^{\frac{1}{u}}\right)}^{pu-1}{\left(1-M\right)}^{q-1}\frac{\text{d}M}{U{G}^{u-1}}}\\ ={\left[B\left(p,q\right)\right]}^{-1}{\displaystyle {\int }_0^1\frac{M^{p-\frac{1}{u}}{\left(1-M\right)}^{q-1}}{M^{1-\frac{1}{u}}}\text{d}M}\\ ={\left[B\left(p,q\right)\right]}^{-1}{\displaystyle {\int }_0^1{M}^{p-1}{\left(1-M\right)}^{q-1}\text{d}M}\end{array}$

${\displaystyle {\int }_0^1{M}^{p-1}{\left(1-M\right)}^{q-1}\text{d}M}=B\left(p,q\right)$

Therefore, $g_{ERLD}\left(z\right)=\frac{B\left(p,q\right)}{B\left(p,q\right)}=1$ .

Hence, the $g_{ERL}$ distribution has a true continuous probability density function.

3. Moments and Generating Function

In this section, we derive and obtain mgf of the distribution $m\left(t\right)=E\left({\text{e}}^{tz}\right)$ and the general rth moment of a beta generated distribution defined by Hosking [11]

${\mu }_r^1=B{\left(p,q\right)}^{-1}{\displaystyle {\int }_0^1{\left[F^{-1}\left(z\right)\right]}^r{z}^{p-1}{\left[1-z\right]}^{q-1}\text{d}z}$ (12)

Cordeiro et al. [12] also discussed another mgf for generated beta distribution.

$m\left(t\right)=B{\left(p,q\right)}^{-1}{\displaystyle {\sum }_{j=0}^{\infty }{\left(-1\right)}^j\left(\begin{array}{c}q-1\\ j\end{array}\right)e\left(a,r;pj-1\right)}$ (13)

where,

$e\left(a,r\right)={\displaystyle {\int }_{-\infty }^{\infty }{\text{e}}^{tz}{\left[F\left(z\right)\right]}^{m}f\left(z\right)\text{d}z}$

then,

$M_z\left(t\right)=B{\left(p,p\right)}^{-1}{\displaystyle {\sum }_{j=0}^{\infty }{\left(-1\right)}^j}\left(\begin{array}{c}q-1\\ j\end{array}\right){\displaystyle {\int }_{-\infty }^{\infty }{\text{e}}^{tz}{\left[F\left(z\right)\right]}^{p\left(j+1\right)-1}f\left(z\right)\text{d}z}$ (14)

Putting pdf and cdf of the Extended Rayleigh Lomax distribution into Equation (14), we get

$M_{ERLD\left(z\right)}\left(t\right)=B{\left(p,q\right)}^{-1}{\displaystyle {\sum }_{j=0}^n{\left(-1\right)}^j}\left(\begin{array}{c}q-1\\ j\end{array}\right){\displaystyle {\int }_{-\infty }^{\infty }{\text{e}}^{tz}{\left[F\left(z\right)\right]}^{p\left(j+1\right)-1}f\left(z\right)}$ (15)

If $p=q=1$ in Equation (14) that becomes the moment generating function of the baseline distribution.

Hence, the rth moment of the ERL distribution is obtained, since the moment generating function of the parent distribution is given by

$\begin{array}{c}{M}_z\left(t\right)={\displaystyle {\sum }_{j=0}^{\infty }\frac{t^j}{j!}}{\displaystyle {\int }_0^{\infty }{z}^{j}f\left(z,\beta ,\lambda ,\theta \right)\text{d}z}{\displaystyle {\sum }_{c=0}^{\infty }\frac{t^k}{k!}}\\ \cdot {\displaystyle {\sum }_{j=0}^{\infty }\left(q{K}_j\right){\theta }^j{\left(\frac{2}{\beta }\right)}^{\frac{1}{2\lambda }}{\left(-\theta \right)}^{q-j}r\left(\frac{j}{2\lambda }+1\right)}\end{array}$ (16)

Equation (16) can be re written as

$\begin{array}{c}{M}_{ERLD\left(z\right)}\left(t\right)=B{\left(p,q\right)}^{-1}{\displaystyle {\sum }_{j=0}^n{\left(-1\right)}^j}\left(\begin{array}{c}q-1\\ j\end{array}\right){\displaystyle {\int }_{-\infty }^{\infty }{\text{e}}^{tz}{\left[F\left(z\right)\right]}^{p\left(j+1\right)-1}{\displaystyle {\sum }_{i=0}^{\infty }\frac{t^i}{i!}}}\\ \cdot {\displaystyle {\sum }_{h=0}^{\infty }\left(i{K}_h\right){\theta }^h{\left(\frac{2}{\beta }\right)}^{\frac{h}{2\lambda }}{\left(-\theta \right)}^{i-h}r\left(\frac{h}{2\lambda }+1\right)}\\ =B{\left(p,q\right)}^{-1}{\displaystyle {\sum }_{i=0}^{\infty }{\displaystyle {\sum }_{j=0}^{\infty }{\displaystyle {\sum }_{h=0}^{\infty }{\left(-1\right)}^j\left(\begin{array}{c}q-1\\ i\end{array}\right)\frac{t^i}{i!}\left(i{K}_h\right){\theta }^h{\left(\frac{2}{\beta }\right)}^{\frac{h}{2\lambda }}}}}\\ \cdot {\left(-\theta \right)}^{i-h}r\left(\frac{h}{2\lambda }+1\right)\cdot {\left[F\left(z\right)\right]}^{p\left(j+1\right)-1}\end{array}$ (17)

and the rth moment of ERL distribution is obtained from Equation (17)

$\begin{array}{c}{\mu }_{ERLD\left(r\right)}^i=E\left(z^r\right)=B{\left(p,q\right)}^{-1}{\displaystyle {\sum }_{j=0}^{\infty }{\displaystyle {\sum }_{h=0}^{\infty }{\left(-1\right)}^j\left(\begin{array}{c}q-1\\ i\end{array}\right)\cdot {\left[F\left(z\right)\right]}^{p\left(j+1\right)-1}}}\\ \cdot \frac{t^i}{r!}\left(\begin{array}{c}r\\ h\end{array}\right){\theta }^h{\left(\frac{2}{\beta }\right)}^{\frac{h}{2\lambda }}{\left(-\theta \right)}^{r-h}r\left(\frac{h}{2\lambda }+1\right)\end{array}$ (18)

Letting $p=q=1$ in (18) gives the rth moment of the baseline distribution by Kawsar et al. [7]

${\mu }_r^i=E\left(z^r\right)={\displaystyle {\sum }_{h=0}^{\infty }\left(\begin{array}{c}r\\ h\end{array}\right){\theta }^h{\left(\frac{2}{\beta }\right)}^{\frac{h}{2\lambda }}{\left(-\theta \right)}^{r-h}r\left(\frac{h}{2\lambda }+1\right)}$

Other measures such as the Skewness (SK_ERLD) (p, q, β, λ, θ) and Kurtosis (KT_ERLD) (p, q, β, λ, θ) are also obtained below:

The rth moment of the ERL distribution is written as:

${\mu }_{ERLD\left(r\right)}^i={\displaystyle {\int }_0^{\infty }{z}^r{F}_{ERLD}\left(z\right)\text{d}z}$

That is,

${\mu }_{ERLD\left(r\right)}^i={\displaystyle {\int }_0^{\infty }{z}^r}\left({\left[B\left(p,q\right)\right]}^{-1}{\left(C\left(z\right)\right)}^{p-1}{\left(1-C\left(z\right)\right)}^{q-1}\text{d}c(\; z\; )\right)$

where, $C\left(z\right)={\left(1-w\left(z\right)\right)}^{-u}$

I.e. $W\left(z\right)={\text{e}}^{-\frac{\beta }{2}\left(\frac{\theta }{\theta +z}\right)}$ and $u=2\lambda$

therefore,

$\begin{array}{c}{\mu }_{ERLD\left(r\right)}^i=\frac{\left(r{K}_h\right){\theta }^h{\left(\frac{2}{\beta }\right)}^{\frac{h}{u}}{\left(-\theta \right)}^{r-h}r\left(\frac{h}{u}+1\right)}{B\left(a,b\right)}\\ \cdot {\displaystyle {\sum }_{j=0}^{\infty }{\displaystyle {\sum }_{h=0}^{\infty }{\left(-1\right)}^j}}\left(\begin{array}{c}q-1\\ i\end{array}\right)\left({\left[{\left(1-w\left(z\right)\right)}^{-u}\right]}^{p\left(j+1\right)-1}\right)\\ =X\left(\left(\begin{array}{c}r\\ h\end{array}\right){\theta }^h{\left(\frac{2}{\beta }\right)}^{\frac{h}{u}}{\left(-\theta \right)}^{r-h}r\left(\frac{h}{u}+1\right)\right)\end{array}$ (19)

where,

$X=\frac{{\displaystyle {\sum }_{j=0}^{\infty }{\displaystyle {\sum }_{h=0}^{\infty }{\left(-1\right)}^j}}\left(\begin{array}{c}q-1\\ i\end{array}\right)\left({\left[{\left(1-w\left(z\right)\right)}^{-u}\right]}^{p\left(j+1\right)-1}\right)}{B\left(p,q\right)}$

At the same time, the first four central moments ${\mu }_r^i=1,2,3,4$ are obtained through (17) as:

Furthermore, the mean and second to fourth moments of the ERL distribution are given as follows:

$\mu ={\mu }_1^{|}$ , ${\mu }_2={\mu }_2^{|}-{\mu }^2$ , ${\mu }_3={\mu }_3^{|}-3\mu {\mu }_2^{|}+2{\mu }^3$ , and ${\mu }_4={\mu }_4^{|}-4\mu {\mu }_3^{|}+6{\mu }^2{\mu }_2^{|}-3{\mu }^4$

${\mu }_1^{|}=X\left(\left(\begin{array}{c}1\\ h\end{array}\right){\theta }^h{\left(\frac{2}{\beta }\right)}^{\frac{h}{u}}{\left(-\theta \right)}^{1-h}r\left(\frac{h}{u}+1\right)\right)$ (20)

${\mu }_2^{|}=X\left(\left(\begin{array}{c}2\\ h\end{array}\right){\theta }^h{\left(\frac{2}{\beta }\right)}^{\frac{h}{u}}{\left(-\theta \right)}^{2-h}r\left(\frac{h}{u}+1\right)\right)$ (21)

${\mu }_3^{|}=X\left(\left(\begin{array}{c}3\\ h\end{array}\right){\theta }^h{\left(\frac{2}{\beta }\right)}^{\frac{h}{u}}{\left(-\theta \right)}^{3-h}r\left(\frac{h}{u}+1\right)\right)$ (22)

${\mu }_4^{|}=X\left(\left(\begin{array}{c}4\\ h\end{array}\right){\theta }^h{\left(\frac{2}{\beta }\right)}^{\frac{h}{u}}{\left(-\theta \right)}^{4-h}r\left(\frac{h}{u}+1\right)\right)$ (23)

Other measures such as skewness, kurtosis and coefficient of variation of the ERL distribution are given below.

3.1. Skewness of the ERL Distribution

The skewness is a means of measuring non symmetry of the distribution. The skewness is given by:

$S{K}_{ERLD}=\frac{{\mu }_3}{{\mu }_2^{1.5}}$ (24)

3.2. Kurtosis of the ERL Distribution

The kurtosis is another measure that measures the peak of the distribution. The kurtosis of the BRL distribution is given as:

$K{T}_{ERLD}=\frac{{\mu }_4}{{\mu }_2^2}-3$ (25)

3.3. Coefficient of Variation of the ERL Distribution

This is also a measure of variability of a probability distribution. The CV of the ERL distribution is given as:

$C{V}_{ERLD}=\frac{\sqrt{{\mu }_2}}{\mu }$ (26)

4. Estimation of Parameter

We made attempt to derive the maximum likelihood estimates (MLEs) of the ERL distribution parameters including: θ, λ, β, a and b which are scale and shape parameters. According to Cordeiro et al. [12] , the log likelihood function is given as:

$\begin{array}{c}L\left(\phi \right)=n\log \left[B\left(p,q\right)\right]+{\displaystyle {\sum }_{i=1}^n\log \left[C\left(z,\delta \right)\right]}\\ +\left(p-1\right){\displaystyle {\sum }_{i=1}^n\log }\left[1+C\left(z,\delta \right)\right]+\left(p-1\right)\left[C\left(z,\delta \right)\right]\end{array}$ (27)

$\phi =\left(p,q,k,\delta \right)$ and $\delta =\left(\theta ,\lambda ,\beta \right)$ are vectors

If c = 1, it becomes Equation (27) which leads to

$\begin{array}{c}L\left(\phi \right)=const-n\log \left[B\left(p,q\right)\right]+{\displaystyle {\sum }_{i=1}^n\log }\left[C\left(z,\delta \right)\right]\\ +\left(p-1\right){\displaystyle {\sum }_{i=1}^n\log }\left[1+C\left(z,\delta \right)\right]+\left(p-1\right)\left[C\left(z,\delta \right)\right]\end{array}$ (28)

$c\left(z,\delta \right)$ and $C\left(z,\delta \right)$ have been stated at the beginning.

The log likelihood function of ERL distribution is given as:

$\begin{array}{c}{L}_{BRLD}\left(\phi \right)=-n\log \left[B\left(p,q\right)\right]+{\displaystyle {\sum }_{i=1}^n\left[c\left(z\right)\right]}+\left(p-1\right){\displaystyle {\sum }_{j=1}^n\log \left[C\left(z\right)\right]}\\ +\left(q-1\right){\displaystyle {\sum }_{i=1}^n\log \left[1-C\left(z\right)\right]}\end{array}$ (29)

Taking the differentiation in respect to p, q, θ, λ and β give the following:

$\frac{\partial L\left(\phi \right)}{\partial p}=-n\frac{r\text{|}\left(p\right)}{r\left(p\right)}+n\frac{r\text{|}\left(p+q\right)}{r\left(p+q\right)}+{\displaystyle {\sum }_{z=1}^n\log \left(1-C\left(Z,\delta \right)\right)}$ (30)

$\frac{\partial L\left(\phi \right)}{\partial q}=-n\frac{r\text{|}\left(q\right)}{r\left(q\right)}+n\frac{r\text{|}\left(p+q\right)}{r\left(p+q\right)}+{\displaystyle {\sum }_{z=1}^n\log \left(C\left(Z,\delta \right)\right)}$ (31)

$\begin{array}{c}\frac{\partial L\left(\phi \right)}{\partial \theta }={\displaystyle {\sum }_{j=1}^n\left[\frac{\frac{\partial \left[C\left(Z,\delta \right)\right]}{\partial \theta }}{C\left(Z,\delta \right)}\right]}+\left(p-1\right){\displaystyle {\sum }_{z=1}^n\left[\frac{\frac{\partial \left[1-C\left(Z,\delta \right)\right]}{\partial \theta }}{1-C\left(Z,\delta \right)}\right]}\\ +\left(q-1\right){\displaystyle {\sum }_{z=1}^n\left[\frac{\frac{\partial \left[C\left(Z,\delta \right)\right]}{\partial \theta }}{\left(Z,\delta \right)}\right]}\end{array}$ (32)

$\begin{array}{c}\frac{\partial L\left(\phi \right)}{\partial \lambda }={\displaystyle {\sum }_{j=1}^n\left[\frac{\frac{\partial \left[C\left(Z,\delta \right)\right]}{\partial \lambda }}{C\left(Z,\delta \right)}\right]}+\left(p-1\right){\displaystyle {\sum }_{z=1}^n\left[\frac{\frac{\partial \left[1-C\left(Z,\delta \right)\right]}{\partial \lambda }}{1-C\left(Z,\delta \right)}\right]}\\ +\left(q-1\right){\displaystyle {\sum }_{z=1}^n\left[\frac{\frac{\partial \left[C\left(Z,\delta \right)\right]}{\partial \lambda }}{\left(Z,\delta \right)}\right]}\end{array}$ (33)

$\begin{array}{c}\frac{\partial L\left(\phi \right)}{\partial \beta }={\displaystyle {\sum }_{j=1}^n\left[\frac{\frac{\partial \left[C\left(Z,\delta \right)\right]}{\partial \beta }}{C\left(Z,\delta \right)}\right]}+\left(p-1\right){\displaystyle {\sum }_{z=1}^n\left[\frac{\frac{\partial \left[1-C\left(Z,\delta \right)\right]}{\partial \beta }}{1-C\left(Z,\delta \right)}\right]}\\ +\left(q-1\right){\displaystyle {\sum }_{z=1}^n\left[\frac{\frac{\partial \left[C\left(Z,\delta \right)\right]}{\partial \beta }}{\left(Z,\delta \right)}\right]}\end{array}$ (34)

4.1. Analysis of Data

The data used for the analysis is a secondary data obtained from Covid-19 situation weekly epidemiological report 39; 5^th-11^th July, 2021 (NCDC website state the website) [13] : Thirty-six (36) States including Federal Capital Territory (FCT) with reported laboratory-confirmed Covid-19 cases, recoveries, deaths, samples tested and active cases (37 data points); and was accessed on Thursday 22^nd July, 2021 put date accesses at reference not here. Only the death cases from all states of the federation are used for the analysis.

4.2. Result and Discussion

The summary of goodness of fit statistics is used to check for normality of the data; skewness, kurtosis, Anderson Darling (AD), Kolmogorov Smirnov (KS) and Cramer-Von-Mises (CVM) shown in Table 1 with their values clearly indicate that the data does not follow normal distribution since p-values less than 5%, skewness greater than 0 (zero) and kurtosis also greater than 3 [14] [15] . While, graphs from Figure 2 show the nature of the data, the scatter, theoretical quantiles, boxplot, histogram, density and empirical cumulative distribution

Figure 2. The scatter, theoretical quantiles, boxplot, histogram, density and distribution plot.

function (ecdf) plot show the data is skewed. For instance, non-linearity by scatter and quantiles plots, outliers by boxplot and skewness by histogram and density plots. The minimum and maximum values in the data set are inclusive.

The results obtained in Table 2 are based on parameter estimates by method of maximum likelihood estimation (MLEs). The standard error values are in

Table 1. Summary of goodness of fit statistics of death cases data.

Table 2. Contains the MLE, standard error (in parenthesis) and model selection criteria.

bracket for all the models. The model ERLD is compared with other six models: ExpLD, LRLD, BRD, RLD, ExpRLD and BLD. Also, model selection criterion is performed on all models considered in the study. From the results, ERLD has the smallest values in all as we can see bold and starred where AIC = 2607.878, CAIC = 2609.813, HQIC = 2610.718 and BIC = 2615.933, which indicates that it is a robust and flexible model.

5. Conclusion

Despite the level of Nigerian Covid-19 death cases data set, the ERL distribution follows the movement of the data and has better representation of the data than any of the other existing distributions. The proposed distribution, being flexible and versatile, can accommodate increasing, decreasing, bathtub and unimodal shape hazard function. It is therefore useful and effective in the analysis of clinical and survival data.

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1]	Badmus, N.I., Olanrewaju, F. and Adeniran, A.T. (2020) Modeling Covid-19 Pandemic Data with Beta Double Exponential Distribution. Asian Journal of Research and Infectious Diseases, 5, 66-79. https://doi.org/10.9734/ajrid/2020/v5i430181
[2]	Dey, S.K., Rahman, M.M., Siddiqi, U.R. and Howlader, A. (2020) Analyzing the Epidemiological Outbreak of COVID19: A Visual Exploratory Data Analysis Approach. Journal of Medical Virology, 92, 632-638. https://doi.org/10.1002/jmv.25743
[3]	World Health Organization (2020) Surveillance Case Definitions for Human Infection with Novel Coronavirus (nCoV).
[4]	Yoo, J.-H. (2020) The Fight against the 2019-nCoV Outbreak: An Arduous March Has Just Begun. Journal of Korean Medical Science, 35, 1598-6357. https://doi.org/10.3346/jkms.2020.35.e56
[5]	Moharraza, P., Richardson, G. and Richardson, G. (2021) A Mathematical Model for Spread of COVID-19 in the World. Journal of Applied Mathematics and Physics, 9, 1890-1895. https://doi.org/10.4236/jamp.2021.98122
[6]	Zheng, H. and Bonasera, A. (2022) Controlling the World Wide Chaotic Spreading of Covid-19 through Vaccinations. Journal of Modern Physics, 13, 1-15. https://doi.org/10.4236/jmp.2022.131001
[7]	Fatima, K., Jan, U. and Ahmad, S.P. (2018) Statistical Properties of Rayleigh Lomax Distribution with Applications in Survival Analysis. Journal of Data Scoence, 16, 531-548. https://doi.org/10.6339/JDS.201807_16(3).0005
[8]	Jones, M.C. (2004) Families of Distributions Arising from Distributions of Order Statistics. Test, 13, 1-43. https://doi.org/10.1007/BF02602999
[9]	El-Bassiouny, A.H., Abdoand, N.F and Shahen, H.S. (2015) Exponential Lomax Distribution. International Journal of Computer Applications, 121, 24-29. https://doi.org/10.5120/21602-4713
[10]	Siddiqui, M.M. (1962) Some Problems Connected with Rayleigh Distributions. Journal of Research of the National Bureau of Standards, 60D, 167-174. https://doi.org/10.6028/jres.066D.020
[11]	Hosking, J.R.M. (1990) L-Moments: Analysis and Estimation of Distributions Using Linear Combinations of Order Statistics. Journal of the Royal Statistical Society: Series B, 52, 105-124. https://doi.org/10.1111/j.2517-6161.1990.tb01775.x
[12]	Cordeiro, G.M. and de Castro, M. (2011) A New Family of Generalized Distributions. Journal of Statistical Computation and Simulation, 81, 883-898. https://doi.org/10.1080/00949650903530745
[13]	NCDC (2021) General Fact Sheet-Data as at 11th July, 2021. States with Reported Laboratory Confirmed Covid-19 Cases, Recoveries, Deaths Samples Tested and Active Cases.
[14]	Yap, B.W and Sim, C.H. (2011) Comparisons of Various Types of Normality Tests. Journal of Statistical Computation and Simulation, 81, 2141-2155. https://doi.org/10.1080/00949655.2010.520163
[15]	Shivam, M (2020) Normality Test with Python in Data Science. https://medium.com/@shivamrkom?source=post_page........5abbefc81fdo