simulating with default initial state and parameters, it looks like what should be the inflow into $E_y$ (force of infection on age group $y$) gets labelled as the inflow into $R_y$:
(indeed, this is verified by manual calculation). mismatched label is also true for inflow reported into $E_y$; it's equal to the value that should be reported for $I_y$.
this may be because the state vector lists R states right after S states and there's an order dependency we were previously unaware of (despite the fact that the state vector is named!):
i imagine this cascades down into the simulated values reported under state. this also affects the main package vignette, which uses the old-and-young model for the first demo. it should not affect other models, which were checked.
either:
- fix the implementation (may already be done with kevin's re-implementation)
- remove this model from the package (it's never really used) and update the vignette to use a new model instead
simulating with default initial state and parameters, it looks like what should be the inflow into$E_y$ (force of infection on age group $y$ ) gets labelled as the inflow into $R_y$ :
(indeed, this is verified by manual calculation). mismatched label is also true for inflow reported into$E_y$ ; it's equal to the value that should be reported for $I_y$ .
this may be because the state vector lists
Rstates right afterSstates and there's an order dependency we were previously unaware of (despite the fact that the state vector is named!):i imagine this cascades down into the simulated values reported under
state. this also affects the main package vignette, which uses theold-and-youngmodel for the first demo. it should not affect other models, which were checked.either: