Next Article in Journal
Valley Fever: Environmental Risk Factors and Exposure Pathways Deduced from Field Measurements in California
Next Article in Special Issue
Preparedness and Response to the COVID-19 Emergency: Experience from the Teaching Hospital of Pisa, Italy
Previous Article in Journal
Promoting Evidence-Based Practice for Improved Occupational Safety and Health at Workplaces in Sweden. Report on a Practice-Based Research Network Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Clustering Approach to Classify Italian Regions and Provinces Based on Prevalence and Trend of SARS-CoV-2 Cases

1
Department of Medical and Surgical Sciences and Advanced Technologies “GF Ingrassia”, University of Catania, 95123 Catania, Italy
2
Azienda Ospedaliero-Universitaria “Policlinico-Vittorio Emanuele”, 95123 Catania, Italy
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2020, 17(15), 5286; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17155286
Submission received: 8 June 2020 / Revised: 10 July 2020 / Accepted: 18 July 2020 / Published: 22 July 2020
(This article belongs to the Special Issue New and Re-emerging Pathogens)

Abstract

:
While several efforts have been made to control the epidemic of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in Italy, differences between and within regions have made it difficult to plan the phase two management after the national lockdown. Here, we propose a simple and immediate clustering approach to categorize Italian regions working on the prevalence and trend of SARS-CoV-2 positive cases prior to the start of phase two on 4 May 2020. Applying both hierarchical and k-means clustering, we identified three regional groups: regions in cluster 1 exhibited higher prevalence and the highest trend of SARS-CoV-2 positive cases; those classified into cluster 2 constituted an intermediate group; those in cluster 3 were regions with a lower prevalence and the lowest trend of SARS-CoV-2 positive cases. At the provincial level, we used a similar approach but working on the prevalence and trend of the total SARS-CoV-2 cases. Notably, provinces in cluster 1 exhibited the highest prevalence and trend of SARS-CoV-2 cases. Provinces in clusters 2 and 3, instead, showed a median prevalence of approximately 11 cases per 10,000 residents. However, provinces in cluster 3 were those with the lowest trend of cases. K-means clustering yielded to an alternative cluster solution in terms of the prevalence and trend of SARS-CoV-2 cases. Our study described a simple and immediate approach to monitor the SARS-CoV-2 epidemic at the regional and provincial level. These findings, at present, offered a snapshot of the epidemic, which could be helpful to outline the hierarchy of needs at the subnational level. However, the integration of our approach with further indicators and characteristics could improve our findings, also allowing the application to different contexts and with additional aims.

1. Introduction

The novel coronavirus (SARS-CoV-2) epidemic began to spread in the Hubei province (China) at the end of 2019, and then to more than 200 countries worldwide [1]. At the end of February 2020, two distinct outbreaks of SARS-CoV-2 occurred in two small Italian areas within the Lombardy and Veneto regions [2]. The epidemic, since then, spread at different times and with different intensities across all the Italian regions [3]. On 10 March 2020, the Italian government has promptly reacted to the epidemic by adopting a first set of national restrictions and recommendations (e.g., travel restrictions, quarantine and contact precautions) [4,5], which were strengthened on 23 March by avoiding non-essential industrial productions and social interactions [4,5]. While these efforts have been made to control the epidemic [4,5], differences between and within regions have made it difficult to plan the phase two management after the national lockdown. As a consequence of the control measures adopted, there has been a 5.3% decrease in SARS-CoV-2 positive patients during the week preceding the beginning of phase two in Italy, 4 May 2020 [3]. This resulted in a total of 100,179 active cases throughout the Italian territory on 3 May 2020 [3]. However, these patients were not evenly distributed across Italian regions, with three regions (i.e., Lombardy, Piedmont, and Emilia-Romagna) that presented more than 60% of the total active cases [3]. Several modelling studies evaluated the efficacy of control measures adopted in Italy and predicted future scenarios to help policymakers in designing the phase two strategy against the SARS-CoV-2 epidemic [6,7,8,9,10]. However, identifying different regional and provincial clusters would be just as crucial to plan strategies for relaxing the restrictions in accordance with specific situations and needs. As for other emerging diseases, useful and innovative tools could be able to provide a snapshot of the epidemic at the regional and provincial level. Here, we used two of the most common clustering methods (i.e., hierarchical clustering and k-means algorithm) [11] to categorize Italian regions and provinces into several groups. At the regional level, this clustering approach works directly on the prevalence and trend of active cases in each region, prior to the start of phase two on 4 May 2020. Specifically, it has been applied to regional data on the prevalence of SARS-CoV-2 positive cases on 3 May 2020 and their trend in the previous week from the 27 April to the 3 May 2020. At the provincial level, we applied a similar approach, but working on the prevalence and trend of total SARS-CoV-2 cases, due to the limited availability of data on positive cases. Thus, our aim was to apply a simple clustering approach in order to classify Italian regions and provinces based on the prevalence and trend of SARS-CoV-2 cases prior to the start of phase two on 4 May 2020.

2. Materials and Methods

We first obtained data on SARS-CoV-2 positive cases and on the number of tests performed, which were collected and released by Italy’s Civil Protection of the Italian Ministry of Health from 27 April to 3 May 2020. These data, deposited on GitHub [12,13], were freely available and provided at the regional and provincial level. On the dedicated website [14], it is also possible to consult the daily prevalence of SARS-CoV-2 cases, their weekly trend and other information on the epidemic in Italy. We also obtained the number of residents of each region and province from the Italian National Institute of Statistics (ISTAT, Istituto Nazionale di Statistica) [15]. Thus, our study did not contain data at the individual level that would require ethical approval. Using the aforementioned data, we calculated the following regional and provincial indicators:
  • The regional prevalence of SARS-CoV-2 positive cases on 3 May 2020 (i.e., expressed as the number of positive cases per 10,000 residents);
  • The weekly regional trend of SARS-CoV-2 positive cases from 27 April to 3 May 2020 (expressed as the percentage of increment/decrement of positive cases);
  • The provincial prevalence of SARS-CoV-2 cases on 3 May 2020 (i.e., expressed as the number of total cases per 10,000 residents);
  • The weekly provincial trend of SARS-CoV-2 cases from 27 April to 3 May 2020 (expressed as percentage of increment of total cases);
  • The number of tests performed per 10,000 residents (i.e., only at the regional level).
We then inspected the structure of the regional dataset by a three-dimensional scatter plot and tested for normality using the Kolmogorov–Smirnov test. According to the skewed distribution of the aforementioned indicators, we tested for correlations using the Spearman’s rank correlation coefficient. We separately tested the linear relationship of the number of tests performed per 10,000 residents (independent variable) with a prevalence or trend of SARS-CoV-2 cases (dependent variables) using linear regression models. Variables included in the regression models were log-transformed before being analyzed and the results were reported as β and its standard error (SE). Since data on the prevalence and trend of SARS-CoV-2 positive cases were related to different testing strategies between Italian regions, we then used the residuals of linear regressions to adjust both indicators for the number of tests performed per 10,000 residents. To account for different scales between the indicators, the adjusted values were further standardized using the Z-Score formula: z = (x − μ)/σ, where x was the log-transformed indicator, μ was its national log-transformed mean, and σ was the standard deviation. Working on these values (i.e., the log-transformed prevalence and trend of SARS-CoV-2 positive cases), we then applied hierarchical clustering using Ward’s criterion. The distance between regions—and hence their similarity—was determined by calculating their squared Euclidean distance in the two-dimensional plane. We also applied the Euclidean distance but obtaining the same cluster solution (data not shown). This clustering was shown by a hierarchical tree—also called dendrogram—in which the height of the branches indicated the distance between the clusters [11]. The optimal number of clusters—k—was determined using the silhouette method. Specifically, we derived four clustering solutions, which differed for the number of clusters k, into a range from two to five clusters (2 ≤ k ≤ 5). For each clustering solution, we calculated the silhouette value that represented how each region was cohered to its own cluster and separated from the others. We chose the optimal clustering solution that maximized the silhouette value [16]. Since slight differences in the clustering solution could be obtained using different methods [11], we consolidated the hierarchical clustering by applying the k-means algorithm. This method partitioned n observations into k clusters by minimizing intra-cluster variation [11]. In contrast to hierarchical clustering, the k-means clustering required a predefined number of clusters, which was determined according to the silhouette method applied to hierarchical clustering. Finally, we evaluated the inter-cluster variability using the Kruskal–Wallis test.
The provincial dataset on the prevalence and trend of total SARS-CoV-2 cases was analyzed in the same way, except for the preliminary adjustment for the number of tests performed. In fact, there were no accurate and complete data on the number of tests performed at provincial level. In brief, the prevalence and trend of total SARS-CoV-2 cases at the provincial level were first log-transformed and standardized using the Z-Score formula. Then, we used hierarchical clustering based on the squared Euclidean distance and Ward’s criterion to derive four clustering solutions (2 ≤ k ≤ 5). The optimal clustering solution was chosen to maximize the silhouette value [16]. Finally, we consolidated the hierarchical clustering solution by applying the k-means algorithm with a predefined number of clusters. All the analyses were performed on the SPSS software (version 23.0, SPSS, Chicago, IL, USA), with a significance level α of 0.05.

3. Results

3.1. Description of Data

The prevalence of SARS-CoV-2 positive cases on 3 May 2020, their trend in the week before, and the number of tests performed in each region are reported in Table 1.
In Figure 1A, we show how Italian regions were distributed in a three-dimensional scatter plot of these indicators. Although this plot showed no clear separation between the regions, we found that six and eight regions were above the national average for the prevalence and trend of SARS-CoV-2 positive cases, respectively. In contrast, 10 regions were below the national average for the number of tests performed per 10,000 residents. We also found a positive correlation between the number of tests performed and the prevalence of SARS-CoV-2 positive cases on 3 May 2020 (Spearman’s correlation coefficient = 0.527; p = 0.015). In contrast, we noted a significant negative correlation between the number of tests performed and the trend of SARS-CoV-2 positive cases in the week from 27 April to 3 May 2020 (Spearman’s correlation coefficient = −0.613; p = 0.003; Figure 1B). Specifically, we found a positive linear relationship between the log-transformed number of tests performed and the log-transformed prevalence of SARS-CoV-2 positive cases (β = 0.786; SE = 0.292; p = 0.015; Figure 2A). By contrast, the log-transformed trend of SARS-CoV-2 positive cases was negatively associated with the log-transformed number of tests performed (β = −0.201; SE = 0.067; p = 0.007; Figure 2B).

3.2. Regional Clusters

We then applied hierarchical clustering, yielding to the dendrogram depicted in Figure 3. After the visual inspection of the hierarchical tree and based on silhouette values, we decided to partition Italian regions into three clusters. Indeed, the clustering solution with a number of three clusters was that with the highest average silhouette width (Table 2).
According to this clustering solution, Lombardy, Piedmont, Liguria and Marche tended to be grouped into cluster 1, Umbria and Valle d’Aosta into cluster 3, and the remaining regions into cluster 2 (Figure 4). Differences between the clusters, in terms of the prevalence and trend of SARS-CoV-2 positive cases, are reported in Table 3. Specifically, the regions in cluster 1 exhibited the highest prevalence and trend of SARS-CoV-2 positive cases. In contrast, the regions in cluster 3 were those with the lowest prevalence and trend of SARS-CoV-2 positive cases. Regions in cluster 2 were those with an intermediate prevalence and trend of SARS-CoV-2 positive cases. We then consolidated this partition by applying a k-means algorithm with a predefined number of three clusters. However, all the regions maintained the same allocation of hierarchical clustering.

3.3. Provincial Clusters

We then applied a similar approach to identify different clusters among 107 Italian provinces, based on the prevalence and trend of total SARS-CoV-2 cases. The hierarchical clustering yielded to the dendrogram shown in Figure 5.
According to the highest silhouette value reported in Table 4, we obtained an optimal clustering solution by partitioning the provincial dataset into three clusters. The hierarchical clusters’ composition is summarized in Table 5.
Specifically, the provinces in cluster 1 exhibited the highest prevalence and trend of SARS-CoV-2 cases; the provinces in cluster 2 had a median prevalence of 11.2 cases per 10,000 residents (IQR = 14.0) and a median increment of 3.9% (IQR = 2.4); and the provinces in cluster 3 were those with the lowest trend and a median prevalence of 11.6 cases per 10,000 residents (IQR = 12.7) (Figure 6A and Table 6). We further consolidated the hierarchical clustering by applying the k-means algorithm with a predefined number of three clusters. This approach yielded to the alternative clustering solution summarized in Table 7. In particular, the provinces in cluster 1 exhibited the highest prevalence of SARS-CoV-2 cases and a median increment of 3.2% (IQR = 1.1); the provinces in cluster 2 had a median prevalence of 42.7 cases per 10,000 residents (IQR = 44.0) and the highest increment; the provinces in cluster 3 were those with the lowest prevalence and trend of SARS-CoV-2 cases (Figure 6B and Table 6).

4. Discussion

Data on the prevalence and trend of total cases could be useful to monitor the SARS-CoV-2 epidemic at the regional level [17]. To the best of our knowledge, our study was the first applying two clustering methods to categorize the Italian regions into different groups, based on the prevalence and trend of SARS-CoV-2 positive cases prior to the starting of phase two. Specifically, we first used hierarchical clustering [11] to identify the groups with a similar prevalence of SARS-CoV-2 positive cases on 3 May and trend from 27 April to 3 May 2020. According to this method, Lombardy and Piedmont—the most hit Italian regions [3]—were grouped with Umbria and Marche into cluster 1, which was characterized by a high prevalence and trend of SARS-CoV-2 positive cases. Looking at the data, the Umbria and Marche were characterized by the values of prevalence and trend that were hybrid between cluster 1 and 2. Most regions, instead, were classified into a more heterogeneous cluster, with an intermediate prevalence and trend of positive cases. Finally, Umbria and Valle d’Aosta—those regions that simultaneously exhibited a low prevalence and trend of active cases [3]—were grouped into cluster 3. We further applied the k-means algorithm, a partitioning clustering used for splitting the dataset into a predefined number of groups [11]. Our combined approach allowed to improve the robustness of our findings by using a mixed algorithm that produced a better clustering solution [11]. In this case, however, all the regions maintained the same allocation of hierarchical clustering.
Our approach did not intend to explore the causality of the observed differences between clusters, nor to account for the spatial distribution of SARS-CoV-2 cases within the Italian territory. Actually, we aimed to provide a simple tool to monitor the SARS-CoV-2 epidemic across Italian regions. Indeed, our findings offered a snapshot of the epidemic, which could be helpful to outline the hierarchy of needs at the regional level. Specifically, some Italian regions—Lombardy and Piedmont, surely, but also Liguria and Marche—required the highest level of attention. For all the other regions, there were hopeful signs that the epidemic was slowly going towards a resolution. Among these, interestingly, Umbria and Valle d’Aosta were those with the best scenario. There were several explanations as to why Italian regions exhibited different prevalences and trends of SARS-CoV-2 cases. First, the epidemic initially started in Lombardy and Veneto and only after spread to the other regions with different times and intensities [2,3]. Thus, temporal and spatial distributions of SARS-CoV-2 cases represented the first reason motivating the differences observed across Italian regions. Beyond this, it is possible to hypothesize other grounds to support our clustering solution, for instance, comparing the epidemic spread between Lombardy and Veneto regions. Indeed, the first outbreaks concomitantly but independently occurred in two small Italian areas within the Lombardy and Veneto regions [2]. However, the scenario observed in Lombardy on the 3 May 2020 was much more serious than that in Veneto. Interestingly, Veneto imposed a wider testing campaign than Lombardy (771 vs. 409 tests per 10,000 residents), which may have contributed to a more rapid resolution of the epidemic [18]. Moreover, it is worth mentioning that performing a lower number of tests might alter the statistics and sustain the epidemic spread through contacts of positive but undocumented SARS-CoV-2 patients [19,20]. Equally important, the variability in the level of compliance with national restrictions and recommendations between and within regions could have influenced their efficacy at the regional level, yielding to the observed differences between clusters [21]. Finally, the demographic structure of each region, as well as environmental factors, could represent more suggestive explanations underlying a part of the variation between regions.
Our study also intended to provide an optimal clustering solution at the provincial level, despite that the data released by Italian authorities were less accurate. Whilst it would have been proper to work on positive cases, we used available data on the prevalence of total cases on 3 May and their trend from 27 April to 3 May 2020. Most provinces, actually, were first included in a cluster with the highest prevalence and trend of total cases. The remaining provinces, instead, were partitioned into two clusters with a similar prevalence of SARS-CoV-2 cases (approximately 11 cases per 10,000 residents) but a different trend from 27 April to 3 May: provinces in cluster 2 exhibited a median increment of 3.9%, while those in cluster 3 had the lowest trend (i.e., 1.5%). An alternative clustering solution obtained by the k-means algorithm categorized Italian provinces in a cluster with the lowest prevalence and trend of SARS-CoV-2 cases and in two clusters with the highest prevalence or the highest trend, respectively. The aforementioned reasons supporting the differences between regions could be equally applied to the partition of Italian provinces.
Our study had several limitations that should be considered when interpreting our results. First, it did not intend to investigate the causality of the observed differences, but rather to provide a simple tool to classify the regions and provinces according to the prevalence and trend of SARS-CoV-2 cases prior to the start of phase two in Italy. In line, our approach did not take into account the temporal and spatial distributions of SARS-CoV-2 cases within the Italian territory. Moreover, it considered only a part of the available data and knowledge, and hence further analyses could include additional indicators, characteristics and epidemic parameters (e.g., demographic structure, environmental factors, adherence to national restrictions, mobility data) at the regional and provincial level that could improve our cluster solutions and motivate the observed partitions [22]. A main limitation was with regard to the accuracy of the data used in our analysis and the differences in ascertainment and reporting capabilities between and within regions. For instance, we did not consider the proportion of unreported cases [20,23,24,25], though previous studies suggested how undocumented SARS-CoV-2 cases could silently sustain the epidemic spread [23]. Indeed, the testing strategy adopted in Italy might have caused an underestimation of infectious individuals, especially those with mild or no symptoms [19]. For this reason, we adjusted the prevalence and trend indicators for the regional number of tests performed per 10,000 residents.

5. Conclusions

In view of these considerations, our study described a simple and immediate approach to monitor the SARS-CoV-2 epidemic at the regional and provincial level. Interestingly, we partitioned both Italian regions and provinces into several groups according to the prevalence and trend of SARS-CoV-2 cases. These findings, at present, could be important to support policymakers in monitoring and planning future strategies against the epidemic. Indeed, they offered a snapshot of the epidemic, which could be helpful to outline the hierarchy of needs at the subnational level. However, the integration of our approach with further indicators and characteristics (i.e., the spatial and temporal distribution of SARS-CoV-2 cases, compliance with national control measures, demographic structure, environmental factors, and mobility data) could improve our findings, also allowing the application to different contexts and with additional aims. To do that, further studies working on additional and more accurate data at the regional and provincial level should be encouraged.

Author Contributions

Conceptualization, A.M. and A.A.; methodology, A.M. and M.B.; software, A.M.; formal analysis, A.M. and M.B.; data curation, A.M.; writing—original draft preparation, A.M. and M.B.; writing—review and editing, A.M., M.B. and A.A.; visualization, A.M.; supervision, A.A.; funding acquisition, A.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Assessorato della Salute, Regione Siciliana—Progetti Obiettivo di Piano Sanitario Nazionale (PSN 2014—4.9.2).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. World Health Organization. Coronavirus Disease (COVID-19) Dashboard. Available online: https://covid19.who.int/ (accessed on 4 May 2020).
  2. Day, M. Covid-19: Italy confirms 11 deaths as cases spread from north. BMJ 2020, 368, m757. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Italian Ministry of Health. Covid-19. Situation Report Update at 4 May 18:00. Available online: http://www.salute.gov.it/portale/nuovocoronavirus/dettaglioNotizieNuovoCoronavirus.jsp?lingua=italiano&menu=notizie&p=dalministero&id=4678 (accessed on 4 May 2020).
  4. Italian Ministry of Health. Novel Coronavirus. Available online: http://www.salute.gov.it/portale/nuovocoronavirus/homeNuovoCoronavirus.jsp?lingua=english (accessed on 4 May 2020).
  5. Signorelli, C.; Scognamiglio, T.; Odone, A. COVID-19 in Italy: Impact of containment measures and prevalence estimates of infection in the general population. Acta Biomed. 2020, 91, 175–179. [Google Scholar] [CrossRef] [PubMed]
  6. Giordano, G.; Blanchini, F.; Bruno, R.; Colaneri, P.; Di Filippo, A.; Di Matteo, A.; Colaneri, M. Modelling the COVID-19 epidemic and implementation of population-wide interventions in Italy. Nat. Med. 2020. [Google Scholar] [CrossRef] [PubMed]
  7. Maugeri, A.; Barchitta, M.; Battiato, S.; Agodi, A. Modeling the novel coronavirus (SARS-CoV-2) outbreak in Sicily, Italy. Int. J. Environ. Res. Public Health 2020, 17, 4964. [Google Scholar] [CrossRef] [PubMed]
  8. Tuite, A.R.; Ng, V.; Rees, E.; Fisman, D. Estimation of COVID-19 outbreak size in Italy. Lancet Infect. Dis 2020, 20, 537. [Google Scholar] [CrossRef] [Green Version]
  9. Gatto, M.; Bertuzzo, E.; Mari, L.; Miccoli, S.; Carraro, L.; Casagrandi, R.; Rinaldo, A. Spread and dynamics of the COVID-19 epidemic in Italy: Effects of emergency containment measures. Proc. Natl. Acad. Sci. USA 2020, 117, 10484–10491. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Remuzzi, A.; Remuzzi, G. COVID-19 and Italy: What next? Lancet 2020, 395, 1225–1228. [Google Scholar] [CrossRef]
  11. Altman, N.; Krzywinski, M. Clustering. Nat. Methods 2017, 14, 545–546. [Google Scholar] [CrossRef]
  12. GitHub. Provincial Data on SARS-CoV-2. Available online: https://github.com/pcm-dpc/COVID-19/tree/master/dati-province (accessed on 4 May 2020).
  13. GitHub. Regional Data on SARS-CoV-2. Available online: https://github.com/pcm-dpc/COVID-19/tree/master/dati-regioni (accessed on 4 May 2020).
  14. Dipartimento della Protezione Civile. Aggiornamento Casi COVID-19. Available online: http://opendatadpc.maps.arcgis.com/apps/opsdashboard/index.html#/b0c68bce2cce478eaac82fe38d4138b1 (accessed on 6 July 2020).
  15. Istituto Nazionale di Statistica, ISTAT. Available online: analysis-and-products/databases (accessed on 4 May 2020).
  16. Lovmar, L.; Ahlford, A.; Jonsson, M.; Syvänen, A.C. Silhouette scores for assessment of SNP genotype clusters. BMC Genom. 2005, 6, 35. [Google Scholar] [CrossRef] [PubMed]
  17. GIMBE. Profili Sanitari Della Cosiddetta Fase due: Strategie Ante e Post Covid-19. Available online: https://coronavirus.gimbe.org/var/uploads/contenuti/allegati/Audizione%20Fondazione%20GIMBE%207%20maggio%202020.pdf (accessed on 4 May 2020).
  18. Sebastiani, G.; Massa, M.; Riboli, E. Covid-19 epidemic in Italy: Evolution, projections and impact of government measures. Eur. J. Epidemiol. 2020, 35. [Google Scholar] [CrossRef] [PubMed]
  19. Onder, G.; Rezza, G.; Brusaferro, S. Case-fatality rate and characteristics of patients dying in relation to COVID-19 in Italy. JAMA 2020. [Google Scholar] [CrossRef] [PubMed]
  20. Maugeri, A.; Barchitta, M.; Battiato, S.; Agodi, A. Estimation of unreported novel coronavirus (SARS-CoV-2) infections from reported deaths: A susceptible-exposed-infectious-recovered-dead model. J. Clin. Med. 2020, 9, 1350. [Google Scholar] [CrossRef] [PubMed]
  21. Kucharski, A.J.; Klepac, P.; Conlan, A.J.K.; Kissler, S.M.; Tang, M.L.; Fry, H.; Gog, J.R.; Edmunds, W.J.; CMMID COVID-19 Working Group. Effectiveness of isolation, testing, contact tracing, and physical distancing on reducing transmission of SARS-CoV-2 in different settings: A mathematical modelling study. Lancet Infect. Dis. 2020. [Google Scholar] [CrossRef]
  22. Ji, Y.; Ma, Z.; Peppelenbosch, M.P.; Pan, Q. Potential association between COVID-19 mortality and health-care resource availability. Lancet Glob. Health 2020, 8, e480. [Google Scholar] [CrossRef] [Green Version]
  23. Li, R.; Pei, S.; Chen, B.; Song, Y.; Zhang, T.; Yang, W.; Shaman, J. Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV2). Science 2020, 368, 489–493. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Wang, H.; Wang, Z.; Dong, Y.; Chang, R.; Xu, C.; Yu, X.; Zhang, S.; Tsamlag, L.; Shang, M.; Huang, J.; et al. Phase-adjusted estimation of the number of Coronavirus Disease 2019 cases in Wuhan, China. Cell Discov. 2020, 6, 10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Zhao, S.; Musa, S.S.; Lin, Q.; Ran, J.; Yang, G.; Wang, W.; Lou, Y.; Yang, L.; Gao, D.; He, D.; et al. Estimating the unreported number of novel coronavirus (2019-nCoV) cases in China in the first half of January 2020: A data-driven modelling analysis of the early outbreak. J. Clin. Med. 2020, 9, 388. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Description and correlation of regional indicators: (A) the distribution of Italian regions by prevalence of SARS-CoV-2 positive cases, their trend, and the number of tests performed; and (B) the correlations between the prevalence of SARS-CoV-2 positive cases, their trend, and the number of tests performed; the results are reported as Spearman’s correlation coefficient from −1 (in red) to 1 (in green).
Figure 1. Description and correlation of regional indicators: (A) the distribution of Italian regions by prevalence of SARS-CoV-2 positive cases, their trend, and the number of tests performed; and (B) the correlations between the prevalence of SARS-CoV-2 positive cases, their trend, and the number of tests performed; the results are reported as Spearman’s correlation coefficient from −1 (in red) to 1 (in green).
Ijerph 17 05286 g001
Figure 2. Scatter plots of the relationships of the number of tests performed with (A) the prevalence of SARS-CoV-2 positive cases and (B) their trend. Indicators are reported as log-transformed values.
Figure 2. Scatter plots of the relationships of the number of tests performed with (A) the prevalence of SARS-CoV-2 positive cases and (B) their trend. Indicators are reported as log-transformed values.
Ijerph 17 05286 g002
Figure 3. Dendrogram of the hierarchical clustering of regions based on Ward’s criterion.
Figure 3. Dendrogram of the hierarchical clustering of regions based on Ward’s criterion.
Ijerph 17 05286 g003
Figure 4. Scatter plot illustrating how the egional clusters were distributed on the prevalence of SARS-CoV-2 positive cases and their trend. Clustering solution was obtained by the hierarchical clustering and consolidated by the k-means algorithm.
Figure 4. Scatter plot illustrating how the egional clusters were distributed on the prevalence of SARS-CoV-2 positive cases and their trend. Clustering solution was obtained by the hierarchical clustering and consolidated by the k-means algorithm.
Ijerph 17 05286 g004
Figure 5. Dendrogram of the hierarchical clustering of the provinces based on Ward’s criterion.
Figure 5. Dendrogram of the hierarchical clustering of the provinces based on Ward’s criterion.
Ijerph 17 05286 g005
Figure 6. Scatter plot illustrating how provincial clusters were distributed on the prevalence of SARS-CoV-2 cases and their trend: (A) the clustering solution obtained by hierarchical clustering; and (B) the clustering solution obtained using the k-means algorithm.
Figure 6. Scatter plot illustrating how provincial clusters were distributed on the prevalence of SARS-CoV-2 cases and their trend: (A) the clustering solution obtained by hierarchical clustering; and (B) the clustering solution obtained using the k-means algorithm.
Ijerph 17 05286 g006
Table 1. Characteristics of the Italian regions.
Table 1. Characteristics of the Italian regions.
RegionsResidentsPrevalence of Positive Cases (Per 10,000 Residents) aTrend of Positive Cases (%) bNumber of Tests (Per 10,000 Residents) a
Abruzzo1,315,19614.2−8.00%309.5
Apulia4,048,2427.31.50%164.1
Basilicata567,1183.4−10.60%250.6
Bolzano c527,75012.6−29.30%838.3
Calabria1,956,6873.6−10.20%198.5
Campania5,826,8604.7−5.30%148.5
Emilia-Romagna4,452,62920.3−26.00%442.6
Friuli Venezia Giulia1,215,5388.9−13.60%616.9
Lazio5,896,6937.4−3.90%255.9
Liguria1,556,98122.8−0.80%350
Lombardy10,036,25836.84.20%409.4
Marche1,531,75320.9−3.40%420.5
Molise308,4935.9−9.50%229.3
Piedmont4,375,86535.70.80%393.5
Sardinia1,648,1764.2−11.20%168.3
Sicily5,026,9894.43.80%170.9
Trento c539,89823.1−27.00%761.2
Tuscany3,736,96814.3−11.00%403.8
Umbria884,6402.1−36.20%438.9
Valle d’Aosta126,2028.6−53.60%641.8
Veneto4,905,03714.9−17.60%771.0
a Data are referred to 3 May 2020. b Weekly trend from 27 April to 3 May 2020. c Autonomous provinces.
Table 2. Average silhouette width for the hierarchical clustering of regions.
Table 2. Average silhouette width for the hierarchical clustering of regions.
Number of ClustersAverage Silhouette Width
(Standard Deviation)
20.566 (0.166)
30.632 (0.099)
40.493 (0.197)
50.530 (0.248)
Table 3. Comparisons of the prevalence of SARS-CoV-2 positive cases, their trend and the number of tests between regional clusters.
Table 3. Comparisons of the prevalence of SARS-CoV-2 positive cases, their trend and the number of tests between regional clusters.
ClustersPrevalence of Positive Cases Per 10,000 Residents aTrend of Positive Cases bNumber of Tests Per 10,000 Residents a
Cluster 129.3 (15.2)0% (6.1)401.5 (56.9)
Cluster 27.4 (9.9)−10.6% (12.3)255.9 (446.0)
Cluster 35.4 (6.5)−44.9% (17.4)540.4 (202.9)
p-Value0.0110.0070.356
Results are reported as the median (interquartile range), with p-values based on the Kruskal–Wallis test. a Data are referred to 3 May 2020 b Weekly trend from 27 April to 3 May 2020.
Table 4. Average silhouette width for the hierarchical clustering of provinces.
Table 4. Average silhouette width for the hierarchical clustering of provinces.
Number of ClustersAverage Silhouette Width
(Standard Deviation)
20.399 (0.169)
30.402 (0.229)
40.368 (0.214)
50.377 (0.225)
Table 5. Hierarchical clusters’ composition of provinces.
Table 5. Hierarchical clusters’ composition of provinces.
ClustersProvinces
Cluster 1Alessandria; Aosta; Asti; Belluno; Bergamo; Biella; Bologna; Brescia; Como; Cremona; Cuneo; Fermo; Firenze; Forlì-Cesena; Genova; Imperia; La Spezia; Lecco; Lodi; Mantova; Massa Carrara; Milano; Modena; Monza e della Brianza; Novara; Parma; Pavia; Pesaro e Urbino; Pescara; Piacenza; Reggio nell’Emilia; Rimini; Savona; Sondrio; Torino; Trento; Trieste; Varese; Venezia; Verbano-Cusio-Ossola; Vercelli; Verona; Vicenza
Cluster 2Arezzo; Avellino; Bari; Barletta-Andria-Trani; Benevento; Brindisi; Cagliari; Caltanissetta; Catania; Chieti; Enna; Ferrara; Foggia; Gorizia; Latina; Lecce; Livorno; Lucca; Macerata; Matera; Messina; Napoli; Nuoro; Palermo; Pistoia; Pordenone; Potenza; Prato; Ragusa; Rieti; Roma; Siracusa; Taranto; Terni; Treviso; Vibo Valentia; Viterbo
Cluster 3Agrigento; Ancona; Ascoli Piceno; Bolzano; Campobasso; Caserta; Catanzaro; Cosenza; Crotone; Frosinone; Grosseto; Isernia; L’Aquila; Oristano; Padova; Perugia; Pisa; Ravenna; Reggio di Calabria; Rovigo; Salerno; Sassari; Siena; Sud Sardegna; Teramo; Trapani; Udine
Table 6. Comparisons of the prevalence of SARS-CoV-2 cases and their trend between clusters of provinces.
Table 6. Comparisons of the prevalence of SARS-CoV-2 cases and their trend between clusters of provinces.
ClustersPrevalence of Total Cases Per 10,000 ResidentsTrend of Total Cases (%)
Hierarchical Clustering
Cluster 161.0 (31.0)5.7% (4.4)
Cluster 211.2 (14.0)3.9% (2.4)
Cluster 311.6 (12.7)1.5% (0.9)
p-Value<0.001<0.001
K-means Clustering
Cluster 153.0 (38.1)3.2% (1.1)
Cluster 242.7 (44.0)6.7% (3.3)
Cluster 37.4 (8.5)2.2% (2.0)
p-Value<0.001<0.001
Results are reported as the median (interquartile range), with p-values based on the Kruskal–Wallis test. Data are referred to 3 May 2020. Weekly trend from 27 April to 3 May 2020.
Table 7. K-means clusters’ composition of provinces.
Table 7. K-means clusters’ composition of provinces.
ClustersProvinces
Cluster 1Ancona; Aosta; Bergamo; Biella; Bolzano; Brescia; Cremona; Enna; Ferrara; Forlì-Cesena; La Spezia; Lecco; Lodi; Lucca; Macerata; Mantova; Massa Carrara; Modena; Padova; Parma; Pesaro e Urbino; Pordenone; Prato; Ravenna; Reggio nell’Emilia; Rieti; Rimini; Sondrio; Treviso; Trieste; Vercelli
Cluster 2Alessandria; Arezzo; Asti; Avellino; Belluno; Bologna; Brindisi; Caltanissetta; Chieti; Como; Cuneo; Fermo; Firenze; Foggia; Genova; Gorizia; Imperia; Matera; Milano; Monza e della Brianza; Novara; Palermo; Pavia; Pescara; Piacenza; Pistoia; Roma; Savona; Terni; Torino; Trento; Varese; Venezia; Verbano-Cusio-Ossola; Verona; Vicenza
Cluster 3Agrigento; Ascoli Piceno; Bari; Barletta-Andria-Trani; Benevento; Cagliari; Campobasso; Caserta; Catania; Catanzaro; Cosenza; Crotone; Frosinone; Grosseto; Isernia; L’Aquila; Latina; Lecce; Livorno; Messina; Napoli; Nuoro; Oristano; Perugia; Pisa; Potenza; Ragusa; Reggio di Calabria; Rovigo; Salerno; Sassari; Siena; Siracusa; Sud Sardegna; Taranto; Teramo; Trapani; Udine; Vibo Valentia; Viterbo

Share and Cite

MDPI and ACS Style

Maugeri, A.; Barchitta, M.; Agodi, A. A Clustering Approach to Classify Italian Regions and Provinces Based on Prevalence and Trend of SARS-CoV-2 Cases. Int. J. Environ. Res. Public Health 2020, 17, 5286. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17155286

AMA Style

Maugeri A, Barchitta M, Agodi A. A Clustering Approach to Classify Italian Regions and Provinces Based on Prevalence and Trend of SARS-CoV-2 Cases. International Journal of Environmental Research and Public Health. 2020; 17(15):5286. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17155286

Chicago/Turabian Style

Maugeri, Andrea, Martina Barchitta, and Antonella Agodi. 2020. "A Clustering Approach to Classify Italian Regions and Provinces Based on Prevalence and Trend of SARS-CoV-2 Cases" International Journal of Environmental Research and Public Health 17, no. 15: 5286. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph17155286

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