 Research
 Open Access
 Published:
Analytical estimation of maximum fraction of infected individuals with oneshot nonpharmaceutical intervention in a hybrid epidemic model
BMC Infectious Diseases volume 22, Article number: 512 (2022)
Abstract
Background
Facing a global epidemic of new infectious diseases such as COVID19, nonpharmaceutical interventions (NPIs), which reduce transmission rates without medical actions, are being implemented around the world to mitigate spreads. One of the problems in assessing the effects of NPIs is that different NPIs have been implemented at different times based on the situation of each country; therefore, few assumptions can be shared about how the introduction of policies affects the patient population. Mathematical models can contribute to further understanding these phenomena by obtaining analytical solutions as well as numerical simulations.
Methods and results
In this study, an NPI was introduced into the SIR model for a conceptual study of infectious diseases under the condition that the transmission rate was reduced to a fixed value only once within a finite time duration, and its effect was analyzed numerically and theoretically. It was analytically shown that the maximum fraction of infected individuals and the final size could be larger if the intervention starts too early. The analytical results also suggested that more individuals may be infected at the peak of the second wave with a stronger intervention.
Conclusions
This study provides quantitative relationship between the strength of a oneshot intervention and the reduction in the number of patients with no approximation. This suggests the importance of the strength and time of NPIs, although detailed studies are necessary for the implementation of NPIs in complicated realworld environments as the model used in this study is based on various simplifications.
Introduction
Because of the global spread of COVID19, our human society is facing a major public health crisis. The COVID19 pandemic is caused by an emerging pathogen, SARSCoV2, for which there is no immunized population, causing an overshooting increase in the number of infected patients and depleting medical resources in many countries. Medical institutions are facing a difficult situation in which they must control second transmissions while treating critically ill patients, and as the number of patients increases, the medical system becomes swiftly tighter. When the number of patients exceeds the capacity, the quality of medical care deteriorates drastically, and the number of medical devices required for life support reaches its limit. This situation further increases the fatality rate of this infectious disease and causes serious damage to our society.
No effective treatment for COVID19 has been established yet as of September, 2020, and only a public health approach can function as a control measure for the epidemic. To mitigate the spread of COVID19, each country is implementing nonpharmaceutical interventions (NPIs) [1] to regulate social activities. NPIs comprise policies such as case isolation, voluntary home quarantine, closure of schools and universities, social distancing, stopping mass gatherings, and border closure. In major European countries, these NPIs were implemented, depending on the epidemic situation, during the first part of spring in 2020, with different timings and intensities [2, 3]. Although the first phase of the epidemic appeared to be suppressed by these mitigation measures, the reepidemic became clearer in almost all countries because of deregulation after the first epidemic.
Because the control of epidemics by NPIs has caused a situation involving the imposition of strong restrictions on human socioeconomic activities, it would be desirable to study in advance the optimal timing, intensity, and duration of interventions that could bring about more promising results with minimal damage to society. Regarding COVID19, for the purpose of expost verification, the influence of NPIs implemented in each country on the effective reproduction number was estimated. In European countries, the correlation between the decay of the effective reproduction number and the implementation of various NPIs has been verified [3]. The impact of travel limitations in China on the spread of infection has been discussed [4, 5]. However, a reliable estimation of the effects of NPIs is difficult, because the differences in NPI strategies employed by each country are strongly related to various background factors, such as the epidemic situations, social structure, legal systems, and culture [2].
Mathematical modelling is an important method for estimating the effect of NPIs. In particular, the global pandemic of COVID19 revealed that a situation in which only NPIs are effective against an emerging infectious disease is possible in societies in the 2020s. Such a situation had been “neglected” as a practical research target, which enhances the importance of theoretical approaches. Recently, compartmental models such as the SIR model [6,7,8,9] have been extended to estimate the effects of NPIs on the number of patients [10,11,12,13]. In other settings, optimal policies where intervention intensity can change continuously over time have been discussed in the context of minimizing objective functions [14,15,16,17]. These solutions are very useful if we can estimate the effect of NPIs on the change in transmission parameters precisely.
The impacts of NPIs should be evaluated as a discontinuous change of the transmission rate in a model to represent the temporal discontinuity of intervention in the real world. The dynamics of the number of infected patients in continuous time in a system that includes discrete parameters, state spaces, and continuoustime dynamics can be modeled as a hybrid dynamical system [18,19,20,21,22]. This framework has been applied to mathematical models of infectious diseases [23,24,25]. In the simplest case, the oneshot intervention model, in which the intervention is implemented only once during the epidemic, can be used to discuss the theoretical dependence of the effects of NPIs on the timing and intensity [26]. As the COVID19 epidemic continues, the accumulation of theoretical research on the effects brought about by NPIs has become even more significant. Recently, compartmental models with intervention have been studied both numerically [27, 28] and using some analytical methods [29, 30] in line with the COVID19 epidemic.
In this study, we provide exact solutions of a simple SIR model with oneshot intervention, represented by a single discrete reduction in the transmission parameter during an epidemic. These solutions describe the dependence of the peak number of infected patients on the reproduction number under consideration of the implementation of NPIs and intervention timing. Theoretical and numerical analyses revealed nontrivial relations among the intensity of suppression of pandemics via NPIs, the number of infected individuals at the peaks, and the final size of infection cases.
The methods and results shown in this study provide basic theoretical understanding in the context of the evaluation of NPIs.
Materials and methods
In this study, we focus on the dependence of the maximum fraction of infected individuals on the timing of the NPIs. Note that this analysis is motivated by the COVID19 epidemic, but we consider a hypothetical epidemic, whose properties are not necessarily the same as that of COVID19.
Model
We here introduce the SIR model, where the time evolution of the fraction of susceptible (s(t)), infected (i(t)), and removed (r(t)) individuals is given by the following ordinary differential equation (ODE):
taking into account the oneshot intervention and the second wave after the intervention. For simplicity, the total population is assumed to be unity, that is, \(s(t)+i(t)+r(t)=1\) holds. The summation of the righthand sides of Eqs. (1)–(3) vanishes, which guarantees conservation of the total population. Therefore, as seen below, only Eqs. (1) and (2) are numerically integrated to obtain the time evolutions of s(t), i(t), and r(t). In this model, the basic reproduction number and the effective reproduction number at time t are given by \(R_0 = \beta /\gamma\) and \(R_t = \beta s(t) / \gamma\), respectively.
Kermack and McKendrick [6] derived the equation for the maximum fraction of infected individuals and showed that the final size is obtained by solving a transcendental equation with a given \(R_0\) in the SIR model without any intervention [7,8,9]. By applying the technique they employed in the derivation, we here provide the relationship between the fractions of susceptible and removed individuals at arbitrary times \(t_0\) and \(t_1\), that is, \(s(t_0)\), \(r(t_0)\), and \(s(t_1)\), \(r(t_1)\) as
See Additional file 1: Section S1 for details of the derivation. All analytical results presented in this paper are derived based on this equality.
Nonpharmaceutical intervention
In the present framework, an NPI in an isolated population is represented by a change in the transmission rate \(\beta\). Let \(\beta = \beta _{\mathrm{off}} (>\gamma )\) be the transmission rate without the intervention, and it is switched to \(\beta _{\mathrm{on}} (<\beta _\mathrm{off})\) when the intervention starts at \(t=t_{\mathrm{on}}\), and restored to \(\beta _{\mathrm{off}}\) at \(t=t_{\mathrm{off}}=t_{\mathrm{on}}+\Delta t\). See Fig. 1 for the schematic diagram of the setting. Let the corresponding basic reproduction numbers be \(R_{0,\mathrm{off}}=\beta _{\mathrm{off}}/\gamma\) and \(R_{0,\mathrm{on}}=\beta _{\mathrm{on}}/\gamma\), respectively. Here, \(\Delta t\) denotes the duration of the intervention. In this paper, we study the effect of oneshot intervention, one of the simplest implementation schemes, where the intervention is implemented only once. The mathematical methods employed in this study can be applied to more complex cases, for example, where multiple interventions are implemented intermittently. If the fraction of susceptible individuals remains large enough and the herd immunity is not achieved after the intervention, a second wave occurs. It is thus necessary to consider this second wave to evaluate the effect of the intervention.
The oneshot intervention setting has been numerically studied by Bootsma and Ferguson [26], showing that the final size, which depends on the timing \(t_{\mathrm{on}}\) and \(t_{\mathrm{off}}\), can be smaller for weaker intervention with larger \(R_{0,\mathrm{on}}\). Di Lauro et al. [28] numerically studied the dependence of the final size and the peak fraction of infected individuals on the timing of the intervention, where the intervention was assumed to start when the fraction of infected and recovered individuals exceeds a threshold value. In particular, they concluded that the onset timing should be chosen so that the two peaks during and after the intervention are comparable. Morris et al. [29] showed that the maximum fraction of infected individuals with oneshot interventions can approach that achieved by the optimal intervention, which requires an unrealistic intervention, such as \(R_{0,\mathrm{on}}=0\). Sadeghi et al. [30] also suggested the existence of the optimal timing of the intervention based on discussion using a linearized equation and numerical simulation.
Analysis
Equations (1) and (2) were numerically integrated using opensource numerical solvers in Python. The codes and datasets generated by them are available in the repository [31]. The time step was set to 0.001. The nonintervention transmission rate and the recovery rate are fixed as \(\beta _{\mathrm{off}}=2/7 \ \mathrm{days}^{1}\) and \(\gamma =1/7 \ \mathrm{day}^{1}\), that is, \(R_{0,\mathrm{off}}=2\). The initial condition \(s(0)=1\epsilon\), \(i(0)=\epsilon\), and \(r(0)=0\), where \(\epsilon =0.001\) is employed in this study. As discussed in detail below, the system behaves qualitatively differently for different values of \(R_{0,\mathrm{on}}\). Here, we primarily focus on the cases for \(R_{0,\mathrm{on}}=1.4\) (Fig. 2) and \(R_{0,\mathrm{on}}=0.7\) (Fig. 3), corresponding to relatively weaker and stronger intervention intensities, respectively. Note that \(R_{0,\mathrm{on}}>1\) in the former case implies that the infection may spread even in the presence of the intervention. The time series without the intervention is shown in Figs. 2A and 3A. They are identical, because \(R_{0,\mathrm{on}}\) does not affect the dynamics without the intervention. Figures 2B–E and 3B–E show the time series with an intervention with a constant intervention duration \(\Delta t=60\) days. The second wave can occur if the herd immunity is not achieved when the intervention ends (Figs. 2C and 3B, C). Note that these time series approximate the time evolution observed in an agentbased model on a random network with the corresponding parameter values and initial conditions. The comparisons of the time series are given in Additional file 1: Section S2.
Maximum fraction of infected individuals
In addition to the numerically obtained time series, we analytically show that peaks of the fraction of infected individuals can appear at the following four timings by applying Eq. (4): after the intervention, during the intervention, at the onset of the intervention, and before the intervention. Note that there can be two peaks of the fraction of infected individuals if a second wave occurs. In such a case, the peak with a larger fraction gives the global maximum. Let us describe the four cases with respect to the timing of the maximum in detail.

(i)
The maximum appears after the intervention (Figs. 2B, C, 3B, C). If the intervention ends before achieving herd immunity, a peak is observed during the second wave after the intervention. The fraction of infected individuals for this peak gives the global maximum if it is higher than the first peak before or during the intervention.

(ii)
The maximum appears during the intervention (Fig. 2D). If the effective reproduction number \(R_t\) declines and crosses unity during the intervention, there is a peak in this timing. The condition \(R_{0, \mathrm{on}}>1\) is necessary for the existence of this peak, because the effective reproduction number has to be larger than unity at the onset of the intervention. Therefore, this peak does not appear for \(R_{0,\mathrm{on}}<1\). This peak is the global maximum if it is larger than the second peak. The peak also appears at this timing in Fig. 2C, but the second peak is slightly higher than this peak.

(iii)
The maximum appears at the onset of the intervention at \(t=t_{\mathrm{on}}\) (Fig. 3D). If \(R_{t_\mathrm{on}}<1\) holds when the intervention starts, the fraction of infected individuals decreases during the intervention, and we observe a local peak at \(t=t_{\mathrm{on}}\). There may be another peak if herd immunity is not achieved during the intervention. If the fraction of infected individuals at the second peak is less than at this peak, the timing of the global maximum is given at \(t_{\mathrm{on}}\). Note that this region also appears for \(R_{0,\mathrm{on}}>1\), although the time series is not shown in Fig. 2B–E.

(iv)
The maximum appears before the onset of the intervention (Figs. 2E and 3E). The fraction of infected individuals reaches its maximum before the intervention. This case implies that the intervention starts too late and fails to mitigate outbreaks in terms of the maximum fraction of infected individuals.
The fraction of infected individuals at the peaks can be calculated analytically. Conditions for the peaks are given in terms of the fraction of susceptible individuals as \(s(t)=1/R_{0,\mathrm{off}}=\gamma /\beta _{\mathrm{off}}\) without the intervention and \(s(t)=1/R_{0,\mathrm{on}}=\gamma /\beta _{\mathrm{on}}\) with the intervention. The maximum fraction of infected individuals for cases (ii), (iii), and (iv) does not depend on the fraction of removed individuals at the offset of the intervention \(r(t_{\mathrm{off}})\), because these peaks appear before \(t_{\mathrm{off}}\). For case (i), the maximum fraction of infected individuals \(i_{\mathrm{max}}\) is given as
which depends on \(\Delta r:= r(t_{\mathrm{off}})r(t_{\mathrm{on}})\), the difference in the fraction of removed individuals between the onset and offset of the intervention. See Additional file 1: Section S3 for the explicit form of \(i_{\mathrm{max}}\) for all cases and its derivation.
The boundaries between regions (i) and (ii), (i) and (iii), (ii) and (iii), and (iii) and (iv) can be obtained analytically with respect to \(r(t_{\mathrm{on}})\) and \(\Delta r\) (Figs. 2F and 3F). On the boundaries, the fraction of infected individuals at two peaks is comparable (Figs. 2C and 3C). See Additional file 1: Section S4 for details of the derivation.
Final size with intervention
The final size also reflects the effect of the intervention. The final size of removed individuals \(r(\infty )\) with the intervention is obtained by solving the equation
in a selfconsistent manner [32]. Specifically, as \(r(\infty )\) appears on both sides, this equation can be solved numerically or using the Lambert W function, except for some special cases. This equation implies that the final size depends on \(\Delta r = r(t_{\mathrm{off}})r(t_{\mathrm{on}})\). Another important implication of this equation is that \(r(t_{\mathrm{off}})\) has the upper bound \({\tilde{r}}\) depending on \(r(t_{\mathrm{on}})\), which is given by
See Additional file 1: Section S5 for the details of the derivation. The final state represented by Eq. (6) is the equilibrium in the presence of the intervention. One can show that the Jacobian matrix of the linearized equation for this final state has zero eigenvalue, which implies that this equilibrium is neutrally stable See Additional file 1: Section S6 for the derivation of the detailed discussion.
Using this equality, one can show that the final size in the presence of the intervention is always smaller than that without the intervention. For
one can achieve \(r(\infty )\approx 11/R_{0,\mathrm{off}}\) by setting \(t_{\mathrm{on}}\) properly, which is the smallest prevalence to achieve herd immunity, with an intervention duration \(\Delta t\) that is large enough [32]. See Supplementary Additional file 1: Section S6 for details of the derivation. Numerical results regarding the final size are summarized in Additional file 1: Section S7.
Results
We report the numerical and analytical results, when the reproduction number under the intervention \(R_{0,\mathrm{on}}\) is large (Fig. 2) and small (Fig. 3), showing qualitatively different behaviors.
In Figs. 2F, G, 3F, G, the dependence of the maximum fraction of infected individuals on the timing of the intervention is plotted. Here, \(i_{\mathrm{max}}\) is normalized by that in the absence of the intervention \(i_{\mathrm{max}}^{0}\) (Figs. 2A and 3A). As the maximum infected fraction drops in the presence of the intervention, \(i_{\mathrm{max}}/i_{\mathrm{max}}^0\) is less than unity and quantifies the effectiveness of the intervention in terms of the maximum fraction of infected individuals. As this ratio decreases, the intervention shows more success in reducing the maximum fraction.
It is difficult to obtain the time series of the SIR model analytically without any approximations, for example, linearization, or the method presented in [6]. Therefore, \(i_{\mathrm{max}}\) is numerically computed with different onset and offset times for the intervention, \(t_{\mathrm{on}}\) and \(t_{\mathrm{off}}\) (Figs. 2G and 3G). However, \(i_{\mathrm{max}}\) can be analytically calculated with respect to the fraction of the recovered individuals at the onset, \(r(t_\mathrm{on})\), and the offset, \(r(t_{\mathrm{off}})\), of the intervention (Figs. 2F and 3F). Note that there exists a onetoone correspondence between \(\{t_\mathrm{on},t_{\mathrm{off}} \}\) (panels (F)) and \(\{ r(t_{\mathrm{on}}), r(t_\mathrm{off}) \}\) (panels (G)) in Figs. 2 and 3.
Weak intervention (large \(R_{0,\mathrm{on}}=1.4\))
In this case, the reproduction number is larger than unity even in the presence of the intervention. Therefore, the fraction of infected individuals may increase during the intervention period. The peaks of infected individuals can be observed during the intervention, and the maximum infected fraction can appear at any of the four timings (i)–(iv) classified above. Figures 2B–E show that \(i_{\mathrm{max}}\) is minimized in the intermediate onset time \(t_{\mathrm{on}}\) near (C), where the peaks during and after the intervention are comparable. This nonmonotonic dependence on \(t_{\mathrm{on}}\) is clearly visualized in Fig. 2G. The peak of infected individuals during the intervention is smaller than that without intervention. This intervention mitigates the second wave.
As shown in Eq. (5) and Fig. 2F, the maximum fraction of infected individuals is linear in \(\Delta r\) if the peak of the second wave is the maximum, that is, case (i). This is verified by the fact that the contours lie horizontally in case (i). To clarify this point, the contours are explicitly shown in Additional file 1: Section S8. If the fraction of the infected individuals reaches the maximum during or at the onset of the intervention (cases (ii) and (iii), respectively), the maximum fraction depends only on \(t_{\mathrm{on}}\), which has onetoone correspondence to \(r(t_{\mathrm{on}})\). This is verified in Fig. 2F, G, and Additional file 1: Section S8. If the maximum appears before the intervention, that is, case (iv), the maximum is independent of both \(t_{\mathrm{on}}\) and \(t_{\mathrm{off}}\).
Onset and offset times for the intervention with constant \(\Delta t\) corresponding to Figs. 2A–E are plotted in Figs. 2F, G. As suggested by the time series, the maximum infected fraction is smallest in the intermediate intervention onset \(t_{\mathrm{on}}\), near (C). It is clear from panel (F) that this point is located close to the boundary between cases (i) and (ii), where peaks during and after the intervention are comparable. In case (ii), the maximum does not depend on \(\Delta r\), which implies that a longer intervention does not reduce the maximum. Note that the relationship between \((t_\mathrm{on},t_{\mathrm{off}})\) and \((r(t_{\mathrm{on}}),r(t_{\mathrm{off}}))\) is nonmonotonous.
Strong intervention (small \(R_{0,\mathrm{on}}=0.7\))
When the transmission rate is small during the intervention, the maximum infected fraction is minimized in the intermediate starting time of the intervention \(t_{\mathrm{on}}\), namely, early implementation of the intervention does not necessarily minimize the infected fraction (Figs. 3F, G). This result is intuitively understood as follows. If the intervention starts too early, the infection does not spread because of the small intervention transmission rate. Therefore, the second wave after the intervention is large, thus the early intervention is not effective. If the timing of the intervention is characterized in terms of \(r(t_{\mathrm{on}})\) and \(r(t_{\mathrm{off}})\) (Fig. 3F), the maximum of the second wave depends only on \(\Delta r\). If the maximum fraction is found during the intervention, its value is independent of \(t_{\mathrm{off}}\) and \(r(t_{\mathrm{off}})\). The peak does not appear during the intervention, that is, case (ii) does not appear because \(R_{t,\mathrm{on}}<1\) holds.
Maximum infected fraction \({\bar{i}}_{\mathrm{max}}\) versus reproduction number under the intervention \(R_{0,\mathrm{on}}\)
For each \(R_{0,\mathrm{on}}\), there exist onset and offset timings of the intervention that minimize the maximum fraction of infected individuals \(i_{\mathrm{max}}\) (Figs. 2F, G and 3F, G). Let this value be \({\bar{i}}_\mathrm{max}(R_{0,\mathrm{on}})\). Figure 4 plots the dependence of \({\bar{i}}_{\mathrm{max}}\) on \(R_{0,\mathrm{on}}\). As seen in the figure, it is minimized at a nontrivial intermediate value of \(R_{0,\mathrm{on}}=R_{0,\mathrm{on}}^{*}\approx 1.23>1\). This implies that the maximum fraction of infected individuals is minimized for a weak intervention under the oneshot condition. For \(R_{0,\mathrm{on}}\ge R_{0,\mathrm{on}}^*\), \({\bar{i}}_{\mathrm{max}}\) is achieved at the boundary between regions (i) and (ii) at \(t_{\mathrm{on}}=0\), where the peaks during and after the intervention are comparable (Fig. 2C). For \(R_{0,\mathrm{on}}^*\ge R_{0,\mathrm{on}} \ge R_{0,\mathrm{on}}^{(1)}\), the boundary between regions (i) and (ii) with \(t_{\mathrm{on}}>0\) gives \({\bar{i}}_{\mathrm{max}}\), where \(R_{0,\mathrm{on}}^{(1)}\approx 1.08\) is the parameter value below which region (ii) does not exist. For strong intervention \(R_{0,\mathrm{on}}\le R_{0,\mathrm{on}}^{(1)}\), \({\bar{i}}_{\mathrm{max}}\) is found at the boundary between regions (i) and (iii), where the peaks of the onset and after the intervention are comparable (Fig. 3(C)), at \(t_{\mathrm{on}}>0\). Conditions for \({\bar{i}}_{\mathrm{max}}\) are analytically derived in all cases, shown by the solid line in Fig. 4. The conditions for \({\bar{i}}_{\mathrm{max}}\) for \(R_{0,\mathrm{on}}\ge R_{0,\mathrm{on}}^{(1)}\) can be explicitly solved. The condition for small \(R_{0,\mathrm{on}}\le R_{0,\mathrm{on}}^{(1)}\) cannot be solved explicitly, and the parametric equations for \({\bar{i}}_{\mathrm{max}}\) and \(R_{0,\mathrm{on}}\) are used to plot the theoretical curve in Fig. 4. See Additional file 1: Section S9 for details of the derivation.
Intervention strategies
It is possile to optimize the onset and offset timings of the NPI by minimizing an objective function under certain constraints. In the present framework, this can be formulated as an optimization problem in a hybrid nonlinear dynamical system. The optimal intervention strategy depends on the objective function. Here, we discuss the following two simple scenarios to minimize \(i_{\mathrm{max}}\). In general, more complex objective functions can be used for optimization. The following discussion may provide some intuition for considering such cases, taking into account the second wave.
Minimizing \(i_{\mathrm{max}}\) with a constraint in the intervention duration
Let us start with a case where the intervention duration \(\Delta t\) is less than a certain value. As \(i_{\mathrm{max}}\) is monotonically decreasing in \(\Delta t\) for a fixed \(t_{\mathrm{on}}\), we can assume that \(\Delta t\) is a constant, and \(t_{\mathrm{on}}\) is varied. This case has been studied in Figs. 2B–E and 3B–E. This scheme fixes the intervention duration, so it is easier to anticipate the economic impact of the intervention, which depends on the duration of the intervention, than in the next scenario based on the fraction of recovered individuals.
As discussed above, \(i_{\mathrm{max}}\) is minimized in the intermediate onset time \(t_{\mathrm{on}}\), when the first peak during the intervention and the second peak after the intervention are comparable. For larger \(R_{0,\mathrm{on}}\), \(t_{\mathrm{on}}\) giving the minimum \(i_{\mathrm{max}}\) converges to zero for large \(\Delta t\) (Fig. 2G), but converges to nonzero for large \(\Delta t\) for smaller \(R_{0,\mathrm{on}}\) (Fig. 3(G)). These results are understood as follows: for a strong intervention with small \(R_{0,\mathrm{on}}\), the intervention immediately suppresses the fraction of infected individuals. Therefore, the early onset of the intervention prevents more individuals from achieving immunity during the intervention and eventually increases the maximum fraction of infected individuals.
Minimizing intervention duration with a constraint in \(i_{\mathrm{max}}\)
Another possible constraint is to minimize the intervention duration \(\Delta t\), keeping the maximum fraction of infected individuals \(i_{\mathrm{max}}\) constant. This corresponds to choosing \(t_{\mathrm{on}}\) along a contour of \(i_{\mathrm{max}}\), which prevents the overcapacity of medical support. In Figs. 2G and 3G, \(\Delta t\) is minimized for an intermediate \(t_{\mathrm{on}}\). Namely, a contour that crosses point (C) (\(\Delta t=60\)) reaches \(\Delta t\approx 20\) if the onset time is later than \(t_{\mathrm{on}}=19.2\) (Fig. 2(C)). As we have already discussed, the maximum fraction depends on \(\Delta r\) in region (i). It takes a shorter time to achieve the same \(\Delta r\) if the intervention starts later in this case. As the maximum fraction of infected individuals depends on \(r(t_{\mathrm{on}})\) but is independent of \(\Delta r\) in cases (ii) and (iii), the optimal \(r(t_{\mathrm{on}})\) is determined by the boundary between regions (i) and (ii) or that between (i) and (iii) in Figs. 2F and 3F. Evidently this timing is an intermediate value of \(t_{\mathrm{on}}> 0\). It may be easier to set a plan for the intervention with respect to \(r(t_{\mathrm{on}})\) and \(r(t_{\mathrm{off}})\), rather than the timing \(t_\mathrm{on}\) and \(t_{\mathrm{off}}\).
Similar optimization problems can be considered using the final size as the objective functions. These problems are discussed in Additional file 1: Section S10.
Conclusion
In the present study, in consideration of the actual COVID19 situation, we studied the situation in which the reproduction number of an infectious disease is temporarily reduced by implementing an NPI once during the epidemic, using a simple mathematical model. The results provide theoretical implications as to how strong NPIs should be introduced during an epidemic of an emerging infectious disease. If the effective reproduction number during the intervention is too small, the fraction of infected individuals at the peak in the second wave may be higher than the first peak. It was also shown numerically and analytically that the fraction of infected individuals can also increase if the intervention is started too early. The upper limit of medical capacity is an essential practical constraint. In particular, for infectious diseases such as COVID19, which is too emergent to expect effective treatments, it is more important to avoid exhausting the medical system. This study suggests that it will be necessary to be alert for a larger second wave that may occur after strong intervention in such cases.
We analytically derived the peak fractions of infected individuals in the SIR model with the oneshot intervention. These analytical results suggest that the peak fractions can be smaller with nontrivial intervention timings. The maximum fraction is smallest for \(R_{0,\mathrm{on}}>1\), that is, the intervention reproduction number is not in the diseasefree regime. In the literature, Bootsma and Ferguson [26] showed that the final size can be minimized for nontrivial intervention timing and \(R_{0,\mathrm{on}}\) numerically. Di Lauro et al. [28] numerically showed that the peak fraction of infected individuals also depends on the timing of the intervention. We obtained analytical expressions for these quantities in this study. The analytical results for the final size are presented in Additional file 1: Section S6. It should be noted that the formulae for the peak fraction depend on the timing of the peak, resulting in various cases compared with the final size. Some analytical results are available for this system. Sadeghi et al. [30] explained these nontrivial effects based on solutions of the linearized equation, which exponentially grows and decays without and with the intervention, respectively. Linearization is one of the simplest approximations and is applicable in this case; in particular, this approximation is useful in discussions regarding timing. Conversely, linearization cannot be used to discuss the important case where \(R_{0,\mathrm{on}}>1\), as the linearized equation cannot explain the declining number of infected individuals. The proposed method in this study provides a unified framework, including the cases where linearization is not feasible. Morris et al. derived [29] an equation for the peaks of the fraction of infected individuals in terms of \(s(t_{\mathrm{on}})\), \(s(t_{\mathrm{off}})\), \(i(t_{\mathrm{on}})\), and \(i(t_{\mathrm{off}})\). In this work, we further show the peak fractions in terms of the two parameters \(r(t_{\mathrm{on}})\) and \(r(t_{\mathrm{off}})\) for a general initial condition.
Although concrete measures and guidelines for COVID19 are required, it is emphasized again that the results in this study were derived using a simplified model with many assumptions. First, the results presented in this study are based on the simplified SIR model, which takes into account neither the realistic pathology of COVID19 nor societal response. This model assumes uniform and random contact within a group and does not consider interactions between different subgroups in the population. Recovered patients are assumed to have complete immunity in this model. These assumptions are not applicable to COVID19, where there are still many unclear factors regarding heterogeneity in contact networks and the immune response of patients after recovery. It has been shown that heterogeneous contacts can affect the infection dynamics [33]. The presented method could be applied in such a case, if the mean field approximation describes the epidemic process well. If a society develops other public health measures during the intervention, the basic reproduction numbers before and after the NPI may be different. Next, this study is a oneshot intervention model, in which the transmission coefficient returns to the original value after a single intervention. Practically, each repandemic requires multiple intermittent interventions [34], making the intervention process much more complex [35]. Change in the parameter \(\beta\) between \(\beta _{\mathrm{off}}\) and \(\beta _{\mathrm{on}}\) occurs discontinuously in the present model. However, it should take a finite time in realistic case reflecting the time for responses of people to an NPI [36]. The finite time interval for the parameter change gives a correction to the results presented here. Detailed discussions are given in Additional file 1: Section S11. Furthermore, the model is based on a deterministic dynamical system of an isolated population. Therefore, important NPI measures such as border control cannot be estimated in the present framework. The deterministic nature of the model assumes that the infectious disease cannot be eradicated, as the number of infected individuals remains nonzero for a finite time. If the population is small enough and no imported cases are assumed, strong and early intervention, which is not necessarily recommended in this study, may eradicate the disease, and a second wave does not occur. Strong intervention would be necessary in other cases, for example, when the number of infected patients approaches the capacity of the medical system. These complex situations may be analyzed in detail by extending the methods presented in this study to e.g., stochastic epidemic models [37], which may lead to different conclusions from those reached in this study. Another important situation would be cases where the endemic state exists [38, 39]. Such systems behave qualitatively differently, e.g., existence of the Lyapunov function may lead to the global stability of the steady state. It would be important to discuss an effective intervention strategies to such systems.
It is also important to consider dynamics using a different criterion, such as the fraction of new cases, as a trigger of an NPI. Here, let us briefly discuss the possibility to use i(t) for onset/offset criteria for the intervention. If we can map i(t) to r(t), then we can apply the present framework. The conservation of the total population \(i(t) = 1  r(t)  s(t) = 1  r(t)  \exp ( R_{0,\mathrm{off}} \, r(t))\) enables to determine i(t) from r(t) at the onset of the intervention. However, there can be multiple r(t) values corresponding to a given i(t). Therefore, if we specify the branch of r(t), we can apply the present framework using i(t) as the trigger of the intervention. The detailed formalization would be an imporntant future work.
In this study, the effect of NPIs was modeled as a hybrid dynamical system, which may further enable us to approach more refined models in future investigations. The influence of NPIs with respect to political decisions and behavioral changes of people can be expressed more accurately by introducing a hybrid dynamical system. In recent years, dynamical systems theory and control theory have been developed, and phenomena specific to hybrid systems such as Zeno solutions and sliding motions have been discussed [18,19,20,21,22]. Some studies, such as [40] and [41], have proposed control of infectious diseases using sliding mode control. The optimal policy under NPIs can be discussed by modeling the effect of economic damages associated with execution and then minimizing the cost function. Discrete changes in parameters such as the transmission coefficient in NPI implementation are strongly linked to the intensity of measures, suppression of economic activity, and changes in human mobility. For example, the correlation between the decrease in human mobility with NPIs and the effective reproduction number for the COVID19 pandemic has been studied [42, 43]. Such studies can contribute to modeling the costs of NPIs. The mathematical model of the epidemic suppression effect can be constructed using a hybrid dynamical system, taking into account the negative socioeconomic impact of NPIs.
Availability of data and materials
The datasets generated and analysed during the current study are available in the repository [31].
References
Aledort JE, Lurie N, Wasserman J, Bozzette SA. Nonpharmaceutical public health interventions for pandemic influenza: an evaluation of the evidence base. BMC Public Health. 2007;7(1):208.
Cohen J, Kupferschmidt K. Countries test tactics in war against COVID19. Science. 2020;367(6484):1287–8.
Flaxman S, Mishra S, Gandy A, Unwin H, Coupland H, Mellan T, Zhu H, Berah T, Eaton J, Perez Guzman P, et al. Report 13: estimating the number of infections and the impact of nonpharmaceutical interventions on COVID19 in 11 European countries 2020. https://doi.org/10.25561/77731. Accessed 1 Mar 2021.
Chinazzi M, Davis JT, Ajelli M, Gioannini C, Litvinova M, Merler S, Pastore y Piontti A, Mu K, Rossi L, Sun K, Viboud C, Xiong X, Yu H, Halloran ME, Longini IM, Vespignani A. The effect of travel restrictions on the spread of the 2019 novel coronavirus (COVID19) outbreak. Science. 2020;368(6489):395–400.
Kraemer MU, Yang CH, Gutierrez B, Wu CH, Klein B, Pigott DM, Du Plessis L, Faria NR, Li R, Hanage WP, et al. The effect of human mobility and control measures on the COVID19 epidemic in china. Science. 2020;368(6490):493–7.
Kermack WO, McKendrick AG. A contribution to the mathematical theory of epidemics. Proc R Soc Lond. 1927;115(772):700–21.
Anderson RM, May RM. Infectious diseases of humans: dynamics and control. Oxford: Oxford University Press; 1992.
Chowell G, Brauer F. The basic reproduction number of infectious diseases: computation and estimation using compartmental epidemic models. In: Mathematical and Statistical Estimation Approaches in Epidemiology, pp. 1–30. Springer, Dordrecht 2009.
Inaba H. Agestructured population dynamics in demography and epidemiology. Singapore: Springer; 2017.
Funk S, Salathé M, Jansen VAA. Modelling the influence of human behaviour on the spread of infectious diseases: a review. J R Soc Interface. 2010;7(50):1247–56.
Merler S, Ajelli M, Pugliese A, Ferguson NM. Determinants of the spatiotemporal dynamics of the 2009 h1n1 pandemic in Europe: implications for realtime modelling. PLOS Comput Biol. 2011;7(9):1–13.
He D, Dushoff J, Day T, Ma J, Earn DJ. Inferring the causes of the three waves of the 1918 influenza pandemic in england and wales. Proc R Soc B Biol Sci. 2013;280(1766):20131345.
Wearing HJ, Rohani P, Keeling MJ. Appropriate models for the management of infectious diseases. PLOS Med. 2005;2(8):320.
Lin F, Muthuraman K, Lawley M. An optimal control theory approach to nonpharmaceutical interventions. BMC Infect Dis. 2010;10(1):32.
Perrings, C., CastilloChavez, C., Chowell, G., Daszak, P., Fenichel, E.P., Finnoff, D., Horan, R.D., Kilpatrick, A.M., Kinzig, A.P., Kuminoff, N.V., et al.: Merging economics and epidemiology to improve the prediction and management of infectious disease. EcoHealth 11(4), 464–475 (2014).
Zhou Y, Yang K, Zhou K, Liang Y. Optimal vaccination policies for an sir model with limited resources. Acta Biotheor. 2014;62(2):171–81.
Ketcheson DI. Optimal control of an sir epidemic through finitetime nonpharmaceutical intervention. 2020. arXiv preprint arXiv:2004.08848. Accessed 1 Mar 2021.
Filippov AF. Differential equations with discontinuous righthand sides: control systems. Dordrecht: Kluwer; 1988.
Lunze J, LamnabhiLagarrigue F. Handbook of hybrid systems control: theory, tools, applications. Cambridge: Cambridge University Press; 2009.
Aihara K, Suzuki H. Theory of hybrid dynamical systems and its applications to biological and medical systems. Philos Trans A Math Phys Eng Sci. 2010;368(1930):4893–914.
di Bernardo M, Hogan SJ. Discontinuityinduced bifurcations of piecewise smooth dynamical systems. Philos Trans A Math Phys Eng Sci. 2010;368(1930):4915–35.
Heemels WPMH, De Schutter B, Lunze J, Lazar M. Stability analysis and controller synthesis for hybrid dynamical systems. Philos Trans A Math Phys Eng Sci. 2010;368(1930):4937–60.
Agur Z, Cojocaru L, Mazor G, Anderson RM, Danon YL. Pulse mass measles vaccination across age cohorts. Proc Natl Acad Sci. 1993;90(24):11698–702.
Wang A, Xiao Y, Cheke RA. Global dynamics of a piecewise epidemic model with switching vaccination strategy. DCDSB. 2014;19(9):2915–40.
Chladná Z, Kopfová J, Rachinskii D, Rouf SC. Global dynamics of sir model with switched transmission rate. J Math Biol. 2020;80(4):1209–33.
Bootsma MC, Ferguson NM. The effect of public health measures on the 1918 influenza pandemic in us cities. Proc Natl Acad Sci. 2007;104(18):7588–93.
Anderson RM, Heesterbeek H, Klinkenberg D, Hollingsworth TD. How will countrybased mitigation measures influence the course of the COVID19 epidemic? Lancet. 2020;395(10228):931–4.
Di Lauro F, Kiss IZ, Miller JC. Optimal timing of oneshot interventions for epidemic control. PLOS Comput Biol. 2021;17(3):1–25.
Morris DH, Rossine FW, Plotkin JB, Levin SA. Optimal, nearoptimal, and robust epidemic control. 2020. arXiv:2004.02209. Accessed 1 Mar 2021.
Sadeghi M, Greene J, Sontag E. Universal features of epidemic models under social distancing guidelines. BioRxiv. 2020. https://doi.org/10.1016/j.arcontrol.2021.04.004.
Fujiwara N, Onaga T, Wada T, Aihara K. Effects of infection control policies on the final size. Seisan Kenkyu. 2020;72(2):141–3 (in Japanese).
Liu QH, Ajelli M, Aleta A, Merler S, Moreno Y, Vespignani A. Measurability of the epidemic reproduction number in datadriven contact networks. Proc Natl Acad Sci. 2018;115(50):12680–5.
Ferguson N, Laydon D, Nedjati Gilani G, Imai N, Ainslie K, Baguelin M, Bhatia S, Boonyasiri A, Cucunuba Perez Z, CuomoDannenburg G, et al. Report 9: impact of nonpharmaceutical interventions (NPIS) to reduce covid19 mortality and healthcare demand 2020. https://doi.org/10.25561/77482. Accessed 1 Mar 2021.
Kissler SM, Tedijanto C, Goldstein E, Grad YH, Lipsitch M. Projecting the transmission dynamics of SARSCOV2 through the postpandemic period. Science. 2020;368(6493):860–8.
Liu QH, Bento AI, Yang K, Zhang H, Yang X, Merler S, Vespignani A, Lv J, Yu H, Zhang W, et al. The COVID19 outbreak in Sichuan, China: epidemiology and impact of interventions. PLoS Comput Biol. 2020;16(12):1008467.
Aslam Noor M, Raza A, Arif MS, Rafiq M, Sooppy Nisar K, Khan I, Abdelwahab SF. Nonstandard computational analysis of the stochastic COVID19 pandemic model: an application of computational biology. Alex Eng J. 2022;61(1):619–30. https://doi.org/10.1016/j.aej.2021.06.039.
Naveed M, Baleanu D, Rafiq M, Raza A, Soori AH, Ahmed N. Dynamical behavior and sensitivity analysis of a delayed coronavirus epidemic model 2020. Comput Mater Continua. 2020;65(1):225–241. https://www.techscience.com/cmc/v65n1/39563.
Akgül A, Ahmed N, Raza A, Iqbal Z, Rafiq M, Baleanu D, Rehman MAu. New applications related to COVID19. Results Phys. 2021;20:103663. https://doi.org/10.1016/j.rinp.2020.103663.
Xiao Y, Xu X, Tang S. Sliding mode control of outbreaks of emerging infectious diseases. Bull Math Biol. 2012;74(10):2403–22.
Khalili Amirabadi R, Heydari A, Zarrabi M. Analysis and control of seir epedemic model via sliding mode control. Adv Model Optim. 2016;18(1):153–62.
Badr HS, Du H, Marshall M, Dong E, Squire MM, Gardner LM. Association between mobility patterns and covid19 transmission in the usa: a mathematical modelling study. Lancet Infect Dis. 2020. https://doi.org/10.1016/S14733099(20)305533.
Yabe T, Tsubouchi K, Fujiwara N, Wada T, Sekimoto Y, Ukkusuri SV. Noncompulsory measures sufficiently reduced human mobility in japan during the COVID19 epidemic. Sci Rep. 2020;10:18053.
Acknowledgements
K.A. and N.F. are grateful to Prof. H. Inaba for their valuable comments.
Funding
N.F and T.N. are supported by Starting Grants for Research toward Resilient Society (SGRRS) of Tohoku University and JST Mirai Program, Japan (Grant Number JPMJMI20B2). N.F. is supported by JSPS KAKENHI Grant Number JP18K11462. T.O. is supported by JSPS KAKENHI Grant Number JP19K14618. T.W. is supoprted by JSPS KAKENHI Grant Number 20H03932. S.T. is supported by JSPS KAKENHI Grant Number 21K10502. K.A. is partially supported by Moonshot R&D Grant Number JPMJMS2021, AMED under Grant Number JP21dm0307009, and Institute of AI and Beyond of UTokyo.
Author information
Authors and Affiliations
Contributions
NF, TW, and KA designed the study. NF derived the analytical results. TO, ST, and KA provided comments from theoretical viewpoints. TO conducted the numerical simulations. TW, JS, and TN discussed the implications of the study with respect to public health. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1.
Detailed discussions carried out in the main text are given.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.
About this article
Cite this article
Fujiwara, N., Onaga, T., Wada, T. et al. Analytical estimation of maximum fraction of infected individuals with oneshot nonpharmaceutical intervention in a hybrid epidemic model. BMC Infect Dis 22, 512 (2022). https://doi.org/10.1186/s12879022074035
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12879022074035
Keywords
 Infectious diseases
 Pandemics
 Nonpharmaceutical interventions
 Hybrid dynamical systems