key: cord-0962846-k3zxudok
authors: Schuster, Peter
title: Ebola—challenge and revival of theoretical epidemiology: Why Extrapolations from early phases of epidemics are problematic
date: 2015-04-28
journal: Complexity
DOI: 10.1002/cplx.21694
sha: 011fc6138b216a77c6b91bf22bad41849efc524f
doc_id: 962846
cord_uid: k3zxudok

nan

model is formulated in the language of ordinary differential equations and can be solved by the standard analytical techniques of calculus dN dt 5 rN 12 N C ; NðtÞ 5 N 0 C N 0 1 ðC 2 N 0 Þ e 2 rt : (1) Herein, N(t) is the number of multiplying individuals or the population size as a function of time, N 0 5 N(0) is the initial value, r is the growth or Malthus parameter, and C is called the carrying capacity and represents a parameter, which determines the maximum number of individuals that can be sustained by the ecosystem (Figure 1 ).

The initial phase of the process (t < 3) yields practically the same curves for all values of C within the line width. It is straightforward to calculate an approximate expression for N(t) at short times by means of series expansions that has the following two first terms:

for small t Equation (2) illustrates well the problem of parameter fitting to data from the initial phase: In cases where N 0 is small and C is large and N 0 /C << 1, the growth parameter r can be determined from the slope of the linearized growth curve, N(t) % N 0 1 N 0 r t, as N 0 /C is only a small correction. The determination of the second parameter C is tricky, in particular for small N 0 /C because the first and the second coefficient of the series expansion are needed and the equations for the fit of C are commonly ill-conditioned. Unfortunately small N 0 and large C, being tantamount to a few infected individuals-in extreme but not uncommon cases we have N 0 5 1-and a comparatively large population of susceptible individuals, are almost always the case in the spread of epidemics. For very good reasons, Aaron King et al. [5] emphasize the requirement of highly reliable data and sophisticated high-level stochastic methods for parameter estimations and potential error predictions from the early phase of an exponentially growing system to a final situation when the process is under control. The problem is very general and not restricted to epidemiology. The politically reasonable but scientifically unacceptable flaw in the predictions of the famous Club of Rome [7] is precisely of the same kind: You cannot extrapolate from the initial phase of an exponentially growing process to the final situation unless additional information is or additional data are available.

The SEIR model of epidemiology is very popular and despite relative simplicity it has interesting dynamical features that are highly relevant for prediction. The model describes infection dynamics in a population of initially susceptible individuals and attributes four different states to the people: S means susceptible to infection, E implies exposed and infected but not yet showing the symptoms and not being infectious, I stands for being infected, showing the symptoms and being infectious, and eventually R stands for the number of individuals being out of the game by either being recovered and immune to further infection or being dead. For population dynamics, the SEIR model can be casted into equations of chemical kinetics:

The rate parameters commonly attributed to the three steps are denoted by b/C, k, and c. As a matter of fact, the epidemiological mechanism (3) is nothing but a chemical reaction consisting of three consecutive steps, and it can be cast into the differential equation

Herein, concentrations or numbers of individuals are denoted by [S] 5 s,

Because of the conservation relation, s 1 i 1 e 1 r 5 C, the model has three independent variables. To take into account intrinsic stochastic effects, the SEIR model can be formulated as a renewal process [4] that is frequently used in theoretical epidemiology [8] . To visualize fluctuations, we shall consider here the expectation values embedded in the one-standard deviation band derived from numerical simulations of the chemical master equation of mechanism (3) [9] . The results are shown in Figure 2 : As expected the solution curves of the deterministic approach lie within the one-standard deviation band but do not exactly coincide with the expectation values E(x), x s,e,i or r-exact coincidence can only be expected for first order or linear processes [10] and, indeed, the step (3a) is corresponding to a bimolecular second order reaction. The SEIR mechanism contains two-step autocatalysis: S 1 I ! E 1 I ! 2 I resulting from (3a) and (3b). One special property of an autocatalytic system concerns the requirement of a non-zero concentration of the autocatalyst as no reaction will take place otherwise. In the language of epidemiology this fact is trivial: There is no outbreak without infectious individuals. Autocatalysis gives also rise to a number of features, which are uncommon in conventional chemical kinetics. One property that is relevant for the problem of prediction will be discussed in detail in the next paragraph. It concerns the longtime result, which in chemistry is commonly a single state that does not depend on initial conditions. In case of autocatalysis, however, the final state can be a one-dimensional (1D) manifold and the particular marginally stable stationary points depend on the initial values of the variables (Figure 3 ; see also a somewhat humoristic but mathematically serious analysis of the equilibria of epimemiological models and their stability see [11] ). In the SEIR and as well in the simpler

Solution curves of the SEIR model (3) . Shown are the numbers of susceptible individuals S(t) (red), exposed individuals E(t) (blue), infectious individuals I(t) (green), and recovered or dead individuals R(t) (black). Expectation values <X(t)> (X S, E, I, R) are given in full color embedded in onestandard deviation error bands, <X(t)> 

The number of individuals, which are not infected by the disease is limited by 0 s < N. Stability analysis of this 1D line of stationary points by means of the evaluation of the eigenvalues of the Jacobian matrix at the steady states yields three eigenvalues, one is zero and two are real:

The zero eigenvalue expresses marginal stability and concerns migration along the line of stationary states. The other two eigenvalues are both negative in the range 0 s < s cr with s cr 5 c N =b and depending on the initial conditions the longtime solution will lie somewhere in this range. An example is shown in Figure 3 , the larger the initial number of infected individuals is, the smaller is s, the number of people that remained unaffected by the epidemic (The parameter values for the calculations are taken from the 1995 Ebola epidemic in Congo [12] ).

In this context, it is illustrative to mention a result of the SIR model in the simplest version, which does not consider exposed individuals-as said above, S 1 I ! 2 I and I ! R with the two rate parameters b and c. For simplicity, we replace the numbers of individuals by normalized variables s ) s/C, i ) i/C, r ) r/C such that s 1 i 1 r 5 1. The analysis of the system reveals again a one-dimensional manifold of stationary points with i 5 0 and s 1 r 5 1: Stability analysis yields stable points in the range 0 s < s cr and unstable points for s cr < s 1, where s cr 5 c=b. In Figure 4 , it is shown how the initial ratio of susceptible to recovered or here better immune individuals determines the final number of people that were affected by the epidemic: The more susceptibles are initially present, the smaller is the number of individuals, which remain spared from the disease. Moreover, in the range 0 s < s cr no epidemic can occur and this has a straightforward consequence for vaccination: If the fraction of immune people, r/(r 1 s), is above 1 2 s cr , the population cannot be invaded by the proliferating disease but, of course, the nonvaccinated individual can still be infected. From c being larger than b follows s cr > 1 and the population is stable in the whole range, 0 s 1, and no outbreak can occur. The fraction of immune people is very hard to know a priori and this implies another uncertainty in modeling.

The parameter b determines the rate of infection and may be controlled by successful intervention strategies. In modeling this can be taken into account using the parameter as a function of time [12] 

The parameter s is the instant when the interventions start and q is

Trajectories of the simple SIR model. Trajectories of individual epidemics, #(t), calculated by the simple SIR model are shown on the unit simplex S 3 , #(t) 5 {s(t), i(t), r(t)} with s 1 i 1 r 5 1. Every trajectory starts at t 5 0 at the right end and progresses with increasing time to the left end that corresponds to the limit t ! 1. The initial conditions are given by s 0 5 a 2 d, i 0 5 d, and r 0 5 1 2 a, with d 5 0.001 and a 5 1.0, 0.9, 0.8, 0.7, 0.6, 0.5, and 0.4, or in other words a population of a susceptible and 1 2 a immune individuals is infected by a small fraction of diseased people. The number of infectious individual increases, passes through a maximum at s cr 5 g=b, and for long times converges to lim t!1 i(t) 5 0, the disease has disappeared. Parameters: b 5 3, c 5 1, and s cr 5 1=3. *The SIR model is distiguished from the SEIR model by the omission of the state ''exposed.'' It contains one-step autocatalysis, S 1 I ! 2 I , and has many features in common with the SEIR model, in particular, the outcome of the SIR process depends also on the initial conditions (see next paragraph). the time constant of the implementation of the countermeasures. This approach is particularly appropriate for Ebola in developing countries since b 1 is there much smaller than b 0 . The quantity, which is the key to modeling the course of an epidemic outbreak, is the basic reproduction number R 0 , which is a measure of the average number of secondary infections caused by a single primary case in a pool of mostly susceptible individuals [12] . For R 0 > 1, the epidemic is spreading whereas R 0 < 1 implies that the infections are under control. The primary goal of politics and healthcare is to drive b(t) below one. The basic reproduction number in a population before intervention is simply given by R 0 5 b 0 /c, and the relation to the critical fraction of susceptible individuals is s cr =N 5 c/b 0 5 R 21 0 . Although being a flexible and appropriate approximation for real situations, Eq. (6) provides an additional problem for the modelers: It contains three new parameters b 1, s, and q that have to be derived from data and further information about the course of the epidemic.

Summarizing the current situation in modeling epidemics in view of the recent Ebola outbreak in West Africa we come to the main conclusions: (i) The epidemic provided a unique occasion for testing epidemiological models. (ii) Extrapolations from the phase of exponential growth to the final outcome-as rigorous and careful they might be [13] -are very risky because they suffer not only from external uncertainties but also from intrinsic instabilities that are common in proc-esses with autocatalytic steps. Elaborate statistics can substantially improve the situation and reduce the problem if enough and reliable data are available [5] . Low robustness of the predictions results also from dynamical systems that are more complex than conventional chemical reaction networks. Simple models can, nevertheless, give important insights and help to understand complex relations.

In Table 1 , finally we present some recent data in the light of the estimates of Ebola cases in early 2015 that were calculated to reach from several hundred thousand to even a million and more. Fortunately enough the severeness of the Ebola epidemic in West Africa so far has been overestimated by orders of magnitude. In addition, the dynamics of the epidemic was and is different in the three countries. At the beginning of March Liberia had reported no new case since 42 days but a few single cases were reported after that, whereas in Guinea and Sierra Leone the disease continues at a slowly down going rate.

The recent Ebola outbreak is so well-documented that it provides a unique opportunity for theoretical epidemiologists to test the existing and to develop new models. The take-home lesson nevertheless is: whenever one cannot avoid extrapolating a process from the growth phase to saturation, special care in needed because the early data contain little information on the parameters that determine the late behavior. Complex dynamics may enlarge the uncertainties. 

Anti-vaccinationists past and present

The age-old struggle against the antivaccinationists

Transmission dynamics and control of Ebola virus disease (EVD): A review

Theoretical perspectives on the infectiousness of Ebola virus disease

Avoidable errors in the modelling of outbreaks of emerging pathogens, with special reference to Ebola

Notice sur la loi que la population pursuit dans son accroisement

The Limits of Growth

Mathematical Epidemiology of Infectious Diseases: Model Building, Analysis and Interpretation

Stochastic simulation of chemical kinetics

Stochastic Processes in Physics and Chemisty

When zombies attack: Mathematical modelling of an outbreak of zombie infection

The basic reproductive number of Ebola and the effect of public health measures: The cases of Congo and Uganda

Infectious Diseases of Humans

Exponential growth in Ebola. Outbreak since May 14

Ebola Epidemiological Update No

World Health Organization (WHO)