Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Census should not start from 0 #187

Closed
cjbayesian opened this issue Mar 22, 2020 · 5 comments
Closed

Census should not start from 0 #187

cjbayesian opened this issue Mar 22, 2020 · 5 comments
Assignees
Labels
bug Something isn't working models Correct/improve the underlying models quickfix STAT We need this immediately!!!

Comments

@cjbayesian
Copy link
Collaborator

Currently, the census is computed by accumulating a proportion of new cases, and discharging them after 'length of stay'. This does not account for the currently hospitalized cases, as can be seen by this example of setting current hospitalized cases to 4000:
image

The initial model assumed the starting point to be very early on when this had much less impact, however this will be increasingly important for localities that are further along in the course of spread.

@cjbayesian cjbayesian added the bug Something isn't working label Mar 22, 2020
@cjbayesian cjbayesian added this to the Release: v0.5.0 milestone Mar 22, 2020
@PhilMiller
Copy link
Collaborator

I'm starting on this. If no one else has done so already and waves me off, I'll try to have a PR posted this evening.

@PhilMiller
Copy link
Collaborator

Ok, I took time looking at #190 and figuring out testing, so no code tonight, just thoughts.

One way to do this is to treat the input "Currently Hospitalized" as 'hidden' admissions on day 0. Since the chart won't show day 0, but the census will include them in the cumsum(), they'll at least be present in the presented count. This is slightly wrong, as a weakened form of the now-invalid very-early-starting-point assumption, since it conflates when they were actually admitted, and doesn't move them with the underlying SIR model values. By starting their stays all at day 0, at least it only over-predicts hospital/icu/vent population.

A more rigorous approach would treat the Infected portion of the population as subpools of outpatient, hosp, icu, and vent, each with their own transition/removal rates.

In the interest of expedience, I'm going to implement the first option.

@mdbecker mdbecker added models Correct/improve the underlying models quickfix STAT We need this immediately!!! labels Mar 23, 2020
@mdbecker
Copy link
Collaborator

Sounds perfect @PhilMiller. We're going to push our a fix to #190 immediately but once you've fixed this we'll push out another release. Thanks!

@quinn-dougherty
Copy link
Collaborator

was piggybacked into #273

@PhilMiller
Copy link
Collaborator

PhilMiller commented Mar 29, 2020

Yeah, what's left is a single patient difference in #352. Closing this

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working models Correct/improve the underlying models quickfix STAT We need this immediately!!!
Projects
None yet
Development

Successfully merging a pull request may close this issue.

5 participants