search

alsnhll / SEIR_COVID19

SEIR model for COVID-19 infection, including different clinical trajectories of infection
194 stars 105 forks source link

Intervention day saved twice #11

Open mpascariu opened 4 years ago

mpascariu commented 4 years ago

Great job, Alison! Thanks a lot for your work.

I am exploring the possibility of modifying the app to allow for multiple interventions of different intensity (like what we see in real world) and also to export the result in a Excel file.

Checking the output of the SimSEIRintB() I can see that the day of the intervention is displayed twice (see below). Same values and I guess plotted like this.

Replicate the results with the code in runIntervention.R

> Tint = 30        # Intervention start time (days)
> simInt = SimSEIRintB(input)
> 
> outInt.df = simInt$out.df
> 
> outInt.df[30:33, ]
   time        S       E0           E1 I0       I1       I2        I3        R         D
30   29 942.5220 22.98966 4.597931e-06  0 15.53003 1.560992 0.3010213 16.99366 0.1026588
31   30 934.5564 26.05645 5.211289e-06  0 17.66609 1.779912 0.3441635 19.47826 0.1187653
32   30 934.5564 26.05645 5.211289e-06  0 17.66609 1.779912 0.3441635 19.47826 0.1187653
33   31 928.2269 27.08032 5.416064e-06  0 19.84364 2.025373 0.3931783 22.29337 0.1371733

Marius