Model error: disch.rand and prog.rand both exit H compartment

[pholck]

It seems to me there is an oversight in the model, in that (in addition to fatalities) there are 2 ways to exit the H compartment. One is disch.rand and the other prog.rand. The problem is that prog.rand ALSO controls transition from I or Q to R. Further issue is that disch.rand transitions into a state ("D"?) that isn't really part of the model (and is certainly distinct from R state).

So...if one wants to have a separate probability and transition for H to R than that occurring for Q or I to R, which seems reasonable, then currently you're out of luck.

Loved the blog and impressive coding...just needs a few little tweaks to really be great! Thanks!