Mathematical Analysis of Epidemic Models with Treatment in Heterogeneous Networks

In this paper, we formulate two different network-based epidemic models to investigate the effect of partly effective treatment on disease dynamics. The first network model represents the individuals with heterogeneous number of contacts in a population as choosing a new partner at each moment, whereas the second one assumes the individuals have fixed or stable neighbors. The basic reproduction number R0 is computed for each model, using the next generation matrix method. In particular, the critical treatment rate is defined for the model, above which the disease can be eliminated through the treatment. The final epidemic size relations are derived, and the solvability of these implicit equations is studied. In particular, a unique solution of the implicit equation for the final epidemic size is determined, and by rewriting the implicit equation as a suitable fixed point problem, it is proved that the iteration of the fixed point problem converges to the unique solution. Stochastic simulations and numerical simulations, including in comparison with the model outputs and the joint influence of network topology and treatment on the final epidemic size, are conducted to illustrate the theoretical results.