Simulated data and code for: The reinfection threshold, revisited

One mode by which infection-derived immunity fails is when recovery leads to a reduced but nonzero risk of reinfection. This type of partial protection is called leaky immunity, with the degree of leakiness quantified by the relative probability a previously infected individual will get infected upon exposure compared to a naively susceptible individual. Previous authors have defined the reinfection threshold, which occurs when the basic reproduction number equals the inverse of the leakiness, however, there has been some debate about whether or not this is a real threshold. Here we show how the reinfection threshold relates to two important occurrences: (1) the point at which the endemic equilibrium changes from being a stable spiral to a stable node, and (2) the point at which the rate of change of the prevalence increases the most relative to leakiness. When the recovery period is short relative to the average lifetime then both occurrences are close to the reinfection threshold. We show how these results are related to the reinfection threshold found in other models of imperfect immunity. To further demonstrate the significance of this threshold in modeling, we conducted a simulation study to evaluate some of the consequences the reinfection threshold might have in parameter estimation and modeling. Using specific parameter values chosen to reflect an acute infection, we found that the basic reproduction number values larger than that of the reinfection threshold value were less identifiable than those below the threshold.

Methods
This is the code and data used for the simulation study in "The reinfection threshold, revisited".