Stochastic dynamics of an SIS epidemic on networks.

We derive a stochastic SIS pairwise model by considering the change of the variables of this system caused by an event. Based on approximations, we construct a low-dimensional deterministic system that can be used to describe the epidemic spread on a regular network. The mathematical treatment of the model yields explicit expressions for the variances of each variable at equilibrium. Then a comparison between the stochastic pairwise model and the stochastic mean-field SIS model is performed to indicate the effect of network structure. We find that the variances of the prevalence of infection for these two models are almost equal when the number of neighbors of every individual is large. Furthermore, approximations for the quasi-stationary distribution of the number of infected individuals and the expected time to extinction starting in quasi-stationary are derived. We analyze the approximations for the critical number of neighbors and the persistence threshold based on the stochastic model. The approximate performance is then examined by numerical and stochastic simulations. Moreover, during the early development phase, the temporal variance of the infection is also obtained. The simulations show that our analytical results are asymptotically accurate and reasonable.