Skip to content

Initialization of recovered compartment fix in ODE SECIR model #1348

@HenrZu

Description

@HenrZu

Bug description

When setting the recovered compartment with real world data, we use the confirmed cases.

However, is should be equal to the secirvvs model

num_rec[i] -= num_InfectedSymptoms[i] / scaling_factor_inf[i];

So, we should take the commulative number of reported infections at t_0 (already in the code) and subtract all compartments except Susceptible in the same region/age group as explained in https://doi.org/10.1371/journal.pcbi.1010054.s002 in Eq.12.

Version

Linux

To reproduce

check the set_cases functionality.

Relevant log output

Add any relevant information, e.g. used compiler, screenshots.

Checklist

  • Attached labels, especially loc:: or model:: labels.
  • Linked to project

Metadata

Metadata

Assignees

Labels

class::bugBugs found in the softwaremodel::odeThis issue concerns any kind of ODE-based model.prio::highThe priority of this task is high.

Type

Projects

Status

Done (Total)

Milestone

No milestone

Relationships

None yet

Development

No branches or pull requests

Issue actions