Adaptive estimation for epidemic renewal and phylogenetic skyline models
File(s)syaa035.pdf (1.73 MB)
Published version
Author(s)
Parag, Kris
Donnelly, Christl
Type
Journal Article
Abstract
Estimating temporal changes in a target population from phylogenetic or count data is an important problem in ecology and epidemiology. Reliable estimates can provide key insights into the climatic and biological drivers influencing the diversity or structure of that population and evidence hypotheses concerning its future growth or decline. In infectious disease applications, the individuals infected across an epidemic form the target population. The renewal model estimates the effective reproduction number, R, of the epidemic from counts of observed incident cases. The skyline model infers the effective population size, N, underlying a phylogeny of sequences sampled from that epidemic. Practically, R measures ongoing epidemic growth while N informs on historical caseload. While both models solve distinct problems, the reliability of their estimates depends on p-dimensional piecewise-constant functions. If p is misspecified, the model might underfit significant changes or overfit noise and promote a spurious understanding of the epidemic, which might misguide intervention policies or misinform forecasts. Surprisingly, no transparent yet principled approach for optimizing p exists. Usually, p is heuristically set, or obscurely controlled via complex algorithms. We present a computable and interpretable p-selection method based on the minimum description length (MDL) formalism of information theory. Unlike many standard model selection techniques, MDL accounts for the additional statistical complexity induced by how parameters interact. As a result, our method optimizes p so that R and N estimates properly and meaningfully adapt to available data. It also outperforms comparable Akaike and Bayesian information criteria on several classification problems, given minimal knowledge of the parameter space, and exposes statistical similarities among renewal, skyline, and other models in biology. Rigorous and interpretable model selection is necessary if trustworthy and justifiable conclusions are to be drawn from piecewise models. [Coalescent processes; epidemiology; information theory; model selection; phylodynamics; renewal models; skyline plots]
Date Issued
2020-11-01
Date Acceptance
2020-04-16
Citation
Systematic Biology, 2020, 69 (6), pp.1163-1179
ISSN
1063-5157
Publisher
Oxford University Press (OUP)
Start Page
1163
End Page
1179
Journal / Book Title
Systematic Biology
Volume
69
Issue
6
Copyright Statement
© The Author(s) 2020. Published by Oxford University Press, on behalf of the Society of Systematic Biologists.
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.
Sponsor
Medical Research Council (MRC)
Grant Number
MR/R015600/1
Subjects
Science & Technology
Life Sciences & Biomedicine
Evolutionary Biology
Coalescent processes
epidemiology
information theory
model selection
phylodynamics
renewal models
skyline plots
DEMOGRAPHIC HISTORY
POPULATION-DYNAMICS
INFERENCE
FRAMEWORK
Classification
Epidemiologic Methods
Models, Biological
Models, Statistical
Phylogeny
Epidemiologic Methods
Models, Statistical
Phylogeny
Models, Biological
Classification
Evolutionary Biology
0603 Evolutionary Biology
0604 Genetics
Publication Status
Published
Date Publish Online
2020-04-25