- Open Access
Cost-effective length and timing of school closure during an influenza pandemic depend on the severity
© Nishiura et al.; licensee BioMed Central Ltd. 2014
- Received: 9 December 2013
- Accepted: 20 January 2014
- Published: 21 January 2014
There has been a variation in published opinions toward the effectiveness of school closure which is implemented reactively when substantial influenza transmissions are seen at schools. Parameterizing an age-structured epidemic model using published estimates of the pandemic H1N1-2009 and accounting for the cost effectiveness, we examined if the timing and length of school closure could be optimized.
Age-structured renewal equation was employed to describe the epidemic dynamics of an influenza pandemic. School closure was assumed to take place only once during the course of the pandemic, abruptly reducing child-to-child transmission for a fixed length of time and also influencing the transmission between children and adults. Public health effectiveness was measured by reduction in the cumulative incidence, and cost effectiveness was also examined by calculating the incremental cost effectiveness ratio and adopting a threshold of 1.0 × 107 Japanese Yen/life-year.
School closure at the epidemic peak appeared to yield the largest reduction in the final size, while the time of epidemic peak was shown to depend on the transmissibility. As the length of school closure was extended, we observed larger reduction in the cumulative incidence. Nevertheless, the cost effectiveness analysis showed that the cost of our school closure scenario with the parameters derived from H1N1-2009 was not justifiable. If the risk of death is three times or greater than that of H1N1-2009, the school closure could be regarded as cost effective.
There is no fixed timing and duration of school closure that can be recommended as universal guideline for different types of influenza viruses. The effectiveness of school closure depends on the transmission dynamics of a particular influenza virus strain, especially the virulence (i.e. the infection fatality risk).
- Cumulative Incidence
- Incremental Cost Effectiveness Ratio
- Influenza Pandemic
- Transmission Dynamic
School closure is one of important non-pharmaceutical countermeasures against influenza pandemic. Among various types of school closure, the so-called "proactive closure", i.e., the closure of schools before observing substantial transmissions among school children, was conducted in Japan during the early stage of H1N1-2009 pandemic, and micro-clade of the viruses that caused the earliest clusters is known to have declined to extinction. Japanese experience demonstrated that the proactive closure as part of concerted effort of containment measure can be very helpful in achieving the local extinction.
However, there has been a variation in published opinions toward another type of closure, "the reactive closure", i.e., the closure of schools when many children, staff or both are experiencing illness, as part of mitigation strategy. Published studies have empirically explored the impact of reactive school closure on an influenza epidemic or pandemic, and some of the studies demonstrated substantial reduction in the rate of transmission among school children during the closure[4–7]. However, others emphasized that the community impact, e.g. reduction in the demand of healthcare service including hospitalization of severe cases, is likely very limited[8–10]. Elucidating the details of school transmission mechanism has been ongoing (see Discussion), and there has been no simple policy (e.g. the timing and duration) to implement the closure in the reactive manner during the course of a pandemic.
If we have a clear quantitative guideline for the reactive closure (e.g. provision of public health conditions at which the closure can be justified and decided), that could greatly benefit public health policymakers. Fundamental insights into the effectiveness of school closure can be gained from a parsimonious mathematical model, exploring possible answers to such key policy questions using simplistic modelling approaches. In the present study, our questions are two-folds. First, we examine when one should close the school during the course of a pandemic. Second, we explore how long the closure should be implemented. Parameterizing the model using published epidemiological estimates of the pandemic H1N1-2009 and accounting for the cost effectiveness of closure, we discuss if the timing and length of school closure could be optimized.
where Na represents the population size of age-group a.
where αi is a relative susceptibility of age group i, Mij represents the number of contacts that an individual in age group i experiences with individuals in age group j per unit time, and Γi is the integral of the survival function of infectious period of age group i (Γi = 1/γi if the infectious period is exponentially distributed with the mean 1/γi days). Hereafter, the matrix M = [Mij] is referred to as the contact matrix.
where ϕ(t) represents the relative increase in the reproduction number between children and young adults due to compensatory behaviour of children with young adults during the closure. For clarity of modelling and due to shortage of scientific evidence, we ignore the influence of compensation on other elements of the contact matrix (e.g. we assume that the contacts within young adults and between children and elderly are not influenced by closure). Further mathematical details of the compensatory contact are described in Appendix.
wa is given as the product of the infection fatality risk (IFR; i.e., the risk of death given infection with influenza virus) and life-expectancy of age group a, assumed to be 65, 45 and 15 years for children, young adults and elderly, respectively. The age-specific estimates of IFR are extracted from empirical study in Hong Kong, assumed as 1, 10 and 500 deaths per 100,000 infections for children, young adults and elderly, respectively. To measure the effectiveness of school closure, we compare the absolute difference of L between two scenarios, i.e. with and without school closure, L1 and L0 (i.e., L0-L1), yielding the life-years saved by the school closure.
where u1 and u0 are the proportions of households in which either father or mother has to take annual leave, with or without an infected child in the household, respectively (i.e., u1 = 0.267 and u0 = 0.104) and N2 is the population size of adults. For simplicity, we ignore small fractions of childless couples and unmarried adults among the total young adults. The cost of u1 is multiplied to the cumulative incidence of children during the closure, ignoring multiple infections in a household (i.e. we ignore brother(s) or sister(s) who are infected at the same time in a household). The population size to multiply u0 is calculated by subtracting the cumulative incidence from N2/2, because it is usually the case that either father or mother is absent from work during closure (and one of them continues to work). According to a white paper of an economic study, c is given as the product of hourly wage, the average working hours (per day) and the cost to be paid for workers during school closure, yielding c = 10,019.52 Japanese Yen (approximately 100 US Dollars) per day. Since we consider the short (and realistic) lengths of closure (e.g. 7, 14 or 21 days), we account for only individual impact of closure and ignore other indirect social cost such as the loss of business opportunities due to extended period of closure or stagnation of overall economic activity.
It should be noted that the quantity Y is theoretically equivalent to the incremental cost effectiveness ratio (ICER) in the cost effectiveness analysis (CEA) studies. We aim to identify reasonable combination of τ and t0 that minimizes our objective function Y. In the United States and United Kingdom, the acceptable threshold of ICER tends to lie around 100,000 US dollars and 30,000 British Pound per life year, respectively. Accordingly, we assume that the corresponding threshold lies in the range from 5.0 × 106 to 1.0 × 107 Japanese Yen (and draw a line for the latter threshold in all associated figures).
We consider an epidemic in a population of 1 million. The size of each age group is assumed as proportional to age-specific population sizes of those aged from 0–19, 20–59 and 60 years and over for the entire Japan, i.e., (N1, N2, N3) = (177207, 500957, 321836). The parameters of contact matrix are derived from previous study with an identical age categorization which essentially assumed that the age-specific contact pattern in Japan is not different from that in England. In our scenario, an epidemic takes place with an introduction of single infected child j1(0) = 1 (while j2(0) = j3(0) = 0), and all residents are assumed as initially susceptible, i.e., (s1(0), s2(0), s3(0)) = (1,1,1). Nevertheless, their relative susceptibility per contact (e.g. the probability of successful transmission per contact) depends on age, i.e., (α1, α2, α3) = (1.000,0.370,0.059) as empirical evidence suggests[11, 19]. Mean generation times of secondary transmission caused by children and young adults or elderly are assumed to be 2.2 days and 2.7 days, respectively, and we assume that the generation time follows an exponential distribution for mathematical convenience (so that the model (1) can also be written as ordinary differential equations using the next-generation matrix parameterized by (6)). The basic reproduction number, R0, is used for scaling the next-generation matrix and is assumed to be 1.4[21, 22].
With respect to the protective effect of intervention, the relative reduction in the reproduction number during school closure (hereafter referred to as the "efficacy" of school closure) has been empirically estimated in limited number of settings[4, 5, 9, 23]. As our baseline, we assume that there is a 70% decline in child-to-child transmission during the closure and thus q = 0.3. In a Japanese survey, no apparent increase in the frequency of child-to-adult contact was observed during the closure, and we set the proportion of child contacts compensated (π; see Appendix) as 0 at the baseline and then vary it from 0 to 0.5 (where 0.5 means that 50% of intervened within-child contacts are alternatively made with young adults).
The above-mentioned parameters correspond to empirically measured results from pandemic H1N1-2009 which is known to have been very mild. Thus, we also measure the sensitivity of ICER to different levels of transmission potential and risks of death, varying R0 from 1.2 to 1.8 and elevating the relative risk of death from 1 to 100 (using H1N1-2009 as the reference).
Epidemic dynamics and school closure
Cost effectiveness analysis of school closure
The present study examined the public health effectiveness and cost effectiveness of school closure which was assumed to be implemented only once during the course of a pandemic. The model was parameterized with reference to empirical estimates of the pandemic H1N1-2009. School closure at the epidemic peak appeared to minimize the cumulative incidence, but the time of epidemic peak was shown to depend on R0. As the duration of school closure was extended, we observed a larger reduction in the cumulative incidence. Strikingly, the cost effectiveness analysis showed that our school closure scenario with parameters derived from the pandemic H1N1-2009 was not cost effective. Nevertheless, if the virulence is three times or greater than that of H1N1-2009, the cost of closure could be justified. These findings were not very sensitive to the efficacy of school closure and compensatory contact behaviour among children.
There are three important learning points from the present study. First, we have shown that there is no fixed timing and duration of school closure that can be recommended as universal for different types of influenza viruses. It is natural that the effectiveness of school closure depends on the underlying transmission dynamics, and the absence of simple universal guideline should be explained and communicated to non-experts. In fact, the dependence of the effectiveness of school closure on the transmission dynamics could explain variations in published opinions toward the epidemiological performance of the school closure as an option of mitigation strategy. Second, school closure during the pandemic H1N1-2009 appeared not to be cost effective even when the ICER took the minimum value. It implies that the cost-effective intervention against mild pandemic strain such as H1N1-2009 is different from that of virulent strains. This echoes the finding by Halder et al.[26, 27] based on a simulation approach. Third, if a particular strain is virulent, school closure could be cost effective. Namely, given that the virulent strain widely spreads in the community, yielding high disease burden, school closure intervention, the cost of which is regarded as small for the high disease burden, should be implemented to reduce the disaster size as much as possible. In addition, since the equation (12) involves the wage of parents, it is important to remember that the cost-effectiveness may also depend on an economic standard of a country (e.g. school closure of pandemic H1N1-2009 may even be justified in a country with much smaller salary than that of Japan).
As a policy implication at minimum, one should remember that optimal school closure depends on the severity of pandemic, characterized by the transmissibility and virulence. Especially, the school closure is likely cost effective for virulent influenza strains. One should also know that the cost-effective interventions of a particular influenza strain are different from those for other strains. In addition, rather than industrialized countries, it may be easier to justify the cost of school closure in developing countries where the parental impact is likely smaller.
Three technical limitations should be noted. First, our model assumes that the transmission from child to child is homogeneous. More rigorous network model has shown that such random mixing assumption could overestimate the effectiveness of school closure. Second, the impact of school closure on social mixing patterns should ideally be based on more realistically socially structured (layered) modelling approach, perhaps classifying transmissions into those occurring in households, schools and community. Third, more precise features of child contact, including weekend contact and the impact of illness on the contact, are recently shown to have a substantial impact on the effectiveness of intervention[31, 32].
Unfortunately, school closure during H1N1-2009 may not be fully justified when it comes to the cost. Nevertheless, the cost effectiveness should be regarded as merely a single aspect of the impact of this intervention to help policymaking. Perhaps, rather than focusing more on the cost, one should carefully reconsider public health objectives of this intervention, e.g. delaying epidemic peak, reducing the height of peak prevalence, or reducing the overall epidemic size, and decide what we would expect from this intervention more in detail. Expecting the effectiveness in all these aspects may not be feasible. Despite the presence of numerous tasks to guide school closure in the next pandemics, our study has at least shown that one can examine the potential performance of school closure using the proposed simplistic modelling approach.
The present study examined the public health effectiveness and cost effectiveness of school closure. The effectiveness of school closure depends on the transmission dynamics of a particular influenza virus strain, especially the virulence. School closure in our scenario with parameters derived from the pandemic H1N1-2009 appeared not to be cost effective. There is no fixed timing and duration of school closure that can be recommended as universal guideline for different types of influenza viruses.
Sum of elements in a single row i or single column i represents the total number of contacts (per day) made by all of those in age-group i. We assume that the sum of each column is decomposed as the product of average contact per person ki and the population size Ni, so that the sum can be rewritten as kiNi. In the presence of compensatory behaviour, we assume that the total number of contact, kiNi, made by age-group i is partially maintained even during the school closure.
It should be noted that the increase in child contact with elderly is ignored, assuming that mostly parents, not elderly, have to take care of children during the school closure. It should also be noted that other contacts were assumed not to have been influenced by the compensation, although, in theory, the maintenance of contacts should influence all other elements (so that the total number of contacts per day remain constant for all age-groups); we ignore this mathematical issue for simplicity.
This study originated from Study Group on Infectious Disease Modelling, supported by the Cooperation with Mathematics Program (2013S01) of the Ministry of Education, Culture, Sports, Science and Technology and The Institute of Statistical Mathematics, Japan. HN wishes to thank Takahashi Industrial and Economic Research Foundation for funding support. KE received scholarship support from the Japan Society for Promotion of Science (JSPS). KE and HI received funding support from the Aihara Project, the FIRST program from JSPS, initiated by CSTP. SN received funding support from Grant-in-Aid for Young Scientists B25871132.
- Cauchemez S, Ferguson NM, Wachtel C, Tegnell A, Saour G, Duncan B, Nicoll A: Closure of schools during an influenza pandemic. Lancet Infect Dis. 2009, 9: 473-481. 10.1016/S1473-3099(09)70176-8.View ArticlePubMedGoogle Scholar
- Nishiura H, Castillo-Chavez C, Safan M, Chowell G: Transmission potential of the new influenza A(H1N1) virus and its age-specificity in Japan. Euro Surveill. 2009, 14: 19227-PubMedGoogle Scholar
- Shiino T, Okabe N, Yasui Y, Sunagawa T, Ujike M, Obuchi M, Kishida N, Xu H, Takashita E, Anraku A, Ito R, Doi T, Ejima M, Sugawara H, Horikawa H, Yamazaki S, Kato Y, Oguchi A, Fujita N, Odagiri T, Tashiro M, Watanabe H: Molecular evolutionary analysis of the influenza A(H1N1)pdm, May-September, 2009: temporal and spatial spreading profile of the viruses in Japan. PLoS One. 2010, 5: e11057-10.1371/journal.pone.0011057.PubMed CentralView ArticlePubMedGoogle Scholar
- Cauchemez S, Valleron AJ, Boëlle PY, Flahault A, Ferguson NM: Estimating the impact of school closure on influenza transmission from Sentinel data. Nature. 2008, 452: 750-754. 10.1038/nature06732.View ArticlePubMedGoogle Scholar
- Wu JT, Cowling BJ, Lau EH, Ip DK, Ho LM, Tsang T, Chuang SK, Leung PY, Lo SV, Liu SH, Riley S: School closure and mitigation of pandemic (H1N1) 2009, Hong Kong. Emerg Infect Dis. 2010, 16: 538-541. 10.3201/eid1603.091216.PubMed CentralView ArticlePubMedGoogle Scholar
- Earn DJ, He D, Loeb MB, Fonseca K, Lee BE, Dushoff J: Effects of school closure on incidence of pandemic influenza in Alberta, Canada. Ann Intern Med. 2012, 156: 173-181. 10.7326/0003-4819-156-3-201202070-00005.View ArticlePubMedGoogle Scholar
- Hens N, Ayele GM, Goeyvaerts N, Aerts M, Mossong J, Edmunds JW, Beutels P: Estimating the impact of school closure on social mixing behaviour and the transmission of close contact infections in eight European countries. BMC Infect Dis. 2009, 9: 187-10.1186/1471-2334-9-187.PubMed CentralView ArticlePubMedGoogle Scholar
- Cowling BJ, Lau EH, Lam CL, Cheng CK, Kovar J, Chan KH, Peiris JS, Leung GM: Effects of school closures, 2008 winter influenza season, Hong Kong. Emerg Infect Dis. 2008, 14: 1660-1662. 10.3201/eid1410.080646.PubMed CentralView ArticlePubMedGoogle Scholar
- Vynnycky E, Edmunds WJ: Analyses of the 1957 (Asian) influenza pandemic in the United Kingdom and the impact of school closures. Epidemiol Infect. 2008, 136: 166-179.PubMed CentralPubMedGoogle Scholar
- House T, Baguelin M, Van Hoek AJ, White PJ, Sadique Z, Eames K, Read JM, Hens N, Melegaro A, Edmunds WJ, Keeling MJ: Modelling the impact of local reactive school closures on critical care provision during an influenza pandemic. Proc R Soc London, Ser B. 2011, 278: 2753-2760. 10.1098/rspb.2010.2688.View ArticleGoogle Scholar
- Nishiura H, Chowell G, Safan M, Castillo-Chavez C: Pros and cons of estimating the reproduction number from early epidemic growth rate of influenza A (H1N1) 2009. Theor Biol Med Model. 2010, 7: 1-10.1186/1742-4682-7-1.PubMed CentralView ArticlePubMedGoogle Scholar
- Ministries associated with Pandemic and Avian Influenza: Guideline of interventions against pandemic and other novel influenza. 2013, Tokyo: Japanese Government, in JapaneseGoogle Scholar
- Charu V, Chowell G, Palacio Mejia LS, Echevarría-Zuno S, Borja-Aburto VH, Simonsen L, Miller MA, Viboud C: Mortality burden of the A/H1N1 pandemic in Mexico: a comparison of deaths and years of life lost to seasonal influenza. Clin Infect Dis. 2011, 53: 985-993. 10.1093/cid/cir644.PubMed CentralView ArticlePubMedGoogle Scholar
- Wong JY, Wu P, Nishiura H, Goldstein E, Lau EH, Yang L, Chuang SK, Tsang T, Peiris JS, Wu JT, Cowling BJ: Infection fatality risk of the pandemic A(H1N1)2009 virus in Hong Kong. Am J Epidemiol. 2013, 177: 834-840. 10.1093/aje/kws314.PubMed CentralView ArticlePubMedGoogle Scholar
- Mizumoto K, Yamamoto T, Nishiura H: Contact behaviour of children and parental employment behaviour during school closures against the pandemic influenza A (H1N1-2009) in Japan. J Int Med Res. 2013, 41: 716-724. 10.1177/0300060513478061.View ArticlePubMedGoogle Scholar
- ITR Corporation: ITR White Paper 2009. Report to consider renovation in working style and Windows 7 following pandemic intervention. 2009, Tokyo: ITR Corporation, in JapaneseGoogle Scholar
- Ejima K, Omori R, Aihara K, Nishiura H: Real-time investigation of measles epidemics with estimate of vaccine efficacy. Int J Biol Sci. 2012, 8: 620-629.PubMed CentralView ArticlePubMedGoogle Scholar
- Mossong J, Hens N, Jit M, Beutels P, Auranen K, Mikolajczyk R, Massari M, Salmaso S, Tomba GS, Wallinga J, Heijne J, Sadkowska-Todys M, Rosinska M, Edmunds WJ: Social contacts and mixing patterns relevant to the spread of infectious diseases. PLoS Med. 2008, 5: e74-10.1371/journal.pmed.0050074.PubMed CentralView ArticlePubMedGoogle Scholar
- Dorigatti I, Cauchemez S, Pugliese A, Ferguson NM: A new approach to characterising infectious disease transmission dynamics from sentinel surveillance: application to the Italian 2009–2010 A/H1N1 influenza pandemic. Epidemics. 2012, 4: 9-21. 10.1016/j.epidem.2011.11.001.PubMed CentralView ArticlePubMedGoogle Scholar
- te Beest DE, Wallinga J, Donker T, van Boven M: Estimating the generation interval of influenza A (H1N1) in a range of social settings. Epidemiology. 2013, 24: 244-250. 10.1097/EDE.0b013e31827f50e8.View ArticlePubMedGoogle Scholar
- Fraser C, Donnelly CA, Cauchemez S, Hanage WP, Van Kerkhove MD, Hollingsworth TD, Griffin J, Baggaley RF, Jenkins HE, Lyons EJ, Jombart T, Hinsley WR, Grassly NC, Balloux F, Ghani AC, Ferguson NM, Rambaut A, Pybus OG, Lopez-Gatell H, Alpuche-Aranda CM, Chapela IB, Zavala EP, Guevara DM, Checchi F, Garcia E, Hugonnet S, Roth C, WHO Rapid Pandemic Assessment Collaboration: Pandemic potential of a strain of influenza A (H1N1): early findings. Science. 2009, 324: 1557-1561. 10.1126/science.1176062.PubMed CentralView ArticlePubMedGoogle Scholar
- Roberts MG, Nishiura H: Early estimation of the reproduction number in the presence of imported cases: pandemic influenza H1N1-2009 in New Zealand. PLoS One. 2011, 6: e17835-10.1371/journal.pone.0017835.PubMed CentralView ArticlePubMedGoogle Scholar
- Omori R, Nishiura H: Theoretical basis to measure the impact of short-lasting control of an infectious disease on the epidemic peak. Theor Biol Med Model. 2011, 8: 2-10.1186/1742-4682-8-2.PubMed CentralView ArticlePubMedGoogle Scholar
- Nishiura H: Case fatality ratio of pandemic influenza. Lancet Infect Dis. 2010, 10: 443-444. 10.1016/S1473-3099(10)70120-1.View ArticlePubMedGoogle Scholar
- Mizumoto K, Yamamoto T, Nishiura H: Age-dependent estimates of the epidemiological impact of pandemic influenza (H1N1-2009) in Japan. Comput Math Methods Med. 2013, 2013: 637064-PubMed CentralView ArticlePubMedGoogle Scholar
- Halder N, Kelso JK, Milne GJ: Developing guidelines for school closure interventions to be used during a future influenza pandemic. BMC Infect Dis. 2010, 10: 221-10.1186/1471-2334-10-221.PubMed CentralView ArticlePubMedGoogle Scholar
- Halder N, Kelso JK, Milne GJ: Cost-effective strategies for mitigating a future influenza pandemic with H1N1 2009 characteristics. PLoS One. 2011, 6: e22087-10.1371/journal.pone.0022087.PubMed CentralView ArticlePubMedGoogle Scholar
- Potter GE, Handcock MS, Longini IM, Halloran ME: Estimating within-school contact networks to understand influenza transmission. Ann Appl Stat. 2012, 6: 1-26. 10.1214/11-AOAS505.PubMed CentralView ArticlePubMedGoogle Scholar
- Cauchemez S, Bhattarai A, Marchbanks TL, Fagan RP, Ostroff S, Ferguson NM, Swerdlow D: Pennsylvania H1N1 working group: Role of social networks in shaping disease transmission during a community outbreak of 2009 H1N1 pandemic influenza. Proc Natl Acad Sci U S A. 2011, 108: 2825-2830. 10.1073/pnas.1008895108.PubMed CentralView ArticlePubMedGoogle Scholar
- Ejima K, Aihara K, Nishiura H: Probabilistic differential diagnosis of Middle East respiratory syndrome (MERS) using the time from immigration to illness onset among imported cases. J Theor Biol. 2014, 346: 47-53. (doi:10.1016/j.jtbi.2013.12.024)View ArticlePubMedGoogle Scholar
- Eames KT, Tilston NL, Edmunds WJ: The impact of school holidays on the social mixing patterns of school children. Epidemics. 2011, 3: 103-108. 10.1016/j.epidem.2011.03.003.View ArticlePubMedGoogle Scholar
- 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: 267-312.View ArticlePubMedGoogle Scholar
- Hollingsworth TD, Klinkenberg D, Heesterbeek H, Anderson RM: Mitigation strategies for pandemic influenza A: balancing conflicting policy objectives. PLoS Comp Biol. 2011, 7: e1001076-10.1371/journal.pcbi.1001076.View ArticleGoogle Scholar
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 cited. 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.