Skip to content

Commit

Permalink
Update README.md
Browse files Browse the repository at this point in the history
  • Loading branch information
kgostic authored Nov 17, 2020
1 parent 9e2fb78 commit f10b909
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions README.md
Original file line number Diff line number Diff line change
@@ -1,6 +1,6 @@
# Example code for deconvolving epidemic curves

Ed Baskerville, Timothy M Pollington & Katie Gostic<br>
Ed Baskerville, Lauren McGough, Timothy M Pollington & Katie Gostic<br>
19 July 2020

This repository contains example code for deconvolving a known delay distribution from observed epidemic time series.
Expand All @@ -10,4 +10,4 @@ None of the code here is ready for prime time, but we’re sharing it in case it
There is code here for two methods, both with attendant warnings:

* The Richardson-Lucy-type deconvolution used in [Goldstein et al. 2009 PNAS](https://doi.org/10.1073/pnas.0902958106), implemented straightforwardly in unoptimized R code. Either a bug in this code or an inherent property of the algorithm makes it extremely sensitive to initial guess. Certainly the observed curve needs to have structure of some kind (i.e. look like a curve) to have a chance of deconvolving into the infection epidemic curve.
* An analogous Bayesian model, implemented straightforwardly in Stan. There are two versions of the model, one (`uncorrelated`) where unobserved states are given identical independent priors, and one (`randomwalk`) where log(unobserved states) follow a Brownian motion. Neither of these implementations have been used for real work, but initial runs seem promising for doing things this way.
* An analogous Bayesian model, implemented straightforwardly in Stan. There are two versions of the model, one (`uncorrelated`) where unobserved states are given identical independent priors, and one (`randomwalk`) where log(unobserved states) follow a Brownian motion. Neither of these implementations have been used for real work, but initial runs seem promising for doing things this way.

0 comments on commit f10b909

Please sign in to comment.