Simulation and prediction of spread of COVID-19 in The Republic of Serbia by SEAIHRDS model of disease transmission

As a response to the pandemic caused by SARS-Cov-2 virus, on 15 March 2020, the Republic of Serbia introduced comprehensive anti-epidemic measures to curb COVID-19. After a slowdown in the epidemic, on 6 May 2020, the regulatory authorities decided to relax the implemented measures. However, the epidemiological situation soon worsened again. As of 7 February 2021, a total of 406,352 cases of SARSCov-2 infection have been reported in Serbia, 4,112 deaths caused by COVID-19. In order to better understand the epidemic dynamics and predict possible outcomes, we have developed an adaptive mathematical model SEAIHRDS (S-susceptible, E-exposed, A-asymptomatic, I-infected, H-hospitalized, R-recovered, d-dead due to COVID-19 infection, S-susceptible). The model can be used to simulate various scenarios of the implemented intervention measures and calculate possible epidemic outcomes, including the necessary hospital capacities. Considering promising results regarding the development of a vaccine against COVID-19, the model is extended to simulate vaccination among different population strata. The findings from various simulation scenarios have shown that, with implementation of strict measures of contact reduction, it is possible to control COVID-19 and reduce number of deaths. The findings also show that limiting effective contacts within the most susceptible population strata merits a special attention. However, the findings also show that the disease has a potential to remain in the population for a long time, likely with a seasonal pattern. If a vaccine, with efficacy equal or higher than 65%, becomes available it could help to significantly slow down or completely stop circulation of the virus in human population. The effects of vaccination depend primarily on: 1. Efficacy of available vaccine(s), 2. Prioritization of the population categories for vaccination, and 3. Overall vaccination coverage of the population, assuming that the vaccine(s) develop solid immunity in vaccinated individuals. With expected basic reproduction number of R(o)=2.46 and vaccine efficacy of 68%, an 87% coverage would be sufficient to stop the virus circulation.