SIS Epidemic Propagation on Hypergraphs Bodó, Ágnes; Katona, Gyula Y.; Simon, Péter L.
Bulletin of mathematical biology,
04/2016, Volume:
78, Issue:
4
Journal Article
Peer reviewed
Open access
Mathematical modelling of epidemic propagation on networks is extended to hypergraphs in order to account for both the community structure and the nonlinear dependence of the infection pressure on ...the number of infected neighbours. The exact master equations of the propagation process are derived for an arbitrary hypergraph given by its incidence matrix. Based on these, moment closure approximation and mean-field models are introduced and compared to individual-based stochastic simulations. The simulation algorithm, developed for networks, is extended to hypergraphs. The effects of hypergraph structure and the model parameters are investigated via individual-based simulation results.
For many diseases (e.g., sexually transmitted infections, STIs), most individuals are aware of the potential risks of becoming infected, but choose not to take action (‘respond’) despite the ...information that aims to raise awareness and to increases the responsiveness or alertness of the population. We propose a simple mathematical model that accounts for the diffusion of health information disseminated as a result of the presence of a disease and an ‘active’ host population that can respond to it by taking measures to avoid infection or if infected by seeking treatment early. In this model, we assume that the whole population is potentially aware of the risk but only a certain proportion chooses to respond appropriately by trying to limit their probability of becoming infectious or seeking treatment early. The model also incorporates a level of responsiveness that decays over time. We show that if the dissemination of information is fast enough, infection can be eradicated. When this is not possible, information transmission has an important effect in reducing the prevalence of the infection. We derive the full characterisation of the global behaviour of the model, and we show that the parameter space can be divided into three parts according to the global attractor of the system which is one of the two disease-free steady states or the endemic equilibrium.
Modelling epidemics on networks represents an important departure from classical compartmental models which assume random mixing. However, the resulting models are high-dimensional and their analysis ...is often out of reach. It turns out that mean-field models, low-dimensional systems of differential equations, whose variables are carefully chosen expected quantities from the exact model provide a good approximation and incorporate explicitly some network properties. Despite the emergence of such mean-field models, there has been limited work on investigating whether these can be used for inference purposes. In this paper, we consider network-based mean-field models and explore the problem of parameter identifiability when observations about an epidemic are available. Making use of the analytical tractability of most network-based mean-field models, e.g. explicit analytical expressions for leading eigenvalue and final epidemic size, we set up the parameter identifiability problem as finding the solution or solutions of a system of coupled equations. More precisely, subject to observing/measuring growth rate and final epidemic size, we seek to identify parameter values leading to these measurements. We are particularly concerned with disentangling transmission rate from the network density. To do this, we give a condition for practical identifiability and we find that except for the simplest model, parameters cannot be uniquely determined, that is, they are practically unidentifiable. This means that there exist multiple solutions (a manifold of infinite measure) which give rise to model output that is close to the data. Identifying, formalising and analytically describing this problem should lead to a better appreciation of the complexity involved in fitting models with many parameters to data.
Several network processes exhibiting transcritical bifurcation have globally stable steady states. Their dynamical behaviour is captured by a simple property of the right hand side of the ...corresponding system of ODEs. Based on this property, a class of dynamical systems is introduced, for which the local stability of the trivial steady state determines the global behaviour of the system. It is shown that this condition is satisfied by three network models, namely the individual-based and degree-based ODE approximations of SIS epidemic propagation on networks and the Hopfield model with non-negative weights. The general result enables us to describe the global behaviour of these systems that was not available for the first and third models and was proved in a significantly more complicated way for the second.
The dynamics of disease transmission strongly depends on the properties of the population contact network. Pair-approximation models and individual-based network simulation have been used extensively ...to model contact networks with non-trivial properties. In this paper, using a continuous time Markov chain, we start from the exact formulation of a simple epidemic model on an arbitrary contact network and rigorously derive and prove some known results that were previously mainly justified based on some biological hypotheses. The main result of the paper is the illustration of the link between graph automorphisms and the process of lumping whereby the number of equations in a system of linear differential equations can be significantly reduced. The main advantage of lumping is that the simplified lumped system is not an approximation of the original system but rather an exact version of this. For a special class of graphs, we show how the lumped system can be obtained by using graph automorphisms. Finally, we discuss the advantages and possible applications of exact epidemic models and lumping.
The dynamical behaviour of a continuous time recurrent neural network model with a special weight matrix is studied. The network contains several identical excitatory neurons and a single inhibitory ...one. This special construction enables us to reduce the dimension of the system and then fully characterize the local and global codimension-one bifurcations. It is shown that besides saddle-node and Andronov-Hopf bifurcations, homoclinic and cycle fold bifurcations may occur. These bifurcation curves divide the plane of weight parameters into nine domains. The phase portraits belonging to these domains are also characterized.
The epidemic threshold is probably the most studied quantity in the modelling of epidemics on networks. For a large class of networks and dynamics, it is well studied and understood. However, it is ...less so for clustered networks where theoretical results are mostly limited to idealised networks. In this paper we focus on a class of models known as pairwise models where, to our knowledge, no analytical result for the epidemic threshold exists. We show that by exploiting the presence of fast variables and using some standard techniques from perturbation theory we are able to obtain the epidemic threshold analytically. We validate this new threshold by comparing it to the threshold based on the numerical solution of the full system. The agreement is found to be excellent over a wide range of values of the clustering coefficient, transmission rate and average degree of the network. Interestingly, we find that the analytical form of the threshold depends on the choice of closure, highlighting the importance of model selection when dealing with real-world epidemics. Nevertheless, we expect that our method will extend to other systems in which fast variables are present.
An adaptive network model using
SIS
epidemic propagation with link-type-dependent link activation and deletion is considered. Bifurcation analysis of the pairwise ODE approximation and the ...network-based stochastic simulation is carried out, showing that three typical behaviours may occur; namely, oscillations can be observed besides disease-free or endemic steady states. The oscillatory behaviour in the stochastic simulations is studied using Fourier analysis, as well as through analysing the exact master equations of the stochastic model. By going beyond simply comparing simulation results to mean-field models, our approach yields deeper insights into the observed phenomena and help better understand and map out the limitations of mean-field models.
The pairwise ODE model for SIS epidemic propagation on an adaptive network with link number preserving rewiring is studied. The model, introduced by Gross et al. (Phys Rev Lett 96:208701,
2006
), ...consists of four ODEs and contains three parameters, the infection rate
τ
, the recovery rate
γ
and the rewiring rate
w
. It is proved that transcritical, saddle-node and Andronov–Hopf bifurcations may occur. These bifurcation curves are determined analytically in the
(
τ
,
w
)
parameter plane by using the parametric representation method, together with the two co-dimensional Takens–Bogdanov bifurcation point. It is shown that this parameter plane is divided into four regions by the above bifurcation curves. The possible behaviours are as follows: (a) globally stable disease-free steady state, (b) stable disease-free steady state with two unstable endemic equilibria and a stable periodic orbit, (c) stable disease-free steady state with a stable and an unstable endemic equilibrium and (d) a globally stable endemic equilibrium. Numerical evidence is shown that homoclinic bifurcation, giving rise to an unstable periodic orbit, and cycle-fold bifurcation also occur.
A key factor in the transmission of infectious diseases is the structure of disease transmitting contacts. In the context of the current COVID-19 pandemic and with some data based on the Hungarian ...population we develop a theoretical epidemic model (susceptible-infected-removed, SIR) on a multilayer network. The layers include the Hungarian household structure, with population divided into children, adults and elderly, as well as schools and workplaces, some spatial embedding and community transmission due to sharing communal spaces, service and public spaces. We investigate the sensitivity of the model (via the time evolution and final size of the epidemic) to the different contact layers and we map out the relation between peak prevalence and final epidemic size. When compared to the classic compartmental model and for the same final epidemic size, we find that epidemics on multilayer network lead to higher peak prevalence meaning that the risk of overwhelming the health care system is higher. Based on our model we found that keeping cliques/bubbles in school as isolated as possible has a major effect while closing workplaces had a mild effect as long as workplaces are of relatively small size.