A network SIRX model for the spreading of COVID-19

Abstract

Infectious diseases, such as the current COVID-19, have a huge economic and societal impact. The ability to model its transmission characteristics is critical to minimize its impact. In fact, predicting how fast an infection is spreading could be a major factor in deciding on the severity, extent and strictness of the applied mitigation measures, such as the recent lockdowns. Even though modelling epidemics is a well studied subject, usually models do not include quarantine or other social measures, such as those imposed in the recent pandemic. The current work builds upon a recent paper by Maier and Brockmann (2020), where a compartmental SIRX model was implemented. That model included social or individual behavioural changes during quarantine, by introducing state X, in which symptomatic quarantined individuals are not transmitting the infection anymore, and described well the transmission in the initial stages of the infection. The results of the model were applied to real data from several provinces in China, quite successfully. In our approach we use a Monte-Carlo simulation model on networks. Individuals are network nodes and the links are their contacts. We use a spreading mechanism from the initially infected nodes to their nearest neighbours, as has been done previously. Initially, we find the values of the rate constants (parameters) the same way as in Maier and Brockmann (2020) for the confirmed cases of a country, on a daily basis, as given by the Johns Hopkins University. We then use different types of networks (random Erdős–Rényi, Small World, and Barabási–Albert Scale-Free) with various characteristics in an effort to find the best fit with the real data for the same geographical regions as reported in Maier and Brockmann (2020). Our simulations show that the best fit comes with the Erdős–Rényi random networks. We then apply this method to several other countries, both for large-size countries, and small size ones. In all cases investigated we find the same result, i.e. best agreement for the evolution of the pandemic with time is for the Erdős–Rényi networks. Furthermore, our results indicate that the best fit occurs for a random network with an average degree of the order of 〈k〉≈ 10–25, for all countries tested. Scale Free and Small World networks fail to fit the real data convincingly.