<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<?covid-19-tdm?>
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Appl. Math. Stat.</journal-id>
<journal-title>Frontiers in Applied Mathematics and Statistics</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Appl. Math. Stat.</abbrev-journal-title>
<issn pub-type="epub">2297-4687</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">786983</article-id>
<article-id pub-id-type="doi">10.3389/fams.2021.786983</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Applied Mathematics and Statistics</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Modeling of COVID-19 Pandemic vis-&#xe0;-vis Some Socio-Economic Factors</article-title>
<alt-title alt-title-type="left-running-head">Oshinubi et&#x20;al.</alt-title>
<alt-title alt-title-type="right-running-head">Modeling of COVID-19 Pandemic</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Oshinubi</surname>
<given-names>Kayode</given-names>
</name>
<uri xlink:href="https://loop.frontiersin.org/people/1500172/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Rachdi</surname>
<given-names>Mustapha</given-names>
</name>
<uri xlink:href="https://loop.frontiersin.org/people/1387995/overview"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Demongeot</surname>
<given-names>Jacques</given-names>
</name>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
</contrib>
</contrib-group>
<aff>
<institution>Laboratory AGEIS EA 7407, Team Tools for e-Gnosis Medical, Faculty of Medicine, University Grenoble Alpes (UGA)</institution>, <addr-line>La Tronche</addr-line>, <country>France</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/767925/overview">Raluca Eftimie</ext-link>, University of Franche-Comt&#xe9;, France</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/844615/overview">Sania Qureshi</ext-link>, Mehran University of Engineering and Technology, Pakistan</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1169351/overview">Mohd Tahir Ismail</ext-link>, Universiti Sains Malaysia (USM), Malaysia</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Jacques Demongeot, <email>Jacques.Demongeot@univ-grenoble-alpes.fr</email>
</corresp>
<fn fn-type="other">
<p>This article was submitted to Mathematical Biology, a section of the journal Frontiers in Applied Mathematics and Statistics</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>05</day>
<month>01</month>
<year>2022</year>
</pub-date>
<pub-date pub-type="collection">
<year>2021</year>
</pub-date>
<volume>7</volume>
<elocation-id>786983</elocation-id>
<history>
<date date-type="received">
<day>30</day>
<month>09</month>
<year>2021</year>
</date>
<date date-type="accepted">
<day>08</day>
<month>11</month>
<year>2021</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2022 Oshinubi, Rachdi and Demongeot.</copyright-statement>
<copyright-year>2022</copyright-year>
<copyright-holder>Oshinubi, Rachdi and Demongeot</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these&#x20;terms.</p>
</license>
</permissions>
<abstract>
<p>The impact of the COVID-19 epidemic on the socio-economic status of countries around the world should not be underestimated, when we consider the role it has played in various countries. Many people were unemployed, many households were careful about their spending, and a greater social divide in the population emerged in 14 different countries from the Organization for Economic Co-operation and Development (OECD) and from Africa (that is, in developed and developing countries) for which we have considered the epidemiological data on the spread of infection during the first and second waves, as well as their socio-economic data. We established a mathematical relationship between Theil and Gini indices, then we investigated the relationship between epidemiological data and socio-economic determinants, using several machine learning and deep learning methods. High correlations were observed between some of the socio-economic and epidemiological parameters and we predicted three of the socio-economic variables in order to validate our results. These results show a clear difference between the first and the second wave of the pandemic, confirming the impact of the real dynamics of the epidemic&#x2019;s spread in several countries and the means by which it was mitigated.</p>
</abstract>
<kwd-group>
<kwd>COVID-19</kwd>
<kwd>regression</kwd>
<kwd>socio-economic factors</kwd>
<kwd>machine learning</kwd>
<kwd>data analysis</kwd>
</kwd-group>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>Modeling of COVID-19 by scientists, epidemiologists, and health experts was considered early on in the pandemic as it began to ravage the world. The socio-economic determinants of this pandemic taken into account in this modeling are important because they condition the severity of an affected country and the way in which it is controlled, which leads to consider the corresponding variables alongside the daily reproduction rates during the period of contagiousness of individuals infected by the pandemic. The aim of this article is to show the joint variations of socio-economic determinants and epidemiological parameters, which can be observed between developed and developing states and between successive epidemic&#x20;waves.</p>
<p>Some researchers have already worked on socio-economic analysis of the COVID-19 pandemic: in [<xref ref-type="bibr" rid="B1">1</xref>], the authors examined the geoclimatic, demographic, and socio-economic determinants of COVID-19 prevalence and have shown that the influence of these determinants varies by comparing the first and second wave of the pandemic. The socio-economic impact of the COVID-19 pandemic in United&#x20;States of America (United&#x20;States) was studied by Barlow and Vodenska [<xref ref-type="bibr" rid="B2">2</xref>], where the authors investigate the systematic risk posted by sector-level industries within the United&#x20;States. Ahmed et&#x20;al. [<xref ref-type="bibr" rid="B3">3</xref>] modeled daily confirmed cases of COVID-19 in different countries across the globe using regression models with predictions for upcoming scenarios. Kong et&#x20;al. [<xref ref-type="bibr" rid="B4">4</xref>] worked on the socio-economic and environmental factors influencing the basic reproduction number of the COVID-19 pandemic by fitting a logistic growth curve to the reported daily cases up to the first peak of the pandemic while Qiu et&#x20;al. [<xref ref-type="bibr" rid="B5">5</xref>] studied the impact of socio-economic factors on the transmission of COVID-19 disease with China as a case study using an empirical model, and the authors conclude that these determinants have rich implications for ongoing efforts in containing the pandemic. The work in this present article is an extension of [<xref ref-type="bibr" rid="B6">6</xref>], which was based on the analysis of the reproduction numbers of COVID-19 based on the Current Health Expenditure as Gross Domestic Product Percentage (CHE/GDP) across several countries using some machine learning tools. The results of this study show that some countries with a high CHE/GDP improved their public health strategy against the virus during the second wave of the pandemic, fighting it all the more effectively against it the more effectively they were. the most affected during the first wave. The difference with the present study lies in the fact that the latter takes into account data from twice as many public sites [<xref ref-type="bibr" rid="B7">7</xref>&#x2013;<xref ref-type="bibr" rid="B17">17</xref>] and that it is more focused on social inequality, quantified for example by the Social Fracture coefficient (SF equal to the ratio between the incomes of the richest 10% and the poorest 10% of a given population), and the Theil and Gini indices. It was shown previously in [<xref ref-type="bibr" rid="B18">18</xref>] that the Gini index was highly correlated (<italic>r</italic>&#x20;&#x3d; 0.93) to another Demo-economic index denoted DI and equal to the quotient (CHE/GDP)/SF, proving that all these indices are closely related and carry part of the causality of inter-country variations in epidemiological parameters.</p>
<p>The main objectives of this article are to establish a relationship between Theil and Gini index, analyze critically some of the socio-economic determinants of the pandemic, correlate them, predict three of the socio-economic variables, and perform some regression analyses. We have also clustered countries according to these parameters with the help of the lasso (least absolute&#x20;shrinkage and selection operator) method, and we were able to select the best variables for the COVID-19 modeling.</p>
<p>The paper is divided into seven sections: after an introductive section, we explain in Section 2 the methodology used in this research, Section 3 deals with the variables used, Section 4 establishes a mathematical relationship between Theil and Gini index, Section 5 is dedicated to the visualization of the results obtained, while we finally give the discussion and conclusion in Sections 6 and 7, respectively.</p>
</sec>
<sec id="s2">
<title>2 Methods</title>
<p>The use of machine learning methods to analyze data has been helpful over the years to get a proper view on how a model behaves. In this research, we used some supervised and unsupervised machine learning methods and we also tried to use one deep learning method for the identification and visualization of clusters. To jointly interpret the socio-economic and epidemiological data, we have chosen these three main classes of the descriptive statistics, which allow us to compare these socio-economic and epidemiological data. Supervised learning is used in its regression function (prediction of a quantitative variable from annotated examples) and unsupervised learning (in which the data is not labeled) in its classification function. As for deep learning, it makes it possible to create a model from large-scale unlabeled&#x20;data.</p>
<p>The supervised machine learning methods we used are first univariate polynomial regression, linear regression, lasso regression, and ridge regression. We also use some of these methods to make prediction by training the model and testing some percentage of the values. Lasso regression helped us to know the best variables to be used in the modeling. After the univariate regressions, we introduced multivariate least square methods, allowing us to test much more complex relations between variables. It can be represented as follows:<disp-formula id="e1">
<mml:math id="m1">
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mi mathvariant="bold-italic">y</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold">&#x3b2;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi mathvariant="bold">x</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold">&#x3b2;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi mathvariant="bold">x</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mo>&#x22ef;</mml:mo>
<mml:mo>&#x2b;</mml:mo>
<mml:mo>&#x2208;</mml:mo>
<mml:mo>,</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
<label>(1)</label>
</disp-formula>Where <inline-formula id="inf1">
<mml:math id="m2">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold">&#x3b2;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold">&#x3b2;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mo>&#x22ef;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> are coefficients or weights, <inline-formula id="inf2">
<mml:math id="m3">
<mml:mo>&#x2208;</mml:mo>
</mml:math>
</inline-formula> is the residual noise, y is the dependent variable, and <inline-formula id="inf3">
<mml:math id="m4">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold">x</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold">x</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mo>&#x22ef;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> are the independent variables. Ridge and lasso regression are simple methods to reduce the model complexity and prevent over-fitting, which may result from linear regression. The cost function for ridge regression is given below:<disp-formula id="e2">
<mml:math id="m5">
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi mathvariant="bold-italic">n</mml:mi>
</mml:msubsup>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">y</mml:mi>
<mml:mi mathvariant="bold-italic">i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">y</mml:mi>
<mml:mo>&#x5e;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mi mathvariant="bold-italic">i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mstyle>
<mml:mo>&#x3d;</mml:mo>
<mml:mstyle displaystyle="true">
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi mathvariant="bold-italic">n</mml:mi>
</mml:msubsup>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">y</mml:mi>
<mml:mi mathvariant="bold-italic">i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mstyle displaystyle="true">
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">j</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mi mathvariant="bold-italic">m</mml:mi>
</mml:msubsup>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold">&#x3b2;</mml:mi>
<mml:mi mathvariant="bold-italic">j</mml:mi>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">x</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">ij</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mstyle>
<mml:mo>&#x2b;</mml:mo>
<mml:mi mathvariant="bold-italic">&#x3b1;</mml:mi>
<mml:mstyle displaystyle="true">
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">j</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mi mathvariant="bold-italic">m</mml:mi>
</mml:msubsup>
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold">&#x3b2;</mml:mi>
<mml:mi mathvariant="bold-italic">j</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>with for some c &#x3e; 0, <inline-formula id="inf4">
<mml:math id="m6">
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">j</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mi mathvariant="bold-italic">m</mml:mi>
</mml:msubsup>
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold">&#x3b2;</mml:mi>
<mml:mi mathvariant="bold-italic">j</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mo>&#x3c;</mml:mo>
<mml:mi mathvariant="bold-italic">c</mml:mi>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
</mml:math>
</inline-formula>, while <inline-formula id="inf5">
<mml:math id="m7">
<mml:mi mathvariant="bold-italic">&#x3b1;</mml:mi>
</mml:math>
</inline-formula> is the penalty term that regularizes the coefficients such that if the coefficients take large values, the optimization function is penalized. Ridge regression puts constraint on the coefficients <inline-formula id="inf6">
<mml:math id="m8">
<mml:mi mathvariant="bold">&#x3b2;</mml:mi>
</mml:math>
</inline-formula>. We define the cost function for lasso regression in the same way, but by replacing the L2 penalty term by an L1 one as:<disp-formula id="e3">
<mml:math id="m9">
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi mathvariant="bold-italic">n</mml:mi>
</mml:msubsup>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">y</mml:mi>
<mml:mi mathvariant="bold-italic">i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">y</mml:mi>
<mml:mo>&#x5e;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mi mathvariant="bold-italic">i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mstyle>
<mml:mo>&#x3d;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi mathvariant="bold-italic">n</mml:mi>
</mml:msubsup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">y</mml:mi>
<mml:mi mathvariant="bold-italic">i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mstyle displaystyle="true">
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">j</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mi mathvariant="bold-italic">m</mml:mi>
</mml:msubsup>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold">&#x3b2;</mml:mi>
<mml:mi mathvariant="bold-italic">j</mml:mi>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">x</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">ij</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:mi mathvariant="bold-italic">&#x3b1;</mml:mi>
<mml:mstyle displaystyle="true">
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">j</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mi mathvariant="bold-italic">m</mml:mi>
</mml:msubsup>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold">&#x3b2;</mml:mi>
<mml:mi mathvariant="bold-italic">j</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>with for some c &#x3e; 0, <inline-formula id="inf7">
<mml:math id="m10">
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi mathvariant="bold-italic">j</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mi mathvariant="bold-italic">m</mml:mi>
</mml:msubsup>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x7c;</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold">&#x3b2;</mml:mi>
<mml:mi mathvariant="bold-italic">j</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo>&#x7c;</mml:mo>
</mml:mrow>
<mml:mo>&#x3c;</mml:mo>
<mml:mi mathvariant="bold-italic">c</mml:mi>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
<p>After the supervised learning methods, we used unsupervised learning approches to cluster variables across countries and the methods we proposed to validate our results were K-means clustering, Hierarchy clustering, and Principal Component Analysis (PCA). We also performed correlation calculations among parameters used in the modeling step and we chose an optimization method called Ordinary Least Square (OLS) for the socio-economic determinants of COVID-19. Eventually, the deep learning methods we used were Neural Network (NN) and Multi-Layer Perceptron (MLP) regressor, which is a class of feedforward Artificial Neural Network (ANN).</p>
</sec>
<sec id="s3">
<title>3 Variables</title>
<sec id="s3-1">
<title>3.1&#x20;Socio-Economic Variables</title>
<p>Socio-economic variables constitute a strong determinant of the spread of the pandemic. We extracted some observed socio-economic variables from [<xref ref-type="bibr" rid="B7">7</xref>&#x2013;<xref ref-type="bibr" rid="B17">17</xref>], while we also calculated other socio-economic variables like the socio-economic fracture index F, equal to the quotient between CHR/GDP and&#x20;SF.</p>
<p>The observed variables used in this research are immigration rate (IR), average life expectancy (LE), Tuberculosis incidence (TB), temperature, percentage of gross domestic product devoted to health expenditure (CHE/GDP) collated from [<xref ref-type="bibr" rid="B16">16</xref>], percentage of 10% lowest (LI) and 10% highest incomes (HI), government response stringency index (SI), sustainable development goal (SDG) index, human development index (HDI), environmental performance index (EPI), consumer confidence index (CCI), stringency index (SI), Theil index (TI), and Gini index (GI). We collated the data on the available countries and most recent years in public databases [<xref ref-type="bibr" rid="B7">7</xref>&#x2013;<xref ref-type="bibr" rid="B17">17</xref>].</p>
<p>The calculated socio-economic variables are as follows:<list list-type="simple">
<list-item>
<p>- Social fracture (SF) index is the ratio between the 10% highest income and the 10% lowest income. In brief, it is expressed by the equation below:</p>
</list-item>
</list>
<disp-formula id="e4">
<mml:math id="m11">
<mml:mrow>
<mml:mi mathvariant="bold-italic">SF</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>10</mml:mn>
<mml:mo>%</mml:mo>
<mml:mi mathvariant="bold-italic">HI</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>10</mml:mn>
<mml:mo>%</mml:mo>
<mml:mi mathvariant="bold-italic">LI</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
<list list-type="simple">
<list-item>
<p>- Demo-economic (DI) index is the ratio between the percentage of GDP devoted to health expenditure and social fracture index. It is expressed by the equation below:</p>
</list-item>
</list>
<disp-formula id="e5">
<mml:math id="m12">
<mml:mrow>
<mml:mi>D</mml:mi>
<mml:mi>I</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
<mml:mi>E</mml:mi>
</mml:mrow>
<mml:mo>/</mml:mo>
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mi>D</mml:mi>
<mml:mi>P</mml:mi>
</mml:mrow>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mi>F</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>
</p>
<p>We give a precise value of all variables in <xref ref-type="table" rid="T1">Table&#x20;1</xref> and in <xref ref-type="sec" rid="s12">Supplementary Table S2&#x2013;S5</xref> in the supplementary material.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Comparison of different regression models for the prediction.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">S/N</th>
<th align="center">Country name</th>
<th align="center">Gini index</th>
<th align="center">Linear regression</th>
<th align="center">Lasso regression</th>
<th align="center">Ridge regression</th>
<th align="center">MLP regression</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">1</td>
<td align="left">Paraguay</td>
<td align="char" char=".">46.2</td>
<td align="char" char=".">46.5</td>
<td align="char" char=".">46.3</td>
<td align="char" char=".">46.5</td>
<td align="char" char=".">46.0</td>
</tr>
<tr>
<td align="left">2</td>
<td align="left">Panama</td>
<td align="char" char=".">49.2</td>
<td align="char" char=".">50.7</td>
<td align="char" char=".">51.2</td>
<td align="char" char=".">51.3</td>
<td align="char" char=".">51.5</td>
</tr>
<tr>
<td align="left">3</td>
<td align="left">Brazil</td>
<td align="char" char=".">53.9</td>
<td align="char" char=".">58.4</td>
<td align="char" char=".">58.9</td>
<td align="char" char=".">59.6</td>
<td align="char" char=".">61.7</td>
</tr>
<tr>
<td align="left">4</td>
<td align="left">Bolivia</td>
<td align="char" char=".">42.2</td>
<td align="char" char=".">42.2</td>
<td align="char" char=".">43.2</td>
<td align="char" char=".">42.7</td>
<td align="char" char=".">42.0</td>
</tr>
<tr>
<td align="left">5</td>
<td align="left">Honduras</td>
<td align="char" char=".">52.1</td>
<td align="char" char=".">56.6</td>
<td align="char" char=".">57.6</td>
<td align="center">&#x2014;</td>
<td align="char" char=".">59.6</td>
</tr>
<tr>
<td align="left">6</td>
<td align="left">Dominican</td>
<td align="char" char=".">43.7</td>
<td align="char" char=".">43.4</td>
<td align="char" char=".">43.5</td>
<td align="char" char=".">43.4</td>
<td align="char" char=".">43.8</td>
</tr>
<tr>
<td align="left">7</td>
<td align="left">Chile</td>
<td align="char" char=".">46.0</td>
<td align="char" char=".">46.9</td>
<td align="center">&#x2014;</td>
<td align="char" char=".">47.0</td>
<td align="char" char=".">48.3</td>
</tr>
<tr>
<td align="left">8</td>
<td align="left">Mexico</td>
<td align="char" char=".">45.4</td>
<td align="center">&#x2014;</td>
<td align="char" char=".">45.8</td>
<td align="char" char=".">45.7</td>
<td align="center">&#x2014;</td>
</tr>
<tr>
<td align="left">9</td>
<td align="left">Columbia</td>
<td align="char" char=".">50.4</td>
<td align="char" char=".">51.3</td>
<td align="char" char=".">51.4</td>
<td align="char" char=".">51.7</td>
<td align="char" char=".">52.8</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s3-2">
<title>3.2 Epidemiologic Variables</title>
<p>We have six epidemiologic variables: first wave maximum <inline-formula id="inf8">
<mml:math id="m13">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mi>o</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, second wave maximum <inline-formula id="inf9">
<mml:math id="m14">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mi>o</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, first wave deterministic <inline-formula id="inf10">
<mml:math id="m15">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mi>o</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, second wave deterministic <inline-formula id="inf11">
<mml:math id="m16">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mi>o</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, and opposite of the initial autocorrelation slope averaged on 6&#xa0;days for both first and second wave of the daily new cases for developed and developing countries. All epidemiologic variables values were taken from [<xref ref-type="bibr" rid="B18">18</xref>] (see also the Appendix in&#x20;[<xref ref-type="bibr" rid="B6">6</xref>]).</p>
<p>The epidemiologic variables were recorded during the exponential phase of the first and second wave of the pandemic. Daily new cases observed during the first 100&#xa0;days were used to calculate the exponential slope for the first and second wave. The opposite of the initial autocorrelation slope was averaged on 6&#xa0;days for the first and second wave. The maximum <inline-formula id="inf12">
<mml:math id="m17">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mi>o</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> was collated from [<xref ref-type="bibr" rid="B6">6</xref>] while observing this value during the first and second waves of countries considered. We also collated from [<xref ref-type="bibr" rid="B6">6</xref>] the deterministic <inline-formula id="inf13">
<mml:math id="m18">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mi>o</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> for the first and second wave of the pandemic taking 6&#xa0;days as length of contagiousness period.</p>
<p>In this present study, we validated our results by performing cross-validation and also training 80% of the data and training&#x20;30%.</p>
</sec>
</sec>
<sec id="s4">
<title>4 The Relationship Between Theil Index and Gini Index</title>
<sec id="s4-1">
<title>4.1 Mathematical Approach</title>
<p>We first show the relationship between Theil index and Gini index mathematically. The Gini index is defined as follows [<xref ref-type="bibr" rid="B16">16</xref>]:<disp-formula id="e6">
<mml:math id="m19">
<mml:mrow>
<mml:mtext>GI</mml:mtext>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>-</mml:mo>
<mml:mstyle displaystyle="true">
<mml:munder>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mtext>k</mml:mtext>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mtext>n</mml:mtext>
</mml:mrow>
</mml:munder>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mtext>Xk</mml:mtext>
<mml:mo>-</mml:mo>
<mml:mtext>Xk</mml:mtext>
<mml:mo>-</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mtext>yk</mml:mtext>
<mml:mo>-</mml:mo>
<mml:mtext>yk</mml:mtext>
<mml:mo>-</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mstyle>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>-</mml:mo>
<mml:mtext>E</mml:mtext>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mi>&#x394;</mml:mi>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mtext>E</mml:mtext>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>-</mml:mo>
<mml:mi>&#x394;</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>where x<sub>k</sub> (respectively y<sub>k</sub>) denotes the <italic>k</italic>th cumulative part of the population (respectively income). If we choose the population increments, d<sub>k</sub> &#x3d; x<sub>k</sub>-x<sub>k-1</sub> are equal to 1/n, and if E(&#x394;) represents the expectation of the increment, &#x394;<sub>k</sub> &#x3d; y<sub>k&#x2212;</sub>y<sub>k&#x2212;1</sub> for the distribution d<sub>k.</sub> Then, the Theil index applied to the percentage y<sub>k</sub> of the total income relative to a percentage x<sub>k</sub> of the total population ([<xref ref-type="bibr" rid="B17">17</xref>]) is defined by the following equation:<disp-formula id="e7">
<mml:math id="m20">
<mml:mrow>
<mml:mtext>TI</mml:mtext>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>-</mml:mo>
<mml:mstyle displaystyle="true">
<mml:munder>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mtext>k</mml:mtext>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mtext>n</mml:mtext>
</mml:mrow>
</mml:munder>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mtext>y</mml:mtext>
<mml:mtext>k</mml:mtext>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:msub>
<mml:mtext>y</mml:mtext>
<mml:mtext>k</mml:mtext>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mstyle>
<mml:mtext>Log</mml:mtext>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mtext>y</mml:mtext>
<mml:mtext>k</mml:mtext>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:msub>
<mml:mtext>y</mml:mtext>
<mml:mtext>k</mml:mtext>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>-</mml:mo>
<mml:mstyle displaystyle="true">
<mml:munder>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mtext>k</mml:mtext>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mtext>n</mml:mtext>
</mml:mrow>
</mml:munder>
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:mtext>kLog</mml:mtext>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi>&#x394;</mml:mi>
<mml:mtext>k</mml:mtext>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>
</p>
<p>If the first increment of y, &#x394;<sub>1</sub> &#x3d; y<sub>1</sub> &#x2264; 1, is close to 1 [which&#x20;corresponds to a square-shaped Lorenz curve, i.e.,&#x20;closed to a left right triangle-shaped income vs. population curve (in red on <xref ref-type="fig" rid="F1">Figure&#x20;1</xref>), or to a high Gini index close to 1], then we have: <sub>&#x2212;</sub>Log(&#x394;<sub>1</sub>) &#x223c; 1<sub>&#x2212;</sub>&#x394;<sub>1</sub> and &#x394;<sub>k</sub> &#x223c; 0, for k &#x3e; 1. Then, we get:<disp-formula id="e8">
<mml:math id="m21">
<mml:mrow>
<mml:mtext>TI</mml:mtext>
<mml:mo>&#x223c;</mml:mo>
<mml:mo>-</mml:mo>
<mml:mtext>Log</mml:mtext>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x394;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>&#x223c;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>-</mml:mo>
<mml:msub>
<mml:mi>&#x394;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x223c;</mml:mo>
<mml:mtext>GI</mml:mtext>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>the equality being available only if the Lorenz curve presents a perfect left right triangle&#x20;shape.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Lorenz curve showing the cumulated part of income vs. cumulated part of a population having this cumulated income. The curve in red represents left right triangle-shaped Lorenz curve.</p>
</caption>
<graphic xlink:href="fams-07-786983-g001.tif"/>
</fig>
</sec>
<sec id="s4-2">
<title>4.2 Statistical Approach and Application to Latino-American Countries</title>
<sec id="s4-2-1">
<title>4.2.1 Correlation</title>
<p>We correlated both Theil and Gini indices with all epidemiologic, demographic, and socio-economic variables, and as it can be seen in <xref ref-type="fig" rid="F2">Figure&#x20;2G</xref>, Theil and Gini indices are highly positively correlated with coefficient&#x20;0.7.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>
<bold>(A)</bold> Linear regression line (in green) with the confidence interval (in gray). <bold>(B)</bold> Residual plot for the linear regression. <bold>(C)</bold> Partial regression plot. <bold>(D)</bold> Fit plot. <bold>(E)</bold> On left-hand side is the cross-validation plot for the linear regression and on right-hand side is the polynomial regression plot. <bold>(F)</bold> Residual plots for polynomial regression. <bold>(G)</bold> Heat map for the correlations between all variables. <bold>(H)</bold> Neural network visualization.</p>
</caption>
<graphic xlink:href="fams-07-786983-g002.tif"/>
</fig>
</sec>
<sec id="s4-2-2">
<title>4.2.2 Neural Network for Theil Index and Gini Index</title>
<p>We used the <italic>neuralnet</italic> package in R in order to visualize the weights of the network and the bias between Theil and Gini index, and as it can be seen in <xref ref-type="fig" rid="F2">Figure&#x20;2H</xref>, the weights are good with low&#x20;bias.</p>
</sec>
<sec id="s4-2-3">
<title>4.2.3 Regression Analysis Between Theil Index and Gini Index</title>
<p>Linear regression models use some historic data concerning independent and dependent variables and consider a linear relationship between both while polynomial regression models use a similar approach but the dependent variable is modeled as a degree m (m &#x3d; 2 in the present study) polynomial in&#x20;x.</p>
<p>Linear regression model is given as:<disp-formula id="e9">
<mml:math id="m22">
<mml:mrow>
<mml:mtext>y</mml:mtext>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold">&#x3b2;</mml:mi>
<mml:mi mathvariant="bold">o</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mstyle displaystyle="true">
<mml:munder>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:munder>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold">&#x3b2;</mml:mi>
<mml:mi mathvariant="bold">i</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi mathvariant="bold">x</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold">i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mo>&#x2208;</mml:mo>
<mml:mi mathvariant="bold-italic">i</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>where <inline-formula id="inf14">
<mml:math id="m23">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold">&#x3b2;</mml:mi>
<mml:mi mathvariant="bold">i</mml:mi>
</mml:msub>
<mml:mi mathvariant="bold-italic">&#x27;s</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> are the weights, <inline-formula id="inf15">
<mml:math id="m24">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold">&#x3b2;</mml:mi>
<mml:mi mathvariant="bold">o</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the intercept and <inline-formula id="inf16">
<mml:math id="m25">
<mml:mo>&#x2208;</mml:mo>
</mml:math>
</inline-formula> is the random error term. The above equation is the linear equation that needs to be obtained with the minimum error. Polynomial regression of order 2 is given below:<disp-formula id="e10">
<mml:math id="m26">
<mml:mrow>
<mml:mtext>y</mml:mtext>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold">&#x3b2;</mml:mi>
<mml:mi mathvariant="bold">o</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold">&#x3b2;</mml:mi>
<mml:mi mathvariant="bold">i</mml:mi>
</mml:msub>
<mml:msub>
<mml:mtext>x</mml:mtext>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mtext>&#xa0;</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mstyle displaystyle="true">
<mml:munder>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:munder>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold">&#x3b2;</mml:mi>
<mml:mi mathvariant="bold">i</mml:mi>
</mml:msub>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">x</mml:mi>
<mml:mi mathvariant="bold-italic">i</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mo>&#x2208;</mml:mo>
<mml:mi mathvariant="bold-italic">i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mstyle>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>
</p>
<p>We present the visualization of the regression results using this approach in <xref ref-type="fig" rid="F2">Figures 2A&#x2013;F</xref>.</p>
<p>For the linear regression as shown in <xref ref-type="fig" rid="F2">Figure&#x20;2A</xref>, the intercept is 31.03, <italic>p</italic>-value is 0.0181, <italic>R</italic>
<sup>2</sup> is 0.4881, residual standard error is 3.116, and all coefficients are significant with <italic>p</italic>&#x20;&#x3c; 0.05 for both the train and test data for linear and polynomial regression. The median of the residual plot in <xref ref-type="fig" rid="F2">Figures 2B,F</xref> are 0.2111 and 0.2566, respectively, for both linear and polynomial regression, which are low values. The normality of the residual was tested using Jarque-Bera and Durbin-Watson tests, which gave a high <italic>p</italic>-value, and we failed to reject the null hypothesis that the skewness and kurtosis of the residuals are statistically equal to zero. In order to know the performance of the linear regression model, we trained 80% of the data and tested the 20% remaining data and also did cross-validation to be sure of the accuracy. The predicted and the observed values are very close to the results presented for the regression models used. For the linear model, we present the cross-validation result in <xref ref-type="fig" rid="F2">Figure&#x20;2E</xref> whose average mean square error for the five portion folds is 11.72794. We observed that the correlation between the tested and the predicted values has high accuracy (<italic>R</italic>
<sup>2</sup> &#x3d; 0.97). The test set <italic>p</italic>-value is 0.02 with a residual standard error of 3.528. For polynomial regression of order 2, the train set has the following results: <italic>R</italic>
<sup>2</sup> &#x3d; 0.6, <italic>p</italic>-value &#x3d; 0.002, and residual standard error &#x3d; 2.935. The test set has the following results: <italic>R</italic>
<sup>2</sup> &#x3d; 0.99, <italic>p</italic>-value &#x3d; 0.008, and residual standard error &#x3d; 0.5639.</p>
</sec>
<sec id="s4-2-4">
<title>4.2.4 Multivariate Analysis for Gini Index and Theil Index Alongside Other Socio-Economic Variables and Epidemiologic Variables</title>
<p>
<xref ref-type="fig" rid="F3">Figure&#x20;3</xref> corresponds to the ordinary multivariate least square methods with <italic>R</italic>
<sup>2</sup> &#x3d; 0.674. <xref ref-type="fig" rid="F3">Figure&#x20;3A</xref> shows Paraguay as outlier not fitting data, <xref ref-type="fig" rid="F3">Figure&#x20;3B</xref> normalizes all countries and does not point any country in the&#x20;plot.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>
<bold>(A)</bold> Influence plot. <bold>(B)</bold> Leverage vs. Normalized residuals squared plot. <bold>(C)</bold> Partial regression plot. <bold>(D)</bold> Component&#x2013;Component plus residual&#x20;plot.</p>
</caption>
<graphic xlink:href="fams-07-786983-g003.tif"/>
</fig>
</sec>
<sec id="s4-2-5">
<title>4.2.5 Prediction of Gini Index Using MLP Regressor, and Linear, Lasso, and Ridge Regression</title>
<p>In this section, we used cross-validation method to choose the best parameter <inline-formula id="inf17">
<mml:math id="m27">
<mml:mi>&#x3b1;</mml:mi>
</mml:math>
</inline-formula> for the modeling as shown in <xref ref-type="fig" rid="F4">Figure&#x20;4C</xref>. For ridge regression, <inline-formula id="inf18">
<mml:math id="m28">
<mml:mrow>
<mml:mi>&#x3b1;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0.142</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> with a mean square error of 1.36 and <inline-formula id="inf19">
<mml:math id="m29">
<mml:mrow>
<mml:mi>&#x3b1;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0.368</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> for lasso regression with a mean square error &#x3d; 5.10. For <xref ref-type="fig" rid="F4">Figure&#x20;4E</xref>, training score &#x3d; 1.000 and test score &#x3d; 0.641; for <xref ref-type="fig" rid="F4">Figure&#x20;4F</xref>, training score &#x3d; 0.992 and test score &#x3d; 0.497; for <xref ref-type="fig" rid="F4">Figure&#x20;4G</xref>, training score &#x3d; 0.99 and test score &#x3d; 0.406; and for <xref ref-type="fig" rid="F4">Figure&#x20;4H</xref>, training score &#x3d; 0.984 and test&#x20;score&#x20;&#x3d; &#x2212;0.077. It is evident from these results that linear regression best predicts Gini index with the highest test score, and predicted values are very close to each other as presented in <xref ref-type="table" rid="T1">Table&#x20;1</xref>. Also, we observed the same pattern of prediction in <xref ref-type="fig" rid="F4">Figures 4E&#x2013;H</xref> showing that all methods used in this section have the same predictive behavior.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>
<bold>(A)</bold> Residual plot for lasso regression. <bold>(B)</bold> Residual plot for ridge regression. <bold>(C)</bold> Lasso regression cross-validation error. <bold>(D)</bold> Prediction error for the lasso regression. <bold>(E)</bold> Linear regression prediction plot. <bold>(F)</bold> Lasso regression and <bold>(G)</bold> ridge regression prediction plots. <bold>(H)</bold> MLP regression prediction&#x20;plot.</p>
</caption>
<graphic xlink:href="fams-07-786983-g004.tif"/>
</fig>
</sec>
<sec id="s4-2-6">
<title>4.2.6 Clustering Analysis of Latino-American Countries for Gini Index and Theil Index Alongside Other Socio-Economic Variables and Epidemiologic Variables</title>
<p>In <xref ref-type="fig" rid="F5">Figure&#x20;5C</xref>, the first two clusters have 14 countries and the third has three countries, which are Uruguay and El Salvador on the same hierarchy while Argentina is on another hierarchy. We only show the cluster dendrogram for the first cluster. In <xref ref-type="fig" rid="F5">Figure&#x20;5F</xref>, the Gini index has the highest positive correlation of 0.44 with the principal component PC 1 and Theil index has only the value 0.34 with PC 1. The main variable causing the separation into three classes is the Gini&#x20;index.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>
<bold>(A)</bold> Scree plot. <bold>(B)</bold> Plot for projection of points for PC1 and PC2. <bold>(C)</bold> Hierarchy clustering dendrogram. <bold>(D)</bold> Parallel coordinates plot for the clusters. <bold>(E,F)</bold> PC&#x2019;s visualization. <bold>(G)</bold> Box plot for the clusters. <bold>(H)</bold> Parallel coordinates plot for the centroids.</p>
</caption>
<graphic xlink:href="fams-07-786983-g005.tif"/>
</fig>
</sec>
</sec>
</sec>
<sec id="s5">
<title>5 Application of the Methods to OECD Countries, African Countries, and Developed and Developing Countries</title>
<sec id="s5-1">
<title>5.1 Developed and Developing Countries</title>
<sec id="s5-1-1">
<title>5.1.1 Regression and Multivariate Analysis for Socio-Economic and Epidemiologic Variables</title>
<p>In <xref ref-type="fig" rid="F6">Figures 6C,D</xref>, we modeled the dependent variable as a degree n (<italic>n</italic>&#x20;&#x3d; 6 in the present study) polynomial in x, an extension of <xref ref-type="disp-formula" rid="e8">Eq. 8</xref>. <xref ref-type="fig" rid="F6">Figure&#x20;6</xref> presents regression analyses with the parameters presented in <xref ref-type="sec" rid="s12">Supplementary Table S6</xref> in the supplementary material.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Linear regression plots for <bold>(A)</bold> the first wave slope vs. the demo-economic index for the developed and developing countries, <bold>(B)</bold> second wave slope vs. immigration rate for developed countries, <bold>(C)</bold> the opposite of the initial auto-correlation slope for first wave vs. social fracture index, <bold>(D)</bold> opposite of initial autocorrelation slope for second wave vs. social fracture index. <bold>(E)</bold> Influence plot. <bold>(F)</bold> Leverage vs. normalized residuals squared plot. <bold>(G)</bold> Partial regression plot. <bold>(H)</bold> Component&#x2013;Component plus residual&#x20;plot.</p>
</caption>
<graphic xlink:href="fams-07-786983-g006.tif"/>
</fig>
<p>
<xref ref-type="fig" rid="F6">Figures 6E&#x2013;H</xref> correspond to the ordinary multivariate least square method with <italic>R</italic>
<sup>2</sup> &#x3d; 0.76. <xref ref-type="fig" rid="F6">Figure&#x20;6E</xref> shows some developing countries as outliers, while Belgium is the only developed country, which does not fit the&#x20;data.</p>
</sec>
<sec id="s5-1-2">
<title>5.1.2 Prediction of Percentage GDP Devoted to Health Expenditure</title>
<p>In this section, we used a cross-validation method to choose the best parameter <inline-formula id="inf20">
<mml:math id="m30">
<mml:mi>&#x3b1;</mml:mi>
</mml:math>
</inline-formula> for the modeling as shown in <xref ref-type="fig" rid="F7">Figure&#x20;7C</xref>. For ridge regression, <inline-formula id="inf21">
<mml:math id="m31">
<mml:mrow>
<mml:mi>&#x3b1;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0.012</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> with a mean square error &#x3d; 2.32 and <inline-formula id="inf22">
<mml:math id="m32">
<mml:mrow>
<mml:mi>&#x3b1;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0.029</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> for the lasso regression with a mean square error &#x3d; 2.21. For <xref ref-type="fig" rid="F7">Figure&#x20;7E</xref>, training score &#x3d; 0.983 and test score &#x3d; 0.607; for <xref ref-type="fig" rid="F7">Figure&#x20;7F</xref>, training score &#x3d; 0.170 and test score &#x3d; 0.021; for <xref ref-type="fig" rid="F7">Figure&#x20;7G</xref>, training score &#x3d; 0.854 and test score &#x3d; 0.115; and for <xref ref-type="fig" rid="F7">Figure&#x20;7H</xref>, training score &#x3d; 0.980 and test score &#x3d; &#x2212;2.386.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>
<bold>(A)</bold> Residual plot for lasso regression. <bold>(B)</bold> Residual plot for ridge regression. <bold>(C)</bold> Lasso regression cross-validation error. <bold>(D)</bold> Prediction error for lasso regression curve (<italic>R</italic>&#x20;&#x3d; 0.8). <bold>(E)</bold> Linear regression prediction plot. <bold>(E)</bold> Linear regression prediction plot. <bold>(F)</bold> Lasso regression prediction plot. <bold>(G)</bold> Ridge regression prediction plot. <bold>(H)</bold> MLP regression prediction&#x20;plot.</p>
</caption>
<graphic xlink:href="fams-07-786983-g007.tif"/>
</fig>
<p>It is evident from the results that linear regression best&#x20;predicts GDP percentage devoted to health expenditure with the highest test score, and all predicted values are very&#x20;close.</p>
</sec>
<sec id="s5-1-3">
<title>5.1.3 Principal Component Analysis and Clustering Results</title>
<p>In <xref ref-type="fig" rid="F8">Figures 8E,F</xref>, the first cluster has 15 countries, the second cluster has 27 countries, while the last cluster has 2 countries, which are Tanzania and Mauritius. We only show the first two cluster dendrograms. With PC1, Gini index GI has the highest positive correlation of 0.52 with PC1 and demo-economic index DI has the second highest negative correlation of &#x2212;0.53, while with PC2, first wave maximum R<sub>0</sub> has the highest positive correlation of 0.70 (<xref ref-type="fig" rid="F8">Figure&#x20;8C</xref>).</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>
<bold>(A)</bold> Heatmap of the parameter&#x2019;s correlations. <bold>(B)</bold> Scree plot. <bold>(C,D)</bold> PCs visualization. </p>
</caption>
<graphic xlink:href="fams-07-786983-g008.tif"/>
</fig>
<p>The first cluster contains a majority of developing countries, and the second cluster contains a majority of developed countries, the main variable causing the separation into two classes being the Gini index in&#x20;PC1.</p>
</sec>
</sec>
<sec id="s5-2">
<title>5.2 African Countries</title>
<sec id="s5-2-1">
<title>5.2.1 Multivariate Analysis for Socio-Economic Variables and Epidemiologic Variables</title>
<p>
<xref ref-type="fig" rid="F9">Figures 9A&#x2013;D</xref> correspond to the ordinary multivariate least square method with <italic>R</italic>
<sup>2</sup> &#x3d; 0.60. <xref ref-type="fig" rid="F9">Figure&#x20;9A</xref> shows Botswana and Tanzania as outliers not fitting the&#x20;data.</p>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>
<bold>(A)</bold> Influence plot. <bold>(B)</bold> Leverage vs. normalized residuals squared plot. <bold>(C)</bold> Partial regression plot of the temperature vs. maximal R<sub>0</sub> of the first wave. <bold>(D)</bold> Fit plot. The first wave maximal R<sub>0</sub> is negatively correlated with the mean temperature of the country ([<xref ref-type="bibr" rid="B19">19</xref>,<xref ref-type="bibr" rid="B20">20</xref>]).</p>
</caption>
<graphic xlink:href="fams-07-786983-g009.tif"/>
</fig>
</sec>
<sec id="s5-2-2">
<title>5.2.2 Prediction of Temperature</title>
<p>In this section, we used the cross-validation method to choose the best parameter <inline-formula id="inf23">
<mml:math id="m33">
<mml:mi>&#x3b1;</mml:mi>
</mml:math>
</inline-formula> for the modeling as shown in <xref ref-type="fig" rid="F10">Figure&#x20;10C</xref>. For ridge regression, <inline-formula id="inf24">
<mml:math id="m34">
<mml:mrow>
<mml:mi>&#x3b1;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1.005</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> with a mean square error of 19.13, and for lasso regression, <inline-formula id="inf25">
<mml:math id="m35">
<mml:mrow>
<mml:mi>&#x3b1;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>6.018</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> with a mean square error &#x3d; 16.93. For <xref ref-type="fig" rid="F10">Figure&#x20;10E</xref>, training score &#x3d; 0.647 and test score &#x3d; &#x2212;2.228; for <xref ref-type="fig" rid="F10">Figure&#x20;10F</xref>, training score &#x3d; 0.316 and test score &#x3d; 0.154; for <xref ref-type="fig" rid="F10">Figure&#x20;10G</xref>, training score &#x3d; 0.573 and test score &#x3d; &#x2212;1.136; and for <xref ref-type="fig" rid="F10">Figure&#x20;10H</xref>, training score &#x3d; &#x2212;6.728 and test score &#x3d; &#x2212;4.714. It is evident from these results that the lasso regression best predicts temperature with the highest test score, and predicted values of temperature for lasso and ridge regression are&#x20;close.</p>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>
<bold>(A)</bold> Residual plot for lasso regression. <bold>(B)</bold> Residual plot for ridge regression. <bold>(C)</bold> Lasso regression cross-validation error. <bold>(D)</bold> Prediction error for lasso regression. </p>
</caption>
<graphic xlink:href="fams-07-786983-g010.tif"/>
</fig>
<p>All the regression methods give about the same result with the maximum accuracy for the ridge regression.</p>
</sec>
<sec id="s5-2-3">
<title>5.2.3 Principal Component Analysis and Clustering Results</title>
<p>In <xref ref-type="fig" rid="F11">Figures 11E,F</xref>, the first cluster has 40 countries, the second&#x20;cluster has 13 countries, while the last cluster has only 1 country, which is Botswana. We only show the two&#x20;cluster dendrograms with many countries. In <xref ref-type="fig" rid="F11">Figure&#x20;11C</xref>, average life expectancy has the highest positive correlation of 0.46 in PC 1 while first wave deterministic <inline-formula id="inf26">
<mml:math id="m36">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold">R</mml:mi>
<mml:mi mathvariant="bold">o</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> (first wave D) has the highest positive correlation in PC 2, equal to&#x20;0.47.</p>
<fig id="F11" position="float">
<label>FIGURE 11</label>
<caption>
<p>
<bold>(A)</bold> Heatmap of the parameter&#x2019;s correlations. <bold>(B)</bold> Scree plot. <bold>(C,D)</bold> PCs visualization. </p>
</caption>
<graphic xlink:href="fams-07-786983-g011.tif"/>
</fig>
<p>The two socio-economic variables explaining the most clustering are the average life expectancy (LE) and the stringency index&#x20;(SI).</p>
</sec>
</sec>
<sec id="s5-3">
<title>5.3 OECD Countries</title>
<sec id="s5-3-1">
<title>5.3.1 Multivariate Analysis for Socio-Economic Variables and Epidemiologic Variables</title>
<p>
<xref ref-type="fig" rid="F12">Figure&#x20;12</xref> corresponds to the ordinary multivariate least square method with <italic>R</italic>
<sup>2</sup> &#x3d; 0.90. <xref ref-type="fig" rid="F12">Figure&#x20;12A</xref> shows Iceland, United&#x20;States, Austria, and Belgium as outliers not fitting the&#x20;data.</p>
<fig id="F12" position="float">
<label>FIGURE 12</label>
<caption>
<p>
<bold>(A)</bold> Influence plot. <bold>(B)</bold> Leverage vs. Normalized residuals squared plot. <bold>(C)</bold> Partial regression plot. <bold>(D)</bold> Component&#x2013;Component plus residual&#x20;plot.</p>
</caption>
<graphic xlink:href="fams-07-786983-g012.tif"/>
</fig>
<p>We see on the partial regression plots of the <xref ref-type="fig" rid="F12">Figure&#x20;12C</xref> that the best correlation observed between parameters is between CHE/GDP and the demo-economic index DI as observed before in&#x20;[<xref ref-type="bibr" rid="B21">21</xref>].</p>
</sec>
<sec id="s5-3-2">
<title>5.3.2 Prediction of Percentage GDP Health Expenditure</title>
<p>In this section, we used the cross-validation method to choose the best parameter <inline-formula id="inf27">
<mml:math id="m37">
<mml:mi>&#x3b1;</mml:mi>
</mml:math>
</inline-formula> for the modeling as shown in <xref ref-type="fig" rid="F13">Figure&#x20;13D</xref>. For ridge regression, <inline-formula id="inf28">
<mml:math id="m38">
<mml:mrow>
<mml:mi>&#x3b1;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0.005</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> with a mean square error of 1.905, and for Lasso regression, <inline-formula id="inf29">
<mml:math id="m39">
<mml:mrow>
<mml:mi>&#x3b1;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mtext>&#xa0;</mml:mtext>
<mml:mn>0.027</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> with a mean square error &#x3d; 1.657. For <xref ref-type="fig" rid="F13">Figure&#x20;13E</xref>, training score &#x3d; 0.993 and test score &#x3d; 0.535; for <xref ref-type="fig" rid="F13">Figure&#x20;13F</xref>, training score &#x3d; 0.898 and test score &#x3d; 0.629; for <xref ref-type="fig" rid="F13">Figure&#x20;13G</xref>, training score &#x3d; 0.983 and test score &#x3d; 0.259; and for <xref ref-type="fig" rid="F13">Figure&#x20;13H</xref>, training score &#x3d; &#x2212;0.072 and test score &#x3d; &#x2212;0.196. It is evident from these results that the lasso regression best predicts percentage of GDP devoted to health expenditure with the highest test score and predicted values are very&#x20;close.</p>
<fig id="F13" position="float">
<label>FIGURE 13</label>
<caption>
<p>
<bold>(A)</bold> Residual plot for lasso regression. <bold>(B)</bold> Residual plot ridge regression. <bold>(C)</bold> Prediction error for lasso regression. <bold>(D)</bold> Lasso regression cross-validation error. </p>
</caption>
<graphic xlink:href="fams-07-786983-g013.tif"/>
</fig>
<p>All the regression methods give about the same result with the maximum accuracy for the ridge regression.</p>
</sec>
<sec id="s5-3-3">
<title>5.3.3 Principal Component Analysis and Clustering Results</title>
<p>In <xref ref-type="fig" rid="F14">Figures 14E,F</xref>, the first cluster has 20 countries, and the&#x20;second has 5 countries, which are United&#x20;States and Bulgaria on the same hierarchy, Mexico and Costa Rica on the same hierarchy, and Chile standing alone. The third cluster has 12 countries. We only show the two highest cluster&#x20;dendrograms. In <xref ref-type="fig" rid="F14">Figure&#x20;14C</xref>, the Gini index and&#x20;social fracture index have the highest positive correlation of 0.45 and 0.46, respectively, in PC 1 while&#x20;the percentage of GDP devoted to health expenditure and demo-economic index have the highest positive correlation in PC 2, whose values equal to 0.65 and 0.41, respectively.</p>
<fig id="F14" position="float">
<label>FIGURE 14</label>
<caption>
<p>
<bold>(A)</bold> Heatmap of the parameter&#x2019;s correlations. <bold>(B)</bold> Scree plot. <bold>(C,D)</bold> PCs visualization. </p>
</caption>
<graphic xlink:href="fams-07-786983-g014.tif"/>
</fig>
<p>The two main clusters correspond both to developed countries, but in the first, countries are more continental, and in the second, countries are more maritime, which could be explained by their difference in consumer confidence index (CCI), which is less important in maritime countries than in continental&#x20;ones.</p>
</sec>
</sec>
</sec>
<sec id="s6">
<title>6 Discussion</title>
<p>We have been able to develop new approaches to the socio-economic determinants for the modeling of the COVID-19 pandemic during the exponential phase. Some of these determinants have shown high correlation with epidemiologic parameters as it can be seen in the heatmap diagrams in <xref ref-type="fig" rid="F2">Figures 2G</xref>, <xref ref-type="fig" rid="F8">8A</xref>, <xref ref-type="fig" rid="F11">11A</xref>, <xref ref-type="fig" rid="F14">14A</xref>, explaining the role of each variable thanks to these correlations.</p>
<p>For developed and developing countries, the lasso regression reduced the correlation between the social fracture index and the 10% highest income, while for OECD countries, the correlation between the Gini index and social fracture index was reduced to zero. Some of our variables were not used in the optimization method-OLS due to multicollinearity observed on results summary. For the two sets of countries, consumer confidence index, opposite of the initial autocorrelation slope averaged on 6&#xa0;days for the first and second wave, 10% lowest income, and 10% highest income were not used in the modeling. The <italic>R</italic>
<sup>2</sup> for OLS results for developed, developing, and OECD countries are 0.76 and 0.90, respectively, which shows a high significance rate (<xref ref-type="fig" rid="F8">Figures 8E</xref>,F,&#x20;<xref ref-type="fig" rid="F12">12</xref>).</p>
<p>The principal component analysis shows high correlation for the numbers of new cases we used in this research. The social fracture index has high correlation in PC1 for both cases, while&#x20;in&#x20;PC2, percentage of GDP devoted to health expenditure was dominant for OECD countries, and maximum <inline-formula id="inf30">
<mml:math id="m40">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold">R</mml:mi>
<mml:mi mathvariant="bold">o</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> for the first wave was dominant for both developed and developing countries (see <xref ref-type="fig" rid="F8">Figures 8C</xref>, <xref ref-type="fig" rid="F14">14C</xref>). We can deduce from all these observations that the socio-economic determinants are a key to the modeling of infectious diseases like COVID-19 as these parameters give high signals on the trend during the spread of the pandemic for various countries ([<xref ref-type="bibr" rid="B22">22</xref>&#x2013;<xref ref-type="bibr" rid="B30">30</xref>]).</p>
</sec>
<sec id="s7">
<title>7 Conclusion</title>
<p>The systematic study of the correlations between socio-economic variables (Gini and Theil indices, percentage of GDP devoted to heath expenditure, etc.) and epidemiological variables (reproduction rate, opposite of the slope of autocorrelation to origin, etc.) shows a disparity between developed and developing countries, as well as between epidemic waves. Developed countries with high indices of social divide, but high health expenditure, did not, for the first wave, react better to the COVID-19 epidemic than developing countries. On the other hand, the rapid implementation of isolation and vaccination measures enabled them to anticipate and reduce the effects of the second wave. In a subsequent work, we will study the evolution of this disparity between developed and developing countries during subsequent waves of SARS CoV-2.</p>
</sec>
</body>
<back>
<sec id="s8">
<title>Data Availability Statement</title>
<p>The original contributions presented in the study are included in the article/<xref ref-type="sec" rid="s12">Supplementary Material</xref>. Further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s9">
<title>Author Contributions</title>
<p>Conceptualization, JD, KO, and MR. Methodology, JD, KO,&#x20;and MR. Software, KO. Validation, JD; KO, and MR. Formal analysis, KO. Investigation, JD and MR. Resources,&#x20;JD. Data curation, KO. Writing&#x2014;original draft preparation, KO.&#x20;Writing&#x2014;review and editing, JD and KO. Visualization, KO. Supervision, JD and MR. Project administration, JD and MR. All authors&#x20;have&#x20;read&#x20;and agreed to the final version of the manuscript.</p>
</sec>
<sec sec-type="COI-statement" id="s10">
<title>Conflict of Interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s11">
<title>Publisher&#x2019;s Note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ack>
<p>The authors wish to acknowledge the Petroleum Technology Development Fund (PTDF) Nigeria doctoral fellowship in collaboration with Campus France Africa&#x20;Unit.</p>
</ack>
<sec id="s12">
<title>Supplementary Material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fams.2021.786983/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fams.2021.786983/full&#x23;supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="DataSheet1.docx" id="SM1" mimetype="application/docx" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<label>1.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Demongeot</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Oshinubi</surname>
<given-names>K</given-names>
</name>
<name>
<surname>Rachdi</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Seligmann</surname>
<given-names>H</given-names>
</name>
</person-group>. <article-title>Geoclimatic, Demographic and Socio-Economic Determinants of the Covid-19 Prevalence</article-title>. <source>EGU Gen Assembly</source> (<year>2021</year>) <fpage>EGU21-7976</fpage>. <pub-id pub-id-type="doi">10.5194/egusphere-egu21-7976</pub-id> </citation>
</ref>
<ref id="B2">
<label>2.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Barlow</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Vodenska</surname>
<given-names>I</given-names>
</name>
</person-group>. <article-title>Socio-Economic Impact of the Covid-19 Pandemic in the U.S</article-title>. <source>Entropy</source> (<year>2021</year>) <volume>23</volume>:<fpage>673</fpage>. <pub-id pub-id-type="doi">10.3390/e23060673</pub-id> </citation>
</ref>
<ref id="B3">
<label>3.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ahmed</surname>
<given-names>HM</given-names>
</name>
<name>
<surname>Elbarkouky</surname>
<given-names>RA</given-names>
</name>
<name>
<surname>Omar</surname>
<given-names>OAM</given-names>
</name>
<name>
<surname>Ragusa</surname>
<given-names>MA</given-names>
</name>
</person-group>. <article-title>Models for COVID-19 Daily Confirmed Cases in Different Countries</article-title>. <source>Mathematics</source> (<year>2021</year>) <volume>9</volume>:<fpage>659</fpage>. <pub-id pub-id-type="doi">10.3390/math9060659</pub-id> </citation>
</ref>
<ref id="B4">
<label>4.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kong</surname>
<given-names>JD</given-names>
</name>
<name>
<surname>Tekwa</surname>
<given-names>EW</given-names>
</name>
<name>
<surname>Gignoux-Wolfsohn</surname>
<given-names>SA</given-names>
</name>
</person-group>. <article-title>Social, Economic, and Environmental Factors Influencing the Basic Reproduction Number of COVID-19 across Countries</article-title>. <source>PLoS ONE</source> (<year>2021</year>) <volume>16</volume>:<fpage>e0252373</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0252373</pub-id> </citation>
</ref>
<ref id="B5">
<label>5.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Qiu</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>X</given-names>
</name>
<name>
<surname>Shi</surname>
<given-names>W</given-names>
</name>
</person-group>. <article-title>Impacts of Social and Economic Factors on&#x20;the&#x20;Transmission of Coronavirus Disease 2019 (COVID-19) in China</article-title>. <source>J&#x20;Popul Econ</source> (<year>2020</year>) <volume>33</volume>:<fpage>1127</fpage>&#x2013;<lpage>72</lpage>. <pub-id pub-id-type="doi">10.1007/s00148-020-00778-2</pub-id> </citation>
</ref>
<ref id="B6">
<label>6.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Oshinubi</surname>
<given-names>K</given-names>
</name>
<name>
<surname>Rachdi</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Demongeot</surname>
<given-names>J</given-names>
</name>
</person-group>. <article-title>Analysis of Reproduction Number R0 of COVID-19 Using Current Health Expenditure as Gross Domestic Product Percentage (CHE/GDP) across Countries</article-title>. <source>Healthcare</source> (<year>2021</year>) <volume>9</volume>:<fpage>1247</fpage>. <pub-id pub-id-type="doi">10.3390/healthcare9101247</pub-id> </citation>
</ref>
<ref id="B7">
<label>7.</label>
<citation citation-type="web">
<collab>World bank</collab>. <article-title>World Bank</article-title> (<year>2021</year>). <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://data.worldbank.org/">https://data.worldbank.org/</ext-link>
</comment>(<comment>accessed on February 12, 2021)</comment>. </citation>
</ref>
<ref id="B8">
<label>8.</label>
<citation citation-type="web">
<collab>OECD</collab>. <article-title>Organisation for Economic Co-operation and Development</article-title> (<year>2021</year>). <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://data.oecd.org/">https://data.oecd.org/</ext-link>
</comment>(<comment>accessed on February 12, 2021)</comment>. </citation>
</ref>
<ref id="B9">
<label>9.</label>
<citation citation-type="web">
<collab>SDSN</collab>. <article-title>2019 Africa SDG Index and Dashboards Report</article-title> (<year>2021</year>). <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://www.sdgindex.org/reports/2019-africa-sdg-index-and-dashboards-report/">https://www.sdgindex.org/reports/2019-africa-sdg-index-and-dashboards-report/</ext-link>
</comment>(<comment>accessed on February 22, 2021)</comment>. </citation>
</ref>
<ref id="B10">
<label>10.</label>
<citation citation-type="web">
<collab>UNDP</collab>. <article-title>Adjusted Net Savings (% of GNI) Dimension: Socio-Economic Sustainability</article-title> (<year>2021</year>). <comment>Available at: <ext-link ext-link-type="uri" xlink:href="http://hdr.undp.org/en/indicators/164406">http://hdr.undp.org/en/indicators/164406</ext-link> (accessed on February 22, 2021)</comment>. </citation>
</ref>
<ref id="B11">
<label>11.</label>
<citation citation-type="web">
<collab>Our World in Data</collab>. <article-title>Covid-Stringency</article-title> (<year>2021</year>). <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://ourworldindata.org/grapher/covid-stringency-index?tab=table">https://ourworldindata.org/grapher/covid-stringency-index?tab&#x3d;table</ext-link> (accessed on February 22, 2021)</comment>. </citation>
</ref>
<ref id="B12">
<label>12.</label>
<citation citation-type="web">
<collab>Wikipedia</collab>. <article-title>List of Countries by Average Yearly Temperature</article-title> (<year>2021</year>). <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://en.wikipedia.org/wiki/List_of_countries_by_average_yearly_temperature">https://en.wikipedia.org/wiki/List_of_countries_by_average_yearly_temperature</ext-link> (accessed on February 22, 2021)</comment>. </citation>
</ref>
<ref id="B13">
<label>13.</label>
<citation citation-type="web">
<collab>Statista</collab>. <article-title>Life Expectancy at Birth in Africa as of 2019, by Country</article-title> (<year>2021</year>). <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://www.statista.com/statistics/1218173/life-expectancy-in-african-countries/">https://www.statista.com/statistics/1218173/life-expectancy-in-african-countries/</ext-link>
</comment>(<comment>accessed on February 22, 2021)</comment>. </citation>
</ref>
<ref id="B14">
<label>14.</label>
<citation citation-type="web">
<collab>EPI</collab>. <article-title>Environmental Performance Index</article-title> (<year>2021</year>). <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://epi.yale.edu/">https://epi.yale.edu/</ext-link>
</comment>(<comment>accessed on February 22, 2021)</comment>. </citation>
</ref>
<ref id="B15">
<label>15.</label>
<citation citation-type="web">
<collab>Wikipedia</collab>. <article-title>List of Countries by Net Migration Rate</article-title> (<year>2021</year>). <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://en.wikipedia.org/wiki/List_of_countries_by_net_migration_rate">https://en.wikipedia.org/wiki/List_of_countries_by_net_migration_rate</ext-link> (accessed on February 22, 2021)</comment>. </citation>
</ref>
<ref id="B16">
<label>16.</label>
<citation citation-type="web">
<collab>Wikipedia</collab>. <article-title>Coefficient de Gini</article-title> (<year>2021</year>). <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://fr.wikipedia.org/wiki/Coefficient_de_Gini">https://fr.wikipedia.org/wiki/Coefficient_de_Gini</ext-link> (accessed on August 22, 2021)</comment>. </citation>
</ref>
<ref id="B17">
<label>17.</label>
<citation citation-type="web">
<collab>Wikipedia</collab>. <article-title>Indice de Theil</article-title> (<year>2021</year>). <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://fr.wikipedia.org/wiki/Indice_de_Theil">https://fr.wikipedia.org/wiki/Indice_de_Theil</ext-link> (accessed on August 22, 2021)</comment>. </citation>
</ref>
<ref id="B18">
<label>18.</label>
<citation citation-type="web">
<collab>Institute of Global Health</collab>. <article-title>COVID-19 Daily Epidemic Forecasting</article-title> (<year>2021</year>). <comment>Available at: <ext-link ext-link-type="uri" xlink:href="https://renkulab.shinyapps.io/COVID-19-Epidemic-Forecasting/_w_f1413dad/?tab=jhu_pred&amp;country=">https://renkulab.shinyapps.io/COVID-19-Epidemic-Forecasting/_w_f1413dad/?tab&#x3d;jhu_pred&#x26;country&#x3d;</ext-link> (accessed on September 22, 2021)</comment>. </citation>
</ref>
<ref id="B19">
<label>19.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Demongeot</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Flet-Berliac</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Seligmann</surname>
<given-names>H</given-names>
</name>
</person-group>. <article-title>Temperature Decreases Spread Parameters of the New Covid-19 Case Dynamics</article-title>. <source>Biology</source> (<year>2020</year>) <volume>9</volume>:<fpage>94</fpage>. <pub-id pub-id-type="doi">10.3390/biology9050094</pub-id> </citation>
</ref>
<ref id="B20">
<label>20.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Seligmann</surname>
<given-names>H</given-names>
</name>
<name>
<surname>Iggui</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Rachdi</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Vuillerme</surname>
<given-names>N</given-names>
</name>
<name>
<surname>Demongeot</surname>
<given-names>J</given-names>
</name>
</person-group>. <article-title>Inverted Covariate Effects for Mutated 2<sup>nd</sup> vs 1<sup>st</sup> Wave Covid-19: High Temperature Spread Biased for Young</article-title>. <source>Biology (Basel)</source> (<year>2020</year>) <volume>9</volume>:<fpage>226</fpage>. <pub-id pub-id-type="doi">10.3390/biology9080226</pub-id> </citation>
</ref>
<ref id="B21">
<label>21.</label>
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Demongeot</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Noury</surname>
<given-names>N</given-names>
</name>
<name>
<surname>Vuillerme</surname>
<given-names>N</given-names>
</name>
</person-group>. <article-title>Data Fusion for Analysis of Persistence in Pervasive Actimetry of Elderly People at home</article-title>. In: <source>IEEE ARES-CISIS&#x27; 08</source>. <publisher-loc>Piscataway</publisher-loc>: <publisher-name>IEEE Proceedings</publisher-name> (<year>2008</year>). p. <fpage>589</fpage>&#x2013;<lpage>94</lpage>. </citation>
</ref>
<ref id="B22">
<label>22.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Demongeot</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Seligmann</surname>
<given-names>H</given-names>
</name>
</person-group>. <article-title>SARS-CoV-2 and miRNA-like Inhibition Power</article-title>. <source>Med Hypotheses</source> (<year>2020</year>) <volume>144</volume>:<fpage>110245</fpage>. <pub-id pub-id-type="doi">10.1016/j.mehy.2020.110245</pub-id> </citation>
</ref>
<ref id="B23">
<label>23.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Demongeot</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Griette</surname>
<given-names>Q</given-names>
</name>
<name>
<surname>Magal</surname>
<given-names>P</given-names>
</name>
</person-group>. <article-title>Computations of the Transmission Rates in SI Epidemic Model Applied to COVID-19 Data in mainland China</article-title>. <source>R Soc Open Sci</source> (<year>2020</year>) <volume>7</volume>:<fpage>201878</fpage>. <pub-id pub-id-type="doi">10.1098/rsos.201878</pub-id> </citation>
</ref>
<ref id="B24">
<label>24.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Soubeyrand</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Demongeot</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Roques</surname>
<given-names>L</given-names>
</name>
</person-group>. <article-title>Towards Unified and Real-Time Analyses of Outbreaks at Country-Level during Pandemics</article-title>. <source>One Health</source> (<year>2020</year>) <volume>11</volume>:<fpage>100187</fpage>. <pub-id pub-id-type="doi">10.1016/j.onehlt.2020.100187</pub-id> </citation>
</ref>
<ref id="B25">
<label>25.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gaudart</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Landier</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Huiart</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Legendre</surname>
<given-names>E</given-names>
</name>
<name>
<surname>Lehot</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Bendiane</surname>
<given-names>MK</given-names>
</name>
<etal/>
</person-group> <article-title>Factors Associated with the Spatial Heterogeneity of the First Wave of COVID-19 in France: a Nationwide Geo-Epidemiological Study</article-title>. <source>The Lancet Public Health</source> (<year>2021</year>) <volume>6</volume>:<fpage>e222</fpage>&#x2013;<lpage>e231</lpage>. <pub-id pub-id-type="doi">10.1016/s2468-2667(21)00006-2</pub-id> </citation>
</ref>
<ref id="B26">
<label>26.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Oshinubi</surname>
<given-names>K</given-names>
</name>
<name>
<surname>Al-Awadhi</surname>
<given-names>F</given-names>
</name>
<name>
<surname>Rachdi</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Demongeot</surname>
<given-names>J</given-names>
</name>
</person-group>. <article-title>Data Analysis and Forecasting of COVID-19 Pandemic in Kuwait</article-title>. <source>Kuwait J. Sci.</source> (<year>2021</year>), <comment>Special Issue</comment> <fpage>1</fpage>&#x2013;<lpage>28</lpage>. <pub-id pub-id-type="doi">10.48129/kjs.splcov.14501</pub-id> </citation>
</ref>
<ref id="B27">
<label>27.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Griette</surname>
<given-names>Q</given-names>
</name>
<name>
<surname>Demongeot</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Magal</surname>
<given-names>P</given-names>
</name>
</person-group>. <article-title>A Robust Phenomenological Approach to Investigate COVID-19 Data for France</article-title>. <source>Math Appl Sci Eng</source> (<year>2021</year>) <volume>2</volume>:<fpage>149</fpage>&#x2013;<lpage>60</lpage>. <pub-id pub-id-type="doi">10.5206/mase/14031</pub-id> </citation>
</ref>
<ref id="B28">
<label>28.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Oshinubi</surname>
<given-names>K</given-names>
</name>
<name>
<surname>Rachdi</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Demongeot</surname>
<given-names>J</given-names>
</name>
</person-group>. <article-title>Functional Data Analysis: Transition from Daily Observation of COVID-19 Prevalence in France to Functional Curves</article-title>. <source>MedRxiv</source> (<year>2021</year>). <pub-id pub-id-type="doi">10.1101/2021.09.25.21264106</pub-id> </citation>
</ref>
<ref id="B29">
<label>29.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Demongeot</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Oshinubi</surname>
<given-names>K</given-names>
</name>
<name>
<surname>Rachdi</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Hobbad</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Alahiane</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Iggui</surname>
<given-names>S</given-names>
</name>
<etal/>
</person-group> <article-title>The Application of ARIMA Model to Analyse Incidence Pattern in Several Countries</article-title>. <source>J&#x20;Math Comput Sci</source> (<year>2021</year>) <volume>26</volume>:<fpage>41</fpage>&#x2013;<lpage>57</lpage>. <pub-id pub-id-type="doi">10.28919/jmcs/6541</pub-id> </citation>
</ref>
<ref id="B30">
<label>30.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Demongeot</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Oshinubi</surname>
<given-names>K</given-names>
</name>
<name>
<surname>Rachdi</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Seligmann</surname>
<given-names>H</given-names>
</name>
<name>
<surname>Thuderoz</surname>
<given-names>F</given-names>
</name>
<name>
<surname>Waku</surname>
<given-names>J</given-names>
</name>
</person-group>. <article-title>Estimation of Daily Reproduction Numbers during the COVID-19 Outbreak</article-title>. <source>Computation</source> (<year>2021</year>) <volume>9</volume>:<fpage>109</fpage>. <pub-id pub-id-type="doi">10.3390/computation9100109</pub-id> </citation>
</ref>
</ref-list>
</back>
</article>