Qualitative study of a quarantine/isolation model with multiple disease stages

Recent studies suggest that, for disease transmission models with latent and infectious periods, the use of gamma distribution assumption seems to provide a better fit for the associated epidemiological data in comparison to the use of exponential distribution assumption. The objective of this study is to carry out a rigorous mathematical analysis of a communicable disease transmission model with quarantine (of latent cases) and isolation (of symptomatic cases), in which the waiting periods in the infected classes are assumed to have gamma distributions. Rigorous analysis of the model reveals that it has a globally-asymptotically stable disease-free equilibrium whenever its associated reproduction number is less than unity. The model has a unique endemic equilibrium when the threshold quantity exceeds unity. The endemic equilibrium is shown to be locally and globally-asymptotically stable for special cases. Numerical simulations, using data related to the 2003 SARS outbreaks, show that the cumulative number of disease-related mortality increases with increasing number of disease stages. Furthermore, the cumulative number of new cases is higher if the asymptomatic period is distributed such that most of the period is spent in the early stages of the asymptomatic compartments in comparison to the cases where the average time period is equally distributed among the associated stages or if most of the time period is spent in the later (final) stages of the asymptomatic compartments. Finally, it is shown that distributing the average sojourn time in the infectious (asymptomatic) classes equally or unequally does not effect the cumulative number of new cases.