Dataset Information

Using information theory to optimise epidemic models for real-time prediction and estimation.

ABSTRACT: The effective reproduction number, Rt, is a key time-varying prognostic for the growth rate of any infectious disease epidemic. Significant changes in Rt can forewarn about new transmissions within a population or predict the efficacy of interventions. Inferring Rt reliably and in real-time from observed time-series of infected (demographic) data is an important problem in population dynamics. The renewal or branching process model is a popular solution that has been applied to Ebola and Zika virus disease outbreaks, among others, and is currently being used to investigate the ongoing COVID-19 pandemic. This model estimates Rt using a heuristically chosen piecewise function. While this facilitates real-time detection of statistically significant Rt changes, inference is highly sensitive to the function choice. Improperly chosen piecewise models might ignore meaningful changes or over-interpret noise-induced ones, yet produce visually reasonable estimates. No principled piecewise selection scheme exists. We develop a practical yet rigorous scheme using the accumulated prediction error (APE) metric from information theory, which deems the model capable of describing the observed data using the fewest bits as most justified. We derive exact posterior prediction distributions for infected population size and integrate these within an APE framework to obtain an exact and reliable method for identifying the piecewise function best supported by available epidemic data. We find that this choice optimises short-term prediction accuracy and can rapidly detect salient fluctuations in Rt, and hence the infected population growth rate, in real-time over the course of an unfolding epidemic. Moreover, we emphasise the need for formal selection by exposing how common heuristic choices, which seem sensible, can be misleading. Our APE-based method is easily computed and broadly applicable to statistically similar models found in phylogenetics and macroevolution, for example. Our results explore the relationships among estimate precision, forecast reliability and model complexity.

SUBMITTER: Parag KV

PROVIDER: S-EPMC7360089 | biostudies-literature | 2020 Jul

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Using information theory to optimise epidemic models for real-time prediction and estimation.

Parag Kris V KV Donnelly Christl A CA

PLoS computational biology 20200701 7

The effective reproduction number, Rt, is a key time-varying prognostic for the growth rate of any infectious disease epidemic. Significant changes in Rt can forewarn about new transmissions within a population or predict the efficacy of interventions. Inferring Rt reliably and in real-time from observed time-series of infected (demographic) data is an important problem in population dynamics. The renewal or branching process model is a popular solution that has been applied to Ebola and Zika vi ...[more]

PMID: 32609732

Similar Datasets

Project description:Radiofrequency ablation is the most common minimally invasive therapy used in the United States to treat hepatocellular carcinoma and liver metastases. The ability to perform real-time temperature imaging while a patient is undergoing ablation therapy may help reduce the high recurrence rates following ablation therapy. Ultrasound echo signals undergo time shifts with increasing temperature due to sound speed and thermal expansion, which are tracked using both 1D cross correlation and 2D block matching based speckle tracking methods. In this paper, we present a quantitative evaluation of the accuracy and precision of temperature estimation using the above algorithms on both simulated and experimental data. A finite element analysis simulation of radiofrequency ablation of hepatic tissue was developed. Finite element analysis provides a method to obtain the exact temperature distribution along with a mapping of the tissue displacement due to thermal expansion. These local displacement maps were combined with the displacement due to speed of sound changes and utilized to generate ultrasound radiofrequency frames at specified time increments over the entire ablation procedure. These echo signals provide an ideal test-bed to evaluate the performance of both speckle tracking methods, since the estimated temperature results can be compared directly to the exact finite element solution. Our results indicate that the 1D cross-correlation (CC) method underestimates the cumulative displacement by 0.20mm, while the underestimation with 2D block matching (BM) is about 0.14 mm after 360 s of ablation. The 1D method also overestimates the size of the ablated region by 5.4% when compared to 2.4% with the 2D method after 720 s of ablation. Hence 2D block matching provides better tracking of temperature variations when compared to the 1D cross-correlation method over the entire duration of the ablation procedure. In addition, results obtained using 1D cross-correlation diverge from the ideal finite element results after 7 min of ablation and for temperatures greater than 65 degrees C. In a similar manner, experimental results presented using a tissue-mimicking phantom also demonstrate that the maximum percent difference with 2D block matching was 5%, when compared to 31% with the 1D method over the 700 s heating duration on the phantom.

Project description:We present Delila-genome, a software system for identification, visualization and analysis of protein binding sites in complete genome sequences. Binding sites are predicted by scanning genomic sequences with information theory-based (or user-defined) weight matrices. Matrices are refined by adding experimentally-defined binding sites to published binding sites. Delila-Genome was used to examine the accuracy of individual information contents of binding sites detected with refined matrices as a measure of the strengths of the corresponding protein-nucleic acid interactions. The software can then be used to predict novel sites by rescanning the genome with the refined matrices.Parameters for genome scans are entered using a Java-based GUI interface and backend scripts in Perl. Multi-processor CPU load-sharing minimized the average response time for scans of different chromosomes. Scans of human genome assemblies required 4-6 hours for transcription factor binding sites and 10-19 hours for splice sites, respectively, on 24- and 3-node Mosix and Beowulf clusters. Individual binding sites are displayed either as high-resolution sequence walkers or in low-resolution custom tracks in the UCSC genome browser. For large datasets, we applied a data reduction strategy that limited displays of binding sites exceeding a threshold information content to specific chromosomal regions within or adjacent to genes. An HTML document is produced listing binding sites ranked by binding site strength or chromosomal location hyperlinked to the UCSC custom track, other annotation databases and binding site sequences. Post-genome scan tools parse binding site annotations of selected chromosome intervals and compare the results of genome scans using different weight matrices. Comparisons of multiple genome scans can display binding sites that are unique to each scan and identify sites with significantly altered binding strengths.Delila-Genome was used to scan the human genome sequence with information weight matrices of transcription factor binding sites, including PXR/RXRalpha, AHR and NF-kappaB p50/p65, and matrices for RNA binding sites including splice donor, acceptor, and SC35 recognition sites. Comparisons of genome scans with the original and refined PXR/RXRalpha information weight matrices indicate that the refined model more accurately predicts the strengths of known binding sites and is more sensitive for detection of novel binding sites.

Dataset Information

Using information theory to optimise epidemic models for real-time prediction and estimation.

Publications

Using information theory to optimise epidemic models for real-time prediction and estimation.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets