- Open Access
Theoretical basis to measure the impact of short-lasting control of an infectious disease on the epidemic peak
Theoretical Biology and Medical Modellingvolume 8, Article number: 2 (2011)
While many pandemic preparedness plans have promoted disease control effort to lower and delay an epidemic peak, analytical methods for determining the required control effort and making statistical inferences have yet to be sought. As a first step to address this issue, we present a theoretical basis on which to assess the impact of an early intervention on the epidemic peak, employing a simple epidemic model.
We focus on estimating the impact of an early control effort (e.g. unsuccessful containment), assuming that the transmission rate abruptly increases when control is discontinued. We provide analytical expressions for magnitude and time of the epidemic peak, employing approximate logistic and logarithmic-form solutions for the latter. Empirical influenza data (H1N1-2009) in Japan are analyzed to estimate the effect of the summer holiday period in lowering and delaying the peak in 2009.
Our model estimates that the epidemic peak of the 2009 pandemic was delayed for 21 days due to summer holiday. Decline in peak appears to be a nonlinear function of control-associated reduction in the reproduction number. Peak delay is shown to critically depend on the fraction of initially immune individuals.
The proposed modeling approaches offer methodological avenues to assess empirical data and to objectively estimate required control effort to lower and delay an epidemic peak. Analytical findings support a critical need to conduct population-wide serological survey as a prior requirement for estimating the time of peak.
The influenza A (H1N1-2009) pandemic began in early 2009, and rapidly spread worldwide. Mathematical epidemiologists characterized the epidemic and provided key insights into its dynamics from the earliest stages of the pandemic . The transmission potential was quantified shortly after the declaration of emergence [2–6], while statistical estimation and relevant discussion of epidemiological determinants were underway before substantial numbers of cases were reported in many countries .
Prior to the pandemic, many countries issued the original pandemic preparedness plans and guidelines, aiming to instruct the public and to advocate community mitigation. The goals of the mitigation have been threefold; (a) to delay epidemic peak, (b) to reduce peak burden on hospitals and infrastructure (by lowering the height of peak) and (c) to diminish overall morbidity impacts . To assess these aspects under different intervention scenarios, various modeling studies have been conducted (e.g. [8–10]), most notably, by simulating the detailed influenza transmission dynamics.
Although simulations have aided our understanding of expected dynamics in realistic situations and in different scenarios, analytical methods that objectively determine the required control effort and that make statistical inference (e.g. evaluation of empirically observed delay) have yet to be developed. Focus on epidemic peak (relating to mitigation goals (a) and (b) above) has been particularly understudied. Goal (c), on the other hand, is readily formulated in terms of the so-called final epidemic size. The time delay of a major epidemic (such as that resulting from international border control) has been explored using simplistic modeling approaches [11, 12]; however, the height and time of an epidemic peak involve nonlinear dynamics, rendering analytical approaches difficult. Despite the mathematical complexity, goals (a) and (b) can be more readily understood from empirical data during early epidemic phase than can goal (c), because an explicit understanding of goal (c) in the presence of interventions requires knowledge of the full epidemiological dynamics over the entire epidemic period.
In the present study, we present a theoretical basis from which the impact of an early intervention on the height and time of epidemic peak may be assessed. As a special case, we consider a scenario in which intervention is implemented only briefly during the early epidemic phase (e.g. unsuccessful containment). We employ a parsimonious epidemic model with homogeneously mixing assumption, because nonlinear epidemic dynamics involve a number of analytical complexities. As a first step towards understanding epidemiological factors that influence the epidemic peak, leading to the eventual statistical inference of relevant effects, we seek fundamental analytical strategies to evaluate the impact of short-lasting control on epidemic peak using the simplest epidemic model . For our model to become fully applicable and to more closely match empirical data, a number of extensions are required. We discuss ways by which these extensions can be practically realized.
We first present our study motivation. During the early epidemic phase of the 2009 pandemic, many countries initially enforced strict countermeasures to locally contain the epidemic. Early intervention includes, but is not limited to, quarantine, isolation, contact tracing and school closure. Nevertheless, once it was realized that a major epidemic was unavoidable, regions and countries across the world were compelled to downgrade control policy from containment to mitigation. Although mitigation also involves various countermeasures (and indeed, mitigation originally intends to achieve the above mentioned goals (a)-(c)), one desires to know the effectiveness of the unsuccessful containment effort. Among its many outcomes, the present study focuses on the height and time of epidemic peak.
The applicability of our theoretical arguments is not restricted to the switch of control policy. In many Northern hemisphere countries, the start of the major epidemic of H1N1-2009 (which may or may not have been preceded by early stochastic phase) corresponds to the summer school holiday period. Adults also take vacation over a part of this period. In addition to strategic school closure as an early countermeasure against influenza [14, 15], school holiday is known to suppress the transmission of influenza , mainly because transmission tends to be maintained by school-age children [2, 17–19]. Following this trend, a decline in instantaneous reproduction number has been empirically observed during the summer holiday period of the 2009 pandemic . Transmission resumes once a new semester starts. The effectiveness of the summer holiday period in lowering and delaying the epidemic peak is, therefore, a matter of great interest.
Both questions are addressed by considering time-dependent increase in the transmission rate. Let β be the transmission rate per unit time in the absence of an intervention of interest (or during the mitigation phase in the case of our first question). Due to intervention (or school holiday) in the early epidemic phase, β is initially reduced by a factor α (0 ≤ α ≤ 1) until time t1 (Figure 1A). Though transmission rate abruptly increases at time t1 when the control policy is eased or when the new school semester starts, we observe a reduced height of, and a time delay in, the epidemic peak compared to the hypothetical situation in which no intervention takes place (Figure 1B). More realistic situations may be envisaged (e.g. a more complex step function or seasonality of transmission), but we restrict ourselves to the simplest scenario in the present study.
Here we consider the simplest form of Kermack and McKendrick epidemic model , formulated in terms of ordinary differential equations. The following assumptions are made: (i) the population is homogeneously mixing, (ii) the epidemic occurs in a population in which the majority of individuals are susceptible, (iii) the time scale of the epidemic is sufficiently shorter than the average life expectancy at birth of the host, and we ignore the background demographic dynamics, (iv) the epidemic occurs in a closed constant population without immigration and emigration again justified based on time scale, and (v) once an infected individual recovers, he/she becomes completely and permanently immune against further infections. Let the numbers of susceptible, infectious and recovered individuals at calendar time t be S(t), I(t) and U(t), respectively. We use the notation U(t) for recovered individuals to avoid confusion with the instantaneous reproduction number at calendar time t, R(t). The population size N remains constant over time (N = S(t) + I(t) + U(t)). The so-called SIR (susceptible-infected-recovered) model is written as
where R(t) is the instantaneous reproduction number (i.e., the average number of secondary cases generated by a single primary case at calendar time t) and γ is the rate of recovery. Given time-dependent transmission rate β(t) and susceptible population size S(t) at time t, R(t) is assumed to be given by
Although β(t) will be dealt with as a simple step function in the following analysis, we use the general notation to motivate future analysis of more complex time-dependent dynamics. We assume that an epidemic starts at time 0 with an initial condition (S(0), I(0), U(0)) = (S0, I0,U0) where I0 = 1 and U0/N ≈ 0, i.e. an epidemic occurs in a population in which the majority of individuals are susceptible at t = 0. Under this initial condition, we consider two different scenarios for R(t). First, a hypothetical scenario in which no intervention takes place, i.e.
which is hereafter referred to as the baseline scenario. Second, we consider an observed scenario in which an intervention takes place during the early stage of the epidemic. Let t1 and tm,0be calendar times at which the intervention terminates, and at which a peak prevalence of infectious individuals is observed in the absence of intervention, respectively. As mentioned above, we assume that the intervention reduces the reproduction number by a factor α (0 ≤ α ≤ 1) for 0 ≤ t < t1. For t ≥ t1, we assume that the transmission rate is recovered to β as in (3).
We assume t1 < tm,0, i.e. we consider a scenario in which transmission rate recovers before the time at which peak prevalence is observed in baseline scenario. We further assume that R(t) > 1 for t < t1. That is, the efficacy α of an intervention effort (or summer holiday) is by itself not sufficient to contain the epidemic.
To illustrate our modeling approaches, we consider the transmission dynamics of pandemic influenza (H1N1-2009), ignoring the detailed epidemiological characteristics (e.g. pre-existing immunity, realistic distribution of generation time and the presence of asymptomatic infection). The initial reproduction number in the absence of interventions R(0) is assumed to be 1.4 . Given that expected values of empirically estimated serial interval ranged from 1.9 to 3.6 days [2, 5, 21–23], the mean generation time 1/γ is assumed to be 3 days [24, 25].
Our study questions are twofold. First, we aim to quantify the decline in peak prevalence (I(t)/N) due to a short-lasting intervention. The peak prevalence of the intervention scenario is always smaller than that of baseline scenario (see below), and we show that this difference can be analytically expressed. Second, we are interested in the time delay in observing peak prevalence in the presence of intervention. We develop an approximate strategy to quantify the difference in times of peak between baseline and intervention scenarios.
Difference in peak prevalence
We move on to consider estimates of peak prevalence in two scenarios. For mathematical convenience, we use the prevalence of infectious individuals (I(t)/N) to consider the epidemic peak. The peak prevalence of infectious individuals is preceded by peak incidence (γR(t)I(t)/N) by approximately the mean infectious period of 1/γ days. As was realized elsewhere , analysis of prevalence is easier than that of incidence. Beginning with two sub-equations of system (1), we have
Note that R(t) is a function of S(t). Integrating (5) in baseline scenario, we obtain 
A theoretical condition for the observation of peak prevalence at time tm,0is dI(tm,0)/dt = 0, or equivalently, R(tm,0) = 1. As evident from equation (2), this condition satisfies S(tm,0) = γ/β. The peak prevalence I(tm,0)/N is then given by 
Note that S0/R(0)N represents the proportion yet to be infected and S0 ln R(0)/R(0)N is the proportion removed at time tm,0. Equation (7) indicates that the peak prevalence of SIR model is determined by the initial condition and the transmission potential R(0). It should be noted that S0/R(0) can be replaced by γ/β, and thus, I(tm,0) is independent of initial condition for U0 = 0 (a special case).
In the intervention scenario, equation (6) with replacement of β by αβ applies for t < t1.
which provides another initial condition at time t = t1 for t ≥ t1. That is, we can also employ (6) to compute peak prevalence for t ≥ t1 with initial condition (S(t1), I(t1),U(t1)). Again, a condition to observe peak prevalence at time tm,1is R(tm,1) = 1, which gives S(tm,1) = γ/β. The peak prevalence I(tm,1)/N of the intervention scenario is given by
Note that R(0) in the above equation refers to βS0/γ (i.e. we use R(0) in our baseline scenario to permit an explicit comparison between the two scenarios). Inserting right-hand side of (8) into (9), we obtain
Consequently, relative reduction in peak prevalence due to intervention within time t1 is ∈ α = (I(tm,0) - I(tm,1))/N, which can be parameterized as
Equation (11) indicates that the difference of peak prevalence between the two scenarios is determined by four different factors; the relative reduction in reproduction number α due to the intervention, initial condition at time 0, transmission potential R(0), and fraction of susceptible individuals at time t1 under the intervention. If the initial condition, the transmission dynamics in the absence of interventions (i.e. R(0), β and γ) and t1 are known, an estimate of α gives S(t1), yielding an estimate of ∈ α .
Delay in epidemic peak
The time to observe peak prevalence is analytically more challenging than the height of peak prevalence, because even an approximate estimate requires an analytical solution to the model (1). We propose a parsimonious approximation strategy which leads to more convenient solutions than those discussed in past studies (e.g. ). Substituting I(t) in the first sub-equation of (1) by (1/γ)(dU(t)/dt), we have
For the baseline scenario (i.e. β(t) = β), integrating (12) from time 0 to t,
Because U(0)/N ≈ 0,
Subsequently, the third sub-equation of (1) is rewritten as
Here we impose another key approximation. Because the quantity βU(t)/γ (≈ R(0)U(t)/N) for influenza (e.g. R(0) = 1.4) tends to be smaller than 1 (especially, before observing epidemic peak), we use a Taylor series expansion, i.e.,
If the quadratic approximation is inadequate for large R(0), a higher order Taylor polynomial function can be used. Inserting the quadratic approximation into (15), and imposing a further approximation (i.e. S0≈ N), we obtain
which appears to be a logistic equation. We use this logistic-form solution instead of the more commonly employed hyperbolic-form solution [29, 30], to illustrate a simpler approximate solution and to demonstrate the problem underlying both solutions. Later, we use a more formal solution (of logarithmic-form) in the intervention scenario, which is numerically identical to the classical hyperbolic-form solution (see below). Assuming that U(0) = U0> 0, the analytical solution of (17) is
The derivative of (18) is dU(t)/dt = γI(t). It follows that
Further differentiation of (19) gives dI(t)/dt, and letting dI(t)/dt = 0, the time to observe peak prevalence is analytically derived. For the logistic equation, the corresponding time has been referred to as the inflection point of the cumulative curve in equation (18) . The inflection point tm,0to observe peak prevalence is
which depends on initial condition and transmission characteristics. In the intervention scenario (in which intervention is short-lasting), an identical approach can be taken for t < t1, replacing β by αβ (or by replacing R(0) by αR(0)). Subsequently, the epidemic peak occurs at tm,1(> t1). We take a similar approach to that used in (15) with a computed initial condition (S(t1), I(t1), U(t1)) using (18) and (19). For t ≥ t1,
Now we apply an approximation
It should be noted that, in the above approximation, we include the term exp(βU(t1)/γ), because U(t) - U(t1) better satisfies the Taylor series approximation than expanding U(t) alone. Let constants A, B and C be
Given these constants, we consider
where z = t - t1 and V (z) = U(z + t1) for t ≥ t1. The initial condition V (0) is V0 = U(t1). Writing (24) in integral form, we have 
Past studies have typically assumed hyperbolic-form functions for the analytical solution of (25) [29, 30]. However, we express the solution in logarithmic-form [31, 32], because logarithmic functions are compatible with spreadsheet programs. The logarithmic-form solution reads
Differentiating (26) with respect to z and taking dV (z)/dz = 0, we find the inflection point zm,1to be
Replacing z by t, we obtain
as an approximate solution of the epidemic peak in intervention scenario. The time delay of this peak, imposed by the intervention in the early epidemic phase, τ α is subsequently calculated as
using (20) and (30) for the right-hand side. The delay depends on initial condition U0, the length of intervention t1 (both of which are apparent from (20) and (30)) and on the efficacy of intervention α (since this quantity influences the initial condition U(t1) in (30)).
Application and illustration
Empirical analysis of influenza A (H1N1-2009)
Here, we apply the above described theory to empirical influenza A (H1N1-2009) data. Figure 2 shows the estimated number of influenza cases based on national sentinel surveillance in Japan from week 31 (week ending 2 August) 2009 to week 13 (week ending 28 March) 2010. The estimates follow an extrapolation of the notified number of cases from a total of 4800 randomly sampled sentinel hospitals to the actual total number of medical facilities in Japan. The cases represent patients who sought medical attendance and who have met the following criteria, (a) acute course of illness (sudden onset), (b) fever greater than 38.0°C, (c) cough, sputum or breathlessness (symptoms of upper respiratory tract infection) and (d) general fatigue, or who were strongly suspected of the disease undertaking laboratory diagnosis (e.g. rapid diagnostic testing). Although the estimates of sentinel surveillance data involve various epidemiological biases and errors, we ignore these issues in the present study. Prior to week 31, the number of cases was small and the dynamics in the early stochastic phase have been examined elsewhere . We arbitrarily assume that the major epidemic starts in week 31.
It is interesting to observe that the period A in Figure 2 corresponds to that of summer school holiday. Due to reporting delay of approximately 1 week , we assume that weeks 31 to 36 inclusive (the latter of which ends on 6 September) reflect the transmission dynamics during the summer school holiday. Subsequently, school opens in September with an epidemic peak in late November (period B), followed by abrupt decline during the winter holiday (period C) and start of winter semester (period D). Among these periods, we focus on the impact of summer holiday (period A), relative to period B, in lowering epidemic peak and delaying the time to observe the peak. More specifically, we estimate the reproduction number R(0) and its reduction α from the data set encompassing weeks 31 to 42. To permit an explicit estimation, we assume that linear approximation holds, as was similarly assumed elsewhere . We assume that the reproduction number is reduced by a factor α from week 31 to 36 due to summer holiday, while the reproduction number recovers to R(0) from week 37 to 42.
Let r0 be the exponential growth rate of cases per day in the absence of summer holiday. Because our SIR model approximates the generation time by an exponential distribution with mean 1/γ days, the estimator of R(0) is [33, 34]
Throughout the summer holiday, we assume that the reproduction number is reduced to R A = αR(0). That is, the growth rate during the summer holiday, r1, is defined by
Let the weekly incidence be J k in week k. During summer holiday period, the conditionally expected value of Jk+1given J k is
where Δt is the length of reporting (i.e. 7 days). In week 37, the conditional expectation is
because, with an initial incidence i k in week k, we have
Although adding the information of test negative individuals could potentially yield a less biased estimate of r0, we do not have access to this data and so we disregard this issue for now. Assuming that the observed counts of cases are Poisson distributed within each period, the likelihood function to estimate r0 and α is
The maximum likelihood estimates of r0 and α are obtained by minimizing the negative logarithm of (39), and the 95% confidence intervals (CI) are computed by profile likelihood. From the maximum likelihood estimates, we compute the differences in peak prevalence and times to observe the peak between baseline and second scenarios using the SIR model (1). For simplicity, we adopt (S0, I0, U0) = (99998, 1, 1) for the numerical computation and t1 is assumed to be 50 days (roughly corresponding to the length of period A plus 8 days).
Although our model (1) adopts exponentially distributed generation time, we can partially address the uncertainty of r0 and α in the parametric assumption of the generation time. That is, we adopt constant generation time (i.e. delta function) as an alternative assumption, which is known to yield a theoretical maximum reproduction number, given identical r0 and mean generation time . Given r0, the estimator of R(0) with constant generation time of 1/γ days reads
As mentioned above, the reproduction number during summer holiday reduces to R A = αR(0). Accordingly, the growth rate during the summer holiday, r1 is written as
Using the above mentioned likelihood (39) and replacing r1 of exponential assumption by (41), we estimate r0 and α. It should be noted that the difference between (33) and (41) indicates that the estimates of both r0 and α depend on the realistic distribution of the generation time .
In addition to the analysis of empirical data, we also examine sensitivity of the height and time of peak prevalence to different values of α and t1 by numerical simulation. As mentioned above, influenza is our case study, and the default value of R(0) of baseline scenario is 1.4, but we also consider R(0) of 1.2 and 1.6. These ranges are adopted, additionally, because we impose an approximation (16). When obtaining the numerical solutions, initial condition is fixed at (S0, I0, U0) = (99998, 1, 1) for clarity. U0 = 1 is adopted to prevent U0 = 0 in (18) and in later equations, and also to select a positive integer value closest to 0 such that U0/N ≈ 0.
Influenza A (H1N1-2009)
Figure 3A compares observed and predicted numbers of influenza cases in Japan from week 31 to 42. Grey bars represent conditionally expected values during summer holiday, and white bars represent the expected values during autumn semester. The estimated growth rate in the absence of summer holiday, r0 is 0.048 (95% CI: 0.029, 0.066) per day. Thus, the estimated reproduction number R(0) is 1.14 (95% CI: 1.09, 1.20) which is likely an underestimate (see below).
The estimated relative proportion of the reproduction number under mitigation conditions α is 0.948 (95% CI: 0.842, 1.053). Although the confidence limits of α include 1, our argument adopts linear dynamics for longer than 10 weeks (note that the largest number of notifications is seen in week 48), and thus, the reproduction number is conservatively estimated (i.e. over the time period that we examine, the transmission dynamics may become nonlinear); therefore, R(0) is likely underestimated due to the linear approximation adopted in our quantitative illustration). Thus, we believe it is appropriate to regard the reduction in the reproduction number during summer holiday as marginally significant. The estimated reproduction number during school holiday is 1.08.
Even adopting constant generation time of 3 days, r0 is of a similar order, i.e. 0.048 (95% CI: 0.029, 0.066) per day. The reproduction number R(0), however, is slightly greater (1.15 (95% CI: 1.09, 1.22)) due to its estimator, exp(r0/γ). α is estimated at 0.942 (95% CI: 0.832, 1.061).
Figure 3B illustrates the number of infectious individuals in a hypothetical population with 100,000 individuals using the estimated α and R(0). In the absence of summer holiday, the epidemic peak would have been observed at Day 171 with I(tm,0) = 822 cases. In the presence of summer holiday from time 0 to t1, the peak is delayed to Day 192 with I(tm,1) = 820 cases. Thus, the estimated α and R(0) do not significantly alter the height of peak prevalence when the effects of summer holiday are included (a difference between the two scenarios of only 2 cases), but the delay effect between the scenarios is as long as 21 days. We do not use our approximate formula for the estimation of time delay in this empirical case study (for reasons explained below). Because R(0) = 1.14 in the absence of intervention, a major epidemic can occur when the initial condition U0/N is smaller than 1 - 1/α R(0) = 7.4% of the population. Assuming a fixed I0 = 1, and varying S0 and U0 from N - 1 to 0.926N and from 0 to 0.074N, respectively, only slight variations in the reduction of peak prevalence (not greater than 1 case) are observed, but the time delay varies greatly; under a scenario with αR(0) ≈ 1 or with S0 = 0.926N and U0 = 0.074N, a possible maximum delay of t1 = 50 days can be readily envisaged.
Differential peak prevalence
Figure 4A examines the sensitivity of relative peak prevalence to α (i.e. reduction in R(0)) for assumed R(0) of 1.2, 1.4 and 1.6. Because αR(0) < 1 prevents major epidemic during the early epidemic phase, possible ranges of α satisfying αR(0) ≥ 1 vary with R(0). It is worth noting that the relative reduction in peak prevalence is a nonlinear function of α. Largest reduction occurs when α lies within the range 0.90 to 0.95, rather than when α is minimum. Figure 4B examines the sensitivity of relative reduction in the peak prevalence as a function of the time length of intervention t1 (e.g. the time required to switch control policy from containment to mitigation). Again, to satisfy t1 < tm,0, possible ranges of t1 vary with R(0). The interpretation of Figure 4B is more straightforward than that of Figure 4A. Essentially, the longer the time period of intervention, the larger the potential reduction of prevalence.
The nonlinear relationship observed in Figure 4A can be explored by combining (11) with approximate solutions. First, because S(t1) ≈ S(0) exp (-α R(0)U(t1)/S0) in the presence of intervention, equation (11) can be expressed in the form
Second, using an approximate solution of U(t1) based on logistic equation (18),
which is a nonlinear function of α. Although the calculation is not shown here due to its mathematical complexity and lack of practical interpretation, the derivative of (43) with respect to α reveals an optimal α yielding the longest delay in Figure 4A.
Delay in epidemic peak
Figure 5A compares epidemic curves of infectious individuals in the absence of intervention between explicit numerical and approximate solutions (i.e. solutions to (1) and (19), respectively). The height of epidemic peak is approximated well for smaller R(0), reflecting the fact that Taylor series expansion is a good approximation to the exponential function. The relationship between R(0) and approximation of epidemic peak height is also analytically expressed. Inserting (20) into (19), the approximate peak prevalence is
It should be noted that S0/R(0) can be replaced by γ/β, and thus, the peak prevalence in the approximated logistic form is independent of the initial condition (indeed, the derivative of a logistic equation is known not to depend on initial condition but rather on the carrying capacity ). On the other hand, an explicit solution can depend on initial condition for U(0) > 0, i.e.,
These two quantities are identical for R(0) = 1, N ≈ S0 and U0/N ≈ 0. Otherwise, numerical observation shows that I(tm,0) > I approx (tm,0) for R(0) > 1 and U0/N ≈ 0. Thus, the smaller the R(0), the better the approximate height of epidemic curve.
While interpretation of the height of peak prevalence is overall straightforward, the time to observe peak prevalence is better captured for larger R(0) (Figure 5A). Clearly, the logistic equation applied to R(0) = 1.2 yields a considerably biased (delayed) time to observation of epidemic peak (with bias longer than 20 days), and thus, we did not apply our approximate solution to the above mentioned case study of influenza A (H1N1-2009). It must be noted that the relationship between R(0) and approximation of epidemic peak in Figure 5A is regulated not only by R(0) but also by the initial condition U0 in (20); that is, the good agreement of the time of peak between two solutions for R(0) = 1.6 is not only due to R(0) but also to U0 = 1 in our simulation setting. For suitable values of U0, good approximations to the time of peak are obtained even for R(0) = 1.2 or smaller (results not shown).
Figure 5B compares the estimated time delay in epidemic peak, induced by an early intervention, between explicit numerical and approximate solutions. For all three R(0) that we investigate, approximation methods result in underestimation of the delay (31). For R(0) = 1.6, the approximate estimate of delay is crudely realized, and its sensitivity to α is close to that of the explicit numerical solution. The approximation is worst for R(0) = 1.2, for which a negative result was yielded for large α. These findings are in concordance with Figure 5A. It should be remembered that approximation of time of epidemic peak can vary with initial condition U0, indicating that the estimation requires knowledge of U0 in addition to R(0) and α (also, as we have seen, good approximation depends on judicious choice of R(0) and U0).
In the present study, we have presented fundamental ideas to assess the impact of a short-lasting intervention of an infectious disease on the epidemic peak. As a first step towards explicit evaluation of control effort in lowering and delaying the epidemic peak, we comprehensively described analytical expressions for the difference in the height of, and the time delay in, the epidemic peak gained by intervention, employing a parsimonious homogenous mixing epidemic model. We restricted our focus to a simple step function (Figure 1A) which adequately illustrated the role of summer holiday in lowering and delaying epidemic peak during the influenza (H1N1-2009) pandemic. Our methods show that both the height and the time of epidemic peak can be readily controlled by varying initial conditions at a given point of time at which transmission rate abruptly changes. The proposed method can be extended in future to encompass more realistic multiple steps for the transmission rate. Analytical solution of the simplest form of Kermack and McKendrick model has been undertaken multiple times, and is documented in many key references [27, 29–32, 38, 39] including the original study in 1927 . However, to our knowledge, the present study is the first to offer a theoretical basis on which to assess the impact of an early countermeasure on the epidemic peak employing logistic and logarithmic-form solutions, with a goal to making statistical inferences in the future. In addition, our analytical expressions not only promote epidemiological understanding of epidemic peak traits, but can be worked backwards to determine required control effort to achieve desired goals of height and delay of the epidemic peak. Although manual adjustment of the efficacy of intervention is not practically feasible, our demonstration of a nonlinear relationship between α and decline in the height of epidemic peak should be considered a key issue for optimal intervention management during an early epidemic phase.
Although public health guidelines have tended to advocate control policy switches from containment to mitigation at some point in time, the likely impact of unsuccessful containment to epidemic dynamics has been seldom discussed. Indeed, common illustration of mitigation (embodied in the three distinct goals (a)-(c) described in the Background section) do not account for the time-dependent transmission rate (rather, a guideline illustrates only differential peaks with various reproduction numbers R(0) ). Statistical modeling studies tend to focus only on time-dependent decreases in transmission potential, with a focus on the instantaneous reproduction number R(t) < 1 to demonstrate successful control of an infectious disease (e.g. [40–42]). Motivated by these problems, we considered the impact of upward change in the transmission rate during the course of an epidemic. In our case study of influenza H1N1-2009, it was shown that summer vacation did not appreciably lower the height of an epidemic peak, but imposed a substantial time delay. Although our approximate estimation of the delay in epidemic peak was biased towards certain choices of R(0) and the initial number of immune individuals U0, this does not imply failure of the model. Rather, by means of approximate analytical computation, we have shown in (20) and (30) that the time of epidemic peak critically depends on the fraction of initially immune individuals prior to an epidemic. Although a preceding study instead emphasized the dependence of the time of peak on initial number of infectious individuals , other studies with hyperbolic-form solutions emphasize dependency of the time of peak on U0[27, 31]. We believe that U0 is more easily quantified than the initial number of infectious individuals I0 in practical settings. The critical importance of the initial number of immune individuals U0 is especially highlighted in the influenza A (H1N1-2009) pandemic because of pre-existing immunity [43–46]. Our analytical undertakings indicate that statistical inference of the time delay gained by early control effort will greatly benefit from population-wide seroepidemiological survey. Depending on the quality of approximation for a given combination of R(0) and U0, one can then decide whether the estimation of delay should be based on analytical or numerical solution.
Despite our motivation to eventually offer a method to estimate the impact of an early intervention, it should be noted that the present study does not account for uncertainty (e.g. confidence interval). Our arguments are based solely on deterministic models, whereas an explicit derivation of the confidence interval requires use of a stochastic Markov jump process . Moreover, potential model extensions are numerous. Relevant factors include, but are not limited to, mobility of host, spatial dynamics, class-age structure (e.g. infection-age dependency), chronological age-structure, social contact patterns, seasonality, strain specificity and immunological dynamics. All of these features would increase the realism, but would greatly complicate analytical inspection, of the model. To illustrate the way forward, we discuss the simplest example of a class-age structured model in the Appendix.
Despite many future tasks to be completed, and our realization that the epidemic peak is vulnerable to heterogeneous patterns of transmission, relevant statistical assessment (including the estimation of R(0), generation time, and incubation period) always starts with a homogeneous modeling assumption [11, 21, 23, 33, 34, 48, 49]. This is particularly true during the early epidemic phase of a pandemic. In this sense, we believe that the present study has successfully offered a methodological avenue to statistically assess empirical data and to assess required control effort to lower and delay epidemic peak.
This study has presented a theoretical basis on which to assess the impact of short-lasting intervention on the epidemic peak of an infectious disease. Employing a homogeneously mixing epidemic model, we derived analytical expressions for the decline in the height of epidemic peak and for the time delay of the peak. Empirical influenza A (H1N1-2009) data were analyzed using a simplistic but practically accessible model, which estimated that the epidemic peak was delayed for 21 days by the summer holiday period in 2009. Approximate logarithmic form solution of the time of epidemic peak appeared to critically depend on initial condition of immune individuals, supporting a need to conduct population-wide serological survey. Despite obvious needs to address various types of heterogeneity, our framework offers a successful methodological avenue to assess relevant empirical data and to advocate required control effort to lower and delay epidemic peak.
Appendix: A way forward
In realistic situations, there is a time delay for newly infected individuals to acquire infectiousness, known as the latent period. This delay is captured by employing the so-called SEIR (susceptible-exposed-infected-recovered) model. In addition to S(t), I(t) and U(t), we consider infected but non-infectious individuals E(t). Assuming that the mean latent period is 1/δ days, the model is written as
The Lyapunov function is known to yield a constant solution (for any t), thus, dW/dt = 0. This is confirmed by
At an epidemic peak at time t m (where we have dE/dt = dI/dt = 0), two obvious conditions apply,
Using Lyapunov function in (47), we have
Writing both sides of (51) in terms of (47), and taking logarithm of both sides, we obtain
Rearranging (52) leads us to
which is very close to (6). By varying initial conditions at each time point when a change in transmission rate occurs, equation (53) permits us to measure the impact of public health intervention on the epidemic peak, under the SEIR assumption. Regarding the time of epidemic peak, a straightforward extension of (15) applies, although the analytical solution may be rather complex or may not exist. We have
S(t) can be replaced as for (14). Expressing E(t) as a function of U(t) requires strong mathematical supports. The problem may be addressed in different ways, but we first consider an analytical solution of dE/dt, i.e.,
where in the right-hand side is equivalent to dS/dσ for which the derivative of approximation (15) can be used. We have yet to derive a simple analytical solution of the time of epidemic peak from (54), but the above discussion demonstrates that it is at least possible to compute the height of epidemic peak with more E compartments (as was discussed in ), using our proposed Lyapunov approach. Nevertheless, it should be remembered that the presence of infection-age dependency in the infectiousness profile is likely to complicate the computation of epidemic peak . The incorporation of other realistic features, however, is greatly aided by the Lyapunov approach. In further studies, we will address epidemic peak with seasonality and age-dependent heterogeneity and in a multi-strain system.
World Health Organization: Mathematical modelling of the pandemic H1N1 2009. Wkly Epidemiol Rec. 2009, 84: 341-348.
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-1661. 10.1126/science.1176062.
Boelle PY, Bernillon P, Desenclos JC: A preliminary estimation of the reproduction ratio for new influenza A(H1N1) from the outbreak in Mexico, March-April 2009. Euro Surveill. 2009, 14: pii:19205
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: pii:19227
McBryde E, Bergeri I, van Gemert C, Rotty J, Headley E, Simpson K, Lester R, Hellard M, Fielding J: Early transmission characteristics of influenza A(H1N1)v in Australia: Victorian state, 16 May - 3 June 2009. Euro Surveill. 2009, 14: pii:19363
White LF, Wallinga J, Finelli L, Reed C, Riley S, Lipsitch M, M P: Estimation of the reproductive number and the serial interval in early phase of the 2009 influenza A/H1N1 pandemic in the USA. Influenza Other Respi Viruses. 2009, 3: 267-276. 10.1111/j.1750-2659.2009.00106.x.
Department of Health and Human Services, Centers for Disease Control and Prevention: Interim Pre-pandemic Planning Guidance: Community Strategy for Pandemic Influenza Mitigation in the United States. Early, Targeted, Layered Use of Nonpharmaceutical Interventions. 2007, Washington, D.C.: Department of Health and Human Services
Ferguson NM, Cummings DA, Fraser C, Cajka JC, Cooley PC, Burke DS: Strategies for mitigating an influenza pandemic. Nature. 2006, 442: 448-452. 10.1038/nature04795.
Halder N, Kelso JK, Milne GJ: Analysis of the effectiveness of interventions used during the 2009 A/H1N1 influenza pandemi. BMC Public Health. 2010, 10: 168-10.1186/1471-2458-10-168.
Lipsitch M, Cohen T, Murray M, Levin BR: Antiviral resistance and the control of pandemic influenza. PLoS Med. 2007, 4: e15-10.1371/journal.pmed.0040015.
Scalia-Tomba G, Wallinga J: A simple explanation for the low impact of border control as a countermeasure to the spread of an infectious disease. Math Biosci. 2008, 214: 70-72. 10.1016/j.mbs.2008.02.009.
Cowling BJ, Lau LL, Wu P, Wong HW, Fang VJ, Riley S, Nishiura H: Entry screening to delay local transmission of 2009 pandemic influenza A (H1N1). BMC Infect Dis. 2010, 10: 82-10.1186/1471-2334-10-82.
Kermack WO, McKendrick AG: Contributions to the mathematical theory of epidemics. I. Proc R Soc Ser A. 1927, 115: 700-721. 10.1098/rspa.1927.0118. (reprinted in Bull Math Biol 1991, 115: 33-55)
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.
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.
Cauchemez S, Valleron AJ, Boelle 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.
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.
Nishiura H: Travel and age of influenza A (H1N1) 2009 virus infection. J Travel Med. 2010, 17: 269-270. 10.1111/j.1708-8305.2010.00418.x.
Nishiura H, Cook AR, Cowling BJ: Assortativity and the probability of epidemic extinction: A case study of pandemic influenza A (H1N1-2009). Interdiscip Perspect Infect Dis. 2011, 2011:Article ID 194507
Cowling BJ, Lau MS, Ho LM, Chuang SK, Tsang T, Liu SH, Leung PY, Lo SV, Lau EH: The effective reproduction number of pandemic influenza: prospective estimation. Epidemiology. 2010, 21: 842-846. 10.1097/EDE.0b013e3181f20977.
Lessler J, Reich NG, Cummings DA, New York City Department of Health and Mental Hygiene Swine Influenza Investigation Team, Nair HP, Jordan HT, Thompson N: Outbreak of 2009 pandemic influenza A (H1N1) at a New York City school. N Engl J Med. 2009, 361: 2628-2636. 10.1056/NEJMoa0906089.
Cauchemez S, Donnelly CA, Reed C, Ghani AC, Fraser C, Kent CK, Finelli L, Ferguson NM: Household transmission of 2009 pandemic influenza A (H1N1) virus in the United States. N Engl J Med. 2009, 361: 2619-2627. 10.1056/NEJMoa0905498.
Cowling BJ, Fang VJ, Riley S, Malik Peiris JS, Leung G M: Estimation of the serial interval of influenza. Epidemiology. 2009, 20: 344-347. 10.1097/EDE.0b013e31819d1092.
Nishiura H: Time variations in the transmissibility of pandemic influenza in Prussia, Germany, from 1918-19. Theor Biol Med Model. 2007, 4: 20-10.1186/1742-4682-4-20.
Nishiura H: Time variations in the generation time of an infectious disease: Implications for sampling to appropriately quantify transmission potential. Math Biosci Eng. 2010, 7: 851-869. 10.3934/mbe.2010.7.851.
Pitzer VE, Lipsitch M: Exploring the relationship between incidence and the average age of infection during seasonal epidemics. J Theor Biol. 2009, 260: 175-185. 10.1016/j.jtbi.2009.06.008.
Bailey NTJ: Macro-modelling and prediction of epidemic spread at community level. Math Modelling. 1986, 7: 689-717. 10.1016/0270-0255(86)90128-4.
Brauer F, Castillo-Chavez C: Mathematical Models in Population Biology and Epidemiology. 2000, New York: Springer
Bailey NTJ: The Mathematical Theory of Infectious Diseases and Its Applications. 1975, London: Griffin, 2
Kendall DG: Deterministic and stochastic epidemics in closed populations. Proceedings of the Third Berkeley Symposium on Mathematical Statistics and Probability. Edited by: Neyman J. 1956, Berkeley: University of California Press, 4: 149-165.
Banks RB: Growth and Diffusion Phenomena: Mathematical Frameworks and Applications. 1991, New York: Springer
Risch H: An approximate solution for the standard deterministic epidemic model. Math Biosci. 1983, 63: 1-8. 10.1016/0025-5564(83)90047-0.
Wallinga J, Lipsitch M: How generation intervals shape the relationship between growth rates and reproductive numbers. Proc R Soc Lond Ser B. 2007, 274: 599-604. 10.1098/rspb.2006.3754.
Roberts MG, Heesterbeek JA: Model-consistent estimation of the basic reproduction number from the incidence of an emerging infection. J Math Biol. 2007, 55: 803-816. 10.1007/s00285-007-0112-8.
Nishiura H, Chowell G, Heesterbeek H, Wallinga J: The ideal reporting interval for an epidemic to objectively interpret the epidemiological time course. J R Soc Interface. 2010, 7: 297-307. 10.1098/rsif.2009.0153.
Lessler J, Santos TD, Aguilera X, Brookmeyer R, PAHO Influenza Technical Working Group, Cummings DA: H1N1pdm in the Americas. Epidemics. 2010, 2: 132-138. 10.1016/j.epidem.2010.07.001.
Nishiura H: Joint quantification of transmission dynamics and diagnostic accuracy applied to influenza. Math Biosci Eng. 2011, 8: 49-64. 10.3934/mbe.2011.8.49.
Nishiura H: Lessons from previous predictions of HIV/AIDS in the United States and Japan: epidemiologic models and policy formulation. Epidemiol Perspect Innov. 2007, 4: 3-10.1186/1742-5573-4-3.
Dietz K: Epidemics: the fitting of the first dynamic models to data. J Contemp Math Anal. 2009, 44: 97-104. 10.3103/S1068362309020034.
Nishiura H, Satou K: Potential effectiveness of public health interventions during the equine influenza outbreak in racehorse facilities in Japan, 2007. Transbound Emerg Dis. 2010, 57: 162-170. 10.1111/j.1865-1682.2010.01134.x.
Heijne JC, Teunis P, Morroy G, Wijkmans C, Oostveen S, Duizer E, Kretzschmar M, Wallinga J: Enhanced hygiene measures and norovirus transmission during an outbreak. Emerg Infect Dis. 2009, 15: 24-30. 10.3201/1501.080299.
Nishiura H, Omori R: An epidemiological analysis of the foot-and-mouth disease epidemic in Miyazaki, Japan, 2010. Transbound Emerg Dis. 2010, 57: 396-403. 10.1111/j.1865-1682.2010.01162.x.
Reichert T, Chowell G, Nishiura H, Christensen RA, McCullers JA: Does Glycosylation as a modifier of Original Antigenic Sin explain the case age distribution and unusual toxicity in pandemic novel H1N1 influenza?. BMC Infect Dis. 2010, 10: 5-10.1186/1471-2334-10-5.
Itoh Y, Shinya K, Kiso M, Watanabe T, Sakoda Y, Hatta M, Muramoto Y, Tamura D, Sakai-Tagawa Y, Noda T, Sakabe S, Imai M, Hatta Y, Watanabe S, Li C, Yamada S, Fujii K, Murakami S, Imai H, Kakugawa S, Ito M, Takano R, Iwatsuki-Horimoto K, Shimojima M, Horimoto T, Goto H, Takahashi K, Makino A, Ishigaki H, Nakayama M, Okamatsu M, Takahashi K, Warshauer D, Shult PA, Saito R, Suzuki H, Furuta Y, Yamashita M, Mitamura K, Nakano K, Nakamura M, Brockman-Schneider R, Mitamura H, Yamazaki M, Sugaya N, Suresh M, Ozawa M, Neumann G, Gern J, Kida H, Ogasawara K, Kawaoka Y: In vitro and in vivo characterization of new swine-origin H1N1 influenza viruses. Nature. 2009, 460: 1021-1025.
Greenbaum JA, Kotturi MF, Kim Y, Oseroff C, Vaughan K, Salimi N, Vita R, Ponomarenko J, Scheuermann RH, Sette A, Peters B: Pre-existing immunity against swine-origin H1N1 influenza viruses in the general human population. Proc Natl Acad Sci USA. 2009, 106: 20365-20370. 10.1073/pnas.0911580106.
Igarashi M, Ito K, Yoshida R, Tomabechi D, Kida H, Takada A: Predicting the antigenic structure of the pandemic (H1N1) 2009 influenza virus hemagglutinin. PLoS One. 2010, 5: e8553-10.1371/journal.pone.0008553.
Andersson H, Britton T: Stochastic Epidemic Models and Their Statistical Analysis. 2000, New York: Springer
Vynnycky E, Trindall A, Mangtani P: Estimates of the reproduction numbers of Spanish influenza using morbidity data. Int J Epidemiol. 2007, 36: 881-889. 10.1093/ije/dym071.
Nishiura H, Inaba H: Estimation of the incubation period of influenza A (H1N1-2009) among imported cases: Addressing censoring using outbreak data at the origin of importation. J Theor Biol. 2011, 272: 123-130. 10.1016/j.jtbi.2010.12.017.
Vincent TL, Grantham WJ: Nonlinear and Optimal Control Systems. 1997, New York: Wiley-Interscience
Feng Z: Final and peak epidemic sizes for SEIR models with quarantine and isolation. Math Biosci Eng. 2007, 4: 675-686.
Wearing HJ, Rohani P, Keeling MJ: Appropriate models for the management of infectious diseases. PLoS Med. 2005, 2: e174-10.1371/journal.pmed.0020174.
Rost G: SEIR epidemiological model with varying infectivity and infinite delay. Math Biosci Eng. 2008, 5: 389-402.
We would like to thank three anonymous reviewers for helpful comments on earlier draft of this paper. HN is supported by the Japan Science and Technology Agency PRESTO program. RO is financially supported by Research Fellowship of Japan Society for the Promotion of Science.
The authors declare that they have no competing interests.
HN conceived of the study, developed methodological ideas and implemented mathematical and statistical analyses. RO revised analytical results. HN drafted the manuscript and RO and HN discussed and revised the manuscript. All authors read and approved the final manuscript.