  Compartmental models are an important quantitative tool in epidemiology,
enabling us to forecast the course of a communicable disease. However, the
model parameters, such as the infectivity rate of the disease, are riddled with
uncertainties, which has motivated the development and use of stochastic
compartmental models. Here, we first show that a common stochastic model, which
treats the uncertainties as white noise, is fundamentally flawed since it
erroneously implies that greater parameter uncertainties will lead to the
eradication of the disease. Then, we present a principled modeling of the
uncertainties based on reasonable assumptions on the contacts of each
individual. Using the central limit theorem and Doob's theorem on Gaussian
Markov processes, we prove that the correlated Ornstein-Uhlenbeck process is
the appropriate tool for modeling uncertainties in the infectivity rate. We
demonstrate our results using a compartmental model of the COVID-19 pandemic
and the available US data from the Johns Hopkins University COVID-19 database.
In particular, we show that the white noise stochastic model systematically
underestimates the severity of the Omicron variant of COVID-19, whereas the
Ornstein-Uhlenbeck model correctly forecasts the course of this variant.
Moreover, using an SIS model of sexually transmitted disease, we derive an
exact closed-form solution for the asymptotic distribution of infected
individuals. This analytic result shows that the white noise model
underestimates the severity of the pandemic because of unrealistic
noise-induced transitions. Our results strongly support the need for temporal
correlations in modeling of uncertainties in compartmental models of infectious
disease.
