A test for constant fatality rate of an emerging epidemic: with applications to severe acute respiratory syndrome in Hong Kong and Beijing

Summary .   The etiology, pathogenesis, and prognosis for a newly emerging disease are generally unknown to clinicians. Effective interventions and treatments at the earliest possible times are warranted to suppress the fatality of the disease to a minimum, and inappropriate treatments should be abolished. In this situation, the ability to extract most information out of the data available is critical so that important decisions can be made. Ineffectiveness of the treatment can be reflected by a constant fatality over time while effective treatment normally leads to a decreasing fatality rate. A statistical test for constant fatality over time is proposed in this article. The proposed statistic is shown to converge to a Brownian motion asymptotically under the null hypothesis. With the special features of the Brownian motion, we are able to analyze the first passage time distribution based on a sequential tests approach. This allows the null hypothesis of constant fatality rate to be rejected at the earliest possible time when adequate statistical evidence accumulates. Simulation studies show that the performance of the proposed test is good and it is extremely sensitive in picking up decreasing fatality rate. The proposed test is applied to the severe acute respiratory syndrome data in Hong Kong and Beijing.     

Case Fatality Proportion

A precise definition of case fatality proportion for compartmental disease transmission models with disease induced mortality rate is given. This is applied in classical epidemic modeling frameworks to models with multiple infectious stages, with multi-groups, with spatial patches, and with age of infection. It is shown that the case fatality proportion is the sum over all stages of the product of the probability of dying from the disease at a given stage and the probability of surviving to that stage. The derived expressions for case fatality can be used to estimate the disease induced death rates from more readily available data.     

Generation interval contraction and epidemic data analysis

The generation interval is the time between the infection time of an infected person and the infection time of his or her infector. Probability density functions for generation intervals have been an important input for epidemic models and epidemic data analysis. In this paper, we specify a general stochastic SIR epidemic model and prove that the mean generation interval decreases when susceptible persons are at risk of infectious contact from multiple sources. The intuition behind this is that when a susceptible person has multiple potential infectors, there is a "race" to infect him or her in which only the first infectious contact leads to infection. In an epidemic, the mean generation interval contracts as the prevalence of infection increases. We call this global competition among potential infectors. When there is rapid transmission within clusters of contacts, generation interval contraction can be caused by a high local prevalence of infection even when the global prevalence is low. We call this local competition among potential infectors. Using simulations, we illustrate both types of competition. Finally, we show that hazards of infectious contact can be used instead of generation intervals to estimate the time course of the effective reproductive number in an epidemic. This approach leads naturally to partial likelihoods for epidemic data that are very similar to those that arise in survival analysis, opening a promising avenue of methodological research in infectious disease epidemiology.     

Impulsive Vaccination of an SEIRS Model with Time Delay and Varying Total Population Size

Pulse vaccination is an effective and important strategy for the elimination of infectious diseases. A delayed SEIRS epidemic model with pulse vaccination and varying total population size is proposed in this paper. We point out, if R* < 1, the infectious population disappear so the disease dies out, while if R _*; > 1, the infectious population persist. Our results indicate that a long period of pulsing or a small pulse vaccination rate is sufficient condition for the permanence of the model.     

Spontaneous behavioural changes in response to epidemics

We study how spontaneous reduction in the number of contacts could develop, as a defensive response, during an epidemic and affect the course of infection events. A model is proposed which couples an SIR model with selection of behaviours driven by imitation dynamics. Therefore, infection transmission and population behaviour become dynamical variables that influence each other. In particular, time scales of behavioural changes and epidemic transmission can be different. We provide a full qualitative characterization of the solutions when the dynamics of behavioural changes is either much faster or much slower than that of epidemic transmission. The model accounts for multiple outbreaks occurring within the same epidemic episode. Moreover, the model can explain “asymmetric waves”, i.e., infection waves whose rising and decaying phases differ in slope. Finally, we prove that introduction of behavioural dynamics results in the reduction of the final attack rate.     

Risk perception and effectiveness of uncoordinated behavioral responses in an emerging epidemic

Beyond control measures imposed by public authorities, human behavioral changes can be triggered by uncoordinated responses driven by the risk perception of an emerging epidemic. In order to account for spontaneous social distancing, a model based on an evolutionary game theory framework is here proposed. Behavioral changes are modeled through an imitation process in which the convenience of different behaviors depends on the perceived prevalence of infections. Effects of misperception of risk induced by partial or incorrect information concerning the state of the epidemic are considered as well. Our findings highlight that, if the perceived risk associated to an epidemic is sufficiently large, then even a small reduction in the number of potentially infectious contacts (as a response to the epidemic) can remarkably affect the infection spread. In particular, the earlier the warning about the epidemic appears, the larger the possible reduction of the peak prevalence, and of the final epidemic size. Moreover, the epidemic spread is delayed if individuals' perception of risk is based on a memory mechanism and the risk of infection is initially overestimated. In conclusion, this analysis allows noteworthy inferences about the role of risk perception and the effectiveness of spontaneous behavioral changes during an emerging epidemic.     

An SEIS epidemic model with transport-related infection

In this paper, an SEIS epidemic model is proposed to study the effect of transport-related infection on the spread and control of infectious disease. New result implies that traveling of the exposed (means exposed but not yet infectious) individuals can bring disease from one region to other regions even if the infectious individuals are inhibited from traveling among regions. It is shown that transportation among regions will change the disease dynamics and break infection out even if infectious diseases will go to extinction in each isolated region without transport-related infection. In addition, our analysis shows that transport-related infection intensifies the disease spread if infectious diseases break out to cause an endemic situation in each region, in the sense of that both the absolute and relative size of patients increase. This suggests that it is very essential to strengthen restrictions of passengers once we know infectious diseases appeared.     

Equilibria of an Epidemic Game with Piecewise Linear Social Distancing Cost

Around the world, infectious disease epidemics continue to threaten people’s health. When epidemics strike, we often respond by changing our behaviors to reduce our risk of infection. This response is sometimes called “social distancing.” Since behavior changes can be costly, we would like to know the optimal social distancing behavior. But the benefits of changes in behavior depend on the course of the epidemic, which itself depends on our behaviors. Differential population game theory provides a method for resolving this circular dependence. Here, I present the analysis of a special case of the differential SIR epidemic population game with social distancing when the relative infection rate is linear, but bounded below by zero. Equilibrium solutions are constructed in closed-form for an open-ended epidemic. Constructions are also provided for epidemics that are stopped by the deployment of a vaccination that becomes available a fixed-time after the start of the epidemic. This can be used to anticipate a window of opportunity during which mass vaccination can significantly reduce the cost of an epidemic.     

Using Combined Diagnostic Test Results to Hindcast Trends of Infection from Cross-Sectional Data

Infectious disease surveillance is key to limiting the consequences from infectious pathogens and maintaining animal and public health. Following the detection of a disease outbreak, a response in proportion to the severity of the outbreak is required. It is thus critical to obtain accurate information concerning the origin of the outbreak and its forward trajectory. However, there is often a lack of situational awareness that may lead to over- or under-reaction. There is a widening range of tests available for detecting pathogens, with typically different temporal characteristics, e.g. in terms of when peak test response occurs relative to time of exposure. We have developed a statistical framework that combines response level data from multiple diagnostic tests and is able to ‘hindcast’ (infer the historical trend of) an infectious disease epidemic. Assuming diagnostic test data from a cross-sectional sample of individuals infected with a pathogen during an outbreak, we use a Bayesian Markov Chain Monte Carlo (MCMC) approach to estimate time of exposure, and the overall epidemic trend in the population prior to the time of sampling. We evaluate the performance of this statistical framework on simulated data from epidemic trend curves and show that we can recover the parameter values of those trends. We also apply the framework to epidemic trend curves taken from two historical outbreaks: a bluetongue outbreak in cattle, and a whooping cough outbreak in humans. Together, these results show that hindcasting can estimate the time since infection for individuals and provide accurate estimates of epidemic trends, and can be used to distinguish whether an outbreak is increasing or past its peak. We conclude that if temporal characteristics of diagnostics are known, it is possible to recover epidemic trends of both human and animal pathogens from cross-sectional data collected at a single point in time.     

Additive‐Multiplicative Regression Models for Spatio‐Temporal Epidemics

An extension of the stochastic susceptible–infectious–recovered (SIR) model is proposed in order to accommodate a regression context for modelling infectious disease data. The proposal is based on a multivariate counting process specified by conditional intensities, which contain an additive epidemic component and a multiplicative endemic component. This allows the analysis of endemic infectious diseases by quantifying risk factors for infection by external sources in addition to infective contacts. Inference can be performed by considering the full likelihood of the stochastic process with additional parameter restrictions to ensure non‐negative conditional intensities. Simulation from the model can be performed by Ogata's modified thinning algorithm. As an illustrative example, we analyse data provided by the Federal Research Centre for Virus Diseases of Animals, Wusterhausen, Germany, on the incidence of the classical swine fever virus in Germany during 1993–2004.     

Estimating absolute and relative case fatality ratios from infectious disease surveillance data

Knowing which populations are most at risk for severe outcomes from an emerging infectious disease is crucial in deciding the optimal allocation of resources during an outbreak response. The case fatality ratio (CFR) is the fraction of cases that die after contracting a disease. The relative CFR is the factor by which the case fatality in one group is greater or less than that in a second group. Incomplete reporting of the number of infected individuals, both recovered and dead, can lead to biased estimates of the CFR. We define conditions under which the CFR and the relative CFR are identifiable. Furthermore, we propose an estimator for the relative CFR that controls for time-varying reporting rates. We generalize our methods to account for elapsed time between infection and death. To demonstrate the new methodology, we use data from the 1918 influenza pandemic to estimate relative CFRs between counties in Maryland. A simulation study evaluates the performance of the methods in outbreak scenarios. An R software package makes the methods and data presented here freely available. Our work highlights the limitations and challenges associated with estimating absolute and relative CFRs in practice. However, in certain situations, the methods presented here can help identify vulnerable subpopulations early in an outbreak of an emerging pathogen such as pandemic influenza.     

Modelling the effect of urbanization on the transmission of an infectious disease

This paper models the impact of urbanization on infectious disease transmission by integrating a CA land use development model, population projection matrix model and CA epidemic model in S-Plus. The innovative feature of this model lies in both its explicit treatment of spatial land use development, demographic changes, infectious disease transmission and their combination in a dynamic, stochastic model. Heuristically-defined transition rules in cellular automata (CA) were used to capture the processes of both land use development with urban sprawl and infectious disease transmission. A population surface model and dwelling distribution surface were used to bridge the gap between urbanization and infectious disease transmission. A case study is presented involving modelling influenza transmission in Southampton, a dynamically evolving city in the UK. The simulation results for Southampton over a 30-year period show that the pattern of the average number of infection cases per day can depend on land use and demographic changes. The modelling framework presents a useful tool that may be of use in planning applications.     

The role of laboratory diagnostics in emerging viral infections: the example of the Middle East respiratory syndrome epidemic

Rapidly emerging infectious disease outbreaks place a great strain on laboratories to develop and implement sensitive and specific diagnostic tests for patient management and infection control in a timely manner. Furthermore, laboratories also play a role in real-time zoonotic, environmental, and epidemiological investigations to identify the ultimate source of the epidemic, facilitating measures to eventually control the outbreak. Each assay modality has unique pros and cons; therefore, incorporation of a battery of tests using traditional culture-based, molecular and serological diagnostics into diagnostic algorithms is often required. As such, laboratories face challenges in assay development, test evaluation, and subsequent quality assurance. In this review, we describe the different testing modalities available for the ongoing Middle East respiratory syndrome (MERS) epidemic including cell culture, nucleic acid amplification, antigen detection, and antibody detection assays. Applications of such tests in both acute clinical and epidemiological investigation settings are highlighted. Using the MERS epidemic as an example, we illustrate the various challenges faced by laboratories in test development and implementation in the setting of a rapidly emerging infectious disease. Future directions in the diagnosis of MERS and other emerging infectious disease investigations are also highlighted.     

An approximate Bayesian approach for estimation of the instantaneous reproduction number under misreported epidemic data

In epidemic models, the effective reproduction number is of central importance to assess the transmission dynamics of an infectious disease and to orient health intervention strategies. Publicly shared data during an outbreak often suffers from two sources of misreporting (underreporting and delay in reporting) that should not be overlooked when estimating epidemiological parameters. The main statistical challenge in models that intrinsically account for a misreporting process lies in the joint estimation of the time-varying reproduction number and the delay/underreporting parameters. Existing Bayesian approaches typically rely on Markov chain Monte Carlo algorithms that are extremely costly from a computational perspective. We propose a much faster alternative based on Laplacian-P-splines (LPS) that combines Bayesian penalized B-splines for flexible and smooth estimation of the instantaneous reproduction number and Laplace approximations to selected posterior distributions for fast computation. Assuming a known generation interval distribution, the incidence at a given calendar time is governed by the epidemic renewal equation and the delay structure is specified through a composite link framework. Laplace approximations to the conditional posterior of the spline vector are obtained from analytical versions of the gradient and Hessian of the log-likelihood, implying a drastic speed-up in the computation of posterior estimates. Furthermore, the proposed LPS approach can be used to obtain point estimates and approximate credible intervals for the delay and reporting probabilities. Simulation of epidemics with different combinations for the underreporting rate and delay structure (one-day, two-day, and weekend delays) show that the proposed LPS methodology delivers fast and accurate estimates outperforming existing methods that do not take into account underreporting and delay patterns. Finally, LPS is illustrated in two real case studies of epidemic outbreaks.     

Endemic persistence or disease extinction: The effect of separation into sub-communities

Consider an infectious disease which is endemic in a population divided into several large sub-communities that interact. Our aim is to understand how the time to extinction is affected by the level of interaction between communities. We present two approximations of the expected time to extinction in a population consisting of a small number of large sub-communities. These approximations are described for an SIR epidemic model, with focus on diseases with short infectious period in relation to life length, such as childhood diseases. Both approximations are based on Markov jump processes. Simulations indicate that the time to extinction is increasing in the degree of interaction between communities. This behaviour can also be seen in our approximations in relevant regions of the parameter space.     

Detecting Presymptomatic Infection Is Necessary to Forecast Major Epidemics in the Earliest Stages of Infectious Disease Outbreaks

We assess how presymptomatic infection affects predictability of infectious disease epidemics. We focus on whether or not a major outbreak (i.e. an epidemic that will go on to infect a large number of individuals) can be predicted reliably soon after initial cases of disease have appeared within a population. For emerging epidemics, significant time and effort is spent recording symptomatic cases. Scientific attention has often focused on improving statistical methodologies to estimate disease transmission parameters from these data. Here we show that, even if symptomatic cases are recorded perfectly, and disease spread parameters are estimated exactly, it is impossible to estimate the probability of a major outbreak without ambiguity. Our results therefore provide an upper bound on the accuracy of forecasts of major outbreaks that are constructed using data on symptomatic cases alone. Accurate prediction of whether or not an epidemic will occur requires records of symptomatic individuals to be supplemented with data concerning the true infection status of apparently uninfected individuals. To forecast likely future behavior in the earliest stages of an emerging outbreak, it is therefore vital to develop and deploy accurate diagnostic tests that can determine whether asymptomatic individuals are actually uninfected, or instead are infected but just do not yet show detectable symptoms.     

Network frailty and the geometry of herd immunity

The spread of infectious disease through communities depends fundamentally on the underlying patterns of contacts between individuals. Generally, the more contacts one individual has, the more vulnerable they are to infection during an epidemic. Thus, outbreaks disproportionately impact the most highly connected demographics. Epidemics can then lead, through immunization or removal of individuals, to sparser networks that are more resistant to future transmission of a given disease. Using several classes of contact networks-Poisson, scale-free and small-world-we characterize the structural evolution of a network due to an epidemic in terms of frailty (the degree to which highly connected individuals are more vulnerable to infection) and interference (the extent to which the epidemic cuts off connectivity among the susceptible population that remains following an epidemic). The evolution of the susceptible network over the course of an epidemic differs among the classes of networks; frailty, relative to interference, accounts for an increasing component of network evolution on networks with greater variance in contacts. The result is that immunization due to prior epidemics can provide greater community protection than random vaccination on networks with heterogeneous contact patterns, while the reverse is true for highly structured populations.     

Stochastic two-group models with transmission dependent on host infectivity or susceptibility

ABSTRACT

Stochastic epidemic models with two groups are formulated and applied to emerging and re-emerging infectious diseases. In recent emerging diseases, disease spread has been attributed to superspreaders, highly infectious individuals that infect a large number of susceptible individuals. In some re-emerging infectious diseases, disease spread is attributed to waning immunity in susceptible hosts. We apply a continuous-time Markov chain (CTMC) model to study disease emergence or re-emergence from different groups, where the transmission rates depend on either the infectious host or the susceptible host. Multitype branching processes approximate the dynamics of the CTMC model near the disease-free equilibrium and are used to estimate the probability of a minor or a major epidemic. It is shown that the probability of a major epidemic is greater if initiated by an individual from the superspreader group or by an individual from the highly susceptible group. The models are applied to Severe Acute Respiratory Syndrome and measles.     

Modeling and real-time prediction of classical swine fever epidemics

We propose a new method to analyze outbreak data of an infectious disease such as classical swine fever. The underlying model is a two-type branching process. It is used to deduce information concerning the epidemic from detected cases. In particular, the method leads to prediction of the future course of the epidemic and hence can be used as a basis for control policy decisions. We test the model with data from the large 1997-1998 classical swine fever epidemic in The Netherlands. It turns out that our results are in good agreement with the data.