Home | About | Journals | Submit | Contact Us | Français |

**|**BMC Public Health**|**v.11(Suppl 1); 2011**|**PMC3317581

Formats

Article sections

- Abstract
- Background
- Methods
- Results
- Conclusions
- Appendix: Effects of vaccination and population structure on influenza epidemic spread in the presence of two circulating strains
- Comparison with mean field model
- Initial conditions for the mean field model
- Authors’ contributions
- Competing interests
- References

Authors

Related links

BMC Public Health. 2011; 11(Suppl 1): S8.

Published online 2011 February 25. doi: 10.1186/1471-2458-11-S1-S8

PMCID: PMC3317581

Murray E Alexander: Murray.Alexander/at/nrc-cnrc.gc.ca

Mathematical Modelling of Influenza

Jean M Tchuenche and Robert Smith?

1471-2458-11-S1-info.pdfCopyright ©2011 Alexander and Kobes; licensee BioMed Central Ltd.

This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article has been cited by other articles in PMC.

Human influenza is characterized by seasonal epidemics, caused by rapid viral adaptation to population immunity. Vaccination against influenza must be updated annually, following surveillance of newly appearing viral strains. During an influenza season, several strains may be co-circulating, which will influence their individual evolution; furthermore, selective forces acting on the strains will be mediated by the transmission dynamics in the population. Clearly, viral evolution and public health policy are strongly interconnected. Understanding population-level dynamics of coexisting viral influenza infections, would be of great benefit in designing vaccination strategies.

We use a Markov network to extend a previous homogeneous model of two co-circulating influenza viral strains by including vaccination (either prior to or during an outbreak), age structure, and heterogeneity of the contact network. We explore the effects of changes in vaccination rate, cross-immunity, and delay in appearance of the second strain, on the size and timing of infection peaks, attack rates, and disease-induced mortality rate; and compare the outcomes of the network and corresponding homogeneous models.

Pre-vaccination is more effective than vaccination during an outbreak, resulting in lower attack rates for the first strain but higher attack rates for the second strain, until a “threshold” vaccination level of ~30-40% is reached, after which attack rates due to both strains sharply dropped. A small increase in mortality was found for increasing pre-vaccination coverage below about 40%, due to increasing numbers of strain 2 infections. The amount of cross-immunity present determines whether a second wave of infection will occur. Some significant differences were found between the homogeneous and network models, including timing and height of peak infection(s).

Contact and age structure significantly influence the propagation of disease in the population. The present model explores only qualitative behaviour, based on parameters derived for homogeneous influenza models, but may be used for realistic populations through statistical estimates of inter-age contact patterns. This could have significant implications for vaccination strategies in realistic models of populations in which more than one strain is circulating.

Human influenza infection is characterized by seasonal epidemics. This occurs because influenza A is able to maintain its presence in human populations by evolutionary adaptations to population-wide immunity, resulting in mutations that gradually change viral antigens allowing the virus to evade immune detection, a process known as “antigenic drift”. On account of these rapid mutations, vaccination for influenza must be updated annually on a global basis, following surveillance to monitor the appearance of new strains [1]. Antigenic drift also diminishes vaccine efficacy for mutant strains, but may still confer partial immunity to these strains. Therefore, understanding the short-term evolution of influenza virus is crucial to developing seasonal vaccines. Conversely, vaccination of a population may influence the short-term evolution of the virus, for example by decreasing the number of hosts in which the virus may replicate.

In general, during a single influenza season, more than one viral strain is circulating. It is known [2,3] that when suitable invasion conditions are satisfied, stable coexistence of two different strains is possible. The coexistence of two or more strains in a population will influence their individual evolution; and furthermore, the selective forces acting on the strains will be mediated by the transmission dynamics in the population. For example, the infection of hosts by one strain will reduce susceptibility to other strains, thereby limiting their spread in the population [4]. In addition, the time lag in emergence of a second strain following onset of an epidemic by a first strain will be influenced by the strategy and timing of vaccination [5]. It is clear that viral evolution and public health policy are strongly interconnected, and understanding the population-level dynamics of coexisting viral influenza infections, when vaccination of the population is to be undertaken, would be of great benefit in designing such vaccination strategies [6].

In [6], a homogeneous model of two viral strains was developed, incorporating cross-immunity and delay in emergence of the second strain. It was found that for small delay and large cross-immunity, infections with both strains appeared as a single epidemic wave; on the other hand, with sufficient delay, a second epidemic wave is possible. Further, for sufficient delay and high cross-immunity, the population of susceptible hosts may become so depleted as to prevent a second wave. These findings, together with possible impact of vaccination on antigenic drift, suggest that vaccination would be an important factor to include [6].

In large populations, contacts between individuals are not uniform, as assumed in the homogeneous model [6]. Typically, the number of contacts per day per individual is much smaller than the population size, and the structure of the corresponding ‘contact matrix’ plays an important role in the development of the pattern of the disease [7]. The effects of spatial correlations [8], such as occur when community structures are present [9], were illustrated in the spread of drug resistance in a network with mild clustering [10]: the spread of the resistant strain occurred more rapidly, and at significantly lower treatment levels, than was predicted by the homogeneous model.

The present paper extends the model in [6] in a number of ways. The model includes either pre-vaccination or vaccination during the epidemic, of a predetermined part of the population. The contact structure is modelled as a Markov network [11], in which the distribution of degrees of the nodes (i.e., number of contacts for individuals in the population) is specified. In addition, the model allows a distribution of ages in the population by incorporating a prescribed number of age classes. The Markov assumption for the contact network allows the specification of structural parameters such as assortativity [12] and clustering [13-15] that are important characteristics of social groupings. These generalizations enable vaccination to be targeted according to age group and ‘contact number’ (degree of node), which in general respond to the vaccine in different ways. The model inevitably contains many parameters and allows a wide range of network structures to be specified; in addition, initial conditions can be specified in many different ways. Therefore, in this paper, only a simplified network model will be investigated. The structure of the network is comprised of uncorrelated nodes, with degree distribution specified as a truncated scale-free form [7]. Furthermore, for simplicity only one or two age-classes are considered, where, for the latter, the median age is chosen to separate the two classes. While a detailed age distribution, characteristic of a real population, could be specified, the present results are intended to be illustrative only and to allow comparison with the corresponding homogeneous model. The network model can potentially be useful in describing specific populations, such as a small or large city, in which case the network structure and age distribution would need to be determined from statistical analysis of demographic and census data [16].

Section 2 describes the model in broad terms, and lists some of the parameter values used; technical details are given in the Appendix. Section 3 presents the results of simulations, in which the cross-immunity and delay in appearance of the second strain infection are varied. These results are also compared with those produced by the corresponding homogeneous model, to ascertain the importance of structure in the network for determining the time-course and final extent (“total attack rate”) of the disease. Finally, Section 4 discusses these results, some possible extensions of the model, and implications for vaccination strategies in more realistic models based on specific demographic data.

The state flow diagram of the model is given in Figure Figure1,1, which represents either the population counts in various compartments in the homogeneous model or the state of any given node (labelled by infection state, degree- and age-class) in the network model. The model describes the evolution of two concurrent strains of influenza infection, over a duration short compared to the natural lifespan of an individual in the population, and for this reason birth and natural death processes are ignored, and furthermore the number of individuals in each age class remains constant. Since we consider a static contact network, the degree class of each individual is also fixed. Therefore, each individual in the population belongs to a unique class (*k*,*a*), where *k* denotes the number of contacts, and *a* the age class, and the total number of individuals in each (*k,a*) class is constant. *S* denotes the susceptibles and *V* denotes individuals receiving vaccination either prior to the onset of the first infection or after this onset.

Vaccination prior to the onset of infection is specified by the fraction of susceptibles in each age class receiving vaccination. For vaccination occurring during an outbreak, the following model is used: for individuals in any given (*k,a*) class, the rate of vaccination at any given time is (i) proportional to the current number of susceptibles in the class; (ii) an increasing function of the total current (symptomatic) infection in the population as a whole, saturating at a prescribed rate. This was done to attempt to model the social response to an outbreak in the population, in which the greater the number of infected individuals the more likely that susceptible individuals would avail themselves of existing vaccination opportunities. The precise mathematical specification of this response is given in the Appendix.

The baseline transmission rate of infection between a susceptible-infected pair of individuals is denoted by τ. The actual rate will depend on the age-classes that these individuals belong to, and whether the susceptible individual of the pair is seeing infection (by either strain) for the first or second time. These various possibilities are accounted for by expressing the actual transmission rate as τ times a factor, which depends on age classes involved, whether this is the first or second infection, and whether the individual has received prior vaccination. Details are given in the Appendix and in Table Table11.

States labelled with *I* denote symptomatic infection, and those labelled with *A* denote asymptomatic infection. The *P* states describe immunity to one strain but not the other: *P _{j}* is the state with immunity to strain

Subscript ‘*V*’ denotes states of infection (or partial recovery) arising from failure of the vaccine; and as before, labels states with infection due to, or partial recovery from, one of the strains. Following vaccination, infection due to strain *j* occurs with probability (1-σ* _{j}*). In general, for seasonal influenza, the vaccine is targeted against the earlier-occurring strain 1 virus; its efficacy against the later-occurring strain 2 (mutated) virus is expected to be less, i.e., σ

In Figure Figure1,1, the diverging pairs of directed edges are labelled with branching ratios for each strain of infection, with two pairs of such edges emanating from *S* and *V* classes. For example, if *S* is infected with one of the strains, it has a probability *p* of being symptomatically infected, and 1-*p* of being asymptomatically infected. (We assume that *p* is the same for both strains). Since *S* may be infected with either strain, there are two pairs of branches emanating from *S* in Figure Figure1.1. Similarly, there are two branch pairs for *V*, representing infection due to failure of the vaccine.

After recovery from one strain of infection, an individual is still, in general, susceptible to infection by the other strain: individuals in state *P _{j}* (i.e., recovered from infection with strain

The converging directed edges in this Figure are labelled with the recovery rates from infection: either µ (symptomatic infection) or *µ _{A}* (asymptomatic infection), where we assume that these rates are the same for both strains, regardless of whether this is the first or second infection for that individual. The parameter values used in the simulations are given in Table Table11.

For the homogeneous model, we may apply the technique of the next-generation matrix [17] to derive the basic reproductive number *R*_{0}. In general, the second strain appears after infection due to the first strain has begun, so that *R*_{0} can be calculated using a one-strain sub-model. With this assumption, we find

where *S*_{0}, *V*_{0}
denote the initial numbers of susceptible and vaccinated individuals, respectively, in the population, and β is the transmission coefficient. To establish a relationship between β and the baseline transmission rate τ between individuals in contact, we construct a (single-age class) network in which the ‘edge probability’ of randomly choosing an edge, one of whose vertices has degree *k*, is uniform. By relating this to the mean field model we derive (see Appendix)

where *k*_{1}
= vertex degree of population sub-class into which the Strain-1 infection is introduced at time *t* = 0, and *k _{max}* = maximum vertex degree in the finite network (

In keeping with the definition of the two age class model (see Appendix), the estimates of death rates [18,19] arising from symptomatic or asymptomatic infection (*d, d _{A}*, respectively) for the two age-class model correspond to the general population above and below the median of the age distribution

In the model described above, the total number of individuals *N _{k}*

The initial state was specified as follows. For pre-vaccination, a prescribed fraction *V*_{0}
(*a*) of individuals in each age class *a* receive vaccination. Infection by strain 1 is introduced into fraction *ε*_{1} of the remaining susceptibles residing in a single class (*k*_{1},*a*_{1}). After the strain 1 infection has spread through the population for a time *T**, a strain 2 infection is introduced into a fraction ε_{2} of class (*k*_{2},*a*_{2}) individuals. In the simulations, we use ε_{1} = ε_{2} = 0.5; *k*_{1} = 5, *k*_{2} = 10, and for the two age class model, *a*_{1} = 1, *a*_{2} = 2. As previously mentioned, it is assumed that no individual may be infected with both strains simultaneously. The simulations were performed using three models: (1) network model with two age classes; (2) network model with one age class; and (3) the homogeneously-mixing 'mean field' model. For (1) and (2), the structure of the network was chosen to have a scale-free form [7], with the number of individuals (nodes of the contact network) with *k* contacts being proportional to *k*^{-2.5}[21], and 1 *≤ k ≤ k _{max}* = 20. Furthermore, the degrees of the nodes of the network were assumed to be uncorrelated: although real networks show significant correlation structure – e.g., clustering and associativity [12-15] - the purpose of the present simulations is to illustrate the general effects of departure from the homogeneous mixing assumption. A value

Figure Figure22 compares the results from the network model using one age class (top row), two age classes (middle row), and the mean field model (bottom row), for delays of *T** = 10 days (left column) and 60 days (right column) for introducing strain-2 infection into a prescribed sub-population. Solid curves correspond to a large (60%) cross-immunity (with *δ* *δ*_{12} =*δ*_{21} = 0.4) between strains, and dashed curves to a low (10%) cross-immunity (*δ* = 0.9). These results show that the level of cross-immunity has a significant effect on whether a second wave appears: if it is too high (60%, δ = 0.4), then the second wave does not appear in any of the three models. For the 2-age class model (middle row), the (first) peak infection occurs at a slightly lower value for the larger cross-immunity: 0.03 vs. 0.033 fraction of the population; however, the timing of these peaks (~75 days after initial infection) is not sensitive to the level of cross-immunity. Similar conclusions apply to the one age class network model and the mean field model, though the timing of these peaks is different for different models.

As expected, if the second strain is introduced after the strain 1 infection has been largely cleared from the population (as is the case when *T** = 60 days: right column), then the first and second waves behave as distinct, non-interacting one-strain epidemics. However, when *T** is only 10 days (left column), there is still a significant presence of strain 1 infection in the population: the infections in the two age-class model merge into a single broad peak, whereas the other two models show two distinct peaks, with the second peak occurring in both models ~50 days after initial infection.

It is therefore apparent that the two age class network model exhibits a larger delay in peak infection – for both first and second waves – compared to the one age class and mean-field models. This can be accounted for by the reduced transmissibility between classes compared to within one class, as well as reduced transmissibility within the second age class (see the *M* matrix in the Appendix). (Recall that strain 1 and 2 infections are introduced into different age classes in the two age class network model). Such differences in delays between mean-field and structured models have been observed elsewhere [10], and underline the importance of spatial structure in determining the course of an epidemic.

Figure Figure33 shows the effect of different levels of pre-vaccination (*V*_{0} = 0.2 [dashed curves] or 0.4 [solid curves]) on the infection profiles, for each of the three models, assuming a delay *T** = 10 days in appearance of strain 2 infection. For all models, the peak infection is reduced by about 50%, and occurs slightly later, when increasing from 20% to 40% coverage. The peak infections in the two age class model, however, occur at significantly lower values than in the other two models: by a factor of 5 for the mean field model, and a factor of 10 for the one age class network model. Again, this may be accounted for by the reduced transmissibility between different age classes. Notice also that, as in Figure Figure2,2, the two age class model exhibits a single peak of infection for both levels of vaccination.

In Figure Figure4,4, the two age class network model is used to explore the effects of vaccination during an epidemic outbreak, with no vaccination prior to the initial appearance of strain 1 infection, where vaccination rates are determined according to the “social response” to total infection in the population, as described in Section 2 and the Appendix. The resulting infection profiles are remarkably insensitive to (i) level of cross-immunity, (ii) delay *T**, and (iii) the rate of vaccination ω_{0}. Apart from a temporary decrease in the rate of disappearance of the infection when *T** = 10 days and δ = 0.4 (top left in Figure Figure4),4), there is only one peak of infection which occurs consistently at very similar times (50 days after strain 1 infection) and with peak magnitudes between 0.07 and 0.08 of the fraction of population. Comparing these results to the two age class model with pre-vaccination (middle row of Figure Figure2),2), it can be seen that the peak infection occurs about 20-30 days earlier for vaccination during the epidemic than when pre-vaccination is carried out.

Tables Tables2,2, ,3,3, ,44 compare the attack rates when individuals are infected by each strain separately (represented by the final values of the *P*_{1}
and *P*_{2}
compartments in Figure Figure1)1) and by both strains in succession (represented by the *R* compartment), for each of the three models. It is seen that in general, for given *T** and δ, the attack rates predicted by the network models are quite different between the model types. However for pre-vaccination, within each model, strain 1 attack rates increase, while strain 2 attack rates decrease, with the delay *T**. For vaccination during the epidemic, these trends are also seen though for strain 1 infection are less pronounced. This is reasonable, because for longer *T**, the strain 2 infection draws upon a smaller pool of susceptibles (and individuals with failed vaccination) which are depleted by strain 1 infection for a longer time interval before strain 2 appears. Also, as expected, strain 2 attack rates are sensitive to the level of cross-immunity of the strains, decreasing sharply as cross-immunity increases from 10% to 60%. Variations between models are manifestations of the importance of heterogeneity of contact structure. The dependence on age distribution between network models in Tables Tables2,2, ,33 is a consequence of our assumption (see matrix *M* in the Appendix) that transmissibility within age class 1 is greater than that within age class 2 or between age classes.

A question of importance to public health is the dependence of both the attack rate and death rate due to infection. The attack rates are shown in Table Table55 for the two age class network model with *T** = 60 days, and cross-immunity of 10% and 60%. For pre-vaccination, as expected, strain 1 attack rates decrease with increasing level of vaccination, at first slowly then dropping off sharply between *V*_{0} = 0.3 and 0.4 regardless of the level of cross-immunity. Interestingly, the maximum attack rate due to strain 2 infections is reached in this range, and drops off sharply thereafter. Thus, for the parameter values used, it appears that vaccination is most effective around these values, with diminishing returns for higher *V*_{0}. Tables Tables6,6, ,77 show that death rates due to pre-vaccination are lower than predicted for the entire range of vaccination rates during an epidemic; and in both cases the death rates decrease with increasing levels of vaccination. Analogous to the attack rates, there is a small increase in mortality for vaccination coverage around 20%, due to increased numbers of strain 2 infections at the expense of reduced numbers of strain 1 infections; however, the mortality rate drops sharply once the vaccination coverage exceeds about 40%.

We have considered extensions of the two viral strain mean field (homogenous) model introduced in [6], to explore the effects of both local network structure and the division of the population into different age classes. The present study used model parameter values (in particular, *R*_{0}) originally estimated for mean field models; and in order to translate these to the network models derived in this paper, a correspondence was established between the mean field model and a limiting case of the network model (see Appendix). The two age class model assumed the age boundary was located at the median age (about 38 years for Vancouver), with vaccine efficacy of 80% in the lower age group and 40% in the upper age group.

Several notable features were observed when comparing the network models to the corresponding mean field case. Firstly, the amount of cross-immunity present is significant in determining whether a second wave of infection occurs. Due to a lower transmission rate between age classes and within the second age class, compared to within the first age class, infection levels were found to be significantly less for the two age class model than for either the one age class or mean field models. The infections occurred as either a single wave or as two successive waves. A second wave is more likely to occur the longer the delay in introduction of the second strain, since when this delay is short (~10 days) infections due to both strains merge into a single, broad peak. When a second wave does occur, the shapes of the two waves depend on when the second strain infection is introduced. If it occurs well after the first infection has run its course, then the two waves behave as distinct, non-interacting infections. The second infection peak is delayed, and its amplitude reduced, in the network model, compared to the mean field case. This behaviour reflects a longer propagation time in the network model, and has been qualitatively observed in other models, reinforcing the importance of including local network structure in realistic models.

As expected, the amount of cross-immunity between the two strains is important in determining the size of the second-strain outbreak. It was found that its size decreased sharply with increasing cross-immunity. As the level of vaccination increases, strain 1 attack rates decrease, with a sharp drop occurring around 30-40% pre-vaccination coverage; at the same time, strain 2 infections *increase* with increasing vaccination coverage, reaching their maximum somewhere in this range, and drop off sharply for higher coverage levels. This phenomenon is reminiscent of the development of drug resistance, where there is an optimal level of drug treatment (compare: vaccination coverage) that minimizes the overall infection [10]. This could have significant implications for vaccination strategies in realistic models of populations in which more than one strain is circulating.

It was found that increasing either pre-vaccination or vaccination during an outbreak, reduces the disease-induced mortality. Furthermore, pre-vaccination appears to be more effective than vaccination during an outbreak in reducing overall mortality, though this needs further investigation as it may depend critically on how the latter is implemented. This study considered only a simple model in which at any given time vaccination rates during an outbreak were governed by the total infection in the population at that time, and considers only vaccination of the susceptible class *S*, neglecting vaccination of other classes (e.g., *P*_{1} and *P*_{2} and asymptomatic cases).

As mentioned earlier, the particular form of the terms included in the model to incorporate local network structure and the effects of age classes was chosen for illustrative purposes. This approach, though, can be used on a specific population if sufficient data are available to determine realistic estimates of the age classes and network structure present and of the parameters of the model. The main difficulty is in determining the form of the two-point correlations between vertices of the contact network for a realistic particular population, and this must be derived indirectly from estimates of network structure extracted from the data [16]. An intermediate approach is to explore the effects of a few network structure parameters – e.g., clustering, associativity, betweenness, and centrality [7,16], obtaining expressions for the two-point probabilities defining the Markov network directly in terms of these parameters. This is currently under investigation.

The various parameters in the model (Figure (Figure11 of main text) are defined below:

τ = baseline transmission rate between a susceptible-infected pair

*p* = probability of developing symptomatic infection with no prior exposure

*p _{V}*

*σ*_{1}, *σ*_{2} = effectiveness of vaccine to strains 1 and 2, respectively

*δ _{V}*

*p*_{12} = probability of developing symptomatic infection with prior exposure to strain 1

*p _{V}*

*p*_{21} = probability of developing symptomatic infection with prior exposure to strain 2

*p _{V}*

*δ _{A}* = reduction in infectiousness due to asymptomatic infection

*µ, µ _{A}* = recovery rates from symptomatic and asymptomatic infections

*δ*_{12}, *δ*_{21} = level of cross-immunity induced by previous exposure to strain 1 and strain 2, respectively

*d, d _{A}* = disease-induced death rates, assumed to be age-dependent but the same for each type of infection.

Define age classes 1 ≤ *a* ≤ *a _{max}*, and network degree classes 1 ≤

for *U* {*A*_{1}, *A _{V}*

In the special case that the contact network is the same for all age-classes, *P*(*k′,a′|k,a*) = *P _{a}*(

In the general case, . If the node-degrees are uncorrelated, then *P*(*k′|k*) *= P _{e}*(

We extend the homogeneous (mean-field) model in [6] to a Markov network model with age structure, and include vaccination of strain 1 prior to onset of influenza outbreak. Furthermore, we extend this model to allow vaccination during an epidemic, by introducing a flux of susceptibles from the *S* classes to their corresponding *V* classes, according to the prescription:

(i) is a function of the total (symptomatic) infection in the population, *I _{tot}* (summed over all

(ii) is proportional to the population in class *S*(*k,a,t*);

(iii) eventually saturates at a maximum value as *I _{tot}* increases.

A functional form that satisfies these criteria is

where ω_{0} is the (age-class dependent) saturation rate of vaccination, α is the value of *I _{tot}* at half-saturation, and

In order to incorporate death due to infection, we add a set of classes {*D*(*k,a,t*)} to the model, and (similar to the recovery rates) assume that death rates are either *d* (for all symptomatic infections) or *d _{A}* (for all asymptomatic infections).

Based on Figure Figure11 (main text), this gives rise to the following set of 24×(*a _{max}*×

where *S* *S*(*k,a,t*), etc., and the force of infection *Θ _{U}* may be expressed as

where *U* {*A*_{1}, *A _{V}*

where

This shows that the various *Θ*(*k,a,t*)*’s* describe the connectivity of a vertex of degree *k* and age-class *a* to all the infected adjacent vertices.

In order to make this comparison, we need to obtain a “limiting” form of the network model that approximates the mean field model. This will enable us to obtain a relationship between the mean field model transmission rate β and the transmission rate τ between a susceptible-infected pair in the network. We consider a simplified network model consisting of a single age class, and an uncorrelated network (so that *P*(*k’|k*)*=P _{e}*(

Therefore, the equation for *S*(*k,t*) becomes

where the term in parentheses is independent of *k*. Assume that *S*(*k,t*) *= P*(*k*)*S*(*t*), etc.; then this becomes

Summing over *k* from 1 to *k _{max}*, we obtain

whereis the mean degree of nodes in the network.

Comparing this approximation with the Mean Field expression for *dS/dt*, suggests we make the following correspondence:

More realistically, since we introduce Strain-1 infections into the sub-population defined by (*k,a*) *=* (*k*_{1},*a*_{1}), and taking account of the fact that *R*_{0}
is defined in terms of the *first* generation of infection, it would be more accurate to replace by *k*_{1}, so that

Using this approximate relationship enables us to compare the simulation of the behaviours of the network and mean-field models, by relating numerical values of the parameters β and τ through the simplified limiting case of a network in which the probability of drawing an *edge* at random from the network is uniform.

In what follows, it is assumed that the total population (including deaths) is normalized to unity, which is permissible since for this model it is constant. The initial conditions for the mean field model must be consistent with those of the network model. The analysis that follows applies to an arbitrary number of age classes and degree distributions.

In the network model, at *t* = 0 a fraction ε_{1} of Strain-1 infection is introduced into the sub-population in class (*k*_{1},*a*_{1}). Therefore, we specify the initial conditions at *t* = 0 according to

where *P _{ka}*

whereare the changes in the susceptible and vaccinated sub-populations, respectively, and *P*_{ka2} *P*(*k*_{2}
)*P _{a}*(

In order to allow comparison between Mean Field and network models, all age-dependent parameters *δ _{A}*, δ

For the network model, for all age classes we set ε_{1} = ε_{2} = 0.5, *p* = 0.6, *V*_{0}
= 0.2, σ_{1} = 0.8, σ_{2} = 0.4. For the two age-class model, we chose (*k*_{1},*a*_{1}) = (5,1), (*k*_{2},*a*_{2}
) = (10,2); and for the one-age class model *k*_{1} = 5, *k*_{2}
= 10. The (truncated) scale-free distribution *P*(*k*) ~ *k ^{-}*

Substituting these values into the expression for *R*_{0} (Section 2 in main paper), and using *R*_{0} = 1.9, *V*_{0} = 0.2, *k _{max}* = 20, and

MEA wrote the bulk of the manuscript, and designed, implemented and simulated the network model from a homogeneous one constructed with assistance from Dr S. Moghadas. RK implemented and simulated the homogeneous model, and wrote the Conclusions. MEA wrote the Appendix.

The authors declare that they have no competing interests

It is with great sadness that I report the sudden passing of my dear friend and co-author, Dr Randy Kobes, whilst this paper was in revision. He will be sorely missed by his many friends and colleagues. I would like to thank Dr Seyed Moghadas for his helpful comments and assistance in designing the mean field model, and the referee for incisive and useful comments that have greatly improved the paper. This work was supported by a Natural Sciences and Engineering Research Council of Canada Discovery Grant.

This article has been published as part of *BMC Public Health* Volume 11 Supplement 1, 2011: Mathematical Modelling of Influenza. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2458/11?issue=S1.

- Boni MF. Vaccination and antigenic drift in influenza. Vaccine. 2008;26(Suppl 3):C8–C14. doi: 10.1016/j.vaccine.2008.04.011. [PMC free article] [PubMed] [Cross Ref]
- Castillo-Chavez C, Hethcote HW, Andreasen V, Levin SA, Liu WM. Epidemiological models with age structure, proportionate mixing, and cross-immunity. J. 1989;27:233–258. doi: 10.1007/BF00290636. [PubMed] [Cross Ref]
- Dietz K. Epidemiologic interference of virus populations. J. 1979;8:291–300. doi: 10.1007/BF00276314. [PubMed] [Cross Ref]
- Andreassen V, Lin J, Levin SA. The dynamics of cocirculating influenza strains conferring partial cross-immunity. J. 1997;35:825–842. doi: 10.1007/s002850050079. [PubMed] [Cross Ref]
- Bowman CS, Arino J, Moghadas SM. Evaluation of vaccination strategies during pandemic outbreaks. Math. 2011;8:113–122. doi: 10.3934/mbe.2011.8.113. [PubMed] [Cross Ref]
- Moghadas SM, Bowman CS, Arino J. Competitive interference between influenza viral strains. Canadian Appl. Math. Quarterly. in press .
- Newman MEJ. The structure and function of complex networks. SIAM Rev. 2003;45:167–256. doi: 10.1137/S003614450342480. [Cross Ref]
- Keeling M. The effects of local spatial structure on epidemiological invasions. Proc Biol Sci. 1999;266:859–867. doi: 10.1098/rspb.1999.0716. [PMC free article] [PubMed] [Cross Ref]
- Fortunato S. Community detection in graphs. Physics Reports. 2010;486:75–174. doi: 10.1016/j.physrep.2009.11.002. [Cross Ref]
- Alexander ME, Dietrich SM, Hua Y, Moghadas SM. A comparative evaluation of modelling strategies for the effect of treatment and host interactions on the spread of drug resistance. J. 2009;259:253–263. doi: 10.1016/j.jtbi.2009.03.029. [PubMed] [Cross Ref]
- Boguñá M, Pastor-Satorras R, Vespignani A. In: Lecture Notes in Physics. Pastor-Satorras R, Rubi M, Diaz-Guilera A, editor. Vol. 625. New York: Springer; 2003. Epidemic spreading in complex networks with degree correlations; pp. 127–147.
- Newman MEJ. Assortative mixing in networks. Phys. Rev. Letts. 2002;89:208701. doi: 10.1103/PhysRevLett.89.208701. [PubMed] [Cross Ref]
- Park J, Newman MEJ. Solution for the properties of a clustered network. Phys. Rev. E. 2005;72:026136. doi: 10.1103/PhysRevE.72.026136. [PubMed] [Cross Ref]
- Serrano MÁ, Boguñá M. Clustering in complex networks. I. General formalism. Phys. Rev. E. 2006;74:056114. doi: 10.1103/PhysRevE.74.056114. [PubMed] [Cross Ref]
- Serrano MA, Boguñá M. Clustering in complex networks. II. Percolation properties. Phys. Rev. E. 2006;74:056115. doi: 10.1103/PhysRevE.74.056115. [PubMed] [Cross Ref]
- Kolaczyk ED. Statistical Analysis of Network Data. New York: Springer; 2009. [PubMed]
- van den Driessche P, Watmough J. Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math. Biosci. 2002;180:29–48. doi: 10.1016/S0025-5564(02)00108-6. [PubMed] [Cross Ref]
- Alexander ME, Bowman C, Moghadas SM, Summers R, Gumel AB, Sahai BM. A vaccination model for transmission dynamics for influenza. SIAM J. 2004;3:503–524. doi: 10.1137/030600370. [Cross Ref]
- Dushoff J. Mortality due to Influenza in the United States–An Annualized Regression Approach Using Multiple-Cause Mortality Data. American Journal of Epidemiology. 2005;163:181–187. doi: 10.1093/aje/kwj024. [PubMed] [Cross Ref]
- Community Services Group. Vancouver; 1996. Data from the Canada Census, Population by Age. City of Vancouver. Vancouver Local Areas 1996.
- Pastor-Satorras R, Vespignani A. Epidemic spreading in scale-free networks. Phys. 2001;86:3200–3203. doi: 10.1103/PhysRevLett.86.3200. [PubMed] [Cross Ref]
- Weber S, Porto M. Generation of arbitrarily two-point-correlated random networks. Phys. Rev. E. 2003;76:046111. doi: 10.1103/PhysRevE.76.046111. [PubMed] [Cross Ref]
- Arino J, Bowman C, Moghadas SM. Antiviral resistance during pandemic influenza: implications for stockpiling and drug use. BMC Infect. 2009. DOI:10.1186/1471-2334-9-8. [PMC free article] [PubMed]
- Cauchemez S, Carrat F, Viboud C, Valleron AJ, Boelle PY. A Bayesian MCMC approach to study transmission of influenza : application to household longitudinal data. Stat. 2004;23:3469–3487. doi: 10.1002/sim.1912. [PubMed] [Cross Ref]
- Gojovic MZ, Sander B, Fisman D, Krahn MD, Bauch CT. Modelling mitigation strategies for pandemic (H1N1) CMAJ. 2009. DOI:10.1503/cmaj.091641. [PMC free article] [PubMed]

Articles from BMC Public Health are provided here courtesy of **BioMed Central**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's Canada Institute for Scientific and Technical Information in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |