Research Article  Open Access
Epidemic Spreading in Complex Networks with Resilient Nodes: Applications to FMD
Abstract
At the outbreak of the animal epidemic disease, farms that recover quickly from partially infected state can delay or even suppress the wide spreading of the infection over farm networks. In this work, we focus on how the spatial transmission of the infection is affected by both factors, the topology of networks and the internal resilience mechanism of nodes. We first develop an individual farm model to examine the influence of initial number of infected individuals and vaccination rate on the transmission in a single farm. Based on such intrafarm model, the farm network is constructed which reflects disease transmission between farms at various stages. We explore the impact of the farms vaccinated at low rates on the disease transmission into entire farm network and investigate the effect of the control on hub farms on the transmission over the farm network. It is shown that intensive control on the farms vaccinated at low rates and hub farms effectively reduces the potential risk of footandmouth disease (FMD) outbreak on the farm network.
1. Introduction
The epidemic spreading in networks has attracted growing attention in recent years [1–3]. One of the main reasons for studying the spread process in networks is reproducing the actual dynamics of the disease and finding an effective control strategy to eradicate the infection. The spreading of infectious animal disease between farms exhibits the vulnerability of the community structure, which has been mostly investigated in terms of the network topology [4–6]. However, whether or not the infection is suppressed before gaining “momentum” to spread over networks also depends on the nodes’ responses to the epidemic [7]. If each farm at a node is vaccinated at a high level and can recover from accidental and partial infection, such resilient mechanism can delay or even suppress the wide spreading of the infection.
Among animal infectious diseases, footandmouth disease (FMD) is an economically important disease. Mathematical models have been developed to simulate the possible progress of FMD transmission under various scenarios, so that relevant researchers and agents can use them for implementing control policies to prevent and reduce the risk of the disease outbreak. A discretized deterministic SLIR model, where , , , and denote the susceptible, the latent, the infectious, and the recovered states, respectively, was developed to describe FMD transmission in a farm [8]. An SVLI model, where denotes the vaccinated state, was considered to investigate FMD transmission in a farm [9]. The authors investigated the stability of the diseasefree equilibrium and performed simulations to illustrate the impact of vaccination and culling on controlling the disease in a farm. Concerning FMD transmission between farms, Keeling et al., using their developed individual farmbased stochastic model, found that the spatial distribution and size of farms have a certain effect on the pattern and regional variability of FMD outbreaks [10] and investigated the effect of vaccination and culling on control of FMD transmission [11]. Moreover, optimal reactive vaccination strategies were investigated in [12] and various vaccination strategies were studied for stochastic SIR model on a random network of social contacts with household structure [13]. Farmtofarm contacts due to movements of operators and vehicles have long been acknowledged as a relevant factor in disease transmission in livestock systems [14–16]. It is observed that the spread of infections can be tremendously strengthened on such contact networks. Particularly, some sorts of complex networks that connect susceptible/infected nodes are prone to the epidemic spreading with the persistence of infections regardless of the spreading rates [17, 18]. The bond percolation on random networks has been used to derive the basic epidemiological quantities and to predict disease transmission. A brief overview of the compartment SIR model in different type of the contact network can be found in [19, 20] and a comprehensive review for modeling FMD transmission in and between farms is done in [21].
In this work, we want to find how much the mutual dependence of farms weakens their resistance to external infections. With intensive and accumulative livestock farming in the industrialized world, farms are more densely connected even in the distance. Farms that are closely linked and share the common resources are generally more vulnerable to disease than an isolated farm. We are going to derive the minimum vaccination level to protect a farm community from sporadic exposure to external infectious source. One of the key problems in the epidemiology is how to effectively control the transmission of infectious disease by immunization of the population. The most undesirable situation is internalization or localization of the disease even under regular practice of vaccination. This implies the network is constantly “echoing” infections through highly connected nodes, resulting in continual reinfection of a constant fraction of farms. We are especially interested in studying under what conditions the infectious animal disease persists and becomes localized under regular vaccination in the intra and interfarm model. However, besides the important role of the hub farms to stop epidemics in scalefree networks, it is also crucial to know how extensively the adverse effect is created by the lowvaccinated farms.
This work focuses on how the spatial transmission of infection is affected by both factors, the topology of networks and the internal resilience mechanism of nodes. We develop a networkbased model for transmission of infectious diseases of livestock, with a focus on internal structure of each node (farm). The model characterizes two distinct dynamical regimes: intrafarm dynamics where infection spreads fast among homogeneously mixing population following conventional compartment model and interfarm dynamics where the disease transmission occurs rather slowly along farmtofarm contact networks. Throughout this paper, we exclude culling and other control measures, since we focus on investigation of the effect of the vaccination on the transmission of FMD in the farm network. All simulation results are performed by MATLAB.
2. Robustness of Networks with Resilient Nodes
Resilience of networks implies how robust a network is to accidental or intentional attack on its vertices. In many models of epidemic spreading, the close relevance of network topology in the burst of the epidemics has been confirmed [17, 22].
The susceptibleinfectedsusceptible (SIS) model has been studied in a network representing potential transmission of the infection. It turned out that, for a scalefree (SF) network, the threshold for epidemic spreading is null in the infinite network limit. In other words, the epidemics in a sufficiently large SF network may never be eradicated even for low spreading rates. However, this twostate model is not directly applicable to farm networks, as the nodes on such networks have more complicated internal structures. One needs to refine the inter and intrafarm spreading process to deal with multitude of animals in the nodes at diverse infection stages.
In order to develop such hierarchical model of epidemic spreading, we first begin with a simple continuous SIS network model. The model will be elaborated with more realistic stages in the following sections. Suppose there are farms and each farm prevalence, say , , takes a continuous value between and where means absence of infection at the farm and means fully infected state; that is, all individuals in the farm are infected. Let denote an adjacency matrix whose th entry is 1 if the farms and are adjacent and the transmission occurs from the farms to , and 0 otherwise. The degree of th farm is defined as , which means the number of connected neighbors of the th farm. Then, the dynamics of the farm prevalence can be formulated as Here represents the internal epidemic development in the farm and represents the transmission rate function from the farm to the farm . If , the farm is said to be resilient from infection with the resiliency . This implies that the farm, if isolated, can recover from small infectious perturbation. In our approach, we assume that the farmtofarm transmission of the infection is proportional to both of the prevalence at the source farm and the susceptible level at the target farm. Hence the transmission rate function is set as where is the transmission rate. We want to see under what condition farm network (1) with (2) maintains the synchronized diseasefree state. The following theorem shows that a community of resilient farms is robust to infection as long as the network is not too dense and the transmission rate is low.
Theorem 1. Consider the farm network defined by (1) and (2) that consists of resilient farms. Let be the maximum degree of the farms, and let be the minimum resiliency. Then the collective diseasefree state is asymptotically stable for .
The proof of Theorem 1 is found at Appendix. The overall resiliency of a farm network depends on the resiliency of individual farms, the network density, and the transmission rate. Figure 1 depicts the mean farm prevalence at equilibrium of thousand farms in a scalefree network according to the transmission rate . Here the simple internal process function was used to model intensive practice of vaccination which leads to exponential decay of the infection ratio. One can observe that all the farms remain diseasefree where is less than a critical value around . In this parameter regime, even if some of the farms are simultaneously infected, the epidemic disease is soon eradicated from the network before it spreads further. However, such collective robustness does not hold and some of the farms fall in an endemic state if the transmission rate is greater than the critical value.
Although model (1) is a networkbased extension of an SIS model, it is still too simple to deal with multitude of animals of varying epidemic states in farms. In the following sections, we develop a hierarchical system that combines a conventional compartment epidemic model (intrafarm model) with a farm network model (interfarm model). The new model enables us to study how robust a farm network is with respect to the vaccination strategies and to what extent a farm network can endure sporadic exposure to external virus attack.
3. Intrafarm Model for FMD Transmission
Our intrafarm model for FMD transmission includes seven compartmental states; the susceptible, latently infected, infectious, and recovered states are denoted by , , , and , respectively, and the vaccinated susceptible, latent, and infectious states are denoted by , , and , respectively. The rate of vaccination followed by antibody formation is included in the parameter and the rate of transition from to is considered in the parameter , which is the rate at which the vaccine wanes off. The transition from to reflects the fact that the vaccinated individuals in can be infected with a reduced susceptibility factor . Moreover, the individuals in have a reduced infectivity, which is considered with the reduced infectivity factor in . An individual in the recovered state has antibody for averagely days after recovery from the disease and then becomes susceptible to the disease.
The schematic diagram of the model is illustrated in Figure 2.
In Figure 2, the parameters and denote the latent and infectious period, respectively, and the parameter denotes the transmission rate, which depends on the total population of livestock [23]. Moreover, the vaccinated individuals have latent and infectious period, and , respectively, and the reduced infectivity and susceptibility factors and , respectively. It is assumed that and due to the effect of the vaccination. The parameter description and values are given in Table 1.

The transmission dynamics of the model is described by the following ODE system: where .
The basic reproduction number denoted by , which reflects the stability property of the diseasefree equilibrium, can be computed by using the method of the next generation matrix [28].
Theorem 2. The basic reproduction number of (3) is where is the total population of livestock.
For the proof of Theorem 2, refer to Appendix. It is well known that if , then the model has a locally asymptotically stable diseasefree equilibrium and if , then the diseasefree equilibrium is unstable [29]. However, the condition is not necessarily a sufficient condition for disease elimination. Even when , an epidemic model with vaccination may have a backward bifurcation, for which a stable diseasefree equilibrium and a stable endemic equilibrium coexist [30, 31]. Here we investigate the bifurcation dynamics of system (3). Since we are mainly interested in the influence of the vaccination rate parameter on the disease transmission, we focus on the bifurcation dynamics of system (3) according to different values.
Figure 3 depicts a bifurcation diagram of system (3) in a certain range of values. For the simulation, we assume and set the values of the parameters except as follows: , , , , , , and and the vaccination effect has a twoday delay from to and a twoday early recovery from to , so that , This setting will be used throughout the rest of the paper. In Figure 3, one can see that if (or ), the farm has the stable diseasefree equilibrium. If (or ), the farm has the stable endemic equilibrium. In the intermediate regime (or ), the farm may end up either diseasefree or endemic, depending on the initial state. That is, in the range of , the model undergoes a backward bifurcation in the disease transmission.
4. Multiple Connected Farms
Our approach conceives a set of farms in a region as a network with connections between farms in the interfarm model. Each node of the network symbolizes a single farm, of which internal structure is described by the compartment model from the previous section. A link represents a relationship between farms that involves repeated contacts through which infections can spread.
There are two types of routes of betweenfarm transmission: direct infection is due to movement of infected animals and indirect infection occurs by sharing contaminated equipment, vehicles, workers, and veterinarians [32, 33]. We extend the assumption adopted in the intrafarm model and claim that the rate of farmtofarm transmission by both direct and indirect contacts is basically proportional to the number of the infected and the susceptible ones. More precisely, the rate at which the infections in farm cause new infections in farm is where and are the number of the susceptible ones in farm and the number of the infected ones in farm , respectively. Here denotes the farmtofarm transmission rate, which should be very small compared to the transmission rate in the intrafarm model. is estimated as the value between and per day [34].
Given a realization of some farm networks, one can simulate the dynamics of weakly coupled compartmental epidemic models on nodes. In the previous section, we observe that an isolated vaccinated farm has a limited resiliency and can return to the diseasefree state for small infections. However, such recovering ability of single farms is generally weakened when connected together, as they can play a role of a temporal repository of disease to one another. Figure 4 illustrates the change in the bifurcation diagram occurring on connection of multiple farms. It is notable that the upper branch representing stable equilibrium states is gradually shifting with the number of connected farms. This implies that a higher rate of vaccination is required to maintain a diseasefree state in a larger farm community. For example, four fully connected farms should have to remain healthy, while a single farm needs
Different from an isolated farm whose steady state may depend on its initial state, a collective state of multiple farms is less sensitive to the initial condition. Figure 5 shows that the backward bifurcation in the single farm model gradually evolves into a transcritical bifurcation as the number of directly connected farms increases. One can see that this change raises the minimum vaccination level to protect the farm community higher.
(a)
(b)
(c)
The dissolution of the backward bifurcation occurs with increase of directly connected neighbors and is also enhanced by networking. The steady states of the collective farms largely depend on the connection topology. In the remaining part of the work, we use scalefree (SF) networks which are considered to be the best suitable models for social networks, including the spread of epidemics [17, 35, 36]. Recent studies have shown that many transit systems found in different areas around the world exhibit scalefree structure [37–39]. This observation partially justifies the use of scalefree networks for the interfarm model, in that one of the major routes of infection is the livestockrelated vehicle movement.
SF networks follow the degree distribution , where is the degree of nodes, and they are characterized by the existence of the small number of nodes with many connections (hubs) and the large number of nodes with few connections. The SF network that we use in this work consists of farms. We set the average degree of the network to be Since every farm in the network is connected to 10 other farms in average, it is understandable that the average infection per farm is similar to that of 10 fully connected farms as in Figure 6. However, the backward bifurcation of the 10 farm models is “smeared out” in the network, raising the minimum vaccination level for diseasefree states from 0.36 to 0.43. This implies that a network, even when it has rather sparse topology, can provide more extensive potential repository of disease.
5. Sporadic Exposure to External Infections
We are interested in how robust overall immunity of vaccinated farms is and how they maintain the collective diseasefree state from the continual external infections. We assume the Poisson infection from the environment following the Poisson distribution with mean [40, 41]. This implies that each farm, besides horizontal infections from other farms in the network, is also exposed to the sporadic external infection at the averaged rate in a day.
In this section, the described model is used to run a set of simulations of epidemic for 10 years. We assume that the farmtofarm transmission rate is [34]. Each farm is assumed to be exposed to an external infection at the Poisson rate of average a day. Although this estimate for the parameter may not reflect real environmental settings, we rather focus on the relative ratio between the inter and intraparameters and study qualitative reactions of the farm networks to external infections.
In a recent outbreak of FMD in Korea, it has been reported that some farms have antibody formation rate less than 20% [42], which is far below the expected rate (about 97.5% for cattle and about 70% for pig) formed by Korea governmental standard vaccination policy. This low rate of antibody formation followed by the negligent vaccination is presumed to be one of the causes of the sudden outbreak and transmission of the disease between farms.
To investigate how the farms vaccinated at a low rate (“bad farms”) affect the disease transmission in the farm network, we assume two groups of farms; one is the group A of farms vaccinated at a high rate () and the other is the group B of farms vaccinated at a low rate (). Figure 7 shows that if the ratio of group B is small, less than 10%, then a massive outbreak of the disease does not occur in the farm network, but if it is big, more than 20%, the outbreak in the farm network occurs and the number of infected individuals increases in proportion to the ratio of the group B. This implies that even if most of the farms are wellvaccinated, that is, if most farms comply with the governmental standard vaccination policy, a small portion of bad farms seriously weakens the immunity of the entire farms. They can easily trigger the outbreak of the disease and moreover spread it long enough until other healthy farms in the network are infected. This may lead to a continual recurrent infection state. Thus, active surveillance and strong control about bad farms are required to prevent such permanent localization of the disease.
The phasetransitionlike jump from zero to a linearly increasing line in Figure 7 can be understood in terms of the hub infection: if the ratio of bad farms is set to be high, there is more chance for them to include one of the hubs. In spreading of disease over networks, it is well known that contamination of the hubs substantially damages the resilience of networks [43–45]. We further investigate how the control of hub farms is effective for preventing the longlasting spread of FMD in the farm network. A hub farm means a farm with heavy traffic frequency to other farms. We classify top 5% farms with highest degree as hub farms and observe their influence on the immunity of the whole networks according to their values of . Vaccination of the other 95% of nonhub farms is set to be relatively low at . In Figure 8, we see that a sufficient vaccination (more than ) for only top 5% hub farms is very effective for preventing the transmission between farms, although other nonhub farms have a low rate of vaccination. The influence of the hub vaccination is clear, especially compared to the random vaccination in which 5% of randomly chosen farms are vaccinated. This implies that a very efficient way to reduce the potential risk of FMD spread is to execute an urgent and intensive control on hub farms.
6. Conclusion
In this paper, we investigated how the spatial transmission of the FMD infection is affected by the two factors; the network topology and the internal resilience mechanism of nodes. We first developed a mathematical model for FMD transmission in farm networks and investigated the effect of vaccination on the control of outbreak of FMD.
While previous works on disease spreading in networks have mostly focused on the effect of the transmission rate, we take a closer look on the role of the internal mechanism of the nodes, especially their recovering ability induced by vaccination. We derived the vaccination level that is required to protect a farm community from sporadic exposure to external infectious source. Since most parameters used in the model are based on the previous studies and actual data, we believe that the result may suggest some practical implications in the FMD control in regional scale.
We found that even in the case that the number of farms with low vaccination rates or bad farms is small compared to the entire number of farms, the bad farms not only trigger an epidemic outbreak, but also hold FMD virus long enough to contaminate other healthy farms. This suggests that active and preemptive surveillance should be implemented for bad farms to prevent internalization of the disease. Moreover, we found that increasing the vaccination rate for entire farms by more than an appropriate level (0.35 or more) can prevent transmission of FMD between farms in a long time scale (10year). It is also observed that a moderate increase in vaccination rate for hub farms reduces very effectively the potential risk of FMD transmission to other farms. This finding suggests that when there is a sudden outbreak of FMD and an urgent control is needed to prevent massive outbreaks of the disease over farms, intensive and prompt vaccination for hub farms is an efficient and proper control measure.
Previous researches pointed out that some epidemic models with vaccination can have a backward bifurcation [30, 31, 46]. We observed that such a backward bifurcation gradually evolves into a transcritical bifurcation, raising the minimum vaccination level to protect the farm community higher. The dissolution of the backward bifurcation occurs with increase of directly connected neighbors and is also enhanced by the smearing effect of networking.
In scalefree networks, effectiveness of the vaccination on the hub nodes is known [43]. However, the hub vaccination for the epidemic diseases, especially for infectious human diseases, is difficult to implement since it is a challenging work to identify potential hub nodes for the disease transmission [47]. In case of FMD transmission, it is a possible work to identify hub nodes when the farm information such as farm location and livestock movement is available. In South Korea, actual geographical farm information and livestockrelated vehicle movement data are available from Korea Animal Health Integrated System (KAHIS). As an extension of this work, we plan to construct a networkbased model for describing FMD transmission between farms in Korea, using the farmrelated actual data obtained from KAHIS.
Appendix
A. Proof of Theorem 1
The Jacobian of system (1) at is Note and . Hence is a strictly diagonally dominant matrix for . Since ’s all diagonal elements are negative, from the Gershgorin circle theorem, the real parts of its eigenvalues are negative.
B. Proof of Theorem 2
We derive the basic reproduction number from system (3), using the method of the next generation matrix. We first find the matrices and the next generation matrix where and denote the values of and at the diseasefree equilibrium, respectively. Note that the basic reproduction number for system (3) is the spectral radius of . We compute the spectral radius of as Since and where is the total population of livestock, we obtain
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This research was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (2017R1D1A1B04032921 (P. Kim), 2016R1D1A1B03934427 (C. H. Lee)).
References
 P. Van Mieghem and R. Van De Bovenkamp, “NonMarkovian infection spread dramatically alters the susceptibleinfectedsusceptible epidemic threshold in networks,” Physical Review Letters, vol. 110, no. 10, Article ID 108701, 2013. View at: Publisher Site  Google Scholar
 C. Castellano, S. Fortunato, and V. Loreto, “Statistical physics of social dynamics,” Reviews of Modern Physics, vol. 81, no. 2, pp. 591–646, 2009. View at: Publisher Site  Google Scholar
 M. E. J. Newman, “Spread of epidemic disease on networks,” Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 66, no. 1, Article ID 016128, p. 016128/11, 2002. View at: Publisher Site  Google Scholar
 M. C. Gates and M. E. J. Woolhouse, “Controlling infectious disease through the targeted manipulation of contact network structure,” Epidemics, vol. 12, pp. 11–19, 2015. View at: Publisher Site  Google Scholar
 S. Rautureau, B. Dufour, and B. Durand, “Vulnerability of animal trade networks to the spread of infectious diseases: a methodological approach applied to evaluation and emergency control strategies in Cattle, France, 2005,” Transboundary and Emerging Diseases, vol. 58, no. 2, pp. 110–120, 2011. View at: Publisher Site  Google Scholar
 G. Rossi, R. L. Smith, S. Pongolini, and L. Bolzoni, “Modelling farmtofarm disease transmission through personnel movements: from visits to contacts, and back,” Scientific Reports, vol. 7, no. 1, article 2375, 2017. View at: Publisher Site  Google Scholar
 J. Gao, X. Liu, D. Li, and S. Havlin, “Recent progress on the resilience of complex networks,” Energies, vol. 8, no. 10, pp. 12187–12210, 2015. View at: Publisher Site  Google Scholar
 T. E. Carpenter, M. C. Thurmond, and T. W. Bates, “A simulation model of intraherd transmission of foot and mouth disease with reference to disease spread before and after clinical diagnosis,” Journal of Veterinary Diagnostic Investigation, vol. 16, no. 1, pp. 11–16, 2004. View at: Publisher Site  Google Scholar
 S. Mushayabasa, C. P. Bhunu, and M. Dhlamini, “Impact of culling stray dogs and vaccination on the control of human rabies: a mathematical modeling approach,” World Journal of Vaccines, vol. 1, pp. 156–161, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 M. J. Keeling, M. E. J. Woolhouse, D. J. Shaw et al., “Dynamics of the 2001 UK foot and mouth epidemic: stochastic dispersal in a heterogeneous landscape,” Science, vol. 294, no. 5543, pp. 813–817, 2001. View at: Publisher Site  Google Scholar
 M. J. Keeling, M. E. J. Woolhouse, R. M. May, G. Davies, and B. T. Grenfell, “Modelling vaccination strategies against footandmouth disease,” Nature, vol. 421, no. 6919, pp. 136–142, 2003. View at: Publisher Site  Google Scholar
 M. J. Tildesley, N. J. Savill, D. J. Shaw et al., “Optimal reactive vaccination strategies for a footandmouth outbreak in the UK,” Nature, vol. 440, no. 7080, pp. 83–86, 2006. View at: Publisher Site  Google Scholar
 F. Ball and D. Sirl, “Evaluation of vaccination strategies for SIR epidemics on random networks incorporating household structure,” Journal of Mathematical Biology, vol. 76, no. 12, pp. 483–530, 2018. View at: Publisher Site  Google Scholar  MathSciNet
 M. E. J. Woolhouse, D. J. Shaw, L. Matthews, W.C. Liu, D. J. Mellor, and M. R. Thomas, “Epidemiological implications of the contact network structure for cattle farms and the 2080 rule,” Biology Letters, vol. 1, no. 3, pp. 350–352, 2005. View at: Publisher Site  Google Scholar
 R. R. Kao, L. Danon, D. M. Green, and I. Z. Kiss, “Demographic structure and pathogen dynamics on the network of livestock movements in Great Britain,” Proceedings of the Royal Society of London B: Biological Sciences, vol. 273, no. 1597, pp. 1999–2007, 2006. View at: Publisher Site  Google Scholar
 K. L. VanderWaal, C. Picasso, E. A. Enns et al., “Network analysis of cattle movements in uruguay: quantifying heterogeneity for riskbased disease surveillance and control,” Preventive Veterinary Medicine, vol. 123, pp. 12–22, 2016. View at: Publisher Site  Google Scholar
 R. PastorSatorras and A. Vespignani, “Epidemic spreading in scalefree networks,” Physical Review Letters, vol. 86, no. 14, pp. 3200–3203, 2001. View at: Publisher Site  Google Scholar
 R. PastorSatorras and A. Vespignani, “Epidemic dynamics and endemic states in complex networks,” Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 63, no. 6, Article ID 066117, 2001. View at: Publisher Site  Google Scholar
 M. J. Keeling and K. T. D. Eames, “Networks and epidemic models,” Journal of the Royal Society Interface, vol. 2, no. 4, pp. 295–307, 2005. View at: Publisher Site  Google Scholar
 L. A. Meyers, “Contact network epidemiology: bond percolation applied to infectious disease prediction and control,” Bulletin of the American Mathematical Society, vol. 44, no. 1, pp. 63–86, 2007. View at: Publisher Site  Google Scholar
 T. KostovaVassilevska, “On the use of models to assess footandmouth disease transmission and control,” Tech. Rep. UCRLTR205241, Lawrence Livermore National Laboratory (LLNL), Livermore, CA, USA, 2004. View at: Publisher Site  Google Scholar
 A. L. Lloyd and R. M. May, “How viruses spread among computers and people,” Science, vol. 292, no. 5520, pp. 13161317, 2001. View at: Publisher Site  Google Scholar
 H. Lee, S. Lee, and C. H. Lee, “Stochastic methods for epidemic models: an application to the 2009 H1N1 influenza outbreak in Korea,” Applied Mathematics and Computation, vol. 286, pp. 232–249, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 A. C. Kinsley, G. Patterson, K. L. VanderWaal, M. E. Craft, and A. M. Perez, “Parameter values for epidemiological models of footandmouth disease in swine,” Frontiers in Veterinary Science, vol. 3, 2016. View at: Publisher Site  Google Scholar
 B. C. Verin, “Validation of commercially available NSP ELISA kits for foot and mouth disease surveillance to support the control and eradication programme,” in Proceedings of the IAEATECDOC–1546, International Atomic Energy Agency (IAEA '07), 2007. View at: Google Scholar
 S. P. Chen, T. M. Ellis, M. C. Lee et al., “Comparison of sensitivity and specificity in three commercial footandmouth disease virus nonstructural protein ELISA kits with swine sera in Taiwan,” Veterinary Microbiology, vol. 119, no. 2–4, pp. 164–172, 2007. View at: Publisher Site  Google Scholar
 C. Manore, B. McMahon, J. Fair, J. M. Hyman, M. Brown, and M. Labute, “Disease properties, geography, and mitigation strategies in a simulation spread of rinderpest across the United States,” Veterinary Research, vol. 42, no. 1, article 55, 2011. View at: Publisher Site  Google Scholar
 F. Brauer and C. CastilloChávez, Mathematical Models in Population Biology and Epidemiology, vol. 40 of Texts in Applied Mathematics, Springer, New York, NY, USA, 2001. View at: Publisher Site  MathSciNet
 P. van den Driessche and J. Watmough, “Reproduction numbers and subthreshold endemic equilibria for compartmental models of disease transmission,” Mathematical Biosciences, vol. 180, pp. 29–48, 2002. View at: Publisher Site  Google Scholar  MathSciNet
 F. Brauer, “Backward bifurcations in simple vaccination models,” Journal of Mathematical Analysis and Applications, vol. 298, no. 2, pp. 418–431, 2004. View at: Publisher Site  Google Scholar  MathSciNet
 A. B. Gumel, “Causes of backward bifurcations in some epidemiological models,” Journal of Mathematical Analysis and Applications, vol. 395, no. 1, pp. 355–365, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 G. Rossi, G. A. De Leo, S. Pongolini et al., “The potential role of direct and indirect contacts on infection spread in dairy farm networks,” PLoS Computational Biology, vol. 13, no. 1, Article ID e1005301, 2017. View at: Publisher Site  Google Scholar
 L. M. Mansley, A. I. Donaldson, M. V. Thrusfield, and N. Honhold, “Destructive tension: Mathematics versus experience—the progress and control of the 2001 foot and mouth disease epidemic in Great Britain,” Revue Scientifique et Technique de l'OIE, vol. 30, no. 2, pp. 483–498, 2011. View at: Publisher Site  Google Scholar
 H. Yoon, H. Kim, Y. Kim et al., “Study on application of a simulation program for outbreaks of footandmouth disease and classical swine fever,” in Research Report, Korea Animal and Plant Quarantine Agency, 2012, http://www.ndsl.kr/ndsl/search/detail/report/reportSearchResultDetail.do?cn=TRKO201400001119. View at: Google Scholar
 S. H. Strogatz, “Exploring complex networks,” Nature, vol. 410, no. 6825, pp. 268–276, 2001. View at: Publisher Site  Google Scholar
 G. Caldarelli, R. Marchetti, and L. Pietronero, “The fractal properties of Internet,” EPL (Europhysics Letters), vol. 52, no. 4, pp. 386–391, 2000. View at: Publisher Site  Google Scholar
 J. Wu, Z. Gao, H. Sun, and H. Huang, “Urban transit system as a scalefree network,” Modern Physics Letters B, vol. 18, no. 1920, pp. 1043–1049, 2004. View at: Publisher Site  Google Scholar
 J. Sienkiewicz and J. A. Hołyst, “Statistical analysis of 22 public transport networks in Poland,” Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 72, no. 4, Article ID 046127, 2005. View at: Publisher Site  Google Scholar
 A. Chatterjee, M. Manohar, and G. Ramadurai, “Statistical analysis of bus networks in India,” PLoS ONE, vol. 11, no. 12, Article ID e0168478, 2016. View at: Publisher Site  Google Scholar
 O. Diekmann and J. A. P. Heesterbeek, Mathematical Epidemiology of Infectious Diseases, Model Building, Analysis and Interpretation, John Wiley & Sons, 2000. View at: MathSciNet
 K. J. Sharkey, “Deterministic epidemiological models at the individual level,” Journal of Mathematical Biology, vol. 57, no. 3, pp. 311–331, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 Korea Animal and Plant Quarantine Agency, “Report on antibody formation rate test for FMD vaccine,” 2016, https://www.qia.go.kr/viewwebQiaCom.do?id=40921&type=2_8ktyzl. View at: Google Scholar
 R. PastorSatorras and A. Vespignani, “Immunization of complex networks,” Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 65, no. 3, Article ID 036104, 2002. View at: Publisher Site  Google Scholar
 G. Hartvigsen, J. M. Dresch, A. L. Zielinski, A. J. Macula, and C. C. Leary, “Network structure, and vaccination strategy and effort interact to affect the dynamics of influenza epidemics,” Journal of Theoretical Biology, vol. 246, no. 2, pp. 205–213, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 H. Zhang, J. Zhang, C. Zhou, M. Small, and B. Wang, “Hub nodes inhibit the outbreak of epidemic under voluntary vaccination,” New Journal of Physics , vol. 12, Article ID 023015, 2010. View at: Publisher Site  Google Scholar
 E. H. Elbasha, C. N. Podder, and A. B. Gumel, “Analyzing the dynamics of an SIRS vaccination model with waning natural and vaccineinduced immunity,” Nonlinear Analysis: Real World Applications, vol. 12, no. 5, pp. 2692–2705, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 F. Takeuchi and K. Yamamoto, “Effectiveness of realistic vaccination strategies for contact networks of various degree distributions,” Journal of Theoretical Biology, vol. 243, no. 1, pp. 39–47, 2006. View at: Publisher Site  Google Scholar  MathSciNet
Copyright
Copyright © 2018 Pilwon Kim and Chang Hyeong Lee. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.