Delay epidemic models determined by latency, infection, and immunity duration

We propose new single and two-strain epidemic models represented by systems of delay differential equations and based on the number of newly exposed individuals. Transitions between exposed, infectious, recovered, and back to susceptible compartments are determined by the corresponding time delays. Existence and positiveness of solutions are proved. Reduction of delay differential equations to integral equations allows the analysis of stationary solutions and their stability. In the case of two strains, they compete with each other, and the strain with a larger individual basic reproduction number dominates the other one. However, if the basic reproduction number exceeds some critical values, stationary solution loses its stability resulting in periodic time oscillations. In this case, both strains are present and their dynamics is not completely determined by the basic reproduction numbers but also by other parameters. The results of the work are illustrated by comparison with data on seasonal influenza.