Next Article in Journal
Alterations in microRNA Expression during Hematopoietic Stem Cell Mobilization
Next Article in Special Issue
COVID-19 Shuts Doors to Flu but Keeps Them Open to Rhinoviruses
Previous Article in Journal
Effect of Foliar Fertigation of Chitosan Nanoparticles on Cadmium Accumulation and Toxicity in Solanum lycopersicum
Previous Article in Special Issue
Dental Students in Germany throughout the COVID-19 Pandemic: A Psychological Assessment and Cross-Sectional Survey
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Protection Strategy against an Epidemic Disease on Edge-Weighted Graphs Applied to a COVID-19 Case

1
Laboratorio de Investigación Lab[e]saM, Departamento de Matemática y Estadística, Universidad de Playa Ancha, Valparaíso 2340000, Chile
2
Escuela de Enfermería, Universidad de Valparaíso, Viña del Mar 2520000, Chile
3
Escuela de Ingeniería Civil Informática, Universidad de Valparaíso, Valparaíso 2340000, Chile
4
Centro Nacional de Sistemas de Información en Salud, Santiago 8320000, Chile
*
Authors to whom correspondence should be addressed.
Biology 2021, 10(7), 667; https://doi.org/10.3390/biology10070667
Submission received: 17 May 2021 / Revised: 9 June 2021 / Accepted: 25 June 2021 / Published: 15 July 2021
(This article belongs to the Special Issue Coronavirus Disease 2019 (COVID-19))

Abstract

:

Simple Summary

Infectious diseases have been part of human history. Countless epidemics have produced high mortality rates in vulnerable populations. With the understanding of the spread of these types of diseases, population groups have been able to adapt and better cope with infections. Given the COVID-19 pandemic, one of the strategies used is the modeling of infectious diseases with the aim of establishing protection measures for people and stopping the spread of the epidemic. Our study evaluates protection strategies through infectious disease modeling with COVID-19 data in a commune in Chile. The results of the simulations indicate that the model generates important protection for the population by recognizing the super-propagating people (bridge nodes). This type of protection can be key in the fight against COVID-19.

Abstract

Among the diverse and important applications that networks currently have is the modeling of infectious diseases. Immunization, or the process of protecting nodes in the network, plays a key role in stopping diseases from spreading. Hence the importance of having tools or strategies that allow the solving of this challenge. In this paper, we evaluate the effectiveness of the DIL-W α ranking in immunizing nodes in an edge-weighted network with 3866 nodes and 6,841,470 edges. The network is obtained from a real database and the spread of COVID-19 was modeled with the classic SIR model. We apply the protection to the network, according to the importance ranking list produced by DIL-W α , considering different protection budgets. Furthermore, we consider three different values for α ; in this way, we compare how the protection performs according to the value of α .

1. Introduction

Infectious diseases have been the focus of multiple fields of research. In public health and epidemiology, efforts are directed at establishing transmission dynamics, the characteristics of infectious agents, and the populations most affected by pathogens, among others, which are of high importance for science [1]. In recent decades, research on infectious diseases has involved the application of complex theories from mathematics and engineering. In particular, the use of network models has allowed explanations of the spread of diseases from infected people (nodes) and their links with others (edges) [2].
Network models establish the connection between population groups, which is useful not only in the field of public health or epidemiology, but also in engineering and social sciences [3] (see [4,5,6,7,8]). In the health field, being a theoretical approach, the importance of recognizing the complexities of community structures has been discussed in order to understand social dynamics in the spread of infectious diseases [9]. For example, Magelinski et al. developed a model to estimate the role played by certain nodes in community structures [10], while Ghalmane et al. included the dimension of centrality in complex networks [11].
Regarding the COVID-19 pandemic, a series of models have been developed which allow the projection and establishment of the progress of this infectious disease based on the data available from different information sources (see, for instance, [12]). In this sense, Manríquez et al. propose the use of weighted graphs at the edges, giving the network model, from a stochastic approach, the possibility of identifying the most important variables in the spread of COVID-19 with real data in a city in Chile [13]. The same authors, in order to classify the importance of the nodes, propose the generalization of the measure of importance of the line with the use of degree of centrality (DIL-W α ), improving the understanding of the network from a local perspective [14].
On the other hand, the same network models have served to search for protection strategies for populations, generally related to the processes of immunization or isolation of people (see [15,16,17,18,19,20]).
Scientific evidence supports immunization strategies being useful in both homogeneous and heterogeneous networks [21]. The strategies seek, first, to develop a dynamic experimental model that includes the classic compartmental systems in epidemiology (SIR, SEIR, SIS, among others) and then to establish immunization measures across the nodes. Immunization measures can be applied in static and dynamic networks, in a random or targeted immunization [22].
Random immunization refers to random strategies, without determining a particular population in the protection process [23]. In contrast, targeted immunization recognizes nodes with a higher degree of connection with other nodes [24]. Nian et al. conclude, through computational simulations in a free scale graph (Model BA), that targeted immunization is more effective than randomized [21]. In the same direction, Wang et al. conclude that, even if the immunization strategy is imperfect or incomplete, it manages to generate positive impacts on the protection of the network [25]. Another investigation by Xia et al. determined that targeted immunization in two rounds of selection provides a greater protective effect compared to progressive strategies [24].
According to the above, a study conducted by Ghalmane et al. proposed an immunization strategy considering the influence of the nodes, the number of communities and the links between them [26]. Along the same lines, Gupta et al. analyzed the importance of protection strategies using information from community networks [27]. Since edge-weighted graphs have been an important way of understanding epidemic diseases, Manríquez et al. measure the effectiveness of DIL-W α with the recognition of bridge nodes. The effectiveness of DIL-W α is high compared to other proposals on four real networks [28].
For all the above, studies that include immunization measures together with community/local network models are highly relevant for public health, both for establishing measures to mitigate epidemic diseases and for the development of vaccination optimization policies to more quickly achieve herd immunization and, therefore, overcome diseases such as COVID-19. This is supported by Zhao et al., who argue that this framework allows for the optimization of immunization resources [29].
This study aims to analyze the protection effect against COVID-19 using the DIL-W α ranking with real data from a city in Chile (Olmué-City), obtained from the Epidemiological Surveillance System of the Ministry of Health of Chile, from which we obtain an edge-weighted graph, denoted by G E , according to the method proposed in [13]. We apply the protection to the G E network, according to the importance ranking list produced by DIL-W α , considering different protection budgets. For the ranking DIL-W α , we consider three different values for α ; they are 0, 0.5 and 1. In this way, we compare how the protection performs according to the value of α . We use a graph-based SIR model, namely, each individual is represented by a vertex in G E . At time t, each vertex v i is in a state v i t belonging to S = { 0 , 1 , 1 } , where 0 , 1 and 1 represent the three discrete states: Susceptible (S), Infected (I) and Recovered or Removed (R). Five hundred simulations were performed on G E ; the initial population contains one infected node and all the simulations, considering δ = 1 15 (recovered rate).
This paper is organized as follows: Section 2 contains generalities about graph theory, includes a graph from a database, and the DIL-W α ranking is explained. In Section 3, we obtain the graph from a real database from a city in Chile (Olmué-City), and we set the protection strategy. In Section 4, the results of the study are presented. Section 5 provides a discussion of the results and potentialities of the method used. Finally, Section 6 provides the conclusions.

2. Basic Definitions

In this section, we establish the definitions and elements used throughout this paper. We summarize the symbols and notations in Table 1.

2.1. Graphs

The following definitions come from [30,31].
Definition 1.
A graph G is a finite nonempty set V of objects called vertices, together with a possibly empty set E of 2-element subsets of V called edges.
To indicate that a graph G has vertex set V and edge set E, we write G = ( V , E ) . If the set of vertices is V = { v 1 , v 2 , , v n } , then the edge between vertex v i and vertex v j is denoted by e i j .
If e i j is an edge of G, then v i and v j are adjacent vertices. Two adjacent vertices are referred to as neighbors of each other. The set of neighbors of a vertex v is called the open neighborhood of v (or simply the neighborhood of v) and is denoted by N ( v ) . If e i j and e j k are distinct edges in G, then e i j and e j k are adjacent edges.
Definition 2.
The number of vertices in a graph G is the order of G and the number of edges is the size of G.
Definition 3.
The degree of a vertex v in a graph G, denoted by d e g ( v ) , is the number of vertices in G that are adjacent to v. Thus, the degree of v is the number of vertices in its neighborhood N ( v ) .
Definition 4.
Let G be a graph of order n, where V ( G ) = { v 1 , v 2 , , v n } . The adjacency matrix of G is the n × n zero-one matrix A ( G ) = [ a i j ] , or simply A = [ a i j ] , where
a i j = 1 i f e i j E ( G ) 0 i f e i j E ( G ) .
On the other hand, an important generalization of the simple graph consists of the definition of a weighted graph, more specifically an edge-weighted graph. Informally, an edge-weighted graph is a graph whose edges have been assigned a weight.
Definition 5.
An edge-weighted graph is a pair ( G , W ) , where G = ( V , E ) is a graph and W : E R is a weight function. If e i j E then W ( e i j ) = w i j .
Definition 6.
The strength of a vertex v i , denoted by S ( v i ) , is defined as the sum of the weights of all edges incident to it, this is to say,
S ( v i ) = v j N ( v i ) w i j .
The following definition comes from [32].
Definition 7
(Degree centrality [32]). The degree centrality of v i V of an edge-weighted graph ( G , w ) , denoted by C D w α ( v i ) , is defined as
C D w α ( v i ) = d e g ( v i ) ( 1 α ) · S ( v i ) α ,
where α [ 0 , 1 ] .
The parameter α is called the tuning parameter. Notice that, when α = 0 , then C D w α ( v i ) = d e g ( v i ) and, when α = 1 , then C D w α ( v i ) = S ( v i ) .

2.2. DIL-W α Ranking

We briefly describe the DIL-W α ranking in this Section. The DIL ranking is a tool for evaluating the node importance based on degree and the importance of lines (DIL) proposed by Liu et al. in [33] for an undirected and unweighted network. Recently, Manríquez et al. in [14] propose DIL-W α rank. This ranking method of node importance for undirected and edge-weighted is a generalization of the measure of line importance (DIL) based on the centrality degree (Definition 7) proposed by Opsahl in [32].
The following comes from [14].
Let us consider an undirected weighted graph ( G , w ) with G = ( V , E ) and V = { v 1 , v 2 , , v n } .
Definition 8
(Importance edge [14]). The importance of an edge e i j E , denoted by I α ( e i j ) , is defined as
I α ( e i j ) = C D w α ( v i ) p i α · C D w α ( v j ) p j α λ α ,
where, for k { i , j } , p k α = ( p + 1 ) ( 1 α ) · t k α with p being the number of triangles, one edge of the triangle is e i j , t k α is the weight of the sum of the edges incident to v k that form a triangle with e i j and λ α = p ( 1 α ) · t i + t j α 2 + 1 .
In order to illustrate the above Definition, let us consider the edge-weighted graph in Figure 1. Moreover, we consider the edges e 78 and e 75 . Notice that they both have the same weight (three). For this example, we set α = 1 . Applying Definition 7, we get:
C D w 1 ( v 7 ) = d e g ( v 7 ) ( 1 1 ) · S ( v t ) 1 = 1 · 12 = 12 and C D w 1 ( v 5 ) = d e g ( v 5 ) ( 1 1 ) · S ( v t ) 1 = 1 · 8 = 8 .
From Definition 8:
p 7 1 = ( 1 + 1 ) ( 1 1 ) · t 7 α = 2 0 · 2 1 = 2 , p 5 1 = ( 1 + 1 ) ( 1 1 ) · t 5 α = 2 0 · 1 1 = 1 , and λ 1 = 1 ( 1 1 ) · 2 + 1 1 2 + 1 = 5 2 .
Therefore,
I 1 ( e 75 ) = C D w 1 ( v 7 ) p 7 1 · C D w 1 ( v 5 ) p 5 1 λ 1 = ( 12 2 ) · ( 8 1 ) 3 2 + 1 = 28 .
In the same way with the edge e 78 , we obtain
I 1 ( e 78 ) = 96 .
In conclusion, edge e 68 is more important than edge e 75 .
The latter is reasonable because the edge e 78 is a bridging edge of the graph.
Definition 9
(Contribution [14]). The contribution that v i V makes to the importance of the edge e i j , denoted by W α ( e i j ) , is defined as
W α ( e i j ) = I α ( e i j ) · C D w α ( v i ) w i j α C D w α ( v i ) + C D w α ( v j ) 2 w i j α ,
where w i j is the weight of e i j .
We have calculated the importance of the edge e 78 of the graph in Figure 1. The contribution that v 7 makes to it is given by Definition 9:
W 1 ( e 78 ) = I 1 ( e 78 ) · C D w 1 ( v 7 ) w 78 1 C D w 1 ( v 7 ) + C D w 1 ( v 8 ) 2 w 78 1 = 96 · 12 3 12 + 8 2 · 3 = 432 7 .
In the same way, the contribution that v 8 makes to I 1 ( e 78 ) is:
W 1 ( e 87 ) = I 1 ( e 78 ) · C D w 1 ( v 8 ) w 78 1 C D w 1 ( v 7 ) + C D w 1 ( v 8 ) 2 w 78 1 = 96 · 8 3 12 + 8 2 · 3 = 240 7 .
The above means that the node v 7 contributes more to the edge e 78 than node v 8 .
Definition 10
(Importance of vertex DIL-W α [14]). The importance of a vertex v i V , denoted by L α ( v i ) , is defined as
L α ( v i ) = C D w α ( v i ) + v j N ( v i ) W α ( e i j ) .
Remark 1.
From the definition of Degree centrality (Definition 7) proposed by Opsahl in [32], we can see that, when the tuning parameter α is 0, the Definitions 8–10 are the same than the proposed by Liu et al.
In order to illustrate the above Definition, we compute the importance of v 7 and v 8 in the graph of Figure 1.
L 1 ( v 7 ) = C D w 1 ( v 7 ) + v j N ( v 7 ) W 1 ( e 7 j ) = 12 + W 1 ( e 71 ) + W 1 ( e 72 ) + W 1 ( e 73 ) + W 1 ( e 74 ) + W 1 ( e 75 ) + W 1 ( e 78 ) = 12 + 12 + 12 + 24 + 75 7 + 18 + 432 7 = 78 + 507 7 = 1053 7 ,
and
L 1 ( v 8 ) = C D w 1 ( v 8 ) + v j N ( v 8 ) W 1 ( e 8 j ) = 8 + W 1 ( e 87 ) + W 1 ( e 8 , 10 ) = 8 + 240 7 + 136 5 = 8 + 2152 35 = 2432 35 .
Since L 1 ( v 7 ) < L 1 ( v 8 ) , then node v 7 is more important than node v 8 (according to DIL-W 1 ranking).

2.3. Graph from a Database

The authors of [13] provide a way to obtain an edge-weighted graph from a database, which we briefly detail.
Let V = v 1 , v 2 , , v N be a set of people registered in a database, denoted by E , with K different variables, denoted by X k . These variables are separated into two categories: the characteristic variables (CHAR) and the relationship variables (REL) (which are those that allow us to assume that some person meets another). Let us denote by K 1 the number of relationship variables and E P I ( i , k ) the response of the person v i to the variable X k .
Definition 11.
We will say that a person v i is related to a person v j if and only if there exists X k R E L for k { 1 , 2 , , K 1 } such that E P I ( i , k ) = E P I ( j , k ) and i j .
To define the weight of each link between two persons, we assume that each X R E L has an associated inherent weight, this is to say, it is possible to discriminate some hierarchical order between the variables. Let p k be the weight associated to the variable X k R E L for k = 1 , , K 1 .
Definition 12.
We will say that for X j , X t R E L , X j is related to X t , denoted by X j R X t , if and only if p j = p t .
Definition 13.
Let A 1 , A 2 , , A c be the different classes that are defined by the different weights p 1 , p 2 , , p c and α 1 , α 2 , , α c and its respective cardinalities. Hence,
p j = α j K 1 ,
for all j { 1 , 2 , , c } .
We denoted by h i , j the number of times that one person is related to another (or the number of variables that matches between them).
Definition 14.
Let v i , v j V be such that v i is related to v j and p k r is the weight of the variable in which v i and v j match, for r = 1 , , h i , j . We will say that
w ˜ i j = r = 1 h i , j p k r ,
is the weight of the link between v i and v j .
Finally, the weighted adjacency matrix, which defines the graph obtained from the database, is the n × n matrix A ( G ) = [ a i j ] , where
a i j = w ˜ i j if e i j E ( G ) 0 if e i j E ( G ) .
Example 1.
In the following example, Table 2 simulates a database with 20 registered people. The data hosted correspond to the city in which they live (City), the workplace (considering school and university as a workplace), gender (Gen.), age, extracurricular activity (EC activity), address, whether they drink alcohol (Drin.), whether they are smokers (Sm.) and marital status (MS). Let us consider A and B as two different cities, and x , y , z , w , u , v , r , s , q , t , p , k , d , g and h as different people’s addresses. Moreover, in the table, Y = Yes, N = No, IC = in couple, M = married, S = single, W = widower.
From Table 2, we have that E P I = { X 1 , X 2 , X 3 , X 4 , X 5 , X 6 , X 7 , X 8 , X 9 } , where X 1 = City, X 2 = Workplace, X 3 = E.P. activity, X 4 = Address, X 5 = Sm., X 6 = Dri., X 7 = Gen., X 8 = M.S. and X 9 = Age. Then, we obtain the sets:
1. 
R E L = { X 1 , X 2 , X 3 , X 4 } and
2. 
C H A R = { X 5 , X 6 , X 7 , X 8 , X 9 } .
In our criteria, the hierarchical order of the variables X 1 , X 2 , X 3 , X 4 in descending form is X 4 , X 2 , X 3 , and X 1 . Moreover, we consider that the variables X 4 and X 2 have the same weight. Hence, A 1 = { X 2 , X 4 } , A 2 = { X 3 } , and A 3 = { X 1 } are the different classes that are defined by the different weights. Hence, by Definition 13
p 1 = 1 2 , p 2 = 1 4 , p 3 = 1 4 .
To construct the graph, we must resort to Definition 11. For instance, person 17 is related to all the people who live in city A or who work at Workplace 8 or who have music as an extra curricular activity or whose address is k. With respect to the weights of the edges, Equation (6) in Definition 14 gives us the answer. For instance, person 6 matches person 11 in the answers of the variables X 1 and X 2 , this is to say, both people live in city A and have the same workplace. Then, the edge v 6 v 11 has weight w 6 11 = 0.5 + 0.25 = 0.75 . Figure 2 shows the obtained graph.

3. Method

The data that are modeled correspond to the city of Olmué (Valparaíso region, Chile) and were obtained from the database of the Epidemiological Surveillance System of the Ministry of Health of Chile, which included the notified cases (positive or negative) and their contacts from 3 March 2020 to 15 January 2021 with a total of 3866 registered persons.
We denote by E p i the database of the Epidemiological Surveillance System of the Ministry of Health of Chile. From the total of variables included in E p i ( K = 279 ) 7 of them are relationship variables ( K 1 = 7 ). They are: full address ( X 1 ); the street where the people live ( X 2 ); town ( X 3 ); place of work ( X 4 ); workplace section ( X 5 ); health facility where they were treated ( X 6 ) and the region of the country where the test was taken to confirm, or not, the contagion ( X 7 ).
In our criteria, the hierarchical order of the seven variables in descending form is X 1 , X 2 , X 3 , X 4 , X 5 , X 6 , X 7 . Moreover, we consider that the variables X 1 , X 2 and X 3 have the same weight. In the same way, we also consider the variables X 4 and X 5 with equal weight. Hence, A 1 = { X 1 , X 2 , X 3 } , A 2 = { X 4 , X 5 } , A 3 = { X 6 } and A 4 = { X 7 } are the different classes that are defined by the different weights. Hence, by Definition 13.
p 1 = 3 7 , p 2 = 2 7 , p 3 = 1 7 , p 4 = 1 7 .
Figure 3 shows the obtained graph.
Let us denote by G E the graph obtained from database E p i .

Strategy Protection

In this section, we provide definitions of the protection of a graph when disease spreads on it. Moreover, we state the protection strategy used in the graph G E obtained in the previous Section. The following definitions come from [28].
Definition 15.
Protecting a vertex means removing all of its corresponding edges. (See Figure 4).
It is also possible to find in the literature that protecting a vertex means removing the vertex from the graph. See, for instance, [20].
Definition 16.
The number of vertices that are allowed to protect is called the protection budget, denoted by k.
Definition 17.
We will say that the survival rate, denoted by σ, is the ratio of vertices that remain uninfected at the end of the disease over the total numbers of vertices.
Therefore, our problem is: given a graph G = ( E , V ) , SIR model, and a protection budget k, the goal is to find a set of vertices S V , such that
θ = arg max S V σ ,
with | S | = k . However, the problem (4) is NP-Hard (see [34]).
Our chosen protection strategy corresponds to the DIL-W α ranking (see [14]). It is well known that an index to measure the connection of a graph is the efficiency of the networks (see [35]). High connectivity of the graph indicates high efficiency. In [14], the authors show that the DIL-W α ranking provides good results regarding the rate of decline in network efficiency (for more detail see [36]), when it comes to eliminating the best positioned nodes by this ranking. One of the good qualities of the DIL-W α ranking is that it recognizes the importance of bridge nodes (see more in [37]). This quality is inherited from the version of the DIL ranking for graphs not weighted at the edges (see [33]). Furthermore, [28] evaluated the effectiveness of the DIL-W α ranking in the immunization of nodes that are attacked by an infectious disease that spreads on an edge-weighted graph using a graph-based SIR model.
Finally, in order to illustrate in a simple way why the DIL-W α ranking has been chosen, let us consider the graph of Figure 5 with 16 edges, 15 nodes, and the respective weights on the edges. When we apply the DIL-W 1 ranking, the first 3 places are occupied by nodes 3, 5 and 1, respectively. These nodes are precise bridge nodes and, when protecting them, according to Definition 15, the graph loses connectivity (see Figure 6). If we apply the Strength ranking, the first three places are occupied by nodes 3, 1 and 4, respectively. Note that the order in which it positions the nodes and the importance it gives to node 4 makes the loss of network connectivity lower than the loss when applying DIL-W 1 (see Figure 6).
In summary, we apply the protection to the G E network, according to the importance ranking list produced by DIL-W α , considering different protection budgets. For the ranking DIL-W α , we consider three different values for α ; they are 0, 0.5 and 1. In this way, we compare how the protection performs according to the value of α .

4. Results

In this paper, we use a graph-based SIR model in the same way as in [13,28], namely, each individual is represented by a vertex in G E . At time t, each vertex v i is in a state v i t belonging to S = { 0 , 1 , 1 } , where 0 , 1 and 1 represent the three discrete states: Susceptible (S), Infected (I) and Recovered or Removed (R). We set
N I ( v i ) = v N ( v i ) : v I .
At time t + Δ t , the vertex v i will change state according to probabilistic rules:
  • The probability ( P I ( v i ) ) that a susceptible vertex v i is infected by one of its neighbors is given by
    P I ( v i ) = v j N I ( v i ) ρ Δ t · w i j ,
    where ρ is a purely biological factor and representative of the disease and w i j is the weight of the edge e i j .
  • The probability ( P R ( v i ) ) that an infected vertex v i at time t will recover is given by
    P R ( v i ) = δ Δ t ,
    where δ is the recovery rate.
Moreover, we assume that the disease is present for a certain period of time and that, when individuals recover, they are immune, that is, reinfection is not considered.
The initial population contains one infected node and all the simulations that consider δ = 1 15 . Five hundred simulations were performed on G E with ρ = 0.00121 . Figure 7 shows the average infected curve and the real infected data in E p i . Moreover, it shows a curve fitted to the data following the SIR model; for this, we used the classic method of least squares to compare with our proposal.
The graph G E was protected with different protection budgets according to the importance of the DIL-W 0 , DIL-W 0.5 and DIL-W 1 rankings. Protection is carried out in week 1, this is to say, at the beginning of the spread of the disease. Figure 8 shows the results.
We can see the survival rate in Figure 9.
Figure 10 shows the relationship between the real infected (450 people) and those immunized according to our proposal.
We can see that 80% of the real infected are located in 60% of the top ranked according to DIL-W α . We think that this is a way to recognize those who will get sick; however, it is not the solution.
Another element that we have considered investigating is the time at which the protection takes place. We modified the protection in the graph as the weeks advanced. In Figure 11, we can see the different infected curves, considering the 10% protection according to the DIL-W α ranking.
Figure 12 shows the relationship between the survival rate and the week in which the protection is carried out with our proposal.
The survival rate is clearly decreasing.

5. Discussion

The results of the present investigation are directed towards the analysis of the effectiveness of immunization using the DIL-W α ranking with real COVID-19 data from the city of Olmué-in Chile. Depending on the importance of the rankings, the immunization results were similar, despite the percentage of protection proposed in the simulations. Our method, therefore, goes in the direction of finding new optimization algorithms in network protection strategies [38].
At the level of protection, it is evident that when the percentage of initial coverage is higher, the epidemic ends with a smaller fraction of people affected by COVID-19. This event is related not only to the random increase in immunization, but also to the possibility, in the model, of recognizing the bridge nodes to increase the effectiveness of vaccination. This is consistent with other investigations that indicate that the recognition of central nodes or high-risk individuals improves the efficiency of immunization strategies in real networks, a situation that favors the protection of the network and the best use of vaccine doses [39,40]. The best use of doses is a challenge for the current scenario of vaccine shortages worldwide, mainly in poor nations [41].
On the other hand, the level of effectiveness of the DIL-W α ranking, given the percentage of protection, is established in the recognition of the bridging nodes in a regular vaccination process. The results using the α parameters of DIL-W α indicate a high survival rate. DIL-W 1 achieves better results with 70% protection and is positioned with the best survival rate, but DIL-W 0 and DIL-W 0.5 show good results. The difference between the different values of α is marginal and can be explained by the adequate representation that the DIL-W α model has and by the values of α , which do not generate excessive differences in the ranking. This is similar to the results of the research by Ophsal et al. who, through Freeman’s EIES network, mention that the centrality degree (Definition 7) is relatively stable among the different α parameters [32].
In a real and regular immunization strategy situation, such as the administration of vaccines, determining the population that infects most frequently is relevant since it allows optimization of these processes. Among our findings, it stands out that 80% of the real infected in the Olmué-Chile commune were located in 60% of the top of the DIL-W α ranking. Consequently, our proposal recognizes the heterogeneity of the network, approaching the reality of human interactions and achieving similar results in complex homogeneous networks [40].
Regarding immunization with 10% protection, a decrease in the survival rate is established by 4% from weeks 5 to 45 of protection. Likewise, with the same percentage of protection, the effectiveness of the immunity strategy tends to be important until week 20. After 20 weeks, the fraction of infected is similar with or without protection. Consequently, our model is strongly effective as a measure of rapid recognition of the epidemic outbreak in a given territory.
Therefore, according to the findings of our research, there are two important variables for the success and effectiveness of immunity strategies against COVID-19: (1) Recognition of bridging nodes (people with the highest probability of contagion) to apply measures of protection; and (2) the development time of this strategy.
Regarding the recognition of bridging nodes, there is evidence to support that targeted immunization schemes significantly reduce epidemic outbreaks [42]. This opens the possibility of changing the traditional perspective of immunization by protecting a small proportion of the population over a long period of time [43]. It is important, therefore, not only to direct COVID-19 immunity efforts towards the population most affected by mortality, but also in those population groups that tend to infect with greater force.
The time of development of the immunization strategy continues to be a variable under discussion in the scientific community regarding the slowness worldwide of the vaccination process, which risks not achieving herd immunity [44]. In summary, both at a theoretical and empirical level, the execution time of immunization strategies is important in overcoming the COVID-19 pandemic.
Finally, our model helps to establish a ranking of bridge nodes in a non-homogeneous network, so it is highly replicable with real COVID-19 dissemination data and it is useful to establish more focused strategies given the reduced number of vaccines available.

6. Conclusions

In this paper, we evaluate the effectiveness of the DIL-W α ranking in the immunization of nodes that are attacked by an infectious disease (COVID-19) that spreads on an edge-weighted graph obtained from the database of the Epidemiological Surveillance System of the Chilean Ministry of Health, using a graph-based SIR model.
Considering survival rates, the DIL-W 1 ranking performs better (by a small margin) than DIL-W 0.5 and DIL-W 0 rankings, subject to the protection budget being equal to 10% of the network nodes.
The period in which immunization or protection is given plays a key role in stopping the spread of the disease (see Figure 11) since around week 25 immunization does not generate a great impact and as time progresses the survival rate decreases almost linearly.
An interesting and complex task to solve is to determine which value of α to choose in the network so that the ranking generated is the optimal one. The same value does not always make the performance the best. One way to explore this is to continue with the ideas proposed in [45], where the selection standard of the optimal turning parameters is proposed for the centrality degree, but is not for DIL-W α ranking. However, when considering this method, there are as many rankings as there are numbers between 0 and 1.

Author Contributions

The authors R.M., C.G.-N. and C.T. contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported and funded by Agencia Nacional de Investigación y Desarrollo (ANID), COVID0739 project.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study were available after being requested by research project COVID-ANID to the Chilean Ministry of Health. The data are not publicly available due to legal restrictions. The source codes are accessible through the Github link: https://github.com/RonaldManriquez/Prot-stra-aga-epi-ewgraphs.git, from 8 June 2021.

Conflicts of Interest

The authors declare that they have no conflict of interest. The funders had no role in the design of the study; in the collection, analysis, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Epidemiology is a science of high importance. Nat. Commun. 2018, 9, 1703. [CrossRef] [Green Version]
  2. Lloyd, A.L.; Valeika, S. Network models in epidemiology: An overview. In Complex Population Dynamics: Nonlinear Modeling in Ecology, Epidemiology and Genetics; World Scientific: London, UK, 2007; pp. 189–214. [Google Scholar]
  3. Haslbeck, J.M.; Waldorp, L.J. How well do network models predict observations? On the importance of predictability in network models. Behav. Res. Methods 2018, 50, 853–861. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. An, X.L.; Zhang, L.; Li, Y.Z.; Zhang, J.G. Synchronization analysis of complex networks with multi-weights and its application in public traffic network. Phys. A Stat. Mech. Appl. 2014, 412, 149–156. [Google Scholar] [CrossRef]
  5. Montenegro, E.; Cabrera, E.; González, J.; Manríquez, R. Linear representation of a graph. Bol. Soc. Parana. MatemÁTica 2019, 37, 97–102. [Google Scholar] [CrossRef] [Green Version]
  6. Wang, G.; Cao, Y.; Bao, Z.Y.; Han, Z.X. A novel local-world evolving network model for power grid. Acta Phys. Sin. 2009, 6, 58. [Google Scholar]
  7. Mersch, D.P.; Crespi, A.; Keller, L. Tracking individuals shows spatial fidelity is a key regulator of ant social organization. Science 2013, 340, 1090–1093. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Firth, J.A.; Sheldon, B.C. Experimental manipulation of avian social structure reveals segregation is carried over across contexts. Proc. R. Soc. B Biol. Sci. 2015, 282, 20142350. [Google Scholar] [CrossRef] [Green Version]
  9. Wu, X.; Liu, Z. How community structure influences epidemic spread in social networks. Phys. A Stat. Mech. Appl. 2008, 387, 623–630. [Google Scholar] [CrossRef]
  10. Magelinski, T.; Bartulovic, M.; Carley, K.M. Measuring Node Contribution to Community Structure With Modularity Vitality. IEEE Trans. Netw. Sci. Eng. 2021, 8, 707–723. [Google Scholar] [CrossRef]
  11. Ghalmane, Z.; Cherifi, C.; Cherifi, H.; El Hassouni, M. Centrality in complex networks with overlapping community structure. Sci. Rep. 2019, 9, 1–29. [Google Scholar] [CrossRef] [Green Version]
  12. Guerrero-Nancuante, C.; Manríquez P, R. An epidemiological forecast of COVID-19 in Chile based on the generalized SEIR model and the concept of recovered. Medwave 2020, 20, e7898. [Google Scholar] [CrossRef]
  13. Manríquez, R.; Guerrero-Nancuante, C.; Martínez, F.; Taramasco, C. Spread of Epidemic Disease on Edge-Weighted Graphs from a Database: A Case Study of COVID-19. Int. J. Environ. Res. Public Health 2021, 18, 4432. [Google Scholar] [CrossRef]
  14. Manríquez, R.; Guerrero-Nancuante, C.; Martínez, F.; Taramasco, C. A generalization of the importance of vertices for an undirected weighted graph. Symmetry 2021, 13, 902. [Google Scholar] [CrossRef]
  15. Tang, P.; Song, C.; Ding, W.; Ma, J.; Dong, J.; Huang, L. Research on the node importance of a weighted network based on the k-order propagation number algorithm. Entropy 2020, 22, 364. [Google Scholar] [CrossRef] [Green Version]
  16. Wang, B.; Chen, G.; Fu, L.; Song, L.; Wang, X. Drimux: Dynamic rumor influence minimization with user experience in social networks. IEEE Trans. Knowl. Data Eng. 2017, 29, 2168–2181. [Google Scholar] [CrossRef]
  17. Zhang, Y.; Prakash, B.A. Dava: Distributing vaccines over networks under prior information. In Proceedings of the 2014 SIAM International Conference on Data Mining, Philadelphia, PA, USA, 24–26 April 2014; SIAM: Philadelphia, PA, USA, 2014; pp. 46–54. [Google Scholar]
  18. Hébert-Dufresne, L.; Allard, A.; Young, J.G.; Dubé, L.J. Global efficiency of local immunization on complex networks. Sci. Rep. 2013, 3, 1–8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Wijayanto, A.W.; Murata, T. Effective and scalable methods for graph protection strategies against epidemics on dynamic networks. Appl. Netw. Sci. 2019, 4, 1–31. [Google Scholar] [CrossRef] [Green Version]
  20. Tong, H.; Prakash, B.A.; Tsourakakis, C.; Eliassi-Rad, T.; Faloutsos, C.; Chau, D. On the Vulnerability of Large Graphs. In Proceedings of the 2010 IEEE International Conference on Data Mining, Sydney, NSW, Australia, 13–17 December 2010; pp. 1091–1096. [Google Scholar] [CrossRef]
  21. Nian, F.; Hu, C.; Yao, S.; Wang, L.; Wang, X. An immunization based on node activity. Chaos Solitons Fractals 2018, 107, 228–233. [Google Scholar] [CrossRef]
  22. Gómez-Gardenes, J.; Echenique, P.; Moreno, Y. Immunization of real complex communication networks. Eur. Phys. J. Condens. Matter Complex Syst. 2006, 49, 259–264. [Google Scholar] [CrossRef] [Green Version]
  23. Cohen, R.; Havlin, S.; Ben-Avraham, D. Efficient immunization strategies for computer networks and populations. Phys. Rev. Lett. 2003, 91, 247901. [Google Scholar] [CrossRef] [Green Version]
  24. Xia, L.; Song, Y.R.; Li, C.C.; Jiang, G.P. Improved targeted immunization strategies based on two rounds of selection. Phys. A Stat. Mech. Appl. 2018, 496, 540–547. [Google Scholar] [CrossRef]
  25. Wang, Y.; Xiao, G.; Hu, J.; Cheng, T.; Wang, L. Imperfect targeted immunization in scale-free networks. Phys. A Stat. Mech. Appl. 2009, 388, 2535–2546. [Google Scholar] [CrossRef]
  26. Ghalmane, Z.; El Hassouni, M.; Cherifi, H. Immunization of networks with non-overlapping community structure. Soc. Netw. Anal. Min. 2019, 9, 1–22. [Google Scholar] [CrossRef] [Green Version]
  27. Gupta, N.; Singh, A.; Cherifi, H. Community-based immunization strategies for epidemic control. In Proceedings of the 2015 7th International Conference on Communication Systems and Networks (COMSNETS), Bangalore, India, 6–10 January 2015; pp. 1–6. [Google Scholar]
  28. Manríquez, R.; Guerrero-Nancuante, C.; Taramasco, C. Protection Strategy for Edge-Weighted Graphs in Disease Spread. Appl. Sci. 2021, 11, 5115. [Google Scholar] [CrossRef]
  29. Zhao, D.; Wang, L.; Li, S.; Wang, Z.; Wang, L.; Gao, B. Immunization of epidemics in multiplex networks. PLoS ONE 2014, 9, e112018. [Google Scholar] [CrossRef] [PubMed]
  30. Chartrand, G.; Lesniak, L. Graphs and Digraphs, 1st ed.; CRC Press: London, UK, 1996. [Google Scholar]
  31. West, D.B. Introduction to Graph Theory, 2nd ed.; Prentice Hall: Englewood Cliffs, NJ, USA, 2001. [Google Scholar]
  32. Opsahl, T.; Agneessens, F.; Skvoretz, J. Node centrality in weighted networks: Generalizing degree and shortest paths. Soc. Netw. 2010, 32, 245–251. [Google Scholar] [CrossRef]
  33. Liu, J.; Xiong, Q.; Shi, W.; Shi, X.; Wang, K. Evaluating the importance of nodes in complex networks. Phys. A Stat. Mech. Appl. 2016, 452, 209–219. [Google Scholar] [CrossRef] [Green Version]
  34. Vinterbo, S.A. Privacy: A machine learning view. IEEE Trans. Knowl. Data Eng. 2004, 16, 939–948. [Google Scholar] [CrossRef]
  35. Latora, V.; Marchiori, M. Efficient Behavior of Small-World Networks. Phys. Rev. Lett. 2001, 87, 198701. [Google Scholar] [CrossRef] [Green Version]
  36. Ren, Z.; Shao, F.; Liu, J.; Guo, Q.; Wang, B.H. Node importance measurement based on the degree and clustering coefficient information. Acta Phys. Sin. 2013, 62, 128901. [Google Scholar]
  37. Musiał, K.; Juszczyszyn, K. Properties of bridge nodes in social networks. In Proceedings of the International Conference on Computational Collective Intelligence, Wroclaw, Poland, 5–7 October 2009; Springer: Berlin/Heidelberg, Germany, 2009; pp. 357–364. [Google Scholar]
  38. Shams, B. Using network properties to evaluate targeted immunization algorithms. Netw. Biol. 2014, 4, 74. [Google Scholar]
  39. Chen, Y.; Paul, G.; Havlin, S.; Liljeros, F.; Stanley, H. Finding a better immunization strategy. Phys. Rev. Lett. 2008, 101, 058701. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Nian, F.; Wang, X. Efficient immunization strategies on complex networks. J. Theor. Biol. 2010, 264, 77–83. [Google Scholar] [CrossRef] [PubMed]
  41. Wouters, O.J.; Shadlen, K.; Salcher-Konrad, M.; Pollard, A.J.; Larson, H.; Teerawattananon, Y.; Jit, M. Challenges in ensuring global access to COVID-19 vaccines: Production, affordability, allocation, and deployment. Lancet 2021. [Google Scholar] [CrossRef]
  42. Liu, Z.; Lai, Y.C.; Ye, N. Propagation and immunization of infection on general networks with both homogeneous and heterogeneous components. Phys. Rev. E 2003, 67, 031911. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Pastor-Satorras, R.; Vespignani, A. Immunization of complex networks. Phys. Rev. E 2002, 65, 036104. [Google Scholar] [CrossRef] [Green Version]
  44. Aschwanden, C. Five reasons why COVID herd immunity is probably impossible. Nature 2021, 591, 520–522. [Google Scholar] [CrossRef] [PubMed]
  45. Wei, D.; Zhang, Y.L.; Deng, Y. Degree centrality based on the weighted network. In Proceedings of the 24th Chinese Control and Decision Conference (CCDC), Taiyuan, China, 23–25 May 2012; pp. 3976–3979. [Google Scholar]
Figure 1. Simple graph.
Figure 1. Simple graph.
Biology 10 00667 g001
Figure 2. Graph obtained from E .
Figure 2. Graph obtained from E .
Biology 10 00667 g002
Figure 3. Graph obtained from database of Olmué city, Chile, with 3866 vertices and 6,841,470 edges.
Figure 3. Graph obtained from database of Olmué city, Chile, with 3866 vertices and 6,841,470 edges.
Biology 10 00667 g003
Figure 4. The left side shows a graph without protection or infection. On the right side, we see an infected node (red) and a protected node (blue) in the same graph.
Figure 4. The left side shows a graph without protection or infection. On the right side, we see an infected node (red) and a protected node (blue) in the same graph.
Biology 10 00667 g004
Figure 5. Graph with 15 nodes, 16 edges, and the respective weights.
Figure 5. Graph with 15 nodes, 16 edges, and the respective weights.
Biology 10 00667 g005
Figure 6. In the left, the first 3 protected places generated by the DIL-W 1 ranking. On the right, the first 3 protected places generated by Strength ranking.
Figure 6. In the left, the first 3 protected places generated by the DIL-W 1 ranking. On the right, the first 3 protected places generated by Strength ranking.
Biology 10 00667 g006
Figure 7. Real infected data (black), fitted curve (red), and infected curve obtained in the spread on G E (cyan).
Figure 7. Real infected data (black), fitted curve (red), and infected curve obtained in the spread on G E (cyan).
Biology 10 00667 g007
Figure 8. Infected curves obtained according to different values of k.
Figure 8. Infected curves obtained according to different values of k.
Biology 10 00667 g008
Figure 9. Survival rate of each ranking.
Figure 9. Survival rate of each ranking.
Biology 10 00667 g009
Figure 10. Relationship between the real infected and those immunized according to the different DIL-W α ranking.
Figure 10. Relationship between the real infected and those immunized according to the different DIL-W α ranking.
Biology 10 00667 g010
Figure 11. The different infected curves, considering the 10% protection according to the DIL-Wαranking and carried out in different weeks.
Figure 11. The different infected curves, considering the 10% protection according to the DIL-Wαranking and carried out in different weeks.
Biology 10 00667 g011aBiology 10 00667 g011b
Figure 12. Survival rate of each ranking.
Figure 12. Survival rate of each ranking.
Biology 10 00667 g012
Table 1. Summary of Symbols and Notations.
Table 1. Summary of Symbols and Notations.
NotationsDefinition and Description
GGraph or network.
( G , w ) Edge-weighted graph.
v i Vertex or node.
N ( v i ) Neighborhood of a vertex v.
e i j Edge between vertex v i and vertex v j .
w i j Weight of the edge e i j .
d e g ( v i ) Degree of the vertex v i .
S ( v i ) Strength of the vertex v i .
α Real number. Tuning parameter.
C D w α ( v i ) Degree centrality of v i V of an edge-weighted graph ( G , w ) .
DIL-W α Ranking based on Degree and importance of line.
I α ( e i j ) Importance of edge e i j .
W α ( e i j ) Contribution that v i makes to the importance of the edge e i j .
L α ( v i ) The importance of a vertex v i .
E Database.
X k Variable of a database.
p k Weight of the variable X k .
kProtection budget (the number of nodes in graph G that can be protected).
σ Ratio of surviving nodes.
Table 2. Database E .
Table 2. Database E .
PersonCityWorkplaceE. C. ActivityAddressSm.Dri.Gen.M. SAge
1AWorkplace 1TheateryYYFIC35
2AWorkplace 3CinemayYYMIC35
3BSchool BFootballzNNFS10
4BWorkplace 1PhotographyxNNFM48
5AWorkplace 5Does not haveuYNFW65
6AWorkplace 4Does not havevYYMIC27
7BWorkplace 2Does not havexYNMM46
8AUniversity 1PhotographyvNNMIC29
9AUniversity 2Does not havewYYMIC19
10BSchool BKaratexNNMS10
11AWorkplace 4Ping-pongrYYFM54
12ASchool AFootballsNNMS8
13AWorkplace 5DancerYYFM57
14BSchool AHandballqNNMS11
15AUniversity 1Does not havetNNFS25
16AWorkplace 7SingingpYYFS60
17AWorkplace 8MusickNYFS28
18AWorkplace 3Does not havedNNMS47
19BSchool AMusicgNNFS8
20AWorkplace 6Does not havehYYMS30
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Manríquez, R.; Guerrero-Nancuante, C.; Taramasco, C. Protection Strategy against an Epidemic Disease on Edge-Weighted Graphs Applied to a COVID-19 Case. Biology 2021, 10, 667. https://doi.org/10.3390/biology10070667

AMA Style

Manríquez R, Guerrero-Nancuante C, Taramasco C. Protection Strategy against an Epidemic Disease on Edge-Weighted Graphs Applied to a COVID-19 Case. Biology. 2021; 10(7):667. https://doi.org/10.3390/biology10070667

Chicago/Turabian Style

Manríquez, Ronald, Camilo Guerrero-Nancuante, and Carla Taramasco. 2021. "Protection Strategy against an Epidemic Disease on Edge-Weighted Graphs Applied to a COVID-19 Case" Biology 10, no. 7: 667. https://doi.org/10.3390/biology10070667

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