Next Article in Journal
A Fully Convolutional Network-Based Tube Contour Detection Method Using Multi-Exposure Images
Next Article in Special Issue
Lot-Size Models with Uncertain Demand Considering Its Skewness/Kurtosis and Stochastic Programming Applied to Hospital Pharmacy with Sensor-Related COVID-19 Data
Previous Article in Journal
Smart-Contract Aware Ethereum and Client-Fog-Cloud Healthcare System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Disjoint and Functional Principal Component Analysis for Infected Cases and Deaths Due to COVID-19 in South American Countries with Sensor-Related Data

by
Carlos Martin-Barreiro
1,2,
John A. Ramirez-Figueroa
1,2,
Xavier Cabezas
2,
Víctor Leiva
3,* and
M. Purificación Galindo-Villardón
1
1
Department of Statistics, Universidad de Salamanca, 37008 Salamanca, Spain
2
Faculty of Natural Sciences and Mathematics, Universidad Politécnica ESPOL, Guayaquil 090902, Ecuador
3
School of Industrial Engineering, Pontificia Universidad Católica de Valparaíso, Valparaíso 2362807, Chile
*
Author to whom correspondence should be addressed.
Submission received: 18 May 2021 / Revised: 7 June 2021 / Accepted: 11 June 2021 / Published: 14 June 2021

Abstract

:
In this paper, we group South American countries based on the number of infected cases and deaths due to COVID-19. The countries considered are: Argentina, Bolivia, Brazil, Chile, Colombia, Ecuador, Peru, Paraguay, Uruguay, and Venezuela. The data used are collected from a database of Johns Hopkins University, an institution that is dedicated to sensing and monitoring the evolution of the COVID-19 pandemic. A statistical analysis, based on principal components with modern and recent techniques, is conducted. Initially, utilizing the correlation matrix, standard components and varimax rotations are calculated. Then, by using disjoint components and functional components, the countries are grouped. An algorithm that allows us to keep the principal component analysis updated with a sensor in the data warehouse is designed. As reported in the conclusions, this grouping changes depending on the number of components considered, the type of principal component (standard, disjoint or functional) and the variable to be considered (infected cases or deaths). The results obtained are compared to the k-means technique. The COVID-19 cases and their deaths vary in the different countries due to diverse reasons, as reported in the conclusions.

1. Introduction

The COVID-19 pandemic has deteriorated the usual dynamics that rule the world [1]. Its impact on health and the worldwide economy is evident [2,3]. The number of people affected by the COVID-19 disease, the significant number of human lives that have been lost, and the consequences of the containment measures are urgent concerns that should be analyzed. This must be done in a technically responsible way to guarantee an effective response from the authorities of each country [4].
The Coronavirus (COVID-19) disease is caused by the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and was first identified in Wuhan City, Hubei Province, China, early in December 2019 [5]. COVID-19 arrived in South America at the end of February 2020 [4]. The coronavirus attacks some organs of the human body, but mainly the lungs [6]. In recent months, experts have detected the presence of some strains, as a result of the mutation of the virus, with higher rates of spread. As mentioned, the world has been greatly affected by the presence of this virus and governments are currently engaged in a mass vaccination process with the aim that people will be able to resume their activities in a normal way [4].
To confirm suspected cases of COVID-19, the real-time transcription-polymerase chain reaction (RT-PCR) test is applied, based on the detection of a certain amount of genetic fragments other than the virus in an individual [7]. Nevetheless, other diagnostic tools to detect COVID-19 cases have been proposed [8].
Proper planning of the response to an emergency is marked by the challenges of economic and political will [8]. However, the application and adaptation of statistical methods for data analysis that strengthen the interpretation of the information, taken from different points of interest, are actions that undoubtedly improve the quality of the response to the emergency [4]. These are actions that the citizens expect from their government authorities. As in other regions around the world, in South American countries, the number of infected people and deaths due to COVID-19 are increasing every day.
One popular multivariate statistical method is the principal component analysis (PCA), which allows for classification of variables or individuals [9,10]. The PCA has been used in several studies related to COVID-19 [11,12,13]. The objective of the present investigation is to apply the PCA to classify South American countries with respect to the number of infected and deaths due to COVID-19. We want to determine how the pandemic evolves in South American countries by using different multivariate methods, based on both matrix and functional approaches. These methods allow us to identify differences and similarities among the countries under study to understand the damage that the COVID-19 pandemic has caused them. The countries considered for this classification are: Argentina (ARG), Bolivia (BOL), Brazil (BRA), Chile (CHI), Colombia (COL), Ecuador (ECU), Peru (PER), Paraguay (PRY), Uruguay (URY), and Venezuela (VEN). We use modern techniques based on (i) standard PCA; (ii) disjoint principal component analysis (DPCA); and (iii) functional principal component analysis (FPCA). The full data set used in this investigation can be downloaded from https://ourworldindata.org/coronavirus (accessed on 12 June 2021). The data regarding the numbers (per million inhabitants) of infected cases and deaths are taken from 1 March 2020 to 15 March 2021. We design an algorithm that permits one to summarize the multivariate methods presented in this work to sense changes in the data by using a sensor and having an updated analysis.
The remainder of the paper is organized as follows: In Section 2, the fundamentals of the PCA, DPCA, and FPCA are provided. Section 3 and Section 4 apply the PCA, DPCA, and FPCA described in Section 2 to analyze the data collected. In Section 5, we discuss the results of the study and address the conclusions as well as ideas for future research.

2. Methods

2.1. Principal Component Analysis

The PCA is used in data matrices, with the measurements of a set of variables on a set of individuals being stored. The objective of the PCA is to reduce the original variables to a few latent variables (principal components) and to work in a reduced-dimension space that facilitates the interpretation. The principal components are obtained by solving an optimization problem defined as
max ( Z = w V w ) , subject to w w = 1 ,
where V = ( 1 / n ) X X is the sample matrix of covariances of the data matrix X (centered by columns), and the constraint w w = 1 indicates that the vector w is of unitary norm, with X being a matrix with n individuals and p variables. When applying a dimension reduction to q components ( q < p ), the solution to the optimization problem defined in (1) is the largest eigenvalue λ of V , that is, V w = λ w . The eigenvalues of V are sorted in descending order as λ 1 > > λ q . The first principal component w 1 is the eigenvector of the unitary norm associated with the largest eigenvalue of V , that is, λ 1 . Similarly, the second principal component w 2 is the eigenvector associated with λ 2 , where w 1 and w 2 are orthogonal, that is, w 1 w 2 = 0 , and so on for the rest of the components.
In order to improve the interpretation, it is possible t o rotate the extracted components. Varimax is one of the most popular rotation methods and is applied to the component loadings (correlations between variables and factors); see details of the varimax method on p. 270 of [14]. The varimax orthogonal rotation tries to maximize the variance of the squared loadings in each factor, leading to large loadings for a few variables and the rest of variables with small loadings. Thus, it is possible to obtain a subset of variables that affect each component. Then, it is desirable that only one loading per observed variable is large in absolute value so that the variables are mainly related to exactly one component [9].

2.2. Disjoint Component Analysis

The DPCA is used as an alternative to classify variables and build clusters. The DPCA seeks to obtain disjoint principal components, that is, to say that the principal components should be highly correlated to some variables but not to others. In the PCA, the loading matrix B is orthogonal of p × q dimension. However, in the DPCA, in addition to it being an orthogonal matrix, B = ( b i j ) must be a disjoint matrix, that is,
(i)
For all i, there is only one j, such that b i j 0 , with i { 1 , , p } and j { 1 , , q } .
(ii)
For all j, there exists at least one i, such that b i j 0 , with i { 1 , , p } and j { 1 , , q } .
The conditions (i) and (ii) above imply an increase of the constraints of the PCA optimization problem defined in (1). To obtain the disjoint components, the heuristic procedure presented in [15,16] can be used. Another procedure for computing disjoint components, based on particle swarm optimization, can be seen in [17], whereas an algorithm to compute disjoint components in three-way tables is presented in [18].

2.3. Functional Component Analysis

Unlike the PCA and DPCA, the FPCA considers data as functions, often over time [19,20,21]. The FPCA allows us to capture the variability of the phenomenon under investigation through its evolution over time. In the FPCA, each eigenvalue corresponds to an eigenfunction. The eigenfunctions describe the variability in the data set [22]. In the FPCA, the components are called harmonics.
In the PCA and DPCA, we have an n × p data matrix with n individuals and p original variables. Nevertheless, in the FPCA, we have p functions x i ( t ) , for i = 1 , , p , where the values of the data matrix x i t are converted into functions x i ( t ) , that is, the value of the variable i at time t, so one function is assigned to each original variable of the PCA [23].
In the multivariate context, a linear combination can be written as the dot product w x , where w = ( w 1 , , w p ) is a vector of weights and x = ( x 1 , , x p ) is a vector of variables. In the FPCA, linear combinations are represented as the dot product in a Hilbert space defined as
L 2 ( T ) : w , x = T w ( t ) x ( t ) d t ,
where T used in (2) is the time interval considered. In the FPCA, the sample covariance variance matrix is replaced by a bivariate function established by
v ( s , t ) = 1 p i = 1 p x i ( s ) x i ( t ) ,
where the functions x i considered in (3) are centered.
Note that a functional principal component must satisfy the expression given by
T v ( s , t ) w ( t ) d t = λ w ( t ) ,
where λ is an eigenvalue and w is its corresponding eigenfunction. The equation given in (4) is the functional version of the matrix equation V w = λ w . There are several ways to solve (4); see, for example, [24,25]. Here, we follow the method outlined in [21]. Therefore, if the eigenvalues obtained by solving the equation given in (4) are sorted in descending order, as λ 1 > > λ q , the first functional principal component w 1 is the eigenfunction corresponding to λ 1 . The second functional principal component w 2 is the eigenfunction associated with λ 2 of unitary norm and is orthogonal to w 1 , and so on in such a way that the new functional principal component is orthogonal to the one calculated in (4), that is,
T w i ( t ) w j ( t ) d t = 0 ,
for i < j . Thus, the score of the variable i in the component k, for k = 1 , , q , is stated as
f i k = T w k ( t ) x i ( t ) d t ,
where the formulation defined in (6) represents the dot product between w k and x i . Note that the expression given in (5) ensures that the functional components are orthogonal.

2.4. Algorithm and Computer Configurations

We have designed an algorithm that allows the multivariate methods presented in this work to sense changes in the data. This algorithm should be implemented in a software that, through the use of a sensor, would permit us to have an updated analysis. Algorithm 1 provides a summary of the methodology proposed in this study. The experimental analyses conducted here were carried out using a computer with the following characteristics: (i) OS: Windows 10 for 64 bits; (ii) RAM: 8 Gigabytes; and (iii) Processor: Intel Core i7-4510U 2–2.60 GHZ. The R software [26] was used for the PCA and FPCA, whereas the DPCA was implemented by the authors in C# based on the algorithm proposed in [15,16].
Algorithm 1 Approach for sensing and updating the component analyses
1:
Update the local data warehouse with the data extracted from the website.
2:
Indicate the presence of new data in the local repository when the data warehouse, through a sensor (trigger), sends a signal to the software that performs the component analyses.
3:
Generate the score matrices and the corresponding loading matrices when the software receives the signal from the sensor and processes the DPCA through its C# implementation.
4:
Obtain the corresponding matrices to process the PCA and FPCA when the software communicates with R and sends the request to it.
5:
Produce the corresponding plots that serve as support for the component analyses when the software sends the request to R.
6:
Store the results of the component analyses in the local data warehouse by means of the software.
7:
End showing all the results of the updated component analyses in the graphical user interface using the software.

3. Results I

3.1. PCA and DPCA Results for the Number of COVID-19 Infected Cases

A PCA was performed using the data in Table 1 with two components based on the correlation matrix. In this table, the columns show the countries and the rows show the number of COVID-19 infected cases (per million inhabitants) in each of the months from March 2020 to March 2021. In addition, a loading matrix was obtained but is not interpretable. Then, the varimax rotation was applied and a loading matrix, which is easier to interpret, was obtained; see Table 2. The percentage of variability explained by the model was 73.92 % . Figure 1a shows the plot of the space of countries with the data of COVID-19 infected cases. From this figure, note that ARG, COL, PRY, and VEN are well represented by the first component, while BOL, BRA, CHI, PER, and URY are well represented by the second component. However, ECU has very similar loadings in both components.
To further improve the interpretation, a DPCA was carried out. Before performing the calculation of its components, as part of the preprocessing, a centering and scaling were applied to the matrix with the data of COVID-19 infected cases. Table 2 also reports the loading matrix with two disjoint components for the number of infected cases. The percentage of variability explained by the model was 69.27 % . Of the variability, 4.65 % was lost, but we gained in terms of interpretation as we can conclude that ECU is better represented in the first component. In summary, when using two components for the data on the number of COVID-19 infected cases, the South American countries are classified as: (Group 1) ARG, COL, PRY, VEN, and ECU; and (Group 2) BOL, BRA, CHI, PER, and URY.
Next, we performed a PCA with three components. A variability percentage of 87.07 % was attained and once again the loading matrix obtained is not interpretable. Then, a varimax rotation was applied and no clearer interpretation was achieved, so we decided to use the DPCA again. When calculating three disjoint components, a variability percentage of 76.60 % was obtained and 10.47 % variability was lost, but an interpretable loading matrix was generated as observed in Table 3. In summary, when using three components for the data of the number of COVID-19 infected cases, the South American countries are classified as: (Group 1) ARG, COL, ECU, and VEN; (Group 2) BOL, BRA, CHI, and PER; and (Group 3) PRY and URY. Note that, when going from two to three components, one country came out of each of the two initial groups (PRY and URY, respectively) to form a third group.

3.2. PCA and DPCA Results for the Number of COVID-19 Deaths

Now, a similar analysis to that conducted in Section 3.1 was performed with the data in Table 4, where the columns show the countries and the rows show the number of COVID-19 deaths (per million inhabitants) in each of the months from March 2020 to March 2021.
Next, a PCA with two components was carried out using the correlation matrix. Then, a loading matrix was obtained but, once again, it is not interpretable. Thus, the varimax rotation was applied and a loading matrix that is easier to interpret was obtained in Table 5. The percentage of variability explained by the model was 69.75 % . Figure 1b displays the plot of the space of countries with the data of deaths.
Notice that ARG, COL, PRY, URY, and VEN are represented by the first component, whereas BOL, BRA, CHI, and PER are represented by the second component. Similarly to the analysis of the COVID-19 infected cases, ECU has very similar loadings in both components. Disjoint components were calculated to improve the analysis. Again, centering and scaling were applied to the matrix with the data on deaths. Table 5 also presents the loading matrix with two disjoint components. The percentage of variability explained by the model was 64.75 % , while 5 % of variability was lost; however, we gain in terms of interpretation, concluding that ECU is better represented in the second component. In summary, when using two components for the data on the number of COVID-19 deaths, the countries are classified as: (Group 1) ARG, COL, PRY, URY, and VEN; and (Group 2) BOL, BRA, CHI, ECU, and PER.
A PCA with three components was also carried out and the model attained 85.28 % of explained variability, but the loading matrix is not interpretable, so a varimax rotation was performed but the interpretation did not improve. When three disjoint components were computed, we obtained 74.78 % of explained variability with 10.50 % of variability being lost. Nevertheless, an interpretable loading matrix was achieved as can be seen in Table 6. In summary, when using three components for the data on the number of COVID-19 deaths, we have: (Group 1) ARG, COL, PRY, URY, and VEN; (Group 2) BRA, CHI, and PER; and (Group 3) BOL and ECU. Observe that, when going from two to three components, the first group keeps the same countries but the second group is subdivided into two others, where BOL and ECU come out of the second initial group to form a third group.

3.3. Analysis of the PCA and DPCA Results

For both data matrices reported in Table 1 and Table 4, the plots displayed in Figure 2a,b for the number of COVID-19 infected cases and deaths suggest two components if it is expected to explain about 70 % of the variance, and three components if it is expected to explain about 90 % of the variance. For more details about this criterion, see [27]. Some alternative criteria can be seen in [28,29].
Regarding the number of COVID-19 infected cases with two components, the countries with the highest score in the first component (C1) are ARG, PRY, and COL. The countries with the highest scores in the second component (C2) are BOL, CHI, and BRA; see Figure 1a. When comparing this clustering with Figure 3a, interpreting the components and considering the countries with the highest scores, note that the component C1 groups countries that took time to reach the first peak of the pandemic, while the C2 component groups countries that reached both peaks of the pandemic at similar times.
About the number of COVID-19 deaths with two components, the countries with the highest scores in C1 are PRY, ARG, and VEN. The countries with the highest scores in C2 are PER, BRA, and CHI; see Figure 1b. When comparing this clustering with Figure 3b, interpreting the components and considering the countries with the highest scores, observe that C1 groups countries that had a small number of deaths after the December 2020 holidays, while C2 groups countries that had an increase in the number of deaths after these holidays.
The previous sections showed the results of the classification of countries for the number of COVID-19 infected cases and deaths with two and three components. Despite the fact that two components are suggested in both situations, we decided to also perform the calculations with three components in order to include in this work an analysis with a higher percentage of variability (using a third principal component).
In relation to the number of COVID-19 infected cases, adding a third component, we gain 13.15 % and 7.33 % in explained variability with standard and disjoint components, respectively. Regarding the number of deaths, when going from two to three components, we gain 15.53 % and 10.03 % in variability with standard and disjoint components, respectively. In general, when using three principal components, there is an explained variability close to 90 % for both the number of infected and the number of deaths due to COVID-19. Figure 3a,b represent Table 1 and Table 4 graphically, respectively. Countries colored in blue, red, and green are grouped by components C1, C2 and C3, respectively; see Table 3 and Table 6. When comparing Figure 3a,b with the results obtained using two components, it is hard to detect similarities and differences in the ten South American countries. However, when using three components, we consider that a better classification is obtained.
Regarding the number of COVID-19 infected cases with three components, in C1 we have ARG, COL, VEN, and ECU; in C2, we have BOL, BRA, PER, and CHI; whereas in C3, we have PRY and URY. By observing Figure 3a, we can interpret that C1 groups those countries that have had similar behavior since December 2020. Note that ARG has a peak and a fall similar to that which happens with COL. In addition, ECU and VEN have similar behavior to each other, with small values in the number of infected. We can interpret component C2 as being those countries that had similar behavior in both peaks; even those times at which the peaks occurred, the time that elapsed between them was also similar. In the case of component C3, it groups those countries that had few infected cases at the beginning of the pandemic and that reached a peak after the December 2020 holidays. These similarities between PRY and URY are detected by component C3, but with only two components this would not have been detected.
About the number of COVID-19 deaths with three components, we have PRY, VEN, COL, ARG, and URY in C1; BRA, PER, and CHI in C2; and in C3, BOL and ECU. From Figure 3b, we can note that component C1 groups the countries with the smallest number of deaths, which are PRY, VEN, and URY. In this first component, we also identify ARG and COL with a smaller number of deaths at the beginning and at the end of the study period. We can interpret component C2 as those countries that had similar behavior in both peaks, the first between June 2020 and July 2020, and the second peak that was reached after the December 2020 holidays. Furthermore, observe that C3 groups those countries that reached their first peak in September 2020 and had an increase in the number of deaths after the holidays. The similarities between BOL and ECU are detected when C3 is considered, but not with C2.
Note that, when using the PCA and DPCA, both for the number of COVID-19 infected cases and deaths, the components are based on the number of peaks and the time in which these peaks are reached. In addition, observe that the beginning, duration, and end of the waves are considered in the groupings.
As mentioned, the PCA enables us to classify the countries analyzed in this study, with the components being latent countries. The DPCA was performed to contrast what was obtained with the PCA obtaining results consistent. The grouping of the countries, according to the number of COVID-19 infected cases and deaths (in both cases per million inhabitants), was carried out to comparatively determine similarities and differences among the countries during the time considered in our investigation (13 months).

3.4. k-Means Cluster Analysis for Countries

Next, a k-means analysis was carried out for the South American countries under study. Figure 4a,b show the results with two clusters, whereas Figure 4c,d display the results with three clusters. Regarding the number of COVID-19 infected cases with two clusters, note that in Figure 4a the countries are grouped in a different way when compared with the PCA/DPCA grouping using two components. From Figure 4b, we can conclude that the same occurs with respect to the number of deaths, that is, the results of the clustering with k-means differ from the clustering obtained with the PCA/DPCA. When comparing the three clusters of the k-means analysis with the three components of the PCA/DPCA, for the number of COVID-19 infected cases and their deaths, the results obtained are more similar than when comparing the two clusters of k-means and the two components of the PCA/DPCA. However, there are also differences. In this work, we consider that the results obtained with the PCA are more reliable and better reflect reality than the k-means analysis for the following reasons:
(i) The PCA reduces the dimensionality of the data, grouping the ten original countries into two or three latent countries. These latent countries (components) are linear combinations of the original countries. The PCA allows us to interpret the components obtained considering the original variables (countries) with the highest score in each of the components.
(ii) In the k-means analysis, we must indicate, from the beginning, the number of groups to be constructed, minimizing the squared distance of the variables (countries) from the centroid. Furthermore, k-means is an association method, which forms groups through their similarity. However, the PCA works with the full data set, maximizing the explained variance to obtain the components, one by one, trapping as much information as possible that it finds in the data.

4. Results II

4.1. FPCA Results for the Number of COVID-19 Infected Cases

Different to the PCA and DPCA, the FPCA allows us to capture the variability of the phenomenon under investigation through its evolution over time. In the FPCA, we associate a proper or harmonic function with each eigenvalue, which describes the main modes of variation existing in the data set. In Figure 5a,c, month 0 corresponds to March 2020, month 1 to April 2020, month 2 to May 2020, and so on until month 12, which corresponds to March 2021.
In Figure 5a, the first component (solid line) reflects the general evolution of the number of COVID-19 infected cases reported per month throughout the pandemic in all the countries considered. The greatest influence can be seen from months 5 to 11 (August 2020 to February 2021). The decline at the end of the curve is due to the fact that the data for March 2021 are only considered until the middle of the month. The second principal functional component (solid line) is shown in Figure 5b and it agrees with the variation detected in the first principal component. This second component only contributes 21.3 % of the explained variability and is not as important as the first one.
Table 7 and Figure 6a show the behavior of the countries over the course of the pandemic. Note that ECU, URY, and VEN have similar behavior, which is in contrast to BRA. In addition, the behavior of ARG and CHI differs from the rest. However, ARG and BOL demonstrate the opposite behavior from each other in the evolution of the number of new COVID-19 cases reported. We must consider that the number of new cases depends on the number of PCR tests applied, which might explain why some countries have a low rate of new cases. The position to the left of the origin in Figure 6a can be due to two situations, either: (i) the health policies are managed correctly; or (ii) insufficient detection PCR tests are carried out. A country in the first situation is URY, which is an example of correct management at the beginning of the pandemic that then ceased to be so.

4.2. FPCA Results for the Number of COVID-19 Deaths

Figure 5c shows the first functional component (solid line), reflecting the temporal evolution of the number of deaths per month and per million inhabitants globally (all countries considered simultaneously). The hightest variability is due to the contribution of the months of June to September 2020. Then, we see a decrease and a rise that can be attributed to the holiday season. As in the previous analysis of the number of COVID-19 infected cases, the decline at the end of the curve is due to the fact that the data for March 2021 are only considered until the middle of the month. The behavior of the second principal functional component (solid line) is similar to that of the first component (see Figure 5b), not providing more information on the temporal evolution of the event. Its percentage of variability explained is much less than for the first principal component.
Table 8 reports the loadings for the first two components and Figure 6b shows the relative positions of the countries considered with respect to the number of COVID-19 deaths per month and per million inhabitants throughout the pandemic. Considering the projections on the first main axis, we see that the countries with the smallest number of deaths (per month and per million) are PRY, URY, and VEN, which are the least affected countries with respect to the number of deaths, but ECU is the average. Forming a group, we have BOL, BRA, CHI, and COL, which means that their behavior regarding the number of deaths is similar, with PER being in the position furthest to the right of the axis, indicating that it is the country hardest hit (in percentage terms) by the pandemic. Note that the officially reported number of deaths due to COVID-19 may differ significantly from the true number of deaths, as has been observed in several countries in the region.

4.3. k-Means Cluster Analysis for Time

To complement the component analysis conducted, we also carried out a k-means analysis for time with three clusters. First, we show the results regarding the number of COVID-19 infected cases. Table 9 reports the scores of the countries in each of the clusters and Table 10 shows the membership of each studied month to each cluster.
Table 9 and Table 10 allow us to conclude that, from March to June 2020 (also March 2021), CHI is the country with the most infections, followed by PER and BRA. In addition, from July to November 2020, ARG is the country with most infections, followed by COL, BRA, and PER. From December 2020 to February 2021, BRA is the country with the most infections, followed by COL, ARG, and URY.
The k-means analysis is consistent with the results of the FPCA, since Cluster 1 corresponds to the months in which the two functional components show a growth in the number of COVID-19 infected cases. Cluster 2 is associated with the months in which the functional components show a decrease. Cluster 3 is related to the months in which there is a more pronounced growth than in Cluster 1 (as a consequence of the end of the year holidays).
Now we show the results regarding the number of COVID-19 deaths. Table 11 reports the scores of the countries in each of the clusters, whereas Table 12 presents the membership of each studied month to its cluster. Table 11 and Table 12 allow us to conclude that, from March to May 2020 (also March 2021), and in November and December 2020, BRA is the country that presented most deaths, followed by ARG, PER, ECU, and COL. From June to August 2020, and in January and February 2021, PER had the largest number of deaths, followed by CHI, BRA, COL, and BOL. In September and October 2020, ARG presented more deaths, followed by ECU, BOL, and COL.
The three clusters detected by the k-means method correspond to the months in which the functional components show growth in the number of deaths, a more attenuated growth, and a decrease. Thus, Cluster 2 corresponds to the months in which the main functional components show a marked growth in the number of deaths. Cluster 3 is related to the months in which the functional components indicate a less accelerated growth in the number of deaths compared to Cluster 2. Note that Cluster 1 is associated with the months in which there is a decrease in the number of deaths due to COVID-19.

5. Conclusions

The behavior of the number of infected cases with COVID-19, and the number of deaths due to this disease, can vary in the countries for different reasons: (i) political and economic decisions; (ii) health infrastructure; (iii) the discipline of the people; (iv) the environment; and (v) the spread rates of the different strains of the virus, among other factors.
Principal component analyses based on recent and modern methods were carried out to study the relationship among the aforementioned variables in ten South American countries. The countries considered were: Argentina, Bolivia, Brazil, Chile, Colombia, Ecuador, Peru, Paraguay, Uruguay, and Venezuela. A k-means analysis was conducted to compare our principal component analyses, which resulted in general consistent.
By using a principal component analysis with varimax rotation and a disjoint component analysis, we report the following results: with two components, and considering the number of COVID-19 infected cases, there are two groups of countries, with Argentina, Colombia, Ecuador, and Venezuela in one group, while Bolivia, Brazil, Chile, Peru, and Uruguay are in another group. When increasing the number of components to three, Paraguay and Uruguay moved away from the other countries and formed a third group.
When considering the number of COVID-19 deaths and two components, we established two groups formed by Argentina, Colombia, Paraguay, Uruguay, and Venezuela in one group, while Bolivia, Brazil, Chile, Ecuador, and Peru were in another group. However, with three components, Bolivia and Ecuador moved away from the other countries and formed a third group.
In the case of using functional components, we conclude that both components showed a general evolution of the number of COVID-19 infected cases reported per month throughout the pandemic, in all the countries considered. The largest values were detected in the months of July, August and September 2020 on the one hand, and in the months of January and February 2021 on the other hand. Regarding the relative positions of the countries, Ecuador, Uruguay, and Venezuela had very similar behavior, and Brazil, Chile, and Peru also behaved similarly to each other. Another group was made up of Argentina and Colombia, while Bolivia and Paraguay were far from the rest of the countries and between them.
Considering the number of deaths using functional components, there are peaks in the months of September 2020 on the one hand, and January and February 2021 on the other hand, which is similar to what happened with the number of COVID-19 infected cases. The countries with the smallest number of deaths were Paraguay, Uruguay, and Venezuela, whereas a second group was formed by Bolivia, Brazil, and Chile, with their behavior regarding the number of deaths being similar and greater than the average. In the case of Peru, this country was the most affected by the pandemic, while Ecuador was within the average; Argentina was a little more affected than Ecuador, but less than Chile. It should be borne in mind that the officially reported number of deaths due to COVID-19 may differ significantly from the number of true deaths, as has been observed in several South American countries.
Data related to COVID-19 can change over time. Therefore, the conclusions to be obtained from them, after statistical analyses, are sensitive to these changes. As mentioned, the data were obtained from the repository for the 2019 Novel Coronavirus Visual Dashboard operated by the Johns Hopkins University Center for Systems Science and Engineering (JHU CSSE) [30], which keeps the data updated, particularly the number of COVID-19 infected cases and the number of deaths due to this disease. These data were also used in [31]. Note that the number of confirmed cases is less than the number of true cases because not all people are tested with PCR. This means that the number of confirmed cases depends on how many PCR tests a country applies.
When looking at Figure 3, we can see that some countries, such as Bolivia, Brazil, Chile, and Peru, reacted late in controlling the pandemic because they were the countries to reach a first peak early with respect to the number of COVID-19 infected cases and deaths. Figure 3 also shows that countries such as Uruguay, Paraguay, and Venezuela had better control of the pandemic since the peaks of infected cases and deaths due to COVID-19 were reached some time after the aforementioned countries.
Despite the fact that some countries already knew that they had COVID-19 infected cases in their territory since the beginning of February 2020, it was in mid-March 2020 when the governments made the first decisions to control the spread of the virus. With full knowledge of what was happening in Asia and Europe, some South American countries minimized the impact that COVID-19 would have on the economy [2,3] and on people’s health. Other countries were a little more cautious and took strong measures from mid-March 2020, such as closing borders and confining citizens. However, it must also be said that, regardless of the policies taken by the governments [4], the climate, the limitations of the health systems (infrastructure, personnel), the behavior of citizens, among other factors, also affected the variables of this study. By classifying countries into groups, we can comparatively study the damage caused by the COVID-19 pandemic, analyzing peaks and waves.
We consider that a better interpretation and classification of the countries is obtained using three principal components, leading to the following conclusions:
(i) For the number of COVID-19 infected cases, the components C1, C2, and C3 can be interpreted as: C1 grouped those countries that had similar behavior since December 2020; C2 grouped those countries with similar behavior throughout the entire period of analysis; and C3 grouped those countries that had few infected cases at the beginning of the pandemic and reached a peak after the December 2020 holidays.
(ii) For the number of COVID-19 deaths, the components C1, C2, and C3 can be interpreted as: C1 grouped those countries with the smallest number of deaths, at the beginning and at the end of the study period; C2 grouped those countries with similar behavior in both peaks; and C3 grouped the countries that reached their first peak in September 2020 and reached a second peak after the December 2020 holidays.
For further research, this study can be replicated for other regions of the planet such as the USA, Europe, and Asia. A study can also be conducted to implement alternative criteria for selecting the correct number of components [28,29].

Author Contributions

Data curation, J.A.R.-F., C.M.-B.; formal analysis, C.M.-B., J.A.R.-F., X.C., V.L., M.P.G.-V.; investigation, C.M.-B., J.A.R.-F., X.C.; methodology, C.M.-B., J.A.R.-F., X.C., V.L., M.P.G.-V.; writing—original draft, C.M.-B., J.A.R.-F., X.C.; writing—review and editing, V.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported partially by project grant Fondecyt 1200525 from the National Agency for Research and Development (ANID) of the Chilean government under the Ministry of Science, Technology, Knowledge, and Innovation.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The full data set used in this investigation can be downloaded from https://ourworldindata.org/coronavirus (accessed on 12 June 2021).

Acknowledgments

The authors would also like to thank the Editor and Reviewers for their constructive comments which led to improve the presentation of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Amaral, F.; Casaca, W.; Oishi, C.M.; Cuminato, J.A. Tsowards providing effective data-driven responses to predict the Covid-19 in São Paulo and Brazil. Sensors 2021, 21, 540. [Google Scholar] [CrossRef] [PubMed]
  2. Chahuan-Jimenez, K.; Rubilar, R.; de la Fuente-Mella, H.; Leiva, V. Breakpoint analysis for the COVID-19 pandemic and its effect on the stock markets. Entropy 2021, 23, 100. [Google Scholar] [CrossRef] [PubMed]
  3. Liu, Y.; Mao, C.; Leiva, V.; Liu, S.; Silva Neto, W.A. Asymmetric autoregressive models: Statistical aspects and a financial application under COVID-19 pandemic. J. Appl. Stat. 2021, 1–25. [Google Scholar] [CrossRef]
  4. Jerez-Lillo, N.; Lagos-Álvarez, B.; Muñoz-Gutierrez, J.; Figueroa-Zúñiga, J.I.; Leiva, V. A statistical analysis for the epidemiological surveillance of COVID-19 in Chile. Signa Vitae 2021, in press. [Google Scholar]
  5. Sun, P.; Lu, X.; Xu, C.; Sun, W.; Pan, B. Understanding of COVID-19 based on current evidence. J. Med. Virol. 2021, 92, 548–551. [Google Scholar] [CrossRef] [PubMed]
  6. Zhu, N.; Zhang, D.; Wang, W.; Li, X.; Yang, B.; Song, J.; Zhao, X.; Huang, B.; Shi, W.; Lu, R.; et al. A novel coronavirus from patients with pneumonia in China. N. Engl. J. Med. 2020, 382, 727–733. [Google Scholar] [CrossRef] [PubMed]
  7. Tang, Y.W.; Schmitz, J.E.; Persing, D.H.; Stratton, C.W. Laboratory diagnosis of COVID-19: Current issues and challenge. J. Clin. Microbiol. 2020, 58, 512–520. [Google Scholar] [CrossRef] [Green Version]
  8. Bustos, N.; Tello, M.; Droppelmann, G.; Garcia, N.; Feijoo, F.; Leiva, V. Machine learning techniques as an efficient alternative diagnostic tool for COVID-19 cases. Signa Vitae 2021, in press. [Google Scholar]
  9. Giordani, P. Principal Component Analysis. In Encyclopedia of Social Network Analysis and Mining; Springer: New York, NY, USA, 2018. [Google Scholar]
  10. Leon-Medina, J.X.; Anaya, M.; Pares, N.; Tibaduiza, D.A.; Pozo, F. Structural damage classification in a Jacket-type wind-turbine foundation using principal component analysis and extreme gradient boosting. Sensors 2021, 21, 2748. [Google Scholar] [CrossRef]
  11. Mahmoudi, M.R.; Heydari, M.H.; Qasem, S.N.; Mosavi, A.; Band, S.S. Principal component analysis to study the relations between the spread rates of covid-19 in high risks countries. Alex. Eng. J. 2021, 60, 457–464. [Google Scholar] [CrossRef]
  12. Noy, I.; Doan, N.; Ferrarini, B.; Park, D. Measuring the economic risk of covid-19. Glob. Policy 2020, 11, 413–423. [Google Scholar] [CrossRef]
  13. Ye, W.; Lu, W.; Tang, Y.; Chen, G.; Li, X.; Ji, C.; Hou, M.; Zeng, G.; Lan, X.; Wang, Y.; et al. Identification of covid-19 clinical phenotypes by principal component analysis-based cluster analysis. Front. Med. 2020, 7, 782. [Google Scholar] [CrossRef]
  14. Jolliffe, I.T. Principal Component Analysis; Springer: New York, NY, USA, 2002. [Google Scholar]
  15. Vichi, M.; Saporta, G. Clustering and disjoint principal component analysis. Comput. Stat. Data Anal. 2009, 53, 3194–3208. [Google Scholar] [CrossRef]
  16. Ferrara, C.; Martella, F.; Vichi, M. Dimensions of well-being and their statistical measurements. In Topics in Theoretical and Applied Statistics; Springer: Cham, Switzerland, 2016; pp. 85–99. [Google Scholar]
  17. Ramirez-Figueroa, J.A.; Martin-Barreiro, C.; Nieto-Librero, A.B.; Leiva, V.; Galindo-Villardón, M.P. A new principal component analysis by particle swarm optimization with an environmental application for data science. Stoch. Environ. Res. Risk Assess. 2021. [Google Scholar] [CrossRef]
  18. Martin-Barreiro, C.; Ramirez-Figueroa, J.A.; Nieto-Librero, A.B.; Leiva, V.; Martin-Casado, A.; Galindo-Villardón, M.P. A new algorithm for computing disjoint orthogonal components in the three-way tucker model. Mathematics 2021, 9, 203. [Google Scholar] [CrossRef]
  19. Giraldo, R.; Herrera, L.; Leiva, V. Cokriging prediction using as secondary variable a functional random field with application in environmental pollution. Mathematics 2020, 8, 1305. [Google Scholar] [CrossRef]
  20. Martinez, S.; Giraldo, R.; Leiva, V. Birnbaum-Saunders functional regression models for spatial data. Stoch. Environ. Res. Risk Assess. 2019, 33, 1765–1780. [Google Scholar] [CrossRef]
  21. Ramsay, J.O.; Silverman, B.W. Functional Data Analysis; Springer: New York, NY, USA, 2005. [Google Scholar]
  22. Di Salvo, F.; Ruggieri, M.; Plaia, A. Functional principal component analysis for multivariate multidimensional environmental data. Environ. Ecol. Stat. 2015, 22, 739–757. [Google Scholar] [CrossRef]
  23. Ramsay, J.; Hooker, G.; Graves, S. Introduction to Functional Data Analysis. In Functional Data Analysis with R and MATLAB; Springer: New York, NY, USA, 2009. [Google Scholar]
  24. Kokoszka, P.; Reimherr, M. Introduction to Functional Data Analysis; Chapman and Hall: New York, NY, USA, 2017. [Google Scholar]
  25. Hsing, T.; Eubank, R. Theoretical Foundations of Functional Data Analysis, with an Introduction to Linear Operators; Wiley: New York, NY, USA, 2013. [Google Scholar]
  26. R Core Team. R: A Language and Environment for Statistical Computing; R Core Team: Vienna, Austria, 2020. [Google Scholar]
  27. Jolliffe, I.T.; Cadima, J. Principal component analysis: A review and recent developments. Philos. Trans. R. Soc. A 2016, 374, 20150202. [Google Scholar] [CrossRef]
  28. Cangelosi, R.; Goriely, A. Component retention in principal component analysis with application to cDNA microarray data. Biol. Direct 2007, 2, 2. [Google Scholar] [CrossRef] [Green Version]
  29. Donoho, D.L.; Gavish, M. The optimal hard threshold for singular values is 4/V3. IEEE Trans. Inf. Theory 2013, 60, 5040–5053. [Google Scholar]
  30. Dong, E.; Du, H.; Gardner, L. An interactive web-based dashboard to track COVID-19 in real time. Lancet Infect. Dis. 2020, 20, 533–534. [Google Scholar] [CrossRef]
  31. Ritchie, H.; Ortiz-Ospina, E.; Beltekian, D.; Mathieu, E.; Hasell, J.; Macdonald, B.; Giattino, C.; Appel, C.; Rodes-Guirao, L.; Roser, M. Coronavirus Pandemic (COVID-19). Our World in Data. 2020. Available online: https://ourworldindata.org/coronavirus/country/pitcairn (accessed on 12 June 2021).
Figure 1. Space of countries for data of the number of COVID-19 (a) infected cases and (b) deaths.
Figure 1. Space of countries for data of the number of COVID-19 (a) infected cases and (b) deaths.
Sensors 21 04094 g001
Figure 2. Cumulative variance plots of the number of COVID-19 (a) infected cases and (b) deaths.
Figure 2. Cumulative variance plots of the number of COVID-19 (a) infected cases and (b) deaths.
Sensors 21 04094 g002
Figure 3. Plots of the number of COVID-19 (a) infected cases and (b) deaths for the indicated country and month.
Figure 3. Plots of the number of COVID-19 (a) infected cases and (b) deaths for the indicated country and month.
Sensors 21 04094 g003
Figure 4. Cluster plots with (a,b) k = 2 and (c,d) 3 of the number of COVID-19 (a,c) infected cases and (b,d) deaths.
Figure 4. Cluster plots with (a,b) k = 2 and (c,d) 3 of the number of COVID-19 (a,c) infected cases and (b,d) deaths.
Sensors 21 04094 g004
Figure 5. Harmonic components (a,c) 1 and (b,d) 2 of the FPCA for data on the number of COVID-19 (a,b) infected cases, and (c,d) deaths.
Figure 5. Harmonic components (a,c) 1 and (b,d) 2 of the FPCA for data on the number of COVID-19 (a,b) infected cases, and (c,d) deaths.
Sensors 21 04094 g005
Figure 6. FPCA space of countries for data on the number of COVID-19 (a) infected cases and (b) deaths.
Figure 6. FPCA space of countries for data on the number of COVID-19 (a) infected cases and (b) deaths.
Sensors 21 04094 g006
Table 1. Data matrix of the number of COVID-19 infected cases per million inhabitants for the indicated month and country.
Table 1. Data matrix of the number of COVID-19 infected cases per million inhabitants for the indicated month and country.
MonthARGBOLBRACHICOLECUPERPRYURYVEN
2020-0323.3209.16826.887148.67117.808126.96132.3009.11297.3034.748
2020-0474.65390.807383.281777.248110.0771286.2831089.14228.18187.8026.963
2020-05274.872755.1602011.9635537.081449.583802.8093866.934100.94551.82041.393
2020-061054.9431990.6544173.8538152.4001345.499982.4823661.822173.15032.533151.992
2020-072804.9513732.5375929.8423990.0173884.6461639.3413708.588437.01394.426448.096
2020-085010.0483410.3415860.8892932.5386280.8161610.3787269.0511727.86095.286990.088
2020-097373.8311603.0954246.6392681.7574217.2341319.4974992.1963238.125129.832998.524
2020-109202.697552.2153409.2592472.5064805.2501765.2822681.4403144.328310.331594.001
2020-115699.844252.8013764.9392170.2544768.2451388.2541891.0102697.360786.763365.067
2020-124446.8971320.6516304.5672993.7966406.2621123.7831595.5113576.2923817.801392.603
2021-016675.9564858.2917192.1426179.8868885.2882171.7313733.5503546.4296511.452470.146
2021-023985.4612756.3546334.8315101.2203081.7062002.3175629.7763679.8984679.701428.651
2021-031379.088638.3953062.8882266.347673.508605.6781703.2702038.4072444.060147.877
Table 2. Loading matrices for data of the number of COVID infected cases with two components (C1 and C2) for the indicated country and method.
Table 2. Loading matrices for data of the number of COVID infected cases with two components (C1 and C2) for the indicated country and method.
VarimaxDPCA
CountryC1C2C1C2
ARG0.9430.0680.4840.000
BOL0.3000.8740.0000.519
BRA0.5850.7430.0000.519
CHI−0.1630.8710.0000.415
COL0.8270.4010.4660.000
ECU0.6490.5540.4170.000
PER0.3160.6870.0000.416
PRY0.8830.1210.4380.000
URY0.3810.4460.0000.340
VEN0.7940.2240.4280.000
Table 3. Loading matrix for data on the number of COVID infected cases with three components (C1, C2, and C3) for the indicated country using the DPCA.
Table 3. Loading matrix for data on the number of COVID infected cases with three components (C1, C2, and C3) for the indicated country using the DPCA.
DPCA
CountryC1C2C3
ARG0.5240.0000.000
BOL0.0000.5410.000
BRA0.0000.5250.000
CHI0.0000.4460.000
COL0.5150.0000.000
ECU0.4720.0000.000
PER0.0000.4830.000
PRY0.0000.0000.707
URY0.0000.0000.707
VEN0.4870.0000.000
Table 4. Data matrix of the number of COVID-19 deaths per million inhabitants for the indicated month and country.
Table 4. Data matrix of the number of COVID-19 deaths per million inhabitants for the indicated month and country.
MonthARGBOLBRACHICOLECUPERPRYURYVEN
2020-030.5980.5140.9460.6250.3154.250.9120.420.2880.105
2020-044.2284.80227.30911.2455.44646.75930.970.984.6080.456
2020-057.10321.505109.65243.26112.695139.319104.7850.141.441.46
2020-0616.99269.389142.453242.4147.0766.258156.8330.841.441.3
2020-0749.472158.827154.691197.164133.07266.599283.3944.4812.3043.982
2020-08113.22175.617135.99195.837187.82448.405296.22238.8372.5927.807
2020-09183.136251.691106.18775.957124.52272.006109.42674.4471.1528.511
2020-10311.20265.10874.95476.687104.47774.53461.11176.6912.885.981
2020-11170.99119.87762.27262.93107.14444.83345.85449.355.4723.489
2020-1299.89717.817102.69562.666126.70432.47753.28970.94329.9434.615
2021-01104.633103.999139.04396.464211.66446.757101.48263.93173.4095.663
2021-0288.304108.795143.198110.903113.65753.959159.92464.91249.5155.453
2021-0330.84320.13373.92733.16519.78916.66456.53328.87920.1511.934
Table 5. Loading matrices for data on the number of COVID deaths with two components for the indicated country and method.
Table 5. Loading matrices for data on the number of COVID deaths with two components for the indicated country and method.
CountryVarimaxDPCA
ARG0.849−0.0630.4380.000
BOL0.5020.7160.0000.473
BRA0.2270.8910.0000.512
CHI−0.0980.8470.0000.452
COL0.7270.5170.4760.000
ECU0.3650.260.0000.247
PER0.0110.9350.0000.499
PRY0.964−0.0740.5090.000
URY0.3940.0250.2590.000
VEN0.8800.4180.5050.000
Table 6. Loading matrix for data on the number of COVID deaths with three components (C1, C2, and C3) for the indicated country using the DPCA.
Table 6. Loading matrix for data on the number of COVID deaths with three components (C1, C2, and C3) for the indicated country using the DPCA.
DPCA
CountryC1C2C3
ARG0.4380.0000.000
BOL0.0000.0000.707
BRA0.0000.5940.000
CHI0.0000.5660.000
COL0.4760.0000.000
ECU0.0000.0000.707
PER0.0000.5710.000
PRY0.5090.0000.000
URY0.2590.0000.000
VEN0.5050.0000.000
Table 7. Loading matrix for data on the number of COVID infected cases with two components (C1 and C2) for the indicated country using the FPCA.
Table 7. Loading matrix for data on the number of COVID infected cases with two components (C1 and C2) for the indicated country using the FPCA.
FPCA
CountryC1C2
ARG7707.767−4791.927
BOL−3319.5352075.921
BRA5959.8082033.667
CHI1755.4685492.041
COL5692.081−1991.679
ECU−4878.641−249.683
PER2686.4483222.517
PRY−2023.569−3145.071
URY−5732.866−1702.523
VEN−7846.961−943.263
Table 8. Loading matrix for data on the number of COVID-19 deaths with two components (C1 and C2) for the indicated country using the FPCA.
Table 8. Loading matrix for data on the number of COVID-19 deaths with two components (C1 and C2) for the indicated country using the FPCA.
FPCA
CountryC1C2
ARG64.104266.481
BOL100.787−6.533
BRA91.109−47.293
CHI86.967−101.372
COL90.49061.045
ECU1.62136.494
PER243.162−133.709
PRY−162.39643.453
URY−252.150−55.442
VEN−263.694−63.125
Table 9. Center of the listed cluster for the number of COVID-19 infected cases in the indicated country using the k-means method.
Table 9. Center of the listed cluster for the number of COVID-19 infected cases in the indicated country using the k-means method.
ClusterARGBOLBRACHICOLECUPERPRYURYVEN
1561.375696.8371931.7743376.349519.295760.8432070.694469.959542.70470.595
26018.2741910.1984642.3142849.4144791.2381544.5504108.4572248.937283.328679.155
35036.1052978.4326610.5134758.3016124.4191765.9443652.9463600.8735002.985430.467
Table 10. Membership of the listed cluster for the number of COVID-19 infected cases in the indicated month using the k-means method.
Table 10. Membership of the listed cluster for the number of COVID-19 infected cases in the indicated month using the k-means method.
Month20-0320-0420-0520-0620-0720-0820-0920-1020-1120-1221-0121-0221-03
Cluster1111222223331
Table 11. Center of the listed cluster for the number of COVID-19 deaths in the indicated country using the k-means method.
Table 11. Center of the listed cluster for the number of COVID-19 deaths in the indicated country using the k-means method.
ClusterARGBOLBRACHICOLECUPERPRYURYVEN
1247.169158.40090.57176.322114.499173.27085.26975.5692.0167.246
252.27714.10862.80035.64945.34947.38448.72425.11910.3172.010
374.524123.325143.075148.556138.65756.396199.57134.60025.8524.841
Table 12. Membership of the listed cluster for the number of COVID-19 deaths in the indicated month using the k-means method.
Table 12. Membership of the listed cluster for the number of COVID-19 deaths in the indicated month using the k-means method.
Month20-0320-0420-0520-0620-0720-0820-0920-1020-1120-1221-0121-0221-03
Cluster2223331122332
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Martin-Barreiro, C.; Ramirez-Figueroa, J.A.; Cabezas, X.; Leiva, V.; Galindo-Villardón, M.P. Disjoint and Functional Principal Component Analysis for Infected Cases and Deaths Due to COVID-19 in South American Countries with Sensor-Related Data. Sensors 2021, 21, 4094. https://0-doi-org.brum.beds.ac.uk/10.3390/s21124094

AMA Style

Martin-Barreiro C, Ramirez-Figueroa JA, Cabezas X, Leiva V, Galindo-Villardón MP. Disjoint and Functional Principal Component Analysis for Infected Cases and Deaths Due to COVID-19 in South American Countries with Sensor-Related Data. Sensors. 2021; 21(12):4094. https://0-doi-org.brum.beds.ac.uk/10.3390/s21124094

Chicago/Turabian Style

Martin-Barreiro, Carlos, John A. Ramirez-Figueroa, Xavier Cabezas, Víctor Leiva, and M. Purificación Galindo-Villardón. 2021. "Disjoint and Functional Principal Component Analysis for Infected Cases and Deaths Due to COVID-19 in South American Countries with Sensor-Related Data" Sensors 21, no. 12: 4094. https://0-doi-org.brum.beds.ac.uk/10.3390/s21124094

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