HIV with contact-tracing: a case study in Approximate Bayesian Computation

Missing data is a recurrent issue in epidemiology where the infection process may be partially observed. Approximate Bayesian Computation, an alternative to data imputation methods such as Markov Chain Monte Carlo integration, is proposed for making inference in epidemiological models. It is a likelihood-free method that relies exclusively on numerical simulations. ABC consists in computing a distance between simulated and observed summary statistics and weighting the simulations according to this distance. We propose an original extension of ABC to path-valued summary statistics, corresponding to the cumulated number of detections as a function of time. For a standard compartmental model with Suceptible, Infectious and Recovered individuals (SIR), we show that the posterior distributions obtained with ABC and MCMC are similar. In a refined SIR model well-suited to the HIV contact-tracing data in Cuba, we perform a comparison between ABC with full and binned detection times. For the Cuban data, we evaluate the efficiency of the detection system and predict the evolution of the HIV-AIDS disease. In particular, the percentage of undetected infectious individuals is found to be of the order of 40%

Capturing the time-varying drivers of an epidemic using stochastic dynamical systems

Epidemics are often modelled using non-linear dynamical systems observed through partial and noisy data. In this paper, we consider stochastic extensions in order to capture unknown influences (changing behaviors, public interventions, seasonal effects etc). These models assign diffusion processes to the time-varying parameters, and our inferential procedure is based on a suitably adjusted adaptive particle MCMC algorithm. The performance of the proposed computational methods is validated on simulated data and the adopted model is applied to the 2009 H1N1 pandemic in England. In addition to estimating the effective contact rate trajectories, the methodology is applied in real time to provide evidence in related public health decisions. Diffusion driven SEIR-type models with age structure are also introduced.Comment: 21 pages, 5 figure

Enhancing Bayesian risk prediction for epidemics using contact tracing

Contact tracing data collected from disease outbreaks has received relatively little attention in the epidemic modelling literature because it is thought to be unreliable: infection sources might be wrongly attributed, or data might be missing due to resource contraints in the questionnaire exercise. Nevertheless, these data might provide a rich source of information on disease transmission rate. This paper presents novel methodology for combining contact tracing data with rate-based contact network data to improve posterior precision, and therefore predictive accuracy. We present an advancement in Bayesian inference for epidemics that assimilates these data, and is robust to partial contact tracing. Using a simulation study based on the British poultry industry, we show how the presence of contact tracing data improves posterior predictive accuracy, and can directly inform a more effective control strategy.Comment: 40 pages, 9 figures. Submitted to Biostatistic

A Comparative Analysis of Influenza Vaccination Programs

The threat of avian influenza and the 2004-2005 influenza vaccine supply shortage in the United States has sparked a debate about optimal vaccination strategies to reduce the burden of morbidity and mortality caused by the influenza virus. We present a comparative analysis of two classes of suggested vaccination strategies: mortality-based strategies that target high risk populations and morbidity-based that target high prevalence populations. Applying the methods of contact network epidemiology to a model of disease transmission in a large urban population, we evaluate the efficacy of these strategies across a wide range of viral transmission rates and for two different age-specific mortality distributions. We find that the optimal strategy depends critically on the viral transmission level (reproductive rate) of the virus: morbidity-based strategies outperform mortality-based strategies for moderately transmissible strains, while the reverse is true for highly transmissible strains. These results hold for a range of mortality rates reported for prior influenza epidemics and pandemics. Furthermore, we show that vaccination delays and multiple introductions of disease into the community have a more detrimental impact on morbidity-based strategies than mortality-based strategies. If public health officials have reasonable estimates of the viral transmission rate and the frequency of new introductions into the community prior to an outbreak, then these methods can guide the design of optimal vaccination priorities. When such information is unreliable or not available, as is often the case, this study recommends mortality-based vaccination priorities

Managing and reducing uncertainty in an emerging influenza pandemic

published_or_final_versio

Model-based comprehensive analysis of school closure policies for mitigating influenza epidemics and pandemics

School closure policies are among the non-pharmaceutical measures taken into consideration to mitigate influenza epidemics and pandemics spread. However, a systematic review of the effectiveness of alternative closure policies has yet to emerge. Here we perform a model-based analysis of four types of school closure, ranging from the nationwide closure of all schools at the same time to reactive gradual closure, starting from class-by-class, then grades and finally the whole school. We consider policies based on triggers that are feasible to monitor, such as school absenteeism and national ILI surveillance system. We found that, under specific constraints on the average number of weeks lost per student, reactive school-by-school, gradual, and county-wide closure give comparable outcomes in terms of optimal infection attack rate reduction, peak incidence reduction or peak delay. Optimal implementations generally require short closures of one week each; this duration is long enough to break the transmission chain without leading to unnecessarily long periods of class interruption. Moreover, we found that gradual and county closures may be slightly more easily applicable in practice as they are less sensitive to the value of the excess absenteeism threshold triggering the start of the intervention. These findings suggest that policy makers could consider school closure policies more diffusely as response strategy to influenza epidemics and pandemics, and the fact that some countries already have some experience of gradual or regional closures for seasonal influenza outbreaks demonstrates that logistic and feasibility challenges of school closure strategies can be to some extent overcome

Nonparametric survival analysis of epidemic data

This paper develops nonparametric methods for the survival analysis of epidemic data based on contact intervals. The contact interval from person i to person j is the time between the onset of infectiousness in i and infectious contact from i to j, where we define infectious contact as a contact sufficient to infect a susceptible individual. We show that the Nelson-Aalen estimator produces an unbiased estimate of the contact interval cumulative hazard function when who-infects-whom is observed. When who-infects-whom is not observed, we average the Nelson-Aalen estimates from all transmission networks consistent with the observed data using an EM algorithm. This converges to a nonparametric MLE of the contact interval cumulative hazard function that we call the marginal Nelson-Aalen estimate. We study the behavior of these methods in simulations and use them to analyze household surveillance data from the 2009 influenza A(H1N1) pandemic. In an appendix, we show that these methods extend chain-binomial models to continuous time.Comment: 30 pages, 6 figure

Estimating infectious disease parameters from data on social contacts and serological status

In dynamic models of infectious disease transmission, typically various mixing patterns are imposed on the so-called Who-Acquires-Infection-From-Whom matrix (WAIFW). These imposed mixing patterns are based on prior knowledge of age-related social mixing behavior rather than observations. Alternatively, one can assume that transmission rates for infections transmitted predominantly through non-sexual social contacts, are proportional to rates of conversational contact which can be estimated from a contact survey. In general, however, contacts reported in social contact surveys are proxies of those events by which transmission may occur and there may exist age-specific characteristics related to susceptibility and infectiousness which are not captured by the contact rates. Therefore, in this paper, transmission is modeled as the product of two age-specific variables: the age-specific contact rate and an age-specific proportionality factor, which entails an improvement of fit for the seroprevalence of the varicella-zoster virus (VZV) in Belgium. Furthermore, we address the impact on the estimation of the basic reproduction number, using non-parametric bootstrapping to account for different sources of variability and using multi-model inference to deal with model selection uncertainty. The proposed method makes it possible to obtain important information on transmission dynamics that cannot be inferred from approaches traditionally applied hitherto.Comment: 25 pages, 6 figure

Interpreting Seroepidemiologic Studies of Influenza in a Context of Nonbracketing Sera

Background: In influenza epidemiology, analysis of paired sera collected from people before and after influenza seasons has been used for decades to study the cumulative incidence of influenza virus infections in populations. However, interpretation becomes challenging when sera are collected after the start or before the end of an epidemic, and do not neatly bracket the epidemic. Methods: Serum samples were collected longitudinally in a community-based study. Most participants provided their first serum after the start of circulation of influenza A(H1N1)pdm09 virus in 2009. We developed a Bayesian hierarchical model to correct for nonbracketing sera and estimate the cumulative incidence of infection from the serological data and surveillance data in Hong Kong. Results: We analyzed 4,843 sera from 2,097 unvaccinated participants in the study, collected from April 2009 to December 2010. After accounting for nonbracketing, we estimated that the cumulative incidence of H1N1pdm09 virus infection was 45% (95% credible interval [CI] = 40%, 49%), 17% (95% CI = 13%, 20%), and 11% (95% CI = 6%, 18%) for children ages 0–18 years, adults 19–50 years, and older adults >50 years, respectively. Including all available data substantially increased precision compared with a simpler analysis based only on sera collected at 6-month intervals in a subset of participants. Conclusions: We developed a framework for the analysis of antibody titers that accounted for the timing of sera collection with respect to influenza activity and permitted robust estimation of the cumulative incidence of infection during an epidemic.postprin

The time to extinction for an SIS-household-epidemic model

We analyse a stochastic SIS epidemic amongst a finite population partitioned into households. Since the population is finite, the epidemic will eventually go extinct, i.e., have no more infectives in the population. We study the effects of population size and within household transmission upon the time to extinction. This is done through two approximations. The first approximation is suitable for all levels of within household transmission and is based upon an Ornstein-Uhlenbeck process approximation for the diseases fluctuations about an endemic level relying on a large population. The second approximation is suitable for high levels of within household transmission and approximates the number of infectious households by a simple homogeneously mixing SIS model with the households replaced by individuals. The analysis, supported by a simulation study, shows that the mean time to extinction is minimized by moderate levels of within household transmission