key: cord-0774306-h0z0nyc1
authors: Makinde, Olusola Samuel; Olusola-Makinde, Olubukola Olayemi; Olamide, Emmanuel Idowu; Abiodun, Gbenga Jacob
title: Impact of COVID-19 prevalence and mode of transmission on mortality cases over WHO regions
date: 2020-10-15
journal: Health Inf Sci Syst
DOI: 10.1007/s13755-020-00127-3
sha: c9b464d9a30f1950e143c23aa8169a1fe5dde308
doc_id: 774306
cord_uid: h0z0nyc1

With the current outbreak of coronavirus disease 2019 (COVID-19), countries have been on rising preparedness to detect and isolate any imported and locally transmitted cases of the disease. It is observed that mode of transmission of the disease varies from one country to the other. Recent studies have shown that COVID-19 cases are not influenced by race and weather conditions. In this study, effect of modes of transmission of COVID-19 is considered with respect to prevalence and mortality counts in World Health Organisation (WHO) regions. Also, a negative binomial model is formulated for new death cases in all WHO regions as a function of confirmed cases, confirmed new cases, total deaths and modes of transmission, with the goal of identifying a model that predicts the total new death cases the best. Results from this study show that there is strong linear relationship among the COVID-19 confirmed cases, total new deaths and mode of transmission in all WHO regions. Findings highlight the significant roles of modes of transmission on total new death cases over WHO regions. Mode of transmission based on community transmission and clusters of cases significantly affects the number of new deaths in WHO regions. Vuong test shows that the formulated negative binomial model fits the data better than the null model.

The on-going pandemic; novel coronavirus disease 2019 , is spreading. All the continents in the world have been affected with daily records of new cases. The causative agent of COVID-19 is a positive-sense singlestranded RNA virus called severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) [6] .

Past experiences with disease outbreaks have shown low level of preparedness in Africa, for instance, the disparaging epidemic by Ebola virus from 2014 to 2016 in West Africa [23] . For COVID-19 outbreak, the continent is issuing different schemes to combat continuous spread of the virus; such schemes include (1) screening for the disease on arrivals at airports and some seaports; trainings on the detection of COVID-19 at point-ofentry has been carried out across Africa by designated authorities [24] . (2) Accompanied with the provision of resources, risk countries like Egypt have been educated on prevention and control of infection, community participation and management at health centres with respect to COVID-19 [9, 13, 19] . On 14th February, 2020, the first COVID-19 confirmed case was detected in Egypt in the whole of Africa, as at the date of this write-up (7th May 2020), Africa's total confirmed cases have risen to 52,821 with 2018 deaths [25] .

SARS-CoV-2 infection took inception from Asia, specifically, China in December, 2019 and has since risen to 605,270 total cases and more than 2205 new case on 6th May, 2020 [25] . For history, China Health Authority noted WHO on 31st December, 2019 about numerous cases of pneumonia-like disease now COVID-19 [14] . Even though, symptoms such as livedo or necrosis has been accrued to COVID-19 [4] . SARS-CoV-2 Phylogenetic studies revealed that SARS- CoV-2 genome  sequence are linked to bat-SL-CoVZC45 and bat-SL-CoVZXC21 which are coronaviruses cultured from bats  at eastern China, 2018 (88% identity) , meanwhile, the genetics of SARS-CoV-2 are distinct from MERS-CoV and SARS-CoV with 79% identity similarity [16] . As of 7th May, 2020, North America has reported a total of 1,378,143 COVID-19 cases, whereas, America has 91.66% of the continent's cases [25] . The analyzed data from April, 2020 in America revealed that the age group that is less than 18 years were least affected, with 2572 cases out of the total number of 149,082 cases [5] . The continent has since reported 82,650 total deaths [25] . Of all the continents, Europe is presently the worst hit of COVID-19 outbreak. It has the highest number of total confirmed cases (1,536,901), deaths (146,751) and new cases (more than 11,997) [25] . On 25th January, 2020, the continent discovered its index case in France with the presence of 790,748 active cases till date [21, 25] . Europe recorded as high as 40,800 new confirmed cases on 4th April, 2020. In order of increase, the European countries; France, Russia, UK, Italy and Spain have the highest number of COVID-19 cases in the continent. On 6th May, 2020, Russia recorded more than 88 new deaths out of about 120 new deaths in the continent [25] . Meanwhile, the continent, Oceania has least been affected by COVID-19 outbreak with no new death cases as at 6th May, 2020. The countries; New Caledonia and Papua New Guinea had total cases of 18 and 4 respectively, whereas all these cases have recovered [25] .

The use of mathematical models has been evident in decision-making around diverse non-pharmaceutical mediations in the confinement of the spread of COVID-19 globally. Modeling is very vital in studying the spread of COVID-19, with a particular model addressing a specific problem [20]. Zhang et al. [26] considered a descriptive study and modeling of transmission dynamics of COVID-19 outside Hubei province, China and established demographic trends of cases at the provincial level using a Bayesian procedure. Hirk et al. [8] formulated numerous situations using data on the number of COVID-19 cases from the hospitalized patients, those in intensive care unit, along with the number of deaths that tested positive in Iceland and Austria using both frequentist and Bayesian approaches.

In another development, Benvenuto et al. [3] predicted prevalence and incidence of the COVID-19 cases using an autoregressive integrated moving average (ARIMA) model on the Johns Hopkins epidemiological data. Khan et al. [11] conducted a case series study for the assessment of the impact of SARS-CoV-2 infection on adverse pregnancy outcomes. The study included both healthcare workers and women from the general population.

Liang [15] fitted nonlinear model to the growth rate of COVID-19. Makinde et al. [18] employed generalized estimating equation to analyze the daily COVID-19 mortality rates in African countries, and found that there was a positive weak linear relationship between the daily reported COVID-19 cases and the population of African countries. However, the magnitude of the observed association was particularly small. Also, there were significant monotone trends in the daily COVID-19 incidence and mortality counts of many countries in Africa.

In this paper, the impact of mode of transmission of COVID-19 on the incidence and mortality; as well as drivers of total new deaths over WHO regions are considered. The linear relationships among the total COVID-19 confirmed cases, total confirmed new cases, total deaths and total new deaths in WHO regions are presented. A negative binomial regression model is formulated for total new death cases accounting for the mode of transmission of COVID-19 over WHO regions.

The World Health Organisations (WHO) situation report 187 data, available at https ://www.who.int/emerg encie s/disea ses/novel -coron aviru s-2019/situa tion-repor ts, present the information on total confirmed cases, total confirmed new cases, total deaths, total new deaths and mode of transmission of COVID-19 in countries, regions and territories defined as WHO regions till July 26, 2020. Figure 1 presents the plot of COVID-19 total confirmed cases (in `000) and total deaths in some countries and territories of WHO regions. The figure shows the 30 countries with highest total confirmed cases and 30 countries with highest total mortality cases in the WHO regions. The United States of America recorded the highest COVID-19 total confirmed cases and total deaths. Figure 2 presents pie of pie chart of WHO regions' total number of COVID-19 deaths. It is observed from the figure that 12% of the WHO regions are yet to record COVID-19 deaths while 9% of the regions recorded 1-5 deaths. 43% of the 215 WHO regions experience more than 100 deaths as at July 26, 2020. Figure 3 presents frequency distribution of mode of transmission of COVID-19 pandemic in the WHO regions. It was observed that the mode of transmission in 36.74% of the WHO regions is by clusters of cases while 42.79% is through community transmission. Also, 14.88% of the WHO regions had their mode of transmission of COVID-19 to be sporadic cases, mode of transmission in 0.93% of the WHO regions is pending while 4.65% of the WHO regions are yet to record a COVID-19 case. This implies that the commonest mode of transmission in all countries, territories and regions defined by WHO is community transmission.

Suppose a random variable Y is distributed as negative binomial, the negative binomial distribution [7] of Y is defined by its probability mass function, defined as where µ is the mean incidence rate of y. Suppose a random variable Y follows the negative binomial distribution. Then its conditional expected value is E Y |X 1 , X 2 , . . . , X p = E p j=1 β j X j and the variance is var Y |X 1 , X 2 , . . . , X p = µ + θµ 2 , where θ is the overdispersion parameter and X 1 , X 2 , . . . , X p are predictor variables. Negative binomial regression is used to model over-dispersed count data. It expresses µ in terms of X 1 , X 2 , . . . , X p .

The total number of new deaths in each country, region or territory may be over-dispersed. That is, the mean of total number of new deaths may be much less than its variance. To investigate this, the mean and variance of total number of new deaths in all countries and territories are 25.8837 and 15,829.75 respectively. In this study, negative binomial (NB) regression model [7, 22] A stepwise model selection approach based on Akaike information criterion (AIC) [2] is used to exclude models with highest AIC values in the fitted NB model in (1) . This is executed with the variables log e (N TC ), N LC and N TD in order to identify the variables that have the least importance which were omitted from further analysis [17] .

The variable N TC is transformed to reduce the variability of data, especially in data that include outlying observations. The variance of N TC is 112,325,761,815. The variable N TC contains some outlying observations. This is shown in the first boxplot in Fig. 4 . The use of the natural log transformation is to reduce variance of N TC . The variance of natural logarithm of N TC is 8.48. This is shown in the second boxplot in Fig. 4 .

Linear relationship among total confirmed cases, total confirmed new cases, total deaths and total new deaths Table 1 presents the estimate of Kendall's rank correlation coefficients [10] among variables ( N TC , N LC , N TD and Y ND ), with p value of t test of significance of correlation estimates in parenthesis. The total confirmed cases and total confirmed new cases are significantly, positively and strongly related linearly. This implies that as the total number of confirmed cases increases, the total number of confirmed new cases increases. This is the same for all combinations of all the variables. Also, the total number of deaths recorded increases with increase in the total number of confirmed cases in either of the countries, territories and area under WHO regions.

Kruskal-Wallis test [12] is applied to test if modes of transmission (clusters of cases, community transmission, pending and sporadic cases) in all the countries are equivalent in terms of each of N TC , N LC , N TD and Y ND . The result of the Kruskal-Wallis test for the modes of transmission based on N TC , N LC , N TD and Y ND is shown in Table 2 in terms of p values of the test. The COVID-19 modes of transmission in all the WHO regions are significantly different at 5% level of significance in terms of N TC , N LC , N TD and Y ND .

The mean and variance of the total number of COVID-19 new deaths reported as at July 11, 2020 across countries, territories or areas are 29.163 and 19,596.35 respectively. Fig. 4 Boxplot of a the total number of incidence cases and b the natural logarithm total number of incidence cases in all WHO regions As a result, negative binomial models are formulated for total number of COVID-19 deaths. A negative binomial model was formulated for the total number of COVID-19 deaths as a function of total number of COVID-19 confirmed cases, total number of new confirmed cases, total number of deaths and mode of transmission. A stepwise count model selection was performed on Eq. (1). The AIC value for the model (Eq. 1) is 832.19. The model with the least AIC (AIC = 829.12) excludes total number of deaths and total number of new confirmed cases in the WHO regions. The coefficient of natural logarithms of N TC , as well as modes of transmission (community transmission and clusters of cases), which predicts Y ND is statistically significant at 5% level of significance. Table 3 presents the estimates of coefficients of negative binomial regression model for the total number of new deaths in the WHO regions. The total number of new deaths in the WHO regions increases by a factor of 2.863 for a 1-unit increase in log of total number of confirmed cases when other variables are held constant. The total number of new deaths in the WHO regions decreases by a factor of 0.0001 and 0.0002 for a 1-unit increase in factor of clusters of cases and community transmission respectively when other variables are held constant. It can be inferred from the table that the total number of new deaths increases with increase in the total number of confirmed cases.

The dispersion parameter ( θ = 0.6598 ) in Table 3 is finite, this implies that the mean is less than variance. As θ approaches 0, the distant is the mean from the variance [1] . Consequentially, negative binomial model will fit the data well. In addition, the Vuong test is used to check if a null model fits the total number of new deaths equivalently as the fitted NB model. The test rejects the null hypothesis at 0.05 level of significance (p value < 2.22e- 16) and suggests that NB model fits the data better than the null model. Figure 5 presents the comparison between observed data and fitted data based on negative binomial model for some countries, territories and areas in WHO regions. It can be observed from the figure that the model fits the data well.

This study provides insight into whether mode of transmission in a WHO region affects number of COVID-19 new death cases over a time period. In particular, a negative binomial regression model was formulated for predicting total new death cases in all WHO regions. The findings showed that COVID-19 total number of new deaths is significantly affected by the total number of confirmed cases and modes of transmission of the disease in the WHO regions. The significance of community transmission and clusters of cases among other modes of transmission on mortality was highlighted. This study is limited by dearth of information on population in all the WHO regions. 

A dynamical and zero-inflated negative binomial regression modelling of malaria incidence in Limpopo Province, South Africa

A new look at the statistical model identification

Application of the ARIMA model on the COVID-2019 epidemic dataset

García-Doval I. Classification of the cutaneous manifestations of COVID-19: a rapid prospective nationwide consensus study in Spain with 375 cases

CDC_AA_refVa l=https %3A%2F%2Fwww .cdc.gov%2Fcor onavi rus%2F201 9-ncov%2Fcas es-inus.html. Retrieved

Origin and evolution of pathogenic coronaviruses

Cambridge: Cambridge University Press

Investigating the dark figure of COVID-19 cases in Austria: borrowing from the deCODE genetics study in Iceland

Delays in global disease outbreak responses: lessons from H1N1, Ebola, and Zika

A new measure of rank correlation

Comparison between observed data and fitted data based on negative binomial model for some countries, territories and areas in WHO regions Page

Association of COVID-19 with pregnancy outcomes in health-care workers and general women

Use of ranks in one-criterion variance analysis

EBOLA and FDA: reviewing the response to the 2014 outbreak, to find lessons for the future

Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia

Mathematical model of infection kinetics and its analysis for COVID-19, SARS and MERS

Genomic characterisation and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding

Modelling of malaria incidence in Akure, Nigeria: negative binomial approach

The significance of daily incidence and mortality cases due to COVID-19 in some African countries

Ebola viral disease in West Africa: a threat to global health, economy and political stability

Number of new coronavirus (COVID-19) cases in Europe from

Modern applied statistics with S

catio ns-detai l/manag ement -of-ill-trave llers -at-point s-of-entry -inter natio nal-airpo rts-seapo rts-andgroun d-cross ings-in-the-conte xt-of-covid -19-outbr eak

Novel Coronavirus (2019-nCoV) situation reports 187

Evolving epidemiology and transmission dynamics of coronavirus disease 2019 outside Hubei province, China: a descriptive and modelling study