 Methodology
 Open Access
 Published:
Determining the dynamics of influenza transmission by age
Emerging Themes in Epidemiology volume 11, Article number: 4 (2014)
Abstract
Background
It is widely accepted that influenza transmission dynamics vary by age; however methods to quantify the reproductive number by age group are limited. We introduce a simple method to estimate the reproductive number by modifying the method originally proposed by Wallinga and Teunis and using existing information on contact patterns between age groups. We additionally perform a sensitivity analysis to determine the potential impact of differential healthcare seeking patterns by age. We illustrate this method using data from the 2009 H1N1 Influenza pandemic in Gauteng Province, South Africa.
Results
Our results are consistent with others in showing decreased transmission with age. We show that results can change markedly when we make the account for differential healthcare seeking behaviors by age.
Conclusions
We show substantial heterogeneity in transmission by age group during the Influenza A H1N1 pandemic in South Africa. This information can greatly assist in targeting interventions and implementing social distancing measures.
Background
The importance of the dynamics of influenza transmission between age groups is wellappreciated [1–6]. Several studies have assessed the nonuniformity of the impact of influenza, particularly pandemic influenza, on different age groups [1–4, 7–9]. The overarching interest in these studies is to gather information in order to influence policy to best determine a strategy to impact on the spread of outbreaks. Which age groups carry the greatest disease burden and which groups are responsible for the greatest amount of disease transmission is an important component of this information.
One key aspect of this work is to estimate the extent to which people in different age groups interact with one another and to what degree they are in contact. This information can then be used as a surrogate for transmission probabilities between age groups [10–12]. Several studies have generated matrices with estimated numbers of contacts between various age groups [5, 12, 13]. Additionally, social network models have been used to estimate these contact patterns [7]. A common finding amongst these studies is that children tend to mix mostly with each other, and to a lesser extent with their parents, while adults mix with individuals from a larger range of ages [5, 7]. These matrices have subsequently been used in modeling exercises to better understand the dynamics of disease spread by age.
When determining which groups to target for interventions in an outbreak, one strategy is to target those who potentially carry the greatest burden of disease [14], which has typically been found to be children [1, 5, 7, 8, 13]. For instance, Bansel et al.[4] consider data from influenza pandemics over the past century and show that the burden of disease is highest amongst children during a pandemic and then shifts to adults the following season.
To better understand the utility of targeting the groups with the greatest burden of disease, it is also important to determine when different age groups tend to have their peak incidence of cases. This can also be seen as a surrogate measure for the age group that is driving an outbreak [15–17]. Most recently Schanzer et al. [16] used 10 years of Canadian surveillance data of laboratory confirmed cases of influenza and found that during seasonal influenza the 10–19 and 20–29 yearold age groups tended to peak one week earlier than other groups. During the pandemic in 2009, the peak came earliest for only the 10–19 year old age group. This is inconsistent with the findings of Brownstein et al. [15] who found that children aged 3–4 were consistently the first to peak.
A different tactic for determining which groups should be the target of interventions is to select those groups most responsible for transmission [18]. Studies examining disease transmission by age have consistently shown that children have higher estimated values of the reproductive number than adults. Recent work has focused on the dynamics during the 2009 Influenza pandemic. During the initial phase of the pandemic in Japan, Nishiura et al. [19] report that children were transmitting illness at higher rates than adults. Glass et al. [2] used Japanese data and a novel method to estimate the reproductive number, R_{0}, for adults and children that assume particular forms for a nextgeneration matrix and estimate the parameters of this matrix, leading to outbreak specific estimates of R_{0}. They estimate R_{0} to be between 2.8 and 3.6 for children and between 0.2 and 0.7 for adults, depending on the assumptions made. In a later study, Glass et al. [3] used serosurvey data and estimate R_{0} from the final size of the epidemic to be 1.6 for kids and less than 1 for adults. Wallinga et al. [18] have similarly shown that the rate in change of the reproductive number for a particular group is related to the incidence of infection and force of infection and suggest allocating resources based on examining these two quantities.
In the present study, we present a new approach to estimating the reproductive number by age group by modifying a method initially proposed by Wallinga and Tuenis [20]. We study age dynamics in South Africa during the 2009 Influenza H1N1 pandemic and illustrate the importance of an appropriately estimated measure of the transmission dynamics on final estimates. Finally, we discuss our results and their implications for future studies on how to respond during an emerging outbreak.
Methods
Data
We use deidentified data previously reported in [21] that includes a line list of the 12,543 confirmed cases reported in South Africa during that outbreak. Included in the data are the ages of the individuals, the provinces where the specimens were collected, the sex of the individuals, the dates of onset of symptoms, and the dates of the reporting of specimens. The information on the date of symptom onset was reported for 758 cases (6%). We use multiple imputation techniques to create 500 different datasets with the missing onset times imputed, as predicted by the province and an indicator of whether the specimen was collected on a weekday or weekend, using Poisson regression [22]. We report the averages and ranges over the 500 imputed datasets. Contact tracing information was collected on 100 initial cases, to provide an estimate of the serial interval, as has been previously reported [21]. We only use data from Gauteng province (n = 5579, 44% of cases) to avoid confounding the results with potential spatial variation in transmission. Gauteng province is the most populous, yet smallest geographically, of the nine provinces in South Africa, with over 10 million inhabitants, predominantly in the cities of Johannesburg and Pretoria.
Statistical methods
Wallinga and Teunis [20] (denoted WT method hereafter) proposed a method for the estimation of the effective reproductive number by making use of the epidemic curve, N = {N_{1},…, N_{T}}, where N_{t} is the number of cases at time point t, and an estimate of the serial interval, p_{1},…, p_{k}, where p_{i} describes the probability of a serial interval of length i and the maximum serial interval length is k. We review this method in Appendix 1. The estimator they obtain for the effective reproductive number for individual j on day t’ is
where n_{s} denotes the number with symptom onset on day s and q_{s,t} denotes the relative probability that case s was infected by case t.
Age transmission data
We propose the use of additional structure in this method to describe the probability of an infection event occurring between two cases that incorporates information on their ages by modifying the probability of transmission to be:
where a_{j} is the age group of individual j and ${\mathit{w}}_{{\mathit{a}}_{\mathit{i}}\mathit{\text{'}}{\mathit{a}}_{\mathit{j}}}$ is a measure of the likelihood of transmission between individuals in age group a_{i} and a_{j}. The matrix W = {${\mathit{w}}_{{\mathit{a}}_{\mathit{i}}{\mathit{a}}_{\mathit{j}}}$} does not necessarily have to be symmetric.
This method requires information on the likelihood of infectious contact between different age groups, or the ${\mathit{w}}_{{\mathit{a}}_{\mathit{i}}{\mathit{a}}_{\mathit{j}}}$. Increasingly studies are being conducted to obtain such information by assuming that transmission is directly related to contact patterns. We use the results of two such studies:
The first is a study of 571 randomly selected individuals in a South Africa township performed in 2010 and reported by JohnstoneRobertson et al. [12]. The authors report two matrices with age specific contact patterns in five year intervals up to a 45+ category. The first matrix only considers contacts that involve all close contacts while the second includes information on only those contacts that involve physical touch.
The second set of matrices we use comes from the European based PolyMod study of Mossong et al. [5]. This study includes information on 97,904 contacts amongst 7,290 participants from eight countries in Europe: Belgium, Finland, Great Britain, Germany, Italy, Luxembourg, the Netherlands, and Poland. Contact matrices describe all close contacts, and then separately, close contacts that involve physical touch. The matrices report agespecific values for five year age groups up to 70+. We modify these matrices to match those presented by JohnstoneRobertson et al. [12] and to match the demographics of South Africa’s young population by averaging all values above 45 years of age to create a single 45+ age category. In our results we focus on those obtained using the contact matrices from South Africa, as these matrices would seem more appropriate for the data at hand. We report results from the PolyMod matrices as a sensitivity analysis.
We estimate R_{t} and R_{0} using the 18 matrices described above with the imputed epidemic data from South Africa, and report age specific estimates of these quantities, as well as aggregate estimates across age groups. The reproductive numbers for each age group represent the expected number of infections generated across the population by an individual in that particular age group.
In the Appendices, we further report the results of two sensitivity analyses: First we test the sensitivity of the results to potential errors in the reporting dates by selecting a single imputed dataset and randomly jittering the onset dates of 10% of the individuals, within 30 days of their observed (or imputed) onset date (Appendix 2). We create 50 such datasets and repeat all analyses on these datasets and compare these results to those obtained without jittering the data. The second sensitivity analysis tests the impact of differential healthcare seeking behaviors by age. We smooth the distribution of the proportion of cases that were hospitalized by age group to serve as a surrogate distribution of healthcare seeking behavior and/or reporting patterns by age. This distribution is Ushaped, indicating that the very young and very old are more likely to seek medical care, a finding that has been reported elsewhere [23]. We attach various weights to this distribution and augment 25 of our imputed datasets according to this distribution. We reanalyze this augmented data to determine the potential impact of differential case reporting by age group on the results (Appendix 3).
Results
Figure 1 provides the epidemic curves across all age groups. Here, school age children and young adults tend to have the greatest number of cases initially in the outbreak.
Figure 2 shows the estimated effective reproductive numbers for the two JohnstoneRoberston et al. matrices. The ${\widehat{\mathit{R}}}_{\mathit{t}}$ estimates are similar regardless of the type of contact matrix assumed (close contact versus those involving physical contact). Overall, ${\widehat{\mathit{R}}}_{\mathit{t}}$ is much higher for those in the 15–19 and 20–24 year old groups throughout much of the epidemic, with the 10–14 and 25–29 year old age groups rapidly achieving high values, as well. Those over 45 initially have fairly high estimates of ${\widehat{\mathit{R}}}_{\mathit{t}}$ but these taper off quickly. Estimates of ${\widehat{\mathit{R}}}_{\mathit{t}}$ are not obtainable for those between 5 and 9 and those less than 5 until the outbreak is well under way, due to the paucity of observed cases for those age groups early on in the epidemic.
We obtain estimates of R_{0}, the basic reproductive number, by averaging the estimates of R_{t} during the epidemic period. In reality this can be viewed as a pseudoR_{0} given the prior immunity to this strain of influenza. We will refer to it as R_{0} throughout the text. We assume that the epidemic period corresponds to the point at which transmission was sustained in Gauteng Province until the overall number of cases peaked. This corresponds to the period between 22 June 2009 and 21 August 2009. Figure 3 and Table 1 show the estimates of R_{0} across age groups along with the number of individuals in each age group who were reported infected throughout the epidemic. Regardless of the choice of matrix, supercritical values of R_{0} are obtained for those between the age of 5 and 24, with the highest values being observed for those in the 10–14 agegroup (R_{0} = 1.53 for close contacts).
We contrast these estimates with those obtained using contact matrices from Europe [5]. Figures 3a and 3b shows the estimates of R_{0} across the 10 age groups obtained when using contact patterns from South Africa and the eight European countries in the PolyMod study for all close contacts (Figure 3a) and all contacts involving physical touch (Figure 3b). There are few notable differences between the estimates. In Figure 4, the mean estimates of R_{0} are shown for each age group. We observe a similar overall trend for the estimate of R_{0} across the age groups.
Finally we provide the overall estimate of R_{0} collapsed over all age groups (Table 2). For comparison purposes, we first estimate R_{0} by using a traditional analysis that assumes homogenous mixing among the age groups (${\widehat{\mathit{R}}}_{0}$=1.28, range: 1.261.31). This is similar to that obtained for all the other contact matrices considered. Additionally there are virtually no differences observed between results from the two contact matrices.
Our first sensitivity analysis, which jitters the onset dates of a subset of the population, (Additional file 1: Table S1) provides results that are consistent with the results presented. Not surprisingly, the impact of reassigning onset dates to a portion of the dataset has the impact of flattening the epidemic curve and thus lowering the estimates of R_{0}. However the results remain consistent and, coupled with the imputation variability reported, provide insight on the overall variability of the estimates reported.
In our second sensitivity analysis (Additional file 2: Table S2 and Additional file 3: Figure S1), examining the potential impact of differential reporting by age, we note a dramatic impact on the results. As we assume a greater underreporting of cases among those who are middle aged, we estimate the bulk of transmission being attributable to those who are older and less transmission being attributed to the very young, a finding contrary to the original results we present.
Discussion
We present a novel approach to estimating the effective and basic reproductive number by age group, and have applied this method to data from the 2009 Influenza H1N1pdm in Gauteng Province in South Africa. This method requires some estimate of contact patterns between age groups. We show results for 18 different possible contact matrices and the impact that these matrices have on the estimates. Additionally, as with the original method proposed by Wallinga and Teunis [20], it is necessary to have an estimate of the serial interval and we use an estimate obtained from contact trace data in South Africa.
As has been previously noted, that the burden of disease appears to be greatest amongst the young [21], a finding consistent with other studies [1, 4, 7]. These data argue that aiming interventions at youth would target the group that carries the largest burden and should have the best chance of success in limiting transmission. This finding is consistent with the strategy proposed by Wallinga et al. [18] and provides further information in the form of actual estimates of the reproductive number.
Our results also illustrate the importance of accounting for the age structure when estimating reproduction numbers. Though the overall estimates of R_{0} are unaffected by the incorporation of this information, we obtain much richer information with the ability to obtain agespecific estimates of the reproductive number. This analysis provides greater insight into the dynamics of disease transmission and informs intervention strategies.
The results obtained using information on transmission dynamics from the study based in South Africa [12] as well as that of the PolyMod study [5], appear to corroborate previous results for influenza pandemics [2, 3, 19] which seem to imply that school aged children are responsible for the bulk of disease transmission. Specifically we estimate that ${\widehat{\mathit{R}}}_{0}$ is highest for 10–14 year olds when using South African contact trace matrices (close contacts: ${\widehat{\mathit{R}}}_{0}$ = 1.53, range, 1.491.58; physical touch contacts: ${\widehat{\mathit{R}}}_{0}$ = 1.47, range, 1.441.51). These results are similar to the results from other studies [2, 19] and those obtained using the European based contact matrices.
Interestingly, it does not appear to make a substantial difference which contact pattern matrix we use in our analysis. One would assume that the matrices obtained in South Africa would be most relevant to the outbreak data we are analyzing and indeed, we have chosen to present the majority of our results using these matrices. We note, however, that when we use contact patterns from other European countries, where the demographics, climate, healthcare system, government, overall health, etc. are different from that of South Africa, there are only minor changes in the results. Indeed, the contact patterns observed in the Polymod study and the South African contact study are not substantially different, however they are not identical. This appears to argue that using some form of adjustment is superior to assuming homogenous mixing, but the method we propose is not overly sensitive to the form the adjustments take [24]. This result is similar to that of Glass et al. [2] who experimented with four forms of next generation matrices to estimate the reproductive number for adults and children separately. They found that the estimates of the reproductive numbers were not overly sensitive to the matrix forms that they assumed.
However, one should still take care in the assumptions used when implementing this method, or others like it. Our study is only one instance and it is not clear that the results we obtain would replicate in other settings. For instance, if one were to always use the Polymod study information for studies throughout the world, there is still the potential for errors if contact patterns do differ dramatically from those observed in Europe. While it is impossible to know with certainty if this is the case without detailed contact pattern information for the area of study, one can, at the least conduct a sensitivity analysis to determine the potential impact of the contact matrix on the analysis. Additionally we implicitly assume that contact patterns are directly related to transmission probabilities, an assumption that has yet to be rigorously tested. There is also work to show that contact patterns can change considerably during illness [25, 26]. Our work relies on the contact patterns of healthy individuals. While we suspect, based on our sensitivity analyses, that this will not have a substantial impact on our results, this is important to note.
It is important to note the caveats and limitations of this study. Our results would be impacted if reporting was inconsistent throughout the outbreak among the age groups. For instance if reporting was very good among one age group initially but declined in quality as the outbreak progressed, we can expect that our results would be biased [27]. In general the default assumption is that the pyramid of disease reporting described in [28, 29] is the same for all age groups. Our second sensitivity analysis (Appendix 3) explores the impact of this assumption and shows that if reporting or healthcare seeking behavior is much lower among middle aged groups than the very old and very young, our results will change dramatically. In the extreme case, we see that transmission is mostly attributable to those who are at least 30 years old and that the very young are unable to sustain transmission. While this result is contrary to what has been reported in the scientific literature to date, the potential for reporting inconsistencies that we explore are not unlikely, and have not been recognized and corrected for in other analyses that we are aware of. BrooksPollock et al. [23] report results from a survey conducted during the 2009 H1N1 Influenza pandemic in the UK that showed that the very young and very old were more likely to seek healthcare when ill. The impact of correcting surveillance data to accommodate this phenomenon was to shift the burden of illness from the very young to the middleaged. Further investigation into potential reporting inconsistencies is important to better understand infectious disease dynamics by age similar to what was previously done but not incorporating age [27]. Another reporting issue arises from silent infections, or those who carry infection and have the potential to transmit it, but are asymptomatic. We did not investigate the impact of these individuals, though the issues are similar to those we have just described.
Additional reporting inconsistencies are possible spatially or across other socioeconomic factors. Our analysis was only performed on data from Gauteng province, the most urban province in South Africa. It is possible that reporting would not be as dramatically variable as it would be if we were to make use of data from the entire country. It is also important to note that we chose to limit our analysis to Gauteng province so as to limit the impact of spatial effects and make the assumption of homogenous mixing more reasonable. This could limit generalizability.
We have also assumed that the contact matrices we use are correct and do not allow for any uncertainty in their estimation. These results might be improved upon and made more realistic by allowing for greater stochastic effects and/or flexibility in the transmission matrix. Ideally we would estimate these parameters in our study, but we do not have sufficient data to do so in the present framework. Glass et al. [2] have shown how to do this for a matrix with adults and children, but are limited to two by two matrices that presume a prespecified structure and are unable to consider a larger number of age groups, thus limiting their ability to gain a more thorough and detailed understanding of transmission.
Conclusions
We have applied a novel method to estimate transmission patterns between individuals from different age groups during the 2009 Influenza H1N1pdm in South Africa. We show that assumptions regarding the assumed contact patterns between age groups do not substantially impact the conclusions one draws from the data analyses in our study. Our results are consistent with other studies that show children are much more likely to become ill and transmit disease than adults during a pandemic, if the completeness of the data reported is independent of the age of the patients. These methods can be used to estimate heterogeneity in transmission parameters in real time by using the modification proposed by Cauchemez et al. [30] and thus inform the use of targeted interventions by age group.
Appendix 1
Wallinga and Teunis method
Wallinga and Teunis (20) (denoted WT method hereafter) proposed a method for the estimation of the effective reproductive number by making use of the epidemic curve, N = {N_{1},…, N_{T}}, where N_{t} is the number of cases at time point t, and an estimate of the serial interval, p_{1},…, p_{k}, where p_{i} describes the probability of a serial interval of length i and the maximum serial interval length is k. For ease of presentation, we assume that the time step is a day. Individuals are placed in a network temporally by symptom onset date and the probability of transmission occurring between two individuals in the network is determined by the serial interval. The calculation of R_{t} occurs in three steps. In what follows, we let t_{i} denote the i^{th} individual with symptom onset on day t, where i = 1,…,N_{t}.

1.
For the i^{th} individual with symptom onset on day t, calculate the probabilities of infection by all those with symptom onset on prior days t^{2} (t2 < t.) These probabilities equal the serial interval probability for the distance in time between the potential infector, t_{j}2 , and infectee, t_{i}, P (t_{j}' →t_{i}) = ${\mathit{p}}_{{\mathit{t}}_{\mathit{i}}{\mathit{t}}_{\mathit{j}}^{\text{'}}}$.

2.
Calculate the relative probability that the case t_{i} was infected by the j^{th} case on day P (t_{j}' →t_{i}), denoted by ${\mathrm{q}}_{{\mathrm{t}}_{\mathit{i}},{\mathrm{t}}_{\mathit{j}}\text{'}}$,
$$\begin{array}{ll}{\mathrm{q}}_{{\mathrm{t}}_{\mathit{i}},{\mathrm{t}}_{\mathit{j}}\text{'}}& =\frac{\mathit{P}\left({\mathit{t}}_{\mathit{j}}^{\mathit{\text{'}}}\to {\mathit{t}}_{\mathit{i}}\right)}{{\displaystyle {\sum}_{\mathit{s}=1}^{min\left(\mathit{k},\mathit{ti}1\right)}}{\displaystyle {\sum}_{\mathit{l}=1}^{{\mathit{n}}_{\mathit{s}}}}\mathit{P}\left({\mathit{s}}_{\mathit{l}}\to {\mathit{t}}_{\mathit{i}}\right)}\\ =\frac{\mathit{P}\left({\mathit{t}}_{\mathit{j}}^{\mathit{\text{'}}}\to {\mathit{t}}_{\mathit{i}}\right)}{{\displaystyle {\sum}_{\mathit{s}=1}^{min\left(\mathit{k},\mathit{ti}1\right)}}{\mathit{n}}_{\mathit{s}}\mathit{P}\left(\mathit{s}\to {\mathit{t}}_{\mathit{i}}\right)}.\end{array}$$ 
3.
Calculate the reproductive number. For the j^{th} case on day t', say t_{j}', the reproductive number is calculated as the sum of the expected values of a Bernoulli random variable. The Bernoulli random variable describes the event that t_{j}' infected another individual on day s, s > t', say s_{i}, and has probability q_{s,t'}. Then the effective reproductive number for individual t_{j}′is
$$\begin{array}{ll}{\mathit{R}}_{{\mathrm{t}}_{\mathit{j}}\mathit{\text{'}}}& =\phantom{\rule{0.5em}{0ex}}{\mathit{\Sigma}}_{\mathit{s}={\mathit{t}}^{\mathit{\text{'}}}+1}^{min\left(\mathit{T},{\mathit{t}}^{\mathit{\text{'}}}+\mathit{k}\right)}{\mathit{\Sigma}}_{\mathit{i}=1}^{{\mathit{n}}_{\mathit{s}}}{\mathit{q}}_{{\mathit{s}}_{\mathit{i}},{\mathit{t}}_{\mathit{j}}\mathit{\text{'}}}\\ =\phantom{\rule{0.5em}{0ex}}{\mathit{\Sigma}}_{\mathit{s}={\mathit{t}}^{\mathit{\text{'}}}+1}^{min\left(\mathit{T},{\mathit{t}}^{\mathit{\text{'}}}+\mathit{k}\right)}\phantom{\rule{0.5em}{0ex}}{\mathit{n}}_{\mathit{s}}{\mathit{q}}_{\mathit{s},{\mathit{t}}_{\mathit{j}}\mathit{\text{'}}},\end{array}$$
where n_{s} denotes the number with symptom onset on day s.
Appendix 2
Sensitivity analysis: Impact of errors in reporting dates
In this analysis, we choose a single imputed dataset and randomly jitter the onset dates of 10% of the sample within 30 days of their observed (or imputed) onset date. We create 50 such datasets and repeat all analyses on these datasets and compare them to the results on the nonjittered dataset.
Appendix 3
Sensitivity analysis: Impact of differential reporting by age
We assume that the reporting distribution by age follows a Ushaped distribution, implying that the very young and very old are most likely to seek healthcare and have their cases reported. To obtain a distribution that follows this shape, we use the distribution of hospitalized cases by age in our data, rescale it so that the highest proportion is one, and smooth the distribution using a loess smoother (Additional file 3: Figure S1).
We use 25 of our imputed datasets and augment each dataset using the distribution f (x), where f (x) is a function of the original age distribution observed in the data, g (x), and the reporting distribution shown in Additional file 3: Figure S1, h (x), as follows:
Here λ ranges between 0 and 1. We run analyses for λ = 0.0, 0.25, 0.50, 0.75 and 1.0 (corresponding to the original analysis). Results for all 25 datasets are shown in Additional file 2: Table S2.
References
 1.
Nishiura H, Oshitani H: Household transmission of influenza (H1N12009) in japan: agespecificity and reduction of household transmission risk by zanamivir treatment. J Int Med Res. 2011, 39 (2): 619628. 10.1177/147323001103900231
 2.
Glass K, Mercer GN, Nishiura H, McBryde ES, Becker NG: Estimating reproduction numbers for adults and children from case data. J R Soc Interface. 2011, 8 (62): 12481259. 10.1098/rsif.2010.0679
 3.
Glass K, Kelly H, Mercer GN: Pandemic influenza H1N1: reconciling serosurvey data with estimates of the reproduction number. Epidemiology. 2012, 23 (1): 8694. 10.1097/EDE.0b013e31823a44a5
 4.
Bansal S, Pourbohloul B, Grenfell B, Meyers LA: The shifting demographic landscape of influenza. PLoS Curr. 2009, 1: RRN1047.
 5.
Mossong J, Hens N, Jit M, Beutels P, Auranen K, Mikolajczyk R, Massari M, Salmaso S, Tomba GS, Wallinga J, Heijne J, SadkowskaTodys M, Rosinska M, Edmunds WJ: Social contacts and mixing patterns relevant to the spread of infectious diseases. PLoS Med. 2008, 5 (3): e74. 10.1371/journal.pmed.0050074
 6.
Apolloni A, Poletto C, Colizza V: Agespecific contacts and travel patterns in the spatial spread of 2009 H1N1 influenza pandemic. BMC Infect Dis. 2013, 15: 13176, 233413176.
 7.
Del Valle SY, Hyman JM, Hethcote HW, Eubank SG: Mixing patterns between age groups in social networks. Social Networks. 2007, 29 (4): 53910.1016/j.socnet.2007.04.005. 10.1016/j.socnet.2007.04.005
 8.
Steens A, Waaijenborg S, Teunis PFM, Reimerink JHJ, Meijer A, van der Lubben M, Koopmans M, van der Sande MAB, Wallinga J, van Boven M: Agedependent patterns of infection and severity explaining the low impact of 2009 influenza A (H1N1): evidence from serial serologic surveys in the netherlands. Am J Epidemiol. 2011, 174 (11): 13071315. 10.1093/aje/kwr245
 9.
Jacobs JH, Archer BN, Baker MG, Cowling BJ, Heffernan RT, Mercer G, Uez O, Hanshaoworakul W, Viboud C, Schwartz J, Tchetgen E, Lipsitch M: Searching for sharp drops in the incidence of pandemic A/H1N1 influenza by single year of age. PLoS One. 2012, 7 (8): e42328. 10.1371/journal.pone.0042328
 10.
Read JM, Edmunds WJ, Riley S, Lessler J, Cummings DA: Close encounters of the infectious kind: methods to measure social mixing behaviour. Epidemiol Infect. 2012, 140 (12): 21172130. 10.1017/S0950268812000842
 11.
Horby P, Pham QT, Hens N, Nguyen TT, Le QM, Dang DT, Nguyen ML, Nguyen TH, Alexander N, Edmunds WJ, Tran ND, Fox A, Nguyen TH: Social contact patterns in vietnam and implications for the control of infectious diseases. PLoS One. 2011, 6 (2): e16965. 10.1371/journal.pone.0016965
 12.
JohnstoneRobertson SP, Mark D, Morrow C, Middelkoop K, Chiswell M, Aquino LD, Bekker LG, Wood R: Social mixing patterns within a south african township community: implications for respiratory disease transmission and control. Am J Epidemiol. 2011, 174 (11): 12461255. 10.1093/aje/kwr251
 13.
Wallinga J, Teunis P, Kretzschmar M: Using data on social contacts to estimate agespecific transmission parameters for respiratoryspread infectious agents. Am J Epidemiol. 2006, 164 (10): 936944. 10.1093/aje/kwj317
 14.
Lipsitch M, Finelli L, Heffernan RT, Leung GM, Redd SC: 2009 H1n1 Surveillance Group. Improving the evidence base for decision making during a pandemic: the example of 2009 influenza A/H1N1. Biosecur Bioterror. 2011, 9 (2): 89115.
 15.
Brownstein JS, Kleinman KP, Mandl KD: Identifying pediatric age groups for influenza vaccination using a realtime regional surveillance system. Am J Epidemiol. 2005, 162 (7): 686693. 10.1093/aje/kwi257
 16.
Schanzer D, Vachon J, Pelletier L: Agespecific differences in influenza A epidemic curves: do children drive the spread of influenza epidemics?. Am J Epidemiol. 2011, 174 (1): 109117. 10.1093/aje/kwr037
 17.
Sebastian R, Skowronski DM, Chong M, Dhaliwal J, Brownstein JS: Agerelated trends in the timeliness and prediction of medical visits, hospitalizations and deaths due to pneumonia and influenza, British Columbia, Canada, 1998–2004. Vaccine. 2008, 26 (10): 13971403. 10.1016/j.vaccine.2007.11.090
 18.
Wallinga J, van Boven M, Lipsitch M: Optimizing infectious disease interventions during an emerging epidemic. Proc Natl Acad Sci U S A. 2010, 107 (2): 923928. 10.1073/pnas.0908491107
 19.
Nishiura H, CastilloChavez C, Safan M, Chowell G: Transmission potential of the new influenza A(H1N1) virus and its agespecificity in japan. Euro Surveill. 2009, 14 (22): 19227.
 20.
Wallinga J, Teunis P: Different epidemic curves for severe acute respiratory syndrome reveal similar impacts of control measures. Am J Epidemiol. 2004, 160 (6): 509516. 10.1093/aje/kwh255
 21.
Archer B, Cohen C, Naidoo D, Thomas J, Makunga C, Blumberg L, Venter M, Timothy GA, Puren A, McAnerney JM, Cengimbo A, Schoub BD: Interim report on pandemic H1N1 influenza virus infections in South Africa, April to October 2009: epidemiology and factors associated with fatal cases. Euro Surveill. 2009, 14 (42): 19369.
 22.
Archer BN, Tempia S, White LF, Pagano M, Cohen C: Reproductive number and serial interval of the first wave of influenza A (H1N1) pdm09 virus in South Africa. PlosOne. 2012, 7 (11): e4948210.1371/journal.pone.0049482. doi:10.1371/journal.pone.0049482. 10.1371/journal.pone.0049482
 23.
BrooksPollock E, Tilston N, Edmunds WJ, Eames KT: Using an online survey of healthcareseeking behaviour to estimate the magnitude and severity of the 2009 H1N1v influenza epidemic in England. BMC Infect Dis. 2011, 11: 68doi:10.1186/147123341168
 24.
White LF, Archer B, Pagano M: Estimating the reproductive number in the presence of spatial heterogeneity of transmission patterns. Int J Health Geogr. 2013, 12: 35doi:10.1186/1476072X1235
 25.
Eames KT, Tilston NL, White PJ, Adams E, Edmunds WJ: The impact of illness and the impact of school closure on social contact patterns. Health Technol Assess. 2010, 14 (34): 267312.
 26.
Van Kerckhove K, Hens N, Edmunds WJ, Eames KT: The impact of illness on social networks: Implications for transmission and control of influenza. Am J Epidemiol. 2013, 178 (11): 16551662. 10.1093/aje/kwt196
 27.
White LF, Pagano M: Reporting errors in infectious disease outbreaks, with an application to pandemic influenza A/H1N1. Epidemiol Perspect Innov. 2010, 7: 12. 10.1186/17425573712
 28.
Reed C, Biggerstaff M, Finelli L, Koonin LM, Beauvais D, Uzicanin A, Plummer A, Bresee J, Redd SC, Jernigan DB: Estimates of the prevalence of pandemic (H1N1) 2009, United States, apriljuly 2009. Emerg Infect Dis. 2009, 15 (12): 20042007. 10.3201/eid1512.091413
 29.
Mead PS, Slutsker L, Dietz V, McCaig LF, Bresee JS, Shapiro C, Griffin PM, Tauxe RV: Foodrelated illness and death in the United States. Emerg Infect Dis. 1999, 5 (5): 607625. 10.3201/eid0505.990502
 30.
Cauchemez S, Boelle PY, Thomas G, Valleron AJ: Estimating in real time the efficacy of measures to control emerging communicable diseases. Am J Epidemiol. 2006, 164 (6): 591597. 10.1093/aje/kwj274
Acknowledgments
The project described was supported by Award Number U54GM088558 from the National Institute of General Medical Sciences to L.F.W. and M.P. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institute of General Medical Sciences or the National Institutes of Health.
Author information
Affiliations
Corresponding author
Additional information
Competing interests
The authors report no conflicts of interest.
Authors’ contributions
LFW and MP conceived the project and developed the methods. LFW performed the analyses. BA provided the data and insight on dynamics in South Africa. LFW drafted the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
12982_2013_122_MOESM1_ESM.doc
Additional file 1: Table S1: Results for the sensitivity analysis using the South African based age contact information. Result presented is the estimate obtained from the original dataset and the values in the parentheses represent the range of values obtained over the 50 datasets generated for the sensitivity analysis. (DOC 56 KB)
12982_2013_122_MOESM2_ESM.doc
Additional file 2: Table S2: Estimates of R_{0} using 25 of the 500 imputations described in the original text. Results shown are the mean and range of estimates across the 25 imputed datasets. λ=1.00 corresponds to the results from the original analysis. (DOC 84 KB)
12982_2013_122_MOESM3_ESM.zip
Additional file 3: Figure S1: Smoothed distribution to reflect potential rates of healthcare seeking behavior and/or casereporting by age. (ZIP 1 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
This article is published under license to 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 credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
White, L.F., Archer, B. & Pagano, M. Determining the dynamics of influenza transmission by age. Emerg Themes Epidemiol 11, 4 (2014). https://doi.org/10.1186/17427622114
Received:
Accepted:
Published:
Keywords
 Pandemic influenza H1N1
 Reproductive number
 Infectious disease