Dynamics and Persistence of a Generalized Multi-strain SIS Model

Autonomous differential equation compartmental models hold broad utility in epidemiology and public health. However, these models typically cannot account explicitly for myriad factors that affect the trajectory of infectious diseases, with seasonal variations in host behavior and environmental conditions as noteworthy examples. Fortunately, using non-autonomous differential equation compartmental models can mitigate some of these deficiencies, as the inclusion of time-varying parameters can account for temporally varying factors. The inclusion of these temporally varying factors does come at a cost though, as many analysis techniques, such as the use of Poincaré maps and Floquet theory, on non-autonomous differential equation compartmental models are typically only tractable numerically. Here, we illustrate a rare <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mi>n</mml:mi></mml:math> -strain generalized Susceptible-Infectious-Susceptible (SIS) compartmental model, with a general time-varying recovery rate, which features Floquet exponents that are algebraic expressions. We completely characterize the persistence and stability properties of our <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mi>n</mml:mi></mml:math> -strain generalized SIS model for <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mi>n</mml:mi> <mml:mo>≥</mml:mo> <mml:mn>1</mml:mn></mml:mrow> </mml:math> . We also derive a closed-form solution in terms of elementary functions for the single-strain SIS model, which is capable of incorporating almost any infectious period distribution. Finally, to demonstrate the applicability of our work, we apply it to recent syphilis incidence data from the United States, utilizing Akaike Information Criteria and Forecast Skill Scores to inform on the model's goodness of fit relative to complexity and the model's capacity to predict future trends.