SIRSimple model questions

[jfunction]

Hi - thanks for this package it's very interesting.

I noticed the SIRSample dataset describes the underlying equations of an SIR model. eg, https://github.com/andy-iskauskas/hmer/blob/7a7f459ad63e2ca9baf6cabf2bc87d05f8c6d452/man/SIRSample.Rd#L31

In this model, the force of infection is effectively aSI*R/(S+I+R). Surely this should have I instead of R in the numerator here? As written, it is effectively saying that the infectious reservoir is those who have recovered.

Additionally, the definition of aSR is "Immunisation: transition rate from S to R". In the equations aSR*R is added to S' and subtracted from R' meaning that it models populations moving from R to S rather than moving from S to R - is this intended?

I may have missed it, but it would be most useful if a data generation script was included as described in https://r-pkgs.org/data.html#sec-data-data-raw - that would help to understand exactly how the data were generated.