Next Article in Journal
Potential of Volunteering in Formal and Informal Medical Education—A Theory-Driven Cross-Sectional Study with Example of the COVID-19 Pandemic
Previous Article in Journal
Developing Ecological Thresholds for Nitrogen and Phosphorus in the Haihe River Basin in China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

COVID-19: A Comparative Study of Contagions Peaks in Cities from Europe and the Americas

1
Centro de Investigación y Modelamiento de Fenómenos Aleatorios-Valparaíso, Instituto de Ingeniería Matemática, Facultad de Ingeniería, Universidad de Valparaíso, Valparaíso 2362905, Chile
2
Departamento de Matemáticas, Facultad de Ciencias, Universidad Nacional de Colombia, Bogotá 111321, Colombia
3
Centro de Modelamiento Matemático, Departamento de Ingeniería Matemática, Unité Mixte Internationale, Centre National de la Recherche Scientifique, Universidad de Chile, Santiago 8370456, Chile
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2022, 19(24), 16953; https://doi.org/10.3390/ijerph192416953
Submission received: 2 December 2022 / Accepted: 7 December 2022 / Published: 16 December 2022
(This article belongs to the Section Public Health Statistics and Risk Assessment)

Abstract

:
Coronavirus disease 2019 (COVID-19) is an infectious disease caused by a group of viruses that provoke illnesses ranging from the common cold to more serious illnesses such as pneumonia. COVID-19 started in China and spread rapidly from a single city to an entire country in just 30 days and to the rest of the world in no more than 3 months. Several studies have tried to model the behavior of COVID-19 in diverse regions, based on differential equations of the SIR and stochastic SIR type, and their extensions. In this article, a statistical analysis of daily confirmed COVID-19 cases reported in eleven different cities in Europe and America is conducted. Log-linear models are proposed to model the rise or drop in the number of positive cases reported daily. A classification analysis of the estimated slopes is performed, allowing a comparison of the eleven cities at different epidemic peaks. By rescaling the curves, similar behaviors among rises and drops in different cities are found, independent of socioeconomic conditions, type of quarantine measures taken, whether more or less restrictive. The log-linear model appears to be suitable for modeling the incidence of COVID-19 both in rises and drops.

1. Introduction

Coronavirus disease 2019 (COVID-19) is an infectious disease caused by a group of viruses called coronaviruses that provoke illnesses ranging from the common cold to more serious illnesses such as pneumonia, Middle East respiratory syndrome (MERS) and severe acute respiratory syndrome (SARS). The coronavirus strain (COVID-19) that has caused the outbreak in China is novel and was not previously known; it has been rapidly exported to other countries and has severely affected different aspects of human health worldwide. On 11 march 2020, the World Health Organization declared the COVID-19 outbreak an epidemic. The COVID-19 virus spread rapidly from a single city to an entire country in just 30 days and to the rest of the world in no more than 3 months. The rapid geographic spread and the sudden increase in the number of cases surprised and quickly overwhelmed the health and public health services of several countries. Hence, governments in different countries promoted border closures, social distancing, school closures and other mandatory control measures in order to reduce the spread of the virus.
Several studies have tried to model the behavior of COVID-19 in diverse regions. Different models based on differential equations of the SIR and stochastic SIR type (and their extensions) have been applied (ref. [1]) and statistical comparisons of the COVID-19 dynamic, especially that of the governmental measures behavior, from different countries have been done (see [2,3,4,5]).
Classical susceptible–infectious–removed (SIR) epidemic models have been widely used in the literature to model infectious diseases. These mathematical models are structured according to the number of susceptible (S), infectious (I) and removed (recovered or deceased) individuals (R). The dynamic behavior of these models depends mainly on the basic reproduction number. The disease is controlled and subsequently eliminated if the basic reproduction number is less than one and continues to spread otherwise [6]. However, in practice, infectious diseases may exhibit periodic oscillations and other nonlinear phenomena during the outbreak. Therefore, several studies have attempted to model this nonlinear behavior.
In this article, another direction was taken, by comparing the speed to reach different contagion peaks and, respectively, the decline in different cities. These were chosen because they provided easy access to their data and because of the type of containment measures adopted by their governments. These cities were Bogota, Lima, Lyon, London, Madrid, Mexico City, Marseille, Paris, New York, Rome and Chile Metropolitan Region.
Several articles compare the impact of COVID-19 infections and mortality in different countries. For example, Ref. [7] studied the early impact from February to July 2020 in the Nordic nations of Sweden, Norway, Denmark and Finland using publicly available case data sources. Ref. [8] studied countries from Asia and Europe in terms of cases, death and case fatality rate. Ref. [9] investigated the spreading rate of COVID-19 comparing 10 European countries. In relation to the comparison of cities by COVID-19, Ref. [10] examined the impacts and outcomes of COVID-19 in three cities, New York, London and Tokyo, taking COVID-19-related deaths in these cities and their respective countries as a variable. Some differences were found, possibly explained in part by political, socioeconomic and cultural differences.
In this paper, we wanted to compare 11 cities to each other, in terms of their downward/upward slopes at different contagion peaks. Since these cities had different population sizes, the number of new daily cases were not on the same scale. Moreover, these data were highly nonlinear. Therefore, the data were normalized by applying a natural logarithm function and a log-linear regression model was applied.

2. Material and Methods

This study used different data sources. In Section 2.1, information about the chosen cities and the data studied is presented, while Section 2.2 explains how the data were modeled.

2.1. Dataset

A dataset was created containing daily infected cases of COVID-19 between 1 February 2020 and 30 August 2021 for eleven cities. Note that the data from all the cities were not exactly from the same dates. Table 1 shows the eleven cities and the abbreviations used, as well as the respective repositories.
In the sequel, for each city, i = 1 , , 11 , N i ( t ) denotes the daily number of infected cases of COVID-19 detected in city i at day t, where the series of infected cases detected have been smoothed with a moving average of 7 days to eliminate weekend variations. Figure 1 shows the daily number of infected cases for these cities up to 24 April 2021, which corresponds to the last day in which all cities had information.
As it can be observed in Figure 1 and Figures in Appendix A, depending on the city, the number of infected had one, two or three peaks. Table 2 summarizes the number of observed rises and drops for each city. In fact, most of the cities had at least 2 peaks (2 upward and 2 downward movements) whereas London (LO) and Mexico City (ME) had only 1 peak.
To compare the development of different cities, with diverse control measures, it was necessary to normalize them to the maximum observed at each peak. More specifically, for i = 1 , , 11 and j = 1 , 2 , 3 , A i j was defined as the maximal value of N i for city i and peak j and the day T i j where this maximum was reached. In fact, note that N i ( T i j ) = A i j . For each city i, peak j and day t, denote by ( C N ) i j ( t ) = N i ( t T i j ) / A i j , the centered and normalized data. Figure 2 shows the normalized and centered data for each city and the corresponding peaks. In this way, each curve of cases in different cities and peaks were overlapped on the same graph, so that the peaks coincided at time 0 (time translation). Thus, −50 (days) represented, for each curve, 50 days prior to the corresponding peak. The curves were also normalized in height, so that the peak of each curve was 1. That is, each curve of cases was divided by its own maximum (normalization of cases). In this way, different curves at different locations, which had different sizes, could be compared to see if they responded to a similar mechanism. It can be seen that the normalized curves had similar behaviors even when distinct containment measures were taken.
In Section 2.2, two regression models for the logarithm of these normalized data, centered around their peaks, are proposed. In these models, increasing and decreasing slopes are estimated. Our goal is then to compare the cities in terms of these estimated slopes (see cluster analysis in Section 2.3 and Section 2.4).

2.2. Regression Models

In this section, two regression models are proposed. The first one models the total number of rises and drops. The second one models the rises and drops in a neighborhood of the peak, that is, between the peak and the days when a city exceeds half of the corresponding peak.

2.2.1. The Total Regression Model

A previous step to the definition of the total model is the definition of the beginning of the rises and the end of the drops. Our model for the rises considered the starting point of a rise as the point at which a significant change in the slope was observed with respect to the pseudo-steady state. The same reasoning was applied for the end of a drop. For example, for the city of Bogota, Figure 3 plots the cases between February 2020 and August 2021. Increases are identified with circles and drops with asterisks. Sometimes, an asterisk is superimposed on a circle. Three rises and drops can be isolated. Figure 4 and Figure 5 show the rises and drops of the logarithm of daily reported cases. For the eleven cities, the daily reported cases and the rises and drops can be found in Appendix A.
Denote by R a rise and by D a drop, for the total model (full rise and drop). For a city i and a peak j, T i j R denotes the starting point of a rise and T i j D , the end point of a drop. Recall that T i j has been defined in Section 2.1 and it is the point where the maximum value of the peak is reached. The following model is then considered
log N i ( T i j R + s ) N i ( T i j ) = α i j R + β i j R s + ϵ i , j R ( s ) , s = 0 , , T i j T i j R ;
log N i ( T i j + s ) N i ( T i j ) = α i j D + β i j D s + ϵ i , j D ( s ) , s = 0 , , T i j D T i j ,
where ϵ i , j R ( s ) and ϵ i , j D ( s ) are independent white noises. Models (1) and (2) are log-linear regression models. Estimates for the slope parameters β i j R , β i j D and the intercept parameters α i j R and α i j D were obtained using the classical least square estimators. The estimators of the upward and downward slopes are given in Table 3 together with their goodness of fit.
Note that the starting and end points of the rises and the drops were not entirely well-defined, as can be seen in the pictures of the Appendix A. Therefore, the model could vary significantly due to this starting point. This was one of the reason for considering the 50% model defined in Section 2.2.2, which did not require obtaining these points.

2.2.2. The 50% Regression Model

Denote by r a rise and by d a drop for the 50% model. Moreover, T i j r = max { t < T i j : N i ( t ) 0.5 N i ( T i j ) } is defined as the day that the city i exceeds half of the corresponding peak j. In a similar way, T i j d = min { t > T i j : N i ( t ) 0.5 N i ( T i j ) } is defined as the first day after the peak j that N i falls below half of the peak.
Our goal was to study the relationship between time and the series N i ( T i j r + s ) : s = 0 , ,
T i j T i j r and N i ( T i j + s ) : s = 0 , , T i j d T i j . The following model was considered
log N i ( T i j r + s ) N i ( T i j ) = α i j r + β i j r s + ϵ i , j r ( s ) , s = 0 , , T i j T i j r ;
log N i ( T i j + s ) N i ( T i j ) = α i j d + β i j d s + ϵ i , j d ( s ) , s = 0 , , T i j d T i j ,
where ϵ i , j r and ϵ i , j d are independent white noise processes. In the same way as done previously, least square estimators were calculated for the intercepts and the slopes (see Table 4). For the rises, the highest slope was observed for the first rise of Madrid at 0.11 along | T i j T i j r | = 6 days. The lowest slope was for the second rise of Marseille at 0.00066 during | T i j T i j r | = 82 days, but with an R 2 close to zero. Most of the models had a good R 2 , several of them close to or greater than 0.9 . Only 4 had values of R 2 lower than 0.7 . Those below this value were the second rise of RM, the third rise of MD, the second rise of MR and the first rise of ME. With respect to the drops, the slopes took values between 0.068 (first drop of LY) and 0.0068 (second drop of NY). Only one case had an R 2 < 0.6 , the second drop of MR. The limiting case corresponded to the first drop of Lima with R 2 = 0.56 .
There were some drops, in which the number of reported cases failed to fall below 50% and subsequently restarting a new rise, a phenomena that we called the wadding effect. In the figures given in Appendix A, this can be seen in the following cases: 3 BO goes up (Figure A2); 1 LI goes down (Figure A6); 1 and 2 MR go up (Figure A16); 1 ME goes up (Figure A19); 2 NY goes down (Figure A22); 2 RM goes up (Figure A27); 1 RO goes down (Figure A31).
In the second Marseille descent, two types of descents were observed: initially a slower one and then a steep one. Moreover, the descents were in general more regular than the rises.

2.3. Upward Slopes in the Total and 50% Models

In this section, the slopes are studied in both regression models. Cities were classified according to the upward slopes estimated in the linear regression models. Table 5 gives the estimates for the upward slopes and the number of days | T i j T i j q | with q { r , R } .
From Table 5, it can be seen that in general, the second ascents in the 50% model were slower than or at least equal to the first one with the exception of MD and PA. The cases of RM and NY were distinguished because in the second peak there was a smaller slope in a greater number of days, which could be interpreted as a lower number of infected per day during the second peak than during the first one. Additionally, the slopes in the total model were in general larger than the slopes in the 50% model. The 50th percentile was considered to be a central point where quarantine or other confinement measures began to result in a slope change. Finally, one way to check that the model had a good fit was whether the product Days × Slope was close to log ( 2 ) 0.69 . Notice this quantity should be close to log ( 2 ) , because it represents the logarithmic increase from the 50% to the peak.
In Figure 6, the slopes of the rises in the total and 50% models are plotted. It is clear that the first rises of Madrid (MD) in the total model (0.277) and at 50% (0.115) were atypical (outliers) points in the graph. So the cluster analysis of the upward slopes was performed without the first rises of Madrid.
A cluster Analysis (CA) is one of the data mining techniques that groups the sample observations into classes depending on the essential similarities within a class and the dissimilarities between the different classes found in the dataset. A cluster analysis has been used in the context of the COVID-19 epidemic for comparing countries in terms of mobility and cases (see [11]). In this work, a cluster analysis based on a K-means method and a squared Euclidean distance was performed using Matlab software.
Figure 7 suggests that three was the optimum number of clusters for each model (total and 50%), and Table 6 shows the three clusters and the corresponding means.
The classes were obtained by classifying the different rises by the means from lowest to highest, indicating the speed of growth of the number of infected. A difference was evident when considering the total model and the 50% model slopes. The means of the classes with the totality of the data were higher than the means of the classes with the 50 % data, which could be attributed to a saturation phenomenon. Several rises from Cluster 1 with the full data moved into Clusters 2 or 3 with respect to the 50th percentile data. NY remained in Cluster 3 with the highest data growth rate. With the exception of LO, the Cluster 2 rises with the entire data remained in Cluster 2 of the 50th percentile data.
The 50th percentile model’s rises had three behaviors or classes, namely: slow (Cluster 1), medium (Cluster 2), and fast (Cluster 3). However, in the total model’s rises, there were essentially two behaviors and an outlier rise in NY, probably due to a faulty measurement, which was also evident in the differences of the slopes (see Figure 7).

2.4. Downward Slopes in the Total and 50% Models

In this section, downward slopes are studied in both regression models. Cities were classified according to the downward slopes observed in the log-linear regression models. Table 7 contains the results for the downward slopes.
Observe that the first slopes were larger in absolute value than during the second peak. RM exhibited a similar behavior in the first and second downslope. The same occurred in RM with the number of days. In BO and LI, an increase occurred in the slope of the second peak and the number of days decreased. In the rest of the cities, the slope decreased, and the number of days increased. This increase in the second dip could be interpreted as the measures taken in the second phase not having much effect or being weaker or less respected. The product −Days × Slope in several runs was close to log ( 2 ) , which indicated that the model fitted the data well.
Figure 8 also suggests that three was the optimum number of clusters for both models (total and 50%), and Table 8 shows the three clusters and the corresponding means.
As with the rises, drops were classified (by the K-means method) into three groups with both total and 50% models. Several drops moved from being in Cluster 2 with the total data to Cluster 1 with 50 % of the data. Only MD and LI moved from Cluster 1 with the total data to Cluster 2 with 50 % of the data.

3. Discussion

In this article, a simple statistical analysis of the new coronavirus disease 2019 (COVID-19) outbreak in 11 cities from Europe and America was provided. Using daily infected data in the cities for approximately the months after the first cases were confirmed in each city, downward and rising slopes for each wave were analyzed using a log-linear model approach. These models were consistent with the rises and drops of SIR-type models if one thinks about each of the pandemic waves.
This study was triggered by the question of whether southern hemisphere cities had anticipated measures learned from the experience of northern hemisphere cities that had their first peaks earlier.
We found that the rises and drops were almost independent of the type of quarantine measures taken, whether more or less restrictive. By rescaling the curves, similar behaviors between the rises and drops of the different cities were obtained. For example, there were similarities between cities such as NY and RM, independently of their socioeconomic conditions. The question that remains here is whether this has self-similarity within each city. Of course, the height of the peaks strongly depended on several reasons, namely climatic conditions, confinement measures, testing intensity among many others. However, the evidence provided by the data showed that when the curves of detected cases were rescaled, they had an unexpected regularity and similarity between different cities of the world with disparate socioeconomic conditions, with very dissimilar health systems and where very distinct containment measures were applied.
In this direction, Figure 9 and Figure 10 shows the social distancing and the use of mask in the nine countries considered (http://www.healthdata.org/covid/data-downloads, accessed on 3 February 2022, data from the corresponding cities were used if available; see the explanation of how the data were collected on the web page); in both figures the horizontal axis is the time measured in days (one corresponds to 2 April 2020). Figure 9 shows the evolution of the social distancing, measured through cell phone data, in different geographical areas, showing a higher variability (except for the first peak) than our normalized curves. The zero on the vertical scale represents the social distancing before the pandemic. A negative number represents the reduction in the social distancing. For example, an index of 80 % means a reduction in 80 % . It is interesting to note how this distancing evolved upwards after the first wave, with smaller declines in the periods of the successive peaks. This disparity was also observed in the use of face masks (see Figure 10) in public places for the different geographical areas. The vertical axis represents the percentage of mask use in different countries. To this must be added the high unpredictability of the evolution of COVID-19, as shown in Canals et al. [12], where a comparison of the main Lyapunov exponent was made on the same data series for different countries, obtaining that these varied between 5 and 15, much higher than in other epidemics such as N1H1, where the main coefficient was of the order of 2.
When applying the log-linear regression model to the daily infected cases, the results for the increases in the 11 chosen cities revealed three groups differentiated by the corresponding slopes. For the 50% model, a higher growth rate for group three and a lower one for group one, between 0.001 and 0.007 for the first, 0.017 and 0.026 for the second and 0.029 and 0.036 for the third. When considering the total data, the results showed three groups also differentiated by a higher growth rate for group three and a lower one for group one, between 0.004 and 0.036 for the first one, 0.045 and 0.061 for the second one and the third one containing a single point or outlier with a slope of 0.122.
For the downslope phases of the outbreaks, the results in the 11 chosen cities revealed three groups differentiated by a higher growth rate for group three and lower for group one, between −0.007 and −0.016 for the third, −0.03 and −0.041 for the second and −0.058 and −0.068 for the 50% model. When considering the total data, the results showed three groups also differentiated by a higher growth rate for group three and a lower one for group one, between −0.008 and −0.035 for the third, −0.042 and −0.059 for the second, and the first one containing a single point or outlier with a slope of −0.078 which corresponded to 1 LY.

4. Conclusions

The log-linear regression models studied for the rises and drops showed that different cities had a similar behavior, independent of their socioeconomic situations. On the other hand, for the same city, very different behaviors were observed from one peak to another.
These results show that estimates can only give reasonable indications in the short term and that other variables such as health interventions, the level of distancing and other public policies are necessary in any model that wants to predict in the medium term.
Despite the simplicity of these models, they provide an interesting insight into the COVID-19 statistics in 11 cities in Europe and the Americas, allowing comparisons to be made between them. The log-linear model appears to be suitable for modeling the incidence of COVID-19 both in rises and drops. In addition, the results could be useful for supporting health policy decisions or government interventions. However, they should be used in conjunction with other more complex mathematical and epidemiological models.

Author Contributions

Conceptualization, J.G. and S.T.; Methodology, K.B., J.G., J.S.M. and S.T.; Formal Analysis, J.S.M.; Investigation, K.B., J.G., J.S.M. and S.T.; Writing original draft, K.B., J.G., J.S.M. and S.T. All authors equally contributed to the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

The authors acknowledge the support of Centro de Modelamiento Matemático (CMM) FB210005, BASAL funds for centers of excellence from ANID-Chile. J. Garzón was partially supported by HERMES projects 51145 and 50803. J. San Martín was partially supported by NLHPC project and OPS/OMS CON21-00013967. S. Torres and K. Bertin were partially supported by the Fondecyt regular ANID Grant 1221373. S. Torres was supported in part by MATHAMSUD AMSUD210023 and ECOS210037 projects.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data and models that support the findings of this study are available from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Graphics

In Figure A1, Figure A4, Figure A7, Figure A9, Figure A12, Figure A15, Figure A18, Figure A20, Figure A23, Figure A26 and Figure A29: Increases are identified with circles and drops with asterisks. Daily cases without moving average (red); Daily cases with 7 days moving average (blue).

Appendix A.1. Bogota

Figure A1. Bogota: on the horizontal axis is the date in days and on the vertical axis are the daily cases.
Figure A1. Bogota: on the horizontal axis is the date in days and on the vertical axis are the daily cases.
Ijerph 19 16953 g0a1
Figure A2. Bogota’s rises. Horizontal axis, days; vertical axis, logarithm of daily cases.
Figure A2. Bogota’s rises. Horizontal axis, days; vertical axis, logarithm of daily cases.
Ijerph 19 16953 g0a2
Figure A3. Bogota’s drops. Horizontal axis, days; vertical axis: logarithm of daily cases.
Figure A3. Bogota’s drops. Horizontal axis, days; vertical axis: logarithm of daily cases.
Ijerph 19 16953 g0a3

Appendix A.2. Lima

Figure A4. Lima: on the horizontal axis is the date in days and on the vertical axis are the daily cases.
Figure A4. Lima: on the horizontal axis is the date in days and on the vertical axis are the daily cases.
Ijerph 19 16953 g0a4
Figure A5. Lima’s rises. Horizontal axis, days; vertical axis: logarithm of daily cases.
Figure A5. Lima’s rises. Horizontal axis, days; vertical axis: logarithm of daily cases.
Ijerph 19 16953 g0a5
Figure A6. Lima’s drops. Horizontal axis, days; vertical axis: logarithm of daily cases.
Figure A6. Lima’s drops. Horizontal axis, days; vertical axis: logarithm of daily cases.
Ijerph 19 16953 g0a6

Appendix A.3. Londres

Figure A7. London: on the horizontal axis is the date in days and on the vertical axis are the daily cases.
Figure A7. London: on the horizontal axis is the date in days and on the vertical axis are the daily cases.
Ijerph 19 16953 g0a7
Figure A8. London’s rise and drop. Horizontal axis, days; vertical axis: logarithm of daily cases.
Figure A8. London’s rise and drop. Horizontal axis, days; vertical axis: logarithm of daily cases.
Ijerph 19 16953 g0a8

Appendix A.4. Lyon

Figure A9. Lyon: on the horizontal axis is the date in days and on the vertical axis are the daily cases.
Figure A9. Lyon: on the horizontal axis is the date in days and on the vertical axis are the daily cases.
Ijerph 19 16953 g0a9
Figure A10. Lyon’s rises. Horizontal axis, days; vertical axis: logarithm of daily cases.
Figure A10. Lyon’s rises. Horizontal axis, days; vertical axis: logarithm of daily cases.
Ijerph 19 16953 g0a10
Figure A11. Lyon’s drops. Horizontal axis, days; vertical axis: logarithm of daily cases.
Figure A11. Lyon’s drops. Horizontal axis, days; vertical axis: logarithm of daily cases.
Ijerph 19 16953 g0a11

Appendix A.5. Madrid

Figure A12. Madrid: on the horizontal axis is the date in days and on the vertical axis are the daily cases.
Figure A12. Madrid: on the horizontal axis is the date in days and on the vertical axis are the daily cases.
Ijerph 19 16953 g0a12
Figure A13. Madrid’s rises. Horizontal axis, days; vertical axis: logarithm of daily cases.
Figure A13. Madrid’s rises. Horizontal axis, days; vertical axis: logarithm of daily cases.
Ijerph 19 16953 g0a13
Figure A14. Madrid’s drops. Horizontal axis, days; vertical axis: logarithm of daily cases.
Figure A14. Madrid’s drops. Horizontal axis, days; vertical axis: logarithm of daily cases.
Ijerph 19 16953 g0a14

Appendix A.6. Marseille

Figure A15. Marseille: on the horizontal axis is the date in days and on the vertical axis are the daily cases.
Figure A15. Marseille: on the horizontal axis is the date in days and on the vertical axis are the daily cases.
Ijerph 19 16953 g0a15
Figure A16. Marseille’s rises. Horizontal axis, days; vertical axis: logarithm of daily cases.
Figure A16. Marseille’s rises. Horizontal axis, days; vertical axis: logarithm of daily cases.
Ijerph 19 16953 g0a16
Figure A17. Marseille’s drops. Horizontal axis, days; vertical axis: logarithm of daily cases.
Figure A17. Marseille’s drops. Horizontal axis, days; vertical axis: logarithm of daily cases.
Ijerph 19 16953 g0a17

Appendix A.7. Mexico City

Figure A18. Mexico City: on the horizontal axis is the date in days and on the vertical axis are the daily cases.
Figure A18. Mexico City: on the horizontal axis is the date in days and on the vertical axis are the daily cases.
Ijerph 19 16953 g0a18
Figure A19. Mexico City’s rise and drop. Horizontal axis, days; vertical axis: logarithm of daily cases.
Figure A19. Mexico City’s rise and drop. Horizontal axis, days; vertical axis: logarithm of daily cases.
Ijerph 19 16953 g0a19

Appendix A.8. New York

Figure A20. New York: on the horizontal axis is the date in days and on the vertical axis are the daily cases.
Figure A20. New York: on the horizontal axis is the date in days and on the vertical axis are the daily cases.
Ijerph 19 16953 g0a20
Figure A21. New York’s rises. Horizontal axis, days; vertical axis: logarithm of daily cases.
Figure A21. New York’s rises. Horizontal axis, days; vertical axis: logarithm of daily cases.
Ijerph 19 16953 g0a21
Figure A22. New York’s drops. Horizontal axis, days; vertical axis: logarithm of daily cases.
Figure A22. New York’s drops. Horizontal axis, days; vertical axis: logarithm of daily cases.
Ijerph 19 16953 g0a22

Appendix A.9. Paris

Figure A23. Paris: on the horizontal axis is the date in days and on the vertical axis are the daily cases.
Figure A23. Paris: on the horizontal axis is the date in days and on the vertical axis are the daily cases.
Ijerph 19 16953 g0a23
Figure A24. Paris’s rises. Horizontal axis, days; vertical axis: logarithm of daily cases.
Figure A24. Paris’s rises. Horizontal axis, days; vertical axis: logarithm of daily cases.
Ijerph 19 16953 g0a24
Figure A25. Paris’s drops. Horizontal axis, days; vertical axis: logarithm of daily cases.
Figure A25. Paris’s drops. Horizontal axis, days; vertical axis: logarithm of daily cases.
Ijerph 19 16953 g0a25

Appendix A.10. Chile Metropolitan Region

Figure A26. Chile Metropolitan Region: on the horizontal axis is the date in days and on the vertical axis are the daily cases.
Figure A26. Chile Metropolitan Region: on the horizontal axis is the date in days and on the vertical axis are the daily cases.
Ijerph 19 16953 g0a26
Figure A27. Chile Metropolitan Region’s rises. Horizontal axis, days; vertical axis: logarithm of daily cases.
Figure A27. Chile Metropolitan Region’s rises. Horizontal axis, days; vertical axis: logarithm of daily cases.
Ijerph 19 16953 g0a27
Figure A28. Chile Metropolitan Region’s drops. Horizontal axis, days; vertical axis: logarithm of daily cases.
Figure A28. Chile Metropolitan Region’s drops. Horizontal axis, days; vertical axis: logarithm of daily cases.
Ijerph 19 16953 g0a28

Appendix A.11. Roma

Figure A29. Roma: on the horizontal axis is the date in days and on the vertical axis are the daily cases.
Figure A29. Roma: on the horizontal axis is the date in days and on the vertical axis are the daily cases.
Ijerph 19 16953 g0a29
Figure A30. Roma’s rises. Horizontal axis, days; vertical axis: logarithm of daily cases.
Figure A30. Roma’s rises. Horizontal axis, days; vertical axis: logarithm of daily cases.
Ijerph 19 16953 g0a30
Figure A31. Roma’s drops. Horizontal axis, days; vertical axis: logarithm of daily cases.
Figure A31. Roma’s drops. Horizontal axis, days; vertical axis: logarithm of daily cases.
Ijerph 19 16953 g0a31

References

  1. Nesteruk, I. Statistics-based predictions of coronavirus epidemic spreading in mainland China. Innov. Biosyst. Bioeng. 2020, 4, 13–18. [Google Scholar] [CrossRef] [Green Version]
  2. Chu, J. A statistical analysis of the novel coronavirus (COVID-19) in Italy and Spain. PLoS ONE 2021, 16, e0249037. [Google Scholar] [CrossRef]
  3. Mishra, S.; Scott, J.A.; Laydon, D.J.; Flaxman, S.; Gandy, A.; Mellan, T.A.; Unwin, H.J.T.; Vollmer, M.; Coupl, H.; Ratmann, O.; et al. Comparing the responses of the UK, Sweden and Denmark to COVID-19 using counterfactual modelling. Sci. Rep. 2021, 11, 16342. [Google Scholar] [CrossRef]
  4. Rozanova, L.; Temerev, A.; Flahault, A. Comparing the Scope and Efficacy of COVID-19 Response Strategies in 16 Countries: An Overview. Int. J. Environ. Res. Public Health 2020, 17, 9421. [Google Scholar] [CrossRef]
  5. Zuo, M.; Saima, K.K.; Zubair, A.; Zahra, A. Comparison of COVID-19 Pandemic Dynamics in Asian Countries with Statistical Modeling. Comput. Math. Methods Med. 2020, 2020, 4296806. [Google Scholar] [CrossRef]
  6. Brauer, F.; Castillo-Chavez, C. Mathematical Models in Population Biology and Epidemiology; Springer: New York, NY, USA, 2012; Volume 2. [Google Scholar]
  7. Yarmol-Matusiak, E.A.; Cipriano, L.E.; Stranges, S. A comparison of COVID-19 epidemiological indicators in Sweden, Norway, Denmark, and Finland. Comp. Study Scand. J. Public Health 2021, 49, 69–78. [Google Scholar] [CrossRef]
  8. Dehkordi, A.H.; Alizadeh, M.; Derakhshan, P.; Babazadeh, P.; Jahandideh, A. Understanding epidemic data and statistics: A case study of COVID-19. J. Med. Virol. 2020, 92, 868–882. [Google Scholar] [CrossRef] [Green Version]
  9. De Leo, S.; Maia, G.; Solidoro, L. Analysing and comparing the COVID-19 data: The closed cases of Hubei and South Korea, the dark March in Europe, the beginning of the outbreak in South America. MedRXiv 2020. [Google Scholar] [CrossRef] [Green Version]
  10. Marcotullio, P.J.; Schmeltz, M. COVID-19 in Three Global Cities: Comparing Impacts and Outcomes. J. Extreme Events 2021, 7, 2150013. [Google Scholar] [CrossRef]
  11. Vitello, P.; Fiandrino, C.; Capponi, A.; Klopp, P.; Connors, R.D.; Viti, F. The Impact of SARS-COVID-19 Outbreak on European Cities Urban Mobility. Front. Future Transp. 2021, 2, 666212. [Google Scholar] [CrossRef]
  12. Canals, M.; Maass, A.; Monreal, J. Non-linear analysis of sequences of SARS-CoV-2 reveals high unpredictability. 2021; submitted. [Google Scholar]
Figure 1. Daily changes of coronavirus disease 2019 cases in different cities during 2020–2021.
Figure 1. Daily changes of coronavirus disease 2019 cases in different cities during 2020–2021.
Ijerph 19 16953 g001
Figure 2. Normalization of reported case data.
Figure 2. Normalization of reported case data.
Ijerph 19 16953 g002
Figure 3. Bogota: on the horizontal axis is the date in days and on the vertical axis are the daily cases. Increases are identified with circles and drops with asterisks. Daily cases without moving average (red). Daily cases with 7 days moving average (blue).
Figure 3. Bogota: on the horizontal axis is the date in days and on the vertical axis are the daily cases. Increases are identified with circles and drops with asterisks. Daily cases without moving average (red). Daily cases with 7 days moving average (blue).
Ijerph 19 16953 g003
Figure 4. Bogota’s rises. Horizontal axis, days; vertical axis, logarithm of daily cases.
Figure 4. Bogota’s rises. Horizontal axis, days; vertical axis, logarithm of daily cases.
Ijerph 19 16953 g004
Figure 5. Bogota’s drops. Horizontal axis, days; vertical axis, logarithm of daily cases.
Figure 5. Bogota’s drops. Horizontal axis, days; vertical axis, logarithm of daily cases.
Ijerph 19 16953 g005
Figure 6. Upward slopes: blue: total, red: 50%.
Figure 6. Upward slopes: blue: total, red: 50%.
Ijerph 19 16953 g006
Figure 7. Cluster analysis for rises: blue: total, red: 50%, (without MD).
Figure 7. Cluster analysis for rises: blue: total, red: 50%, (without MD).
Ijerph 19 16953 g007
Figure 8. Downward slope: blue: total, red: 50%.
Figure 8. Downward slope: blue: total, red: 50%.
Ijerph 19 16953 g008
Figure 9. Social distancing in different countries’ cities between 4 April 2020 and 16 December 2021.
Figure 9. Social distancing in different countries’ cities between 4 April 2020 and 16 December 2021.
Ijerph 19 16953 g009
Figure 10. Use of mask in different countries’ cities between 4 April 2020 and 17 December 2021.
Figure 10. Use of mask in different countries’ cities between 4 April 2020 and 17 December 2021.
Ijerph 19 16953 g010
Table 1. List of chosen cities.
Table 2. Numbers of rises and drops for the eleven cities.
Table 2. Numbers of rises and drops for the eleven cities.
CityNo. IncreasesNo. Decreases
BO33
LI32
LO11
LY22
MD32
MR22
ME11
NY22
PA22
RM22
RO22
Table 3. Estimation in the total regression model.
Table 3. Estimation in the total regression model.
CityPeakc = R, D | T ij T ij c | α ij c p-Value β ij c p-Value R 2
BO1R69−2.05 6.2 × 10 59 0.0336 1.1 × 10 46 0.95
BO2R64−1.26 5.5 × 10 67 0.0206 3.6 × 10 53 0.98
BO3R103−2.39 2.7 × 10 71 0.0271 8.8 × 10 55 0.91
LI1R57−2.28 1.5 × 10 36 0.045 5.9 × 10 27 0.88
LI2R36−0.82 1.3 × 10 34 0.0253 6.4 × 10 28 0.97
LI3R47−1.91 5.5 × 10 36 0.047 5.9 × 10 28 0.93
LO1R20−1.95 7.1 × 10 30 0.061 1.1 × 10 23 0.95
LY1R77−2.81 1.2 × 10 66 0.036 7.1 × 10 49 0.95
LY2R46−1.13 3.4 × 10 31 0.023 2.2 × 10 20 0.86
MD1R11−1.81 1.5 × 10 04 0.270.0010.81
MD2R63−2.56 6.9 × 10 38 0.046 2.0 × 10 27 0.86
MD3R45−1.58 6.0 × 10 29 0.035 2.7 × 10 19 0.85
MR1R92−2.25 4.3 × 10 52 0.023 1.9 × 10 31 0.78
MR2R97−0.73 1.0 × 10 31 0.004 5.4 × 10 07 0.23
ME1R114−2.01 2.7 × 10 89 0.017 9.5 × 10 63 0.92
NY1R28−2.55 1.5 × 10 11 0.12 5.5 × 10 09 0.74
NY2R73−1.84 4.8 × 10 62 0.026 2.4 × 10 46 0.95
PA1R23−2.77 6.0 × 10 77 0.032 2.5 × 10 58 0.95
PA2R92−1.69 6.0 × 10 95 0.126 1.0 × 10 61 0.92
RM1R54−2.62 6.7 × 10 39 0.058 1.6 × 10 30 0.92
RM2R112−1.13 9.0 × 10 50 0.011 1.5 × 10 30 0.7
RO1R59−2.92 1.4 × 10 48 0.058 7.2 × 10 39 0.95
RO2R28−0.65 5.1 × 10 26 0.0257 8.0 × 10 21 0.97
BO1D46−0.12 1.5 × 10 06 −0.0225 2.2 × 10 28 0.94
BO2D54−0.16 8.2 × 10 04 −0.043 4.8 × 10 34 0.94
BO3D44−0.26 4.2 × 10 09 −0.045 4.9 × 10 31 0.96
LI1D39−0.09 1.5 × 10 02 −0.025 2.1 × 10 18 0.88
LI2D29−0.13 1.0 × 10 02 −0.057 1.1 × 10 17 0.93
LO1D600.13 1.3 × 10 10 −0.586 3.7 × 10 79 1
LY1D29−0.17 4.2 × 10 04 −0.078 4.3 × 10 22 0.97
LY2D610.08 2.0 × 10 03 −0.0275 3.67 × 10 42 0.96
MD1D48−0.03 1.8 × 10 01 −0.0581 2.0 × 10 52 0.99
MD2D79−0.32 7.0 × 10 16 −0.0183 6.0 × 10 40 0.90
MD3D55−0.07 3.1 × 10 02 −0.033 1.7 × 10 36 0.95
MR1D34−0.05 1.9 × 10 01 −0.053 1.3 × 10 23 0.96
MR2D63−0.26 5.1 × 10 07 −0.035 9.3 × 10 36 0.92
ME1D97−0.11 1.5 × 10 05 −0.0195 4.5 × 10 66 0.96
NY1D640.00 5.5 × 10 107 −0.0419 1.1 × 10 60 0.99
NY2D101−0.13 5.6 × 10 09 −0.0077 3.8 × 10 38 0.81
PA1D41−0.29 5.7 × 10 01 −0.057 4.2 × 10 26 0.95
PA2D690.20 4.3 × 10 10 −0.0268 1.0 × 10 47 0.96
RM1D103−0.35 1.4 × 10 11 −0.0242 8.4 × 10 54 0.91
RM2D580.040.1−0.0423 2.5 × 10 49 0.98
RO1D95−0.22 2.5 × 10 11 −0.0093 1.4 × 10 30 0.76
RO2D107−0.54 2.1 × 10 18 −0.0333 1.2 × 10 65 0.94
Table 4. Estimation in the 50% regression model.
Table 4. Estimation in the 50% regression model.
CityPeakc = r, d | T ij T ij c | α ij c p-Value β ij c p-Value R 2
BO1r29−0.47 7.2 × 10 22 0.017 8.3 × 10 16 0.91
BO2r35−0.58 2.4 × 10 30 0.017 1.6 × 10 22 0.95
BO3r44−0.72 1.4 × 10 25 0.02 1.6 × 10 19 0.86
LI1r28−0.70 4.6 × 10 22 0.023 3.7 × 10 15 0.91
LI2r29−0.61 2.4 × 10 25 0.023 5.8 × 10 20 0.96
LI3r25−0.53 1.2 × 10 18 0.021 5.2 × 10 13 0.9
LO1r20−0.58 1.7 × 10 15 0.031 1.4 × 10 11 0.93
LY1r19−0.60 1.0 × 10 10 0.036 1.1 × 10 07 0.82
LY2r19−0.55 7.3 × 10 14 0.034 7.0 × 10 11 0.92
MD1r6−0.50 1.9 × 10 03 0.110.0070.87
MD2r36−0.85 1.5 × 10 20 0.02 2.1 × 10 11 0.74
MD3r19−0.49 1.2 × 10 07 0.022 6.6 × 10 04 0.50
MR1r15−0.41 1.8 × 10 06 0.031 1.7 × 10 04 0.67
MR2r82−0.50 1.1 × 10 21 0.00060.450.07
ME1r43−0.50 2.8 × 10 15 0.007 5.9 × 10 05 0.33
NY1r18−0.52 2.8 × 10 11 0.036 6.5 × 10 09 0.89
NY2r36−0.76 1.6 × 10 27 0.019 7.8 × 10 18 0.89
PA1r23−0.55 3.8 × 10 16 0.023 9.9 × 10 11 0.87
PA2r36−0.81 1.0 × 10 10 0.029 9.3 × 10 17 0.93
RM1r25−0.48 4.4 × 10 15 0.018 1.3 × 10 09 0.80
RM2r84−0.49 2.9 × 10 30 0.0039 1.1 × 10 09 0.37
RO1r25−0.51 3.5 × 10 15 0.025 8.2 × 10 12 0.87
RO2r28−0.64 5.1 × 10 26 0.026 8.0 × 10 21 0.97
BO1d23−0.01 3.4 × 10 01 −0.033 4.4 × 10 19 0.98
BO2d150.01 4.7 × 10 01 −0.048 2.5 × 10 11 0.97
BO3d260.10 2.8 × 10 04 −0.031 1.0 × 10 15 0.94
LI1d24−0.14 1.4 × 10 02 −0.02 2.9 × 10 05 0.56
LI2d18−0.02 5.6 × 10 01 −0.034 1.0 × 10 08 0.88
LO1d170.02 1.2 × 10 01 −0.045 1.7 × 10 14 0.98
LY1d130.14 5.5 × 10 02 −0.068 1.6 × 10 05 0.83
LY2d31−0.02 4.7 × 10 01 −0.019 6.1 × 10 13 0.84
MD1d13−0.00 8.8 × 10 01 −0.063 2.4 × 10 09 0.97
MD2d17−0.12 2.1 × 10 02 −0.038 9.3 × 10 07 0.81
MD3d190.0 9.9 × 10 01 −0.035 1.2 × 10 08 0.86
MR1d120.05 2.9 × 10 01 −0.062 1.5 × 10 06 0.91
MR2d310.00 9.8 × 10 01 −0.014 1.5 × 10 06 0.56
ME1d240.07 2.8 × 10 02 −0.03 8.8 × 10 12 0.88
NY1d190.03 9.2 × 10 02 −0.04 3.0 × 10 15 0.98
NY2d84−0.16 5.3 × 10 11 −0.0068 6.4 × 10 26 0.74
PA1d130.08 9.0 × 10 02 −0.058 8.9 × 10 07 0.90
PA2d340.03 1.7 × 10 01 −0.016 3.8 × 10 14 0.8
RM1d200.08 2.0 × 10 04 −0.041 9.1 × 10 16 0.97
RM2d190.07 4.5 × 10 04 −0.040 1.0 × 10 14 0.97
RO1d260.08 9.9 × 10 07 −0.03 1.0 × 10 22 0.98
RO2d510.07 1.6 × 10 06 −0.0147 3.8 × 10 35 0.96
Table 5. Upward slopes for total and 50% models.
Table 5. Upward slopes for total and 50% models.
Total50%
CityPeakDaysSlopeDaysSlopeDays × Slope
BO1680.034290.0170.49
BO2630.021350.0170.60
BO31000.027440.020.88
LI1560.045280.0230.64
LI2350.025290.0230.67
LI3460.047250.0210.53
LO1360.061200.0310.62
LY1760.036190.0360.68
LY2450.023190.0340.65
MD180.27760.1150.69
MD2620.046360.020.72
MD3440.035190.0220.42
MR1910.023150.0310.47
MR2960.004820.0010.08
ME11130.017430.0070.30
NY1270.122180.0360.65
NY2720.026360.0190.68
PA1910.032230.0230.53
PA21130.014290.0290.84
RM1530.058250.0180.45
RM21110.011840.0040.34
RO1580.058250.0250.63
RO2270.026280.0260.73
Table 6. Cluster analysis for rises at total and 50%, without the first rise of MD.
Table 6. Cluster analysis for rises at total and 50%, without the first rise of MD.
50%Total
CityPeakClassMeanSlopeClassMeanSlope
RM210.0040.00410.0220.011
MR210.0040.00110.0220.004
ME110.0040.00710.0220.017
RM120.0210.01820.0480.058
NY220.0210.01910.0220.026
BO120.0210.01710.0220.034
BO220.0210.01710.0220.021
BO320.0210.0210.0220.027
RO120.0210.02520.0480.058
RO220.0210.02610.0220.026
LI120.0210.02320.0480.045
LI220.0210.02310.0220.025
LI320.0210.02120.0480.047
MD220.0210.0220.0480.046
MD320.0210.02210.0220.035
PA120.0210.02310.0220.032
NY130.0330.03630.1220.122
LO130.0330.03120.0480.061
LY130.0330.03610.0220.036
LY230.0330.03410.0220.023
MR130.0330.03110.0220.023
PA230.0330.02910.0220.013
Table 7. Downward slopes for total and 50% models.
Table 7. Downward slopes for total and 50% models.
Total50%
CityPeakDaysSlopeDaysSlope−Days × Slope
BO146−0.02323−0.0330.76
BO254−0.04315−0.0480.72
BO344−0.04526−0.0310.81
LI138−0.02524−0.020.48
LI228−0.05718−0.0340.61
LO159−0.05917−0.0450.77
LY128−0.07813−0.0680.88
LY260−0.02731−0.0190.59
MD147−0.05813−0.0630.82
MD278−0.01817−0.0380.65
MD354−0.03319−0.0350.67
MR133−0.05312−0.0620.74
MR262−0.03531−0.0140.43
ME196−0.0224−0.030.72
NY163−0.04219−0.040.76
NY2100−0.00884−0.0070.59
PA140−0.05713−0.0580.75
PA268−0.02734−0.0160.54
RM1102−0.02420−0.0410.82
RM257−0.04219−0.040.76
RO194−0.00926−0.030.78
RO2106−0.03351−0.0150.77
Table 8. Cluster analysis for drops at total and 50%.
Table 8. Cluster analysis for drops at total and 50%.
50%Total
CityPeakClassMeanSlopeClassMeanSlope
LY11−0.063−0.0681−0.078−0.078
MD11−0.063−0.0632−0.051−0.058
MR11−0.063−0.0622−0.051−0.053
PA11−0.063−0.0582−0.051−0.057
RM12−0.037−0.0413−0.023−0.024
RM22−0.037−0.042−0.051−0.042
NY12−0.037−0.042−0.051−0.042
BO12−0.037−0.0333−0.023−0.023
BO22−0.037−0.0482−0.051−0.043
BO32−0.037−0.0312−0.051−0.045
RO12−0.037−0.033−0.023−0.009
LI22−0.037−0.0342−0.051−0.057
LO12−0.037−0.0452−0.051−0.059
MD22−0.037−0.0383−0.023−0.018
MD32−0.037−0.0353−0.023−0.033
ME12−0.037−0.033−0.023−0.02
NY23−0.015−0.0073−0.023−0.008
RO23−0.015−0.0153−0.023−0.033
LI13−0.015−0.023−0.023−0.025
LY23−0.015−0.0193−0.023−0.027
MR23−0.015−0.0143−0.023−0.035
PA23−0.015−0.0163−0.023−0.027
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bertin, K.; Garzón, J.; San Martín, J.; Torres, S. COVID-19: A Comparative Study of Contagions Peaks in Cities from Europe and the Americas. Int. J. Environ. Res. Public Health 2022, 19, 16953. https://doi.org/10.3390/ijerph192416953

AMA Style

Bertin K, Garzón J, San Martín J, Torres S. COVID-19: A Comparative Study of Contagions Peaks in Cities from Europe and the Americas. International Journal of Environmental Research and Public Health. 2022; 19(24):16953. https://doi.org/10.3390/ijerph192416953

Chicago/Turabian Style

Bertin, Karine, Johanna Garzón, Jaime San Martín, and Soledad Torres. 2022. "COVID-19: A Comparative Study of Contagions Peaks in Cities from Europe and the Americas" International Journal of Environmental Research and Public Health 19, no. 24: 16953. https://doi.org/10.3390/ijerph192416953

APA Style

Bertin, K., Garzón, J., San Martín, J., & Torres, S. (2022). COVID-19: A Comparative Study of Contagions Peaks in Cities from Europe and the Americas. International Journal of Environmental Research and Public Health, 19(24), 16953. https://doi.org/10.3390/ijerph192416953

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop