Volume 11, Number 9—September 2005
Potential Impact of Antiviral Drug Use during Influenza Pandemic
The recent spread of highly pathogenic strains of avian influenza has highlighted the threat posed by pandemic influenza. In the early phases of a pandemic, the only treatment available would be neuraminidase inhibitors, which many countries are considering stockpiling for pandemic use. We estimate the effect on hospitalization rates of using different antiviral stockpile sizes to treat infection. We estimate that stockpiles that cover 20%–25% of the population would be sufficient to treat most of the clinical cases and could lead to 50% to 77% reductions in hospitalizations. Substantial reductions in hospitalization could be achieved with smaller antiviral stockpiles if drugs are reserved for persons at high risk.
Recent outbreaks of highly pathogenic avian influenza in poultry in East Asia (H5N1), Canada (H7N3), and the Netherlands (H7N7), and their subsequent transmission to humans, have intensified concern over the emergence of a novel strain of influenza with pandemic potential. Three influenza pandemics occurred during the 20th century, with varying degrees of severity; outcomes ranged from the high levels of illness and death observed during the 1918 Spanish flu pandemic (estimates of deaths range from 20 to 100 million ) to the much lower levels observed during the pandemics of 1957 and 1968 (≈1 million deaths each ). While recognizing that the characteristics of future influenza pandemics are difficult to predict, the World Health Organization (WHO) has recommended that nations prepare pandemic contingency plans (3). Several have been drafted, and some have been published (4–7), although all are subject to continuous refinement. Surveillance, on both a local and global scale, will enable policy makers and practitioners to act during the early phases of a pandemic. However, the likely rapid global spread of a pandemic strain will limit the time available to implement appropriate mitigating strategies, and preemptive contingency planning is needed.
A number of intervention strategies can reduce the impact of influenza pandemics. During interpandemic years, influenza vaccination is used to reduce deaths and disease. However, vaccine is unlikely to be available in time or in sufficient quantities for use during a pandemic (8,9). Other, nontherapeutic, disease control options may be used, such as those used during the outbreak of severe acute respiratory syndrome (10).
However, 2 groups of antiviral drugs are available for the treatment and prophylaxis of influenza. These are the adamantanes (amantadine and rimantadine) and the neuraminidase inhibitors (oseltamivir and zanamivir). The adamantanes may be effective against pandemic strains, but concern exists about adverse reactions and the development of antiviral resistance. Resistance to amantadine has been demonstrated in a number of avian H5 strains (11) and its use for treatment of influenza is not recommended (12).
The neuraminidase inhibitors (NIs) reduce the period of symptomatic illness from both influenza A and B viruses (13) and both are recommended for use in the United Kingdom for treatment of at-risk adults who are able to begin treatment within 48 hours of onset of symptoms. Oseltamivir is also recommended for the treatment of at-risk children >12 months of age (12). The development of antiviral resistance has been reported for NIs, particularly related to oseltamivir use in children (14), although current evidence suggests that resistant strains are pathogenically weakened (15). The use of NIs for treatment of pandemic influenza remains an option since they may improve individual disease outcomes and the effect of the disease in the population.
An influenza pandemic is likely to increase demands on healthcare providers, especially in hospitals. Except in Japan, current levels of NI use are low. Any strategy involving NI use would require stockpiles of these drugs. The potential use of antiviral agents for prophylaxis has been investigated elsewhere and may be of greatest use in the earliest phases of a pandemic to retard the spread of the virus (16,17). Earlier pandemic influenza modeling studies have also focused on the economic effect of vaccination (18) and the use of NI prophylaxis for disease control (19). We assessed the potential effect of using NIs for treatment on the estimated number of influenza-related hospitalizations likely to occur during a pandemic. Unlike in previous studies (20), we have also taken into account the reduction in infectivity that antiviral treatment may have on community transmission.
Our models focused on using NIs to treat different age and risk groups and the potential effects treatment might have on influenza hospitalizations. These effects have been quantified by using the mathematical model described in the Appendix. The length of the latent, noninfectious period was assumed to be 2 days (19), and the infectious period was assumed to be 4 days (19,21). Hospitalization rates for the basline scenario were calculated by using data from interpandemic influenza and are given for different and age risk groups (Table 1).
To be effective, NI treatment must be administered within 48 hours of symptom onset. The efficacy of NI treatment appears to prevent 50% of hospitalizations, mirroring efficacy rates against developing complications; this efficacy rate is approximately the same for oseltamivir and zanamivir (13). Symptoms were also reduced by ≈1.5 days; treatment was assumed to produce the same decrease in the infectious period.
The population was stratified as for seasonal influenza; persons were considered to be either at high risk for severe outcome or at low risk (22). The at-risk group included those with chronic respiratory disease, chronic heart disease, chronic renal failure, diabetes mellitus, and immunosuppression; this group also included all persons living in long-term care facilities, such as nursing homes (23), and all those >65 years of age (24).
Demographic data used in the model were based on age-specific distribution of the UK population (Office for National Statistics. http://www.statistics.gov.uk). The model was used to simulate a number of scenarios, on the basis of contingency plans and previous pandemics, to investigate the effect of targeting NIs to different age and risk groups on the expected number of hospitalizations during a pandemic.
The baseline scenario for this study was that advocated by WHO (3) and was also used previously by Meltzer et al. (18). This scenario assumes a clinical attack rate, in the absence of interventions, of 25% of the population, which occurs during a single wave. Assuming that half of infections are nonclinical or asymptomatic (i.e., a serologic attack rate across the population of 50%) (25), a value for the basic reproduction number, R0, of 1.39 can be calculated. When these parameters are used in the model in the Appendix, the effect of different-sized antiviral stockpiles on the overall clinical attack rate can be estimated.
The outputs from the first set of simulations are shown in Figure 1. The baseline scenario is shown alongside a range of other clinical attack rates (20%–40%) (i.e., varying R0 from 1.28 to 2.0) in the absence of interventions. For these scenarios, antiviral treatment is assumed to be possible within 48 hours of onset for all symptomatic patients until the stockpile is exhausted, with the exception of those <1 year of age, who are not treated at any stage (treatment for this age group is contraindicated ). The points on the curves in Figure 1, where the gradients change from vertical to horizontal, indicate the points at which the stockpile is sufficient to treat all patients; increasing the stockpile size would produce no additional benefit and would therefore result in a surplus of antiviral treatments.
For the baseline scenario, a stockpile large enough to treat 12% of the population (i.e., a 12% stockpile) would be sufficient to treat all patients, even if the clinical attack rate in the absence of treatment is 25%. This difference is due to a reduction in the effective reproduction number of the disease, Rε, caused by shortening the infectious period of those treated by 1.5 days. Across the different attack rates, stockpiles sufficient to treat <1% of the population are unlikely to result in major changes to disease dynamics. Outputs are most sensitive to the clinical attack rate when the reduction in the infection period caused by treatment is sufficient to bring Rε <1. When Rε is <1, the number of secondary cases produced by each person is <1, and incidence, therefore, decreases. The value of Rε can be calculated as
where S is the proportion of the population susceptible. With treatment, this equation can be rewritten as
where It is the decrease in the infectious period due to treatment, Ip the infectious period, and ci the proportion of infections in each of the different population subgroups, i, that are treated. For the scenarios in Figure 1, It = 1.5 days, Ip = 4.0 days and ci = 0.5 for all groups except those <1 year of age, who only constitute 1.1% of the population. Therefore, the term within the brackets for this scenario can be calculated as 0.81. At the start of the pandemic, S is assumed to be 1; therefore, if R0 is <1.23, the outbreak can be controlled by treating all patients. For pandemics in which R0 is >1.23, depletion of susceptible persons through infection is also required before Rε decreases to <1, which is equivalent to S = (0.81R0) – 1.
The effect of different treatment strategies on hospitalization rates was generated from the baseline scenario: treating all patients, only at-risk groups, only children and the elderly (1–14 and >65 years of age), and only the working population (15–64 years of age). These scenarios were of potential interest to public health planners; outputs are shown in Figure 2. Given a large enough stockpile, the best option to minimize hospitalizations would be to treat all patients; for this scenario, a 12% antiviral coverage would reduce hospitalizations by up to 77%. An alternative strategy of treating the whole working population reduces the hospitalization rate by up to 40% but requires a similar antiviral stockpile size, and treating the working population consistently fails to reduce the number of hospitalizations below the number that would be expected if everyone were treated, regardless of stockpile size. This increase is because the hospitalization rate for the working population is less than the average in the population and also because treating a smaller proportion of the population has less effect on the overall transmission rate. For stockpile sizes only large enough to treat <5% of the population, the best strategy would be to treat at-risk groups; this strategy is also best for stockpile sizes up to 7%, with hospitalizations at this level reduced by up to 45%. For stockpile sizes from 7% to 10%, the best strategy is to treat children and the elderly (reducing hospitalizations by up to 48%) and for stockpile sizes >10%, to treat everyone.
The optimum treatment strategy is therefore dependent on treating those at highest risk for hospitalization. The simulations for the baseline scenario were based on a uniform age-specific attack rate and on age- and risk-specific hospitalization rates from interpandemic years because of the uncertainty over the precise characteristics of a future pandemic. Since the age-specific clinical attack rate has varied between pandemics, we repeated the analysis above, as far as possible, using the age-specific attack rates from previous pandemics (26–28) (Table 2) for comparison with the baseline scenario.
The 1957 UK pandemic began with imported cases in July 1957; deaths peaked in November 1957, with a reported overall clinical attack rate of 31% (26). The proportion of infections resulting in clinical illness was calculated from a small serologic survey of general practitioners; only 46% of the general practitioners surveyed with a positive antibody titer actually had symptoms (26). The serologic attack rate was calculated as 67%, which would require R0 = 1.65. The epidemic curve that this figure would generate is shown in Figure 3A, with the curve scaled to fit the 1957 epidemic curve for deaths (26). The only additional change from the baseline scenario is the 1957 hospitalization rate, which was reported to be 188/100,000 population (26). Using the age-specific attack rates for 1957 (Table 2) in the model, we scaled hospitalization rates to achieve an overall hospitalization rate of 188/100,000 (Table 3).
The results (Figure 3B) show that a 20%–25% antiviral stockpile would be sufficient to treat all patients during the first wave, a figure that is larger than that seen for the baseline scenario, as both the clinical and serologic clinical attack rates were higher. However, qualitatively, the results are similar in spite of the differences in attack rates between different age groups. With a stockpile as large as 20%–25%, an estimated reduction in hospitalizations of ≈67% could be expected. As in the baseline scenario, effective targeting of smaller stockpiles to at-risk groups can also be used to produce large reductions in hospitalization rates. For stockpiles <11%, the best strategy is to treat those at risk, which results in a reduction of 36%. For stockpiles sizes from 11% to 17%, the best strategy is to treat the young and elderly, which results in a 39% reduction. The highest reduction from treating the working population is 31% and remains a suboptimal strategy for any stockpile size.
The implications of different treatment strategies on the hospitalization rates with a 10% stockpile are shown in Figure 3C. Strategies with larger proportions of the 10% stockpile had the greatest effect on the epidemic, steadily delaying, but not diminishing, the peak of hospitalizations. Treating only the working population results in a 15% decrease in hospitalizations, treating all patients results in a 22% decrease, and treating children and the elderly a 32% reduction. With each of these strategies, the antiviral stockpile is exhausted before the end of the pandemic, whereas the fourth strategy of treating at-risk groups reduces hospitalizations by 36% and only requires a 5% stockpile. Therefore, treating those at risk is the most efficient strategy, but further targeting may be considered to avoid surplus treatments.
The 1968 pandemic was characterized by 2 waves, the first relatively small, occurring from February to April 1969; the larger wave occurred from November 1969 to January 1970 (27). We predominately considered the second wave. A confounding factor is that a proportion of the population would have been immune because of the first wave. Weighting age-specific clinical attack rates (Table 1) by age-group sizes from census data, we calculated the overall clinical attack rates for the first and second waves to be 6% and 21%, respectively (27; Office for National Statistics [http://www.statistics.gov.uk]). The serologic attack rate was derived by fitting the model to the data for the second wave from the Royal College of General Practitioners (provided by Douglas Fleming; http://www.rcgp.org.uk); we assumed a similar proportion of asymptomatic cases in both waves. The fit of the model to the data is shown in Figure 4A, from which is derived a 15% residual immunity from the first wave and a 65% serologic attack rate for the second wave, which produces an effective reproduction number of 1.85 for the second wave. The overall hospitalization rate for the second wave was reported as 144 per 100,000 (29), and using the age-specific attack rates for 1968 in Table 2, we adjusted the values in Table 1 to fit this value.
The size of the stockpile required to treat all patients is ≈18% (which is relatively small compared to the 1957 pandemic because of the lower clinical attack rate), which leads to fewer patients being treated and less reduction in overall transmission. If all persons whose infections resulted in clinical illness (i.e., patients) were treated, the hospitalization rate would drop by ≈56% (Figure 4B). For the 1968 pandemic, the effects of the different antiviral targeting strategies were different than in the previous scenarios as a result of the different age-specific attack rates, which are shifted more towards the working population (Table 2). Thus, relatively small stockpiles are required to treat either the at-risk group or the young and elderly group (≈3% for each group), since most patients are in the working population and neither of these 2 groups. For stockpiles of up to 12%, treating the at-risk group is marginally better than treating the young and the elderly (37% reduction in hospitalization as opposed to 32%), and for stockpiles >12%, treating all clinical patients would be the best strategy.
The effects of the different treatment strategies with a 10% stockpile are shown in Figure 4C. Hospitalizations would drop by ≈29% if all patients were treated and by 16% if the working population were treated; both treatment strategies would lead to the stockpiles' being exhausted. As above, treating those at risk would reduce hospitalizations by 37%, whereas treating only children and the elderly would reduce hospitalizations by 32% and only require a 3% stockpile per group. Of these 4 strategies, treating the at-risk groups is the most efficient, but given surplus stockpile, further extension of the groups to be targeted may be considered.
The characteristics for the 1918 pandemic differ substantially from the other 2 in that 3 distinct waves occurred; the age-specific attack rates were highest for those in their teens, 20s, and 30s; and the mortality rates were higher (2). In addition, age-specific attack rates and mortality rates differed for each of the 3 waves (28). Modeling based on the 1918 pandemic was therefore considerably less straightforward than for the previous 2 pandemics, and an approach was taken to fit the transmission model to each of the 3 waves, separately. No cross-immunity was assumed between different waves since studies suggested only weak effects; indeed, some studies suggested greater susceptibility in the third wave if a person had had influenza in the first pandemic wave (28). Clinical attack rates were calculated from reported weekly mortality data and clinical case-fatality rates (28). Serologic attack rates were then fitted separately to each of the curves (Figure 5), from which values of R0 = 2.0, 1.55, and 1.7 were derived from each of the respective waves. The estimate for the second wave is lower than other estimates of ≈3 (30) derived from US cities and is probably because our estimates were derived from data from throughout England and Wales, thereby incorporating spatial heterogeneity.
Since hospitalization rates were not available for any of the 3 waves, we considered the effect of antiviral treatment on death. The potential efficacy of antiviral treatments in preventing death between waves may have differed, but it was assumed to provide 50% protection against death. This estimate was based on the assumption that 50% protection from the more serious outcomes of influenza can be translated to equivalent protection from death (20).
A pandemic with the characteristics of that in 1918 would, without antiviral treatment, produce an estimated number of deaths equivalent to ≈0.5% of the population across all 3 waves. However, a 20% stockpile sufficient to treat all patients across the 3 waves would result in ≈53% reduction in deaths. With a smaller stockpile of 10%, the reduction in deaths was only 17% because the stockpile becomes exhausted during the second wave, before most of the deaths occur (Figure 6).
The baseline scenario with an overall clinical attack rate of 25%, as currently advised by WHO (3), is roughly in accordance with data from previous pandemics. The general conclusion from our study is that antiviral treatments for 20% to 25% of the population are likely to be sufficient to treat all patients for pandemics with characteristics that have been observed to date. The size of the stockpile required will depend on the clinical attack rate of the pandemic and the R0 value.
However, with smaller stockpile sizes, substantial reductions in hospitalizations can be achieved through targeting. For the smallest stockpiles, the best strategy was to treat conventional influenza at-risk groups. Treating the young and elderly is only slightly less effective. Treating the working population may have benefits beyond reducing hospitalizations, such as reducing illness-related absenteeism, but it consistently fails to be the best strategy for reducing hospitalizations. For large stockpiles, treating all patients is consistently the best strategy in reducing hospitalization and transmission. When all patients are treated, the marginal effect of treatment on reduced transmission increases with the number of patients treated, until all patients have been treated.
Further studies regarding the effects of antiviral treatments would improve the robustness of the parameter estimates. In particular, better estimates on the efficacy of NI treatment against hospitalization and death rates for different age and risk groups and estimates on the reduction in the infectious period are required. Also, the issue of antiviral resistance needs to be resolved since it could compromise NI effectiveness.
The scenarios above assume that clinical patients were treated within 48 hours of onset of symptoms; however, in reality, some cases will be diagnosed or reported too late, and other patients will be administered drugs mistakenly. To maximize the benefits of antiviral treatment, patients should be strongly encouraged to seek treatment and treatment should be supported by sound clinical judgment and diagnostic capability. If high levels of treatment are not achievable, disproportionately higher hospitalization rates than those calculated here would ensue. In addition, identifying groups with higher transmission rates for targeting treatment would result in greater reductions in transmission than reported here.
Assessments will need to be recalculated in the earliest phases of a pandemic with real-time data to confirm or update the assumptions used and ensure that the model parameters are appropriate. Therefore, were a pandemic to occur, intensive analysis of its dynamics would be required at its start.
Mathematical Model Used To Calculate Outputs
The model used was based on Kermack and McKendrick (31), which has formed the basis of a number of models for both epidemic and pandemic influenza (32–34) and is implemented by using the set of differential equations given in equation 1.
where α = 1/L = 0.5, γ = 1/PP = 0.4, λ = 1/IP = 2/3, and β = R0/(PP + IP). LP represents the length of the latent period; PP represents the length of the nonsymptomatic infectious period; and IP represents the length of the infectious symptomatic period. S represents the total proportion susceptible, E the total proportion incubating, Pi the proportion from the total population in each group i within the first 2.5 days of their infectious period, Ii the proportion of total population in each group i within the final 1.5 days of their infectious period, and R the total proportion recovered and immune or dead. ci is the proportion of infections resulting in clinical cases, and Ti is the proportion of group i receiving treatment. The average number of secondary cases per primary case when the population is entirely susceptible is represented by R0, and the proportion of the population in each group i is given by Ni. The proportion treated in each group each week can be given as Ai(t) as
where t is time in days.
The proportion of the population within each group being hospitalized each week, Hi(t), can be calculated as
where ε is the efficacy of antiviral treatment against hospitalization and hi is the hospitalization rate for each group i.
Supplementary information on the probability of hospitalization in the absence of vaccination is available from the author (see Comments to the Authors).
Dr Gani is a mathematical modeler. His research interests are the impact of pandemic influenza and other emerging and reemerging infectious diseases on human populations and assessments of policy options available to mitigate these impacts.
We thank members of the UK Department of Health Steering Group for their comments and help with setting model parameters.
Financial support for this work was provided by the UK Health Protection Agency. The views expressed in this publication are those of the authors and not necessarily those of the Health Protection Agency.
- Johnson NP, Mueller J. Updating the accounts: global mortality of the 1918–1920 "Spanish" influenza pandemic.Bull Hist Med. 2002;76:105–15.
- Nguyen-Van-Tam JS, Hampson AW. The epidemiology and clinical impact of pandemic influenza.Vaccine. 2003;21:1762–8.
- World Health Organization. Influenza pandemic plan. The role of WHO and guidelines for national and regional planning. Geneva: The Organization; 1999.
- Patriarca PA, Cox NJ. Influenza pandemic preparedness plan for the United States.J Infect Dis. 1997;176(Suppl 1):S4–7.
- UK Department of Health. UK pandemic influenza contingency plan. March 2005 [cited 2005 Mar 1]. Available from http://www.dh.gov.uk/assetRoot/04/10/44/37/04104437.pdf
- Fock R, Bergmann H, Bußmann H, Fell G, Finke E-J, Koch U, Influenza pandemic: preparedness planning in Germany.Eurosurveillance.2002;7:1–5.
- Health Canada. Canadian Pandemic Influenza Plan, Health Canada; 2004. [cited 2005 Mar 1]. Available from http://www.phac-aspc.gc.ca/cpip-pclcpi/
- Webby RJ, Webster RG. Are we ready for pandemic influenza?Science. 2003;302:1519–22.
- Fedson DS. Pandemic influenza and the global vaccine supply.Clin Infect Dis. 2003;36:1552–61.
- Pang X, Zu Z, Xu F, Guo J, Gong X, Liu D, Evaluation of control measures implemented in the severe acute respiratory syndrome outbreak in Beijing, 2003.JAMA. 2003;290:3215–21.
- Wainright PO, Perdue ML, Brugh M, Beard CW. Amantadine resistance among hemagglutinin subtype 5 strains of avian influenza virus.Avian Dis. 1991;35:31–9.
- National Institute for Clinical Excellence. Full guidance on the use of zanamivir, National Institute for Clinical Excellence. Full guidance on the use of zanamivir, oseltamivir and amantadine for the treatment of influenza. Available from http://www.nice.org.uk/pdf/58_Flu_fullguidance.pdf (2005)
- Stiver G. The treatment of influenza with antiviral drugs.CMAJ. 2003;168:49–57.
- Kiso M, Mitamura K, Sakai-Tagawa Y, Shiraishi K, Kawakami C, Kimura K, Resistant influenza A viruses in children treated with oseltamivir: descriptive study.Lancet. 2004;364:759–65.
- McKimm-Breschkin J, Trivedi T, Hampson A, Hay A, Klimov A, Tashiro M, Neuraminidase sequence analysis and susceptibilities of influenza virus clinical isolates to zanamivir and oseltamivir.Antimicrob Agents Chemother. 2003;47:2264–72.
- Balicer RD, Huerta M, Grotto I. Tackling the next influenza pandemic.BMJ. 2004;328:1391–2.
- Nguyen-Van-Tam JS, Leach SA, Cooper B, Gani R, Goddard NJ, Watson JM, Tackling the next influenza pandemic: ring prophylaxis may prove useful early on, but is unlikely to be effective or practical to implement once the pandemic is established. BMJ. B 22 July 2004]. [Cited 2005 Mar 1]. Available from http://bmj.bmjjournals.com/cgi/eletters/328/7453/1391#68042
- Meltzer MI, Cox NJ, Fukuda K. The economic impact of pandemic influenza in the United States: priorities for intervention.Emerg Infect Dis. 1999;5:659–71.
- Longini IM, Halloran ME, Nizam A, Yang Y. Containing pandemic influenza with antiviral agents.Am J Epidemiol. 2004;159:623–33.
- van Genugten ML, Heijnen ML, Jager JC. Pandemic influenza and healthcare demand in the Netherlands: scenario analysis.Emerg Infect Dis. 2003;9:531–8.
- 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 Med. 2004;23:3469–87.
- Fleming D, Charlton J, McCormick A. The population at risk in relation to influenza immunisation policy in England and Wales.Health Trends. 1997;29:42–7.
- Department of Health. Influenza immunisation. CMO's update. 1997. [cited 2005 Mar 1]. Available from http://www.dh.gov.uk/assetRoot/04/01/35/74/04013574.pdf
- Monto AS. Influenza: quantifying morbidity and mortality.Am J Med. 1987;82:20–6.
- Mann PG, Pereira MS, Smith JW, Hart RJ, Williams WO. A five-year study of influenza in families. Joint Public Health Laboratory Service/Royal College of General Practitioners Working Group.J Hyg (Lond). 1981;87:191–200.
- The influenza epidemic in England and Wales 1957–58. Reports on public health and medical subjects no. 100. London: Her Majesty's Stationary Office; 1960.
- Taylor MP. Influenza 1968–1970 incidence in general practice based on a population survey.J R Coll Gen Pract. 1971;21:17–22.
- Reports on public health and medical subjects no. 4: report on the pandemic of influenza 1918–19. London: Her Majesty's Stationary Office; 1920.
- Barker WH, Mullooly JP. Impact of epidemic type A influenza in a defined adult population.Am J Epidemiol. 1980;112:798–811.
- Mills CE, Robins JM, Lipsitch M. Transmissibility of 1918 pandemic influenza.Nature. 2004;432:904–6.
- Kermack WO, McKendrick AG. A contribution to the mathematical theory of epidemics. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences. 1927 Series A;115:700–21.
- Rvachev LA, Longini IM. A mathematical model for the global spread of influenza.Math Biosci. 1985;75:3–22.
- Flahault A, Letrait S, Blin P, Hazout S, Menares J, Valleron AJ. Modelling the 1985 influenza epidemic in France.Stat Med. 1988;7:1147–55.
- Flahault A, Deguen S, Valleron AJ. A mathematical model for the European spread of influenza.Eur J Epidemiol. 1994;10:471–4.