What would be a valid observation model for example for the number of infected people (I_obs
) in a SIR model? In this case study, a negative binomial distribution is used with mean equal to the ODE solution obtained from the SIR model (I_ODE
) and variance that depends on I_ODE
and a noise parameter phi
. However, it places probability mass also for I_obs > population_size
which doesn’t really make sense as it’s not possible that somehow the measured number of infected people is larger than the population size.
Here are visualizations of the prior predictive distribution for the above case study.
Although 0.95 interval seems to be mostly under the population size, there are still draws that are above:
EDIT: the first plot is the same as this plot
in Section 2 of the case study, which shows 90% interval.