/fitting_model_epidemiology

Code to estimate models in epidemiology

Primary LanguageJupyter NotebookApache License 2.0Apache-2.0

The impact of contact tracing: simulation study

I am puzzled by the fact that the current debate around a possible tracking app is focusing on the privacy issue instead of its relevance for containing the virus. If we are ready to give our geolocation data to avoid a traffic jam, I think we would do the same to save lives...

There are quite a few limitations to the current proposals for contact tracing as it is very well explained by Ross Anderson in Contact Tracing in the Real World. But technology will probably play a crucial role if it is shown to be efficient and if public trust it.

Here I try to explain the impact of contact tracing on a simple model through simulations and compare it to a testing startegy.

Here is an LinkedIn article version of the study Impact_tracing_simulation

Code is available here: Impact_tracing_simulation

Estimating models in epidemiology with optimization

This repo describes a method to find parameters for a SIR like model from a partial view of the trajectory (i.e. with hidden variables).

For example, if we have only access to the number of individuals in hospital (dashed green curve above) and the number of deceased individuals (dashed red curve above) as a function of time, we can still estimate the number of infectious individuals as shown in the picture above. The estimated number of infectious individuals is given in blue and the true number of infectious individuals (as given by our synthetic model) is given in dashed blue.

The problem and the method is described in the Synthetic_data notebook and its application to real data is done in the Real data

Here is an html version of Synthettic data

Author: Marc Lelarge

Disclaimer

I am an applied mathematician, not an epidemiologist. The main purpose of this repository is to test ideas but they need to be validated by an expert. There are still problems with the model but I think it can be fixed with your help!

Please use the code in this repository with care.