Source author record

Péter L. Simon

Péter L. Simon appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
13topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2022arXiv

Learning the parameters of a differential equation from its trajectory via the adjoint equation

The paper contributes to strengthening the relation between machine learning and the theory of differential equations. In this context, the inverse problem of fitting the parameters, and the initial condition of a differential equation to some measurements constitutes a key issue. The paper explores an abstraction that can be used to construct a family of loss functions with the aim of fitting the solution of an initial value problem to a set of discrete or continuous measurements. It is shown, that an extension of the adjoint equation can be used to derive the gradient of the loss function as a continuous analogue of backpropagation in machine learning. Numerical evidence is presented that under reasonably controlled circumstances the gradients obtained this way can be used in a gradient descent to fit the solution of an initial value problem to a set of continuous noisy measurements, and a set of discrete noisy measurements that are recorded at uncertain times.

preprint2022arXiv

On parameter identifiability in network-based epidemic models

Many models in mathematical epidemiology are developed with the aim to provide a framework for parameter estimation and then prediction. It is well-known that parameters are not always uniquely identifiable. 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 define strong and weak identifiability and we find that except for the simplest model, parameters cannot be uniquely determined, that is they are weakly identifiable. This means that there exists 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.

preprint2015arXiv

Bounds for the expected value of one-step processes

Mean-field models are often used to approximate Markov processes with large state-spaces. One-step processes, also known as birth-death processes, are an important class of such processes and are processes with state space $\{0,1,\ldots,N\}$ and where each transition is of size one. We derive explicit bounds on the expected value of such a process, bracketing it between the mean-field model and another simple ODE. Our bounds require that the Markov transition rates are density dependent polynomials that satisfy a sign condition. We illustrate the tightness of our bounds on the SIS epidemic process and the voter model.

preprint2015arXiv

SIS epidemic propagation on hypergraphs

Mathematical modeling 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.

preprint2015arXiv

Super compact pairwise model for SIS epidemic on heterogeneous networks

In this paper we provide the derivation of a super compact pairwise model with only 4 equations in the context of describing susceptible-infected-susceptible (SIS) epidemic dynamics on heterogenous networks. The super compact model is based on a new closure relation that involves not only the average degree but also the second and third moments of the degree distribution. Its derivation uses an a priori approximation of the degree distribution of susceptible nodes in terms of the degree distribution of the network. The new closure gives excellent agreement with heterogeneous pairwise models that contain significantly more differential equations.

preprint2014arXiv

Dynamic control of modern, network-based epidemic models

In this paper we make the first steps to bridge the gap between classic control theory and modern, network-based epidemic models. In particular, we apply nonlinear model predictive control (NMPC) to a pairwise ODE model which we use to model a susceptible-infectious-susceptible (SIS) epidemic on non-trivial contact structures. While classic control of epidemics concentrates on aspects such as vaccination, quarantine and fast diagnosis, our novel setup allows us to deliver control by altering the contact network within the population. Moreover, the ideal outcome of control is to eradicate the disease while keeping the network well connected. The paper gives a thorough and detailed numerical investigation of the impact and interaction of system and control parameters on the controllability of the system. The analysis reveals, that for certain set parameters it is possible to identify critical control bounds above which the system is controllable. We foresee, that our approach can be extended to even more realistic or simulation-based models with the aim to apply these to real-world situations.

preprint2014arXiv

Oscillating epidemics in a dynamic network model: stochastic and mean-field analysis

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. A compact pairwise approximation for the dynamic network case is also developed and, for the case of link-type independent rewiring, the outcome of epidemics and changes in network structure are concurrently presented in a single bifurcation diagram. 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.

preprint2014arXiv

PDE approximation of large systems of differential equations

A large system of ordinary differential equations is approximated by a parabolic partial differential equation with dynamic boundary condition and a different one with Robin boundary condition. Using the theory of differential operators with Wentzell boundary conditions and similar theories, we give estimates on the order of approximation. The theory is demonstrated on a voter model where the Fourier method applied to the PDE is of great advantage.

preprint2013arXiv

Exact deterministic representation of Markovian SIR epidemics on networks with and without loops

In a previous paper Sharkey et al. [13] proved the exactness of closures at the level of triples for Markovian SIR (susceptible-infected-removed) dynamics on tree-like networks. This resulted in a deterministic representation of the epidemic dynamics on the network that can be numerically evaluated. In this paper, we extend this modelling framework to certain classes of networks exhibiting loops. We show that closures where the loops are kept intact are exact, and lead to a simplified and numerically solvable system of ODEs (ordinary-differential-equations). The findings of the paper lead us to a generalisation of closures that are based on partitioning the network around nodes that are cut-vertices (i.e. the removal of such a node leads to the network breaking down into at least two disjointed components or subnetworks). Exploiting this structural property of the network yields some natural closures, where the evolution of a particular state can typically be exactly given in terms of the corresponding or projected sates on the subnetworks and the cut-vertex. A byproduct of this analysis is an alternative probabilistic proof of the exactness of the closures for tree-like networks presented in Sharkey et al. [13]. In this paper we also elaborate on how the main result can be applied to more realistic networks, for which we write down the ODEs explicitly and compare output from these to results from simulation. Furthermore, we give a general, recipe-like method of how to apply the reduction by closures technique for arbitrary networks, and give an upper bound on the maximum number of equations needed for an exact representation.

preprint2012arXiv

Identification of criticality in neuronal avalanches: I. A theoretical investigation of the non-driven case

In this paper we study a simple model of a purely excitatory neural network that, by construction, operates at a critical point. This model allows us to consider various markers of criticality and illustrate how they should perform in a finite-size system. By calculating the exact distribution of avalanche sizes we are able to show that, over a limited range of avalanche sizes which we precisely identify, the distribution has scale free properties but is not a power law. This suggests that it would be inappropriate to dismiss a system as not being critical purely based on an inability to rigorously fit a power law distribution as has been recently advocated. In assessing whether a system, especially a finite-size one, is critical it is thus important to consider other possible markers. We illustrate one of these by showing the divergence of susceptibility as the critical point of the system is approached. Finally, we provide evidence that power laws may underlie other observables of the system, that may be more amenable to robust experimental assessment.

preprint2011arXiv

Differential equation approximations of stochastic network processes: an operator semigroup approach

The rigorous linking of exact stochastic models to mean-field approximations is studied. Starting from the differential equation point of view the stochastic model is identified by its Kolmogorov equations, which is a system of linear ODEs that depends on the state space size ($N$) and can be written as $\dot u_N=A_N u_N$. Our results rely on the convergence of the transition matrices $A_N$ to an operator $A$. This convergence also implies that the solutions $u_N$ converge to the solution $u$ of $\dot u=Au$. The limiting ODE can be easily used to derive simpler mean-field-type models such that the moments of the stochastic process will converge uniformly to the solution of appropriately chosen mean-field equations. A bi-product of this method is the proof that the rate of convergence is $\mathcal{O}(1/N)$. In addition, it turns out that the proof holds for cases that are slightly more general than the usual density dependent one. Moreover, for Markov chains where the transition rates satisfy some sign conditions, a new approach for proving convergence to the mean-field limit is proposed. The starting point in this case is the derivation of a countable system of ordinary differential equations for all the moments. This is followed by the proof of a perturbation theorem for this infinite system, which in turn leads to an estimate for the difference between the moments and the corresponding quantities derived from the solution of the mean-field ODE.

preprint2011arXiv

Exact and approximate epidemic models on networks: a new, improved closure relation

Recently, research that focuses on the rigorous understanding of the relation between simulation and/or exact models on graphs and approximate counterparts has gained lots of momentum. This includes revisiting the performance of classic pairwise models with closures at the level of pairs and/or triples as well as effective-degree-type models and those based on the probability generating function formalism. In this paper, for a fully connected graph and the simple $SIS$ (susceptible-infected-susceptible) epidemic model, a novel closure is introduced. This is done via using the equations for the moments of the distribution describing the number of infecteds at all times combined with the empirical observations that this is well described/approximated by a binomial distribution with time dependent parameters. This assumption allows us to express higher order moments in terms of lower order ones and this leads to a new closure. The significant feature of the new closure is that the difference of the exact system, given by the Kolmogorov equations, from the solution of the newly defined approximate system is of order $1/N^2$. This is in contrast with the $\mathcal{O}(1/N)$ difference corresponding to the approximate system obtained via the classic triple closure.