rajeshrinet / pyross

PyRoss: inference, forecasts, and optimised control of epidemiological models in Python. github.com/rajeshrinet/pyross
https://pyross.readthedocs.io
MIT License
169 stars 57 forks source link

Any chance of adding asymptomatic? #6

Closed philip-mach closed 4 years ago

philip-mach commented 4 years ago

This paper models the effect of asymptomatic transmission of dengue.

I am concerned that this is a big factor being missed as there are reports that asymptomatic could be as much as 50%. If you compare Germany with Italy and Spain, Germany has an order of magnitude fewer deaths / case and its ICUs are not being swamped. It has a much higher rate of testing too. The difference cannot only be accounted for by counting more actual cases. My hypothesis: by doing aggressive track and test, the Germans reduced the next wave of contagion from asymptomatic carriers.

rajeshrinet commented 4 years ago

@philip-mach thanks for sharing the paper. We already have the asymptomatic transmission in our model. We set the fraction to zero (alpha=0 in the example)

alpha can be tuned based on the data.

technOslerphile commented 4 years ago

@rajeshrinet Two related questions:

  1. If I add a non-zero 'alpha' parameter, then do the india.txt file need to be modified so that it reflects the number of cases that would include the confirmed cases + asymptomatic cases? I am asking because the current example simulations are all based on the assumption that 'alpha' is zero.

  2. The current covid case data in 'india.txt' is updated only till March 25th. We have two more weeks of data and do you think the model predictions will be more robust if we use the updated data (let's say till April 8th)? Forget about the non-pharmacological interventions (NPIs) that started in India on March 24th and which was basically the focus of your paper. My question is only with regard to whether having the current updated data would make the model forecasts more robust. In this context, I am only interested in forecast and not effects of NPIs

rajeshrinet commented 4 years ago

@technOslerphile

  1. Yes, one needs to get estimates of asymptomatic cases or update india.txt so that true symptomatic and asymptomatic cases are listed. The parameter alpha can be modified but that means the fitting parameter beta has to be reevaluated. The data file should also be reflected in initial number of symptomatic and asymptomatic cases for the simulations.

  2. Yes, it is a good idea. I have not been able to do it myself because right now I am working on developing code for other models with more and more details. In our paper, we had only considered SIR model, but we now have SEIR, and other models. So even for the same data, it provides a way to do the model comparison...

technOslerphile commented 4 years ago

@rajeshrinet

  1. We have some decent estimates of asymptomatic proportion of cases from certain 'randomly' sampled populations in Iceland, Diamond Princess ship and Italian municipality of Vo. https://medium.com/@dennis.robert.nm/how-deadly-is-covid-19-e80c973e2823 We would get better idea with serological testing being done more widely, but as of now about 50% seems to be a reasonable informed alpha. So I was thinking of using a range of alpha and see what happens. So when you say I have to update india.txt, do you mean to update the number of cases as well? I mean, for example, on March 25th, the # is 606. If we assume alpha is 0.5, then this number should be updated to 1212 [606 (symptomatic) + 606 (asymptomatic)] . Is my understanding correct?
  2. Got it. My question was just to clarify my doubt. You have already done quite enough!!! :-)
rajeshrinet commented 4 years ago

@technOslerphile in that sense, you want to partition all reported cases in symptomatic and asymptomatic cases. That should work as well. This has to be reflected in initial number of symptomatic and asymptomatic cases for the simulations.