key: cord-0005445-udoioj4c authors: Nian, Fuzhong; Wang, Ke title: A two-step high-risk immunization based on high-risk immunization date: 2014-07-09 journal: Nonlinear Dyn DOI: 10.1007/s11071-014-1555-2 sha: ecec7746daa9bd284d8a39e02fa351b9a3581c54 doc_id: 5445 cord_uid: udoioj4c In this paper, the two-step high-risk immunization was investigated based on high-risk immunization for the SIRS model in small-world networks and scale-free networks. First, the effects of various immunization schemes are studied and compared. When the number of immune is same, the research result shows that the immune effect of the two-step high-risk immunization strategy is not the best nor the worst. However, the practicability is better compare with others. Furthermore, by changing the proportional of immunization the optimal immune effect can be achieved in the two-step high-risk immunization. Computation results verify that the two-step high-risk immunization is effective, and it is economic and feasible in practice. Epidemic is one of the most important issues related to our real lives. From SARS in 2003, to 2005 H1N1 flu, to today's H7N9, the epidemic problem has attracted a lot of attention. The dynamical behaviors of epidemic diseases have been studied for a long time; it also has been a hot topic [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] . For the problem of epidemic, the networks and its dynamics are some basic issue need to be solved firstly. The traditional study of networks includes lattices, regular tree, and ER random graph [15] and then, many real networks were investigated, such as the Internet, the www, the scientific web, and so on. The small-world phenomenon and scale-free property are two major statistical topological characteristics, and the corresponding network models were named as small-world network [16] [17] [18] and scale-free network [19] [20] [21] , respectively. So many epidemic spreading problems are investigated on the two networks. The epidemic system can be represented as a network where nodes stand for individuals, and edges stand for the relationships between them. Based on this hypothesis, several models have been proposed for epidemic dynamics. The typical Epidemic mode for the disease dynamics includes the susceptible-infected (SI) model [22] , susceptible-infected-recovered (SIR) model [23] , and susceptible-infected-recovered-susceptible (SIRS) model [24] , which represent the development of each individual's disease in a network. In this paper, the SIRS model was selected as our study object. In this model, an infected individual (I) is infective for some period of time and then is recovered (R), and no longer be infected by other infected individuals (I), and then he (she) would lose the immunization and become a susceptible individual again. The purpose of immune strategy is protecting the part of network nodes would not be infected by cutting off the propagation path of virus. It is one of effec-tive methods to prevent the virus from spreading on the population. In fact, an efficient immunization strategy can not only be applied to the computer network, but also have an important significance to the diseases spread and control in the human social networks. Currently, various immunization strategies have been proposed, for example: In 2000, Callway proposed a random immunization [25] which the nodes to be vaccinate are selected randomly, while in order to control the spread of the disease we almost vaccinate to all individuals. Obviously, it is not practical. Lately, Pastor-Satorras and Vespignani proposed targeted immunization schemes [26] , which immune the large degree nodes firstly. Although there are good effects, but it needs to obtain the whole network topology structure information. In response to this issue, Cohen designs an efficient immunization strategy named "acquaintance immunization" [27] . In this scheme, first randomly select some nodes, then randomly immune one of the nodes from their acquaintances. It is evident that the nodes with large degree have more chance to be vaccinated. In recent years, some other immune strategies have been proposed. For example: J. Gómez-Gardenes put forward on top of real communication networks the immune strategy [28] based on covering algorithm, the effect is basically consistent with the targeted immunization strategy. Madar and Kalisky also proposed the double acquaintance immunization [29] that first randomly selected some nodes, and then randomly immune their two acquaintances. The double acquaintance immunization shows a higher fraction of endemic states than the targeted, but its fraction of infected individuals is still lower than that obtained with the acquaintance immunization strategy. And we proposed a high-risk immunization strategy [30] , which the individuals with infected neighbors were looked as highrisk groups, and by immunizing the high-risk groups to controls the spread of the disease. However, the neighbor's neighbors were not considered in our last work. In this paper, the two-step high-risk immunization strategy was designed based on high-risk immunization. In this scheme, not only the high-risk neighbors were immunized, but also the indirect neighbors of infected individual were immunized too. The rest of the paper is organized as follows: in Sect. 2, the two-step high-risk immunization in different WS small-world complex networks is investigated, and the standard SIRS model is built. In Sect. 3, we investigate the behaviors of density of infected indi-viduals on different BA scale-free networks, and theoretical analysis and computer simulations are given. Finally, conclusions are given in last section. Suppose there are N nodes in the WS small-world network, each node is symmetrical connected to the nearest K nodes with probability p randomly re-connected to the each edge in the original network, and then the average degree of the network is k = 2K . There are three node statuses: S represents susceptible state, I represents infected state, and R represents removed state. S state nodes are infected with a certain rate of β, I state nodes are cured with a cure rate γ , while the R state nodes are loss of immunity with a certain probability σ , then the model's effective communication rate is λ = β γ . Fig. 1 is the sketch of SIRS. Like the high-risk immunization [30] , define is the probability that an arbitrary given node is a particular node's neighbor. The probability of an arbitrary node for access to immunization through vaccination is δ; at the same time suppose p is the proportion that immunes the direct neighbors of high-risk people, and q is the proportion that immune the indirect neighbors. In homogeneous network, the two-step high-risk immunization model can be described as follows: Here, S(t), I (t), and R(t), respectively, represent the average density of susceptible population, infected population, and removed population. In Eq. (1), the first S I R γ β σ Fig. 1 Sketch of SIRS term on the right side of the first equation considers the average density of an infected node infection produces new infections. The secondary term on the right side of the first equation considers the recovered nodes loss of immunity probability. Because S(t), I (t), and R(t) obey the normalization condition S(t) + I (t) + R(t) = 1, and the static conditions of no infectious diseases spread on networks are From the Eq. (1) and Eq. (2), we obtain From the secondary equation of Eq. (3), Substitute (4) into Eq. (5), we obtain Because I (t) = 0, the above equation can be simplified as So, the non-zero epidemic threshold is λ c = If the vaccinating rate δ = 0, obtain λ c = 1 k . From the results obtained by the above mathematical model, we can see that the non-zero epidemic threshold is associated with the average degree k of the network, the immune probability δ, the probability of a node with specific nodes which are neighbors to each other, and the numbers of the direct neighbors and indirect neighbors. When the probability of k , δ, and is fixed, we can obtain the diffident non-zero epidemic threshold by change the immune proportion of neighbors and indirect neighbors. Then we can achieve the optimal immune effect. Based on the analysis, the corresponding experiments are performed to further verify it. Our simulations are implemented at fixed N = 2, 000, k = 4, and the rewiring rate is 0.1. Initially, 10 % of the susceptible nodes in the network are infected. Suppose that the chance of being infected by its neighbor is equal for all nodes. In other words, in each iteration, the probability of being infected by its neighbor is the same. Simulations were computed averaging over 30 different starting configurations, performed on 20 different realizations of the network. Figure 2 is the comparison of random immunization, high-risk immunization, and the two-step highrisk immunization. Under the condition that the vaccinating number is same, the effect of the high-risk immu- Fig. 3a , we can see, when the proportion of immune direct neighbor nodes from 0.1 to 0.9 increased, and the proportion of immune indirect neighbor nodes from 0.9 to 0.1 decreased, the epidemic threshold λ increases at begin and then decreases. Figure 3b is changing of total number of immune nodes with p and q. Form Fig. 3b , we can see that the total numbers generally equaled when the immune proportion is changed. Here, p represents the proportion of immune direct neighbor node; q represents the proportion of immune indirect neighbor node. Figure 3 indicates the immune effect when the number of vaccinated nodes is same or similar. For the stronger infectious virus, the more intimate neighbors ( p = 0.9, q = 0.1) or more indirect neighbors ( p = 0.1, q = 0.9) can be vaccinated to achieve the optimal immune effect; for the not strong infectious virus, the spread of virus can be controlled effectively when evenly choose the number of direct and indirect neighbors individuals ( p = 0.5, q = 0.5). The characteristic of WS small-world model is the connectivity distribution can be approximately represented by the Poisson distribution. The distribution has a peak in the average value k , then rapid degradation exponentially, which means that when k k , it almost not exists a node which degree is K . So the network is also known as the Homogeneous network. In recent years, another major discovery is in many complex networks, including Internet, WWW, and the new supersedes the old network, connection degree distribution function has the power-law form. Because in these network the nodes' degree without significantly characteristic length, so called scale-free networks [19] . Barabási and Albert proposed the BA scale-free model based on the characteristic of growth and preferential attachment. We study the SIRS two-step high-risk immunization infectious disease model based on the scale-free network model is described as follows: constructs a BA scale-free network, it satisfies the growth and preferential connection characteristic, after iterating this scheme a sufficient number of times, we obtain a networks composed of N nodes with degree distribution P(k) = k −3 and degree k = 2m. Let S k (t), I k (t), and R k (t) be the densities of susceptible, infected, and recovered vertex of degree k. Then we have the following dynamics model: In Eq. (6), the first term on the rhs of the first equation considers the average density of newly infected nodes generated by each infected node. And the probability is (t) that any given link points to an infected node. The secondary term on the rhs of the first equation represents the nodes losing immunity with rate σ . The first term on the rhs of the last equation of Eq. (6) denotes the infected nodes recovering with unit rate. As we know, in scale-free networks a randomly chosen link is more likely to be connected to an infected node with high connectivity, yielding [31] : Also, S k (t), I k (t), and R k (t) obey the normalization condition for each k While non-epidemic stationary conditions can be given as follows: d S k (t) Combine above equation with Eq. (8), we obtain From the second equation of Eq. (6) Substitute (11) into (10) in the following. We obtain the following equation: From the formula, the solution I k (t) = 0 and (t) = 0 always satisfies the Eq. (12) . In order to get a non-zero static solution, to both sides of the equation (12) represents a function F( ), The range of this function is 0 < ≤ 1. Therefore, a nontrivial solution that exists provided is [30] Must be fulfilled. i.e., d d (14) is the probability of any given node which is a neighbor of some specific node. Therefore, As for BA scale-free network, Substitute Eqs. (15) and (16) into (14), we obtain λ k The value of λ yielding the inequality (17), i.e., From above equation, the critical epidemic threshold λ c , that is given by When σ = 1 and δ = 0, λ c = k k 2 Here, Substitute the above two equations into Eq. (18) In the continuous k approximation, calculate the approximation of k 1 k and k Here, M is the maximum degree, substitute Eqs. (20) and (21) into Eq. (19) 1 When M is big enough (when m is constant, From Eq. (22), we can see In our real life, the longer the time immune in our body, the harder the virus transmit in the crowd, namely there is a higher spreading threshold; while the shorter the time with the immune in our body, the easier the virus transmit in the crowd, which means the spreading threshold is lower. (2) λ c increases with δ increasing. δ is disease vaccination rate. If the vaccination rate is very high namely to have more individual acquired immunity, the disease is not easy to spread and the threshold higher. (3) λ c decreases with N increasing. The more population, the more easily outbreak infectious diseases, this is consistent with our real life. (4) λ c increases with p increasing or increases with q increasing, and exists a optimal proportion p and q of control the spread of the disease, makes λ c maximum. When immunizing the more direct neighbors, the disease is not easy to spread. When immunizing the more indirect neighbors, the disease is not easy to spread too. While when the total immunity number is same or similar, an optimal immune proportion can be achieved by changing the immune proportion of direct neighbors and indirect neighbors. In this case, Fig. 4 The comparison of random immunization, high-risk immunization, and the two-steps high-risk immunization in BA scale-free network (when N = 2, 000, k = 4, σ = 0.2 Infection density varies with propagation probability) the spread of virus can be most effectively controlled, and the propagation threshold reached the maximum. Based on the above analysis of mathematical theory, the corresponding experiments also were implemented to verification. The network was selected as scale-free network with N = 2, 000, average degree k = 4, and the initial node m = 5. The 0.1 * N infected nodes were selected randomly from the network. Suppose the chance of be infected by infected neighbor is equal to all nodes, and the probability of infection is the same. Simulations were computed averaging over 30 different starting configurations, performed on 20 different realizations of the network. Figure 4 is the comparison of random immunization, high-risk immunization, and the two-step highrisk immunization effect. Under the condition that the vaccinating number is same, the results are similar to the small-world network; the effect of the high-risk immunization is best, while the effect of the two-step high-risk immunization is not only slightly inferior to high-risk immunization but also better than random immunization. Figure 5a is the comparison of the immune effect when we change the immune proportion of direct neighbors and indirect numbers in the two-step highrisk immunization. In Fig. 5a , when the proportion of immune direct neighbor nodes from 0.1 to 0.9 increased, and the proportion of immune indirect neighbor nodes from 0.9 to 0.1 decreased, the epidemic threshold λ decreases at begin and then increases. Figure 5b is changing of total number of immune nodes with p and q. Form Fig. 5b , we can see that the total numbers change from 400 to 1,000 when the immune proportion is changed as above. This phenomenon is mainly determined by the characteristics of scale-free network itself. In the construction of scale-free networks there are two important features which are the growth and preferential attachment, so some degree of nodes is relatively large. This leads to a rapidly increase in the number of direct neighbors and indirect neighbors we looked for, which affect the immune quantity. So if a smaller proportion of direct neighbors and a larger proportion of indirect neighbors ( p = 0.1, q = 0.9) are vaccinated, it means that almost all the points are vaccinated in the network, at this time the propagation threshold is large. On the contrary when a larger proportion of direct neighbors and a smaller proportion of indirect neighbors ( p = 0.9, q = 0.1) are vaccinated, the number of immune nodes is relatively small, but the virus transmission route be cut off from the root, and the propagation threshold is large too; once the virus break through the first layer of immune defense, the spread of the virus will unbridle. This is also the reason why the infection density sudden increases when p = 0.9, q = 0.1 in Fig. 5a . Here, p represents the proportion of immune direct neighbor node; q represents the proportion of immune indirect neighbor node. Figure 5b is changing of total number of immune nodes with p and q. In Fig. 5b the total number of immune nodes is so large difference, while in practical application, the immune effect is expected to compare in similarity immune nodes. So immune number controlled as is shown in Fig. 6 . Figure 6a shows the variation of infection density with propagation probability when the proportion of immune direct neighbor nodes from 0.1 to 0.7 increased and the proportion of immune indirect neighbor nodes from 0.55 to 0.25 decreased. From Fig. 6b we can see that the total number changes from 500 to 600 generally not big fluctuations when the immune proportion is changed as above. Figure 6 indicates the immune effect when the number of vaccinated nodes is same or similar. In Fig. 6 , for the stronger infectious virus, more indirect neighbors can be immunized to active the optimal immune effect ( p = 0.1, q = 0.55); for not strong infectious virus, more direct neighbors can be immunized to control the spread of disease effectively ( p = 0.7, q = 0.25). In this paper, the two-step high-risk immunization strategy was proposed by experimenting in the WS small-world network model and BA scale-free network model; further the effects of various immunization schemes are compared; at the same time the application effect of the two-step high-risk immunization is discovered. According to the different strengths of infectious disease, the different immunization proportional can be taken to achieve the optimal immune effect. The results verified that the two-step high-risk immunization not only has its own advantages, but also is effective, economic, and feasible. Epidemic spreading with nonlinear infectivity in weighted scale-free networks Epidemic spreading of an SEIRS model in scale-free networks Evolution of scaling emergence in large-scale spatial epidemic spreading Epidemic spreading on interconnected networks Spreading dynamics and global stability of a generalized epidemic model on complex heterogeneous networks Dynamical interplay between awareness and epidemic spreading in multiplex networks Analysis of epidemic spreading of an SIRS model in complex heterogeneous networks Epidemic spreading on weighted complex networks An SIR model with infection delay and propagation vector in complex networks Analysis of stability and bifurcation for an SEIV epidemic model with vaccination and nonlinear incidence rate High-risk immunization strategies for multiethnic regions Pattern dynamics of an epidemic model with nonlinear incidence rate Pattern formation of an epidemic model with diffusion The existence of codimension-two bifurcation in a discrete SIS epidemic model with standard incidence On the evolution of random graphs Models of the small world Collective dynamics of 'smallworld'networks Renormalization group analysis of the small-world network model Emergence of scaling in random networks Scale-free networks: a decade and beyond Presented at the Proceedings of the 2012 International Conference on Communication, Electronics and Automation Engineering Global stability of an SI epidemic model with feedback controls Exact analytical solutions of the Susceptible-Infected-Recovered (SIR) epidemic model and of the SIR model with equal death and birth rates Global stability of a delayed SEIRS epidemic model with saturation incidence rate Network robustness and fragility: percolation on random graphs Immunization of complex networks Efficient immunization strategies for computer networks and populations Immunization of real complex communication networks Immunization and epidemic dynamics in complex networks Efficient immunization strategies on complex networks Epidemic spreading in scale-free networks